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