da_transform_xtotpw_adj.inc

1 subroutine da_transform_xtotpw_adj(grid)
2 
3    !---------------------------------------------------------------------
4    ! Purpose: TBD
5    !---------------------------------------------------------------------
6 
7    implicit none
8 
9    type (domain),  intent(inout) :: grid
10 
11    integer :: i, j, k, is, js, ie, je
12 
13    real    :: pw, dzpw
14 
15    if (trace_use) call da_trace_entry("da_transform_xtotpw_adj")
16 
17    is = its
18    js = jts
19 
20    ie = ite
21    je = jte
22 
23    if (test_wrfvar) then
24       is = its-1
25       js = jts-1
26 
27       ie = ite+1
28       je = jte+1
29 
30       if (is < ids) is = ids
31       if (js < jds) js = jds
32 
33       if (ie > ide) ie = ide
34       if (je > jde) je = jde
35    end if
36 
37    do j=js, je
38       do i=is, ie
39          pw = 0.1*grid%xa%tpw(i,j)
40 
41          do k=kts, kte
42             dzpw = (grid%xb%hf(i,j,k+1)-grid%xb%hf(i,j,k))*pw
43 
44             grid%xa%  q(i,j,k)=grid%xa%  q(i,j,k)+dzpw*grid%xb%rho(i,j,k)
45             grid%xa%rho(i,j,k)=grid%xa%rho(i,j,k)+dzpw*grid%xb%  q(i,j,k)
46          end do
47       end do
48    end do
49 
50    if (trace_use) call da_trace_exit("da_transform_xtotpw_adj")
51  
52 end subroutine da_transform_xtotpw_adj
53 
54