var/da/da_ssmi/da_transform_xtozrhoq_adj.inc

   1 subroutine da_transform_xtozrhoq_adj(grid, i, j, zh, zf, adj_zrhom)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (domain), intent(inout) :: grid
  10    integer,       intent(in)    :: i, j
  11    real,          intent(in)    :: ADJ_zrhom  ! gridded height weighted moisture
  12    real,          intent(in)    :: zh(mkz)
  13    real,          intent(in)    :: zf(mkz+1)
  14
  15    integer                      :: k
  16
  17    if (trace_use) call da_trace_entry("da_transform_xtozrhoq_adj")
  18
  19    do k = 1,mkz
  20       grid%xa%rho(i,j,k) = (zf(k)-zf(k+1))*zh(k)*grid%xb%q(i,j,k)*ADJ_zrhom   + grid%xa%rho(i,j,k)
  21       grid%xa%q(i,j,k)   = (zf(k)-zf(k+1))*zh(k)*ADJ_zrhom*grid%xb%rho(i,j,k) + grid%xa%q(i,j,k)
  22    end do
  23
  24    if (trace_use) call da_trace_exit("da_transform_xtozrhoq_adj")
  25
  26 end subroutine da_transform_xtozrhoq_adj
  27
  28