wrfv2_fire/var/da/da_physics/da_integrat_dz.inc

   1 subroutine da_integrat_dz(grid)
   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 (domain), intent(inout) :: grid
  23
  24    integer :: i, j, K, ij
  25
  26    real    :: pw
  27
  28    if (trace_use) call da_trace_entry("da_integrat_dz")
  29
  30    ! weighted sum of vertical column
  31
  32    !$OMP PARALLEL DO &
  33    !$OMP PRIVATE (ij, i, j, pw)
  34    do ij = 1, grid%num_tiles
  35
  36    do j=jts, jte
  37       do i=its, ite
  38          pw = 0.0
  39          do k=kts, kte
  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)
  41          end do
  42
  43          grid%xb%tpw(i,j) = 0.1*pw
  44       end do
  45    end do
  46
  47    end do
  48    !$OMP END PARALLEL DO
  49
  50    if (trace_use) call da_trace_exit("da_integrat_dz")
  51
  52 end subroutine da_integrat_dz
  53
  54