1 subroutine da_integrat_dz(grid)
3 !---------------------------------------------------------------------------
4 ! Non-linear PW forward operator.
5 ! ===============================
7 ! Purpose: To calculate the IWV from the model QV and PP, TT.
9 ! Method: IWV = sum {QV * RHO * dZ}
11 ! Unit: Qv (Kg/Kg), RHO(Kg/M^3), dZ(M)
18 !---------------------------------------------------------------------------
22 type (domain), intent(inout) :: grid
24 integer :: i, j, K, ij
28 if (trace_use) call da_trace_entry("da_integrat_dz")
30 ! weighted sum of vertical column
33 !$OMP PRIVATE (ij, i, j, pw)
34 do ij = 1, grid%num_tiles
40 pw = pw + (grid%xb%hf(i,j,k+1)-grid%xb%hf(i,j,k)) * grid%xb%q(i,j,k)*grid%xb%rho(i,j,k)
43 grid%xb%tpw(i,j) = 0.1*pw
50 if (trace_use) call da_trace_exit("da_integrat_dz")
52 end subroutine da_integrat_dz