var/da/da_physics/da_transform_xtotpw.inc

   1 subroutine da_transform_xtotpw(grid)
   2
   3    !---------------------------------------------------------------------
   4    ! Purpose: weighted sum of vertical column
   5    !---------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (domain),  intent(inout) :: grid
  10
  11    integer :: i, j, k
  12
  13    real    :: pw
  14
  15    if (trace_use) call da_trace_entry("da_transform_xtotpw")
  16
  17    do j=jts, jte
  18       do i=its, ite
  19          pw = 0.0
  20          do k=kts, kte
  21             pw = pw + (grid%xb%hf(i,j,k+1)-grid%xb%hf(i,j,k)) &
  22                     * (grid%xa%q(i,j,k)*grid%xb%rho(i,j,k) &
  23                     +  grid%xb%q(i,j,k)*grid%xa%rho(i,j,k))
  24          end do
  25
  26          ! To convert the unit of PW to cm:
  27          grid%xa%tpw(i,j) = 0.1*pw
  28       end do
  29    end do
  30
  31    if (trace_use) call da_trace_exit("da_transform_xtotpw")
  32
  33 end subroutine da_transform_xtotpw
  34
  35