da_transform_xtozrhoq.inc

1 subroutine da_transform_xtozrhoq(xb, i, j, zh, zf, 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    real,           intent(out)     :: 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")
18    
19    zrhom = 0.0
20 
21    do k = 1,mkz
22       zrhom = (zf(k)-zf(k+1))*zh(k)*(xb%q(i,j,k)*xb%rho(i,j,k))+zrhom
23    end do
24 
25    if (trace_use) call da_trace_exit("da_transform_xtozrhoq")
26  
27 end subroutine da_transform_xtozrhoq
28 
29