da_integrat_dz.inc

1 subroutine da_integrat_dz(xb)
2 
3    !---------------------------------------------------------------------------
4    ! Non-linear PW forward operator.
5    ! ===============================
6    !
7    ! Purpose: To calculate the IWV from the model QV and PP, TT.
8    !
9    ! Method:  IWV = sum {QV * RHO * dZ}
10    !
11    !           Unit: Qv (Kg/Kg), RHO(Kg/M^3), dZ(M)
12    !                 PW (cm)
13    !
14    ! input     : QV, PP, TT
15    !
16    ! output    : PW
17    !
18    !---------------------------------------------------------------------------
19 
20    implicit none
21 
22    type (xb_type), intent(inout) :: xb    ! first guess (local).
23 
24    integer :: i, j, K 
25 
26    real    :: pw
27 
28    ! WEIGHTED sum OF VERTICAL COLUMN 
29 
30    do j=xb%jts, xb%jte
31       do i=xb%its, xb%ite
32          pw = 0.0
33          do k=xb%kts, xb%kte
34             pw = pw + (xb%hf(i,j,k+1)-xb%hf(i,j,k)) &
35                  *  xb%q(i,j,k)*xb%rho(i,j,k)
36          end do
37 
38          xb%tpw(i,j) = 0.1*pw
39       end do
40    end do
41 
42 end subroutine da_integrat_dz
43 
44