var/da/da_ssmi/da_transform_xtozrhoq_lin.inc

   1 subroutine da_transform_xtozrhoq_lin(xb, xa, i, j, zh, zf, tgl_zrhom)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    integer,        intent(in)      :: i, j
  10    type (xb_type), intent(in)      :: xb         ! first guess state.
  11    type (x_type) , intent(in)      :: xa         ! increment
  12    real,           intent(out)     :: TGL_zrhom  ! gridded height weighted moisture
  13    real,           intent(in)      :: zh(mkz)
  14    real,           intent(in)      :: zf(mkz+1)
  15
  16    integer                         :: k
  17
  18    if (trace_use) call da_trace_entry("da_transform_xtozrhoq_lin")
  19
  20    TGL_zrhom = 0.0
  21
  22    do k = 1,mkz
  23       TGL_zrhom = (zf(k)-zf(k+1))*zh(k)*(xb%q(i,j,k)*xa%rho(i,j,k) + &
  24                                          xa%q(i,j,k)*xb%rho(i,j,k)   &
  25                                        ) + TGL_zrhom
  26    end do
  27
  28    if (trace_use) call da_trace_exit("da_transform_xtozrhoq_lin")
  29
  30 end subroutine da_transform_xtozrhoq_lin
  31
  32