da_transform_xtotb.inc

1 subroutine da_transform_xtotb(grid)
2 
3    !----------------------------------------------------------------------
4    ! Purpose: TBD
5    !----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (domain), intent(inout)   :: grid
10 
11    integer            :: i,j,k
12 
13    real               :: psfc,ta,gamma,sst,htpw,speed,alw,zcld,tpw,dum1,zrhom
14 
15    if (trace_use) call da_trace_entry("da_transform_xtotb")
16 
17    do j=jts,jte
18       do i=its,ite
19          ! surface pressure (mb) (940 -1030)
20 
21          psfc          = 0.01*grid%xb%psfc(i,j)
22 
23          ! sea surface temperature (k) (273 - 303) (doesnot change) 
24 
25          sst           = grid%xb%tgrn(i,j)
26 
27          ! effective surface air temperature (263 - 303)
28 
29          ta          = grid%xb%tgrn(i,j) + &
30                        (grid%xb%t(i,j,kts)-grid%xb%tgrn(i,j))*log(2.0/0.0001)/ &
31                        log((grid%xb%h(i,j,kts) - grid%xb%terr(i,j))/0.0001)
32 
33          ! gamma is an emperical formula and zcld is given for now
34 
35          gamma   = (ta-270)*0.023 + 5.03  ! effective lapse rate(km) (4.0-6.5)
36 
37          zcld    = 1                           ! effective cloud height (km)
38                                                ! = 1 if no cloud infomation
39          ! total precipitable water in cm
40          ! total precipitable water in (kg/m**2) (0 - 70)
41 
42          tpw          = grid%xb%tpw(i,j)*10.0
43          speed        = grid%xb%speed(i,j)
44 
45          ! Column liquid water (kg/m**2)  (0-0.5) (no data now. So do it later)
46 
47          alw          = 0.0
48 
49          ! Column height weighted moisture density on the grid locally 
50 
51          zrhom = 0.0
52          do k=kts,kte
53             zrhom=zrhom+(grid%xb%hf(i,j,k+1)-grid%xb%hf(i,j,k))*grid%xb%h(i,j,k)*grid%xb%q(i,j,k)* &
54                grid%xb%rho(i,j,k)
55          end do
56 
57          ! Column moisture density on the grid locally
58 
59          htpw          = zrhom/tpw/1000.0
60 
61          call tb(1,53.0,psfc,ta,gamma,sst,tpw,htpw, &
62             speed,alw,zcld,grid%xb%tb19v(i,j),grid%xb%tb19h(i,j))
63          call tb(2,53.0,psfc,ta,gamma,sst,tpw,htpw, &
64             speed,alw,zcld,grid%xb%tb22v(i,j),dum1)
65          call tb(3,53.0,psfc,ta,gamma,sst,tpw,htpw, &
66             speed,alw,zcld,grid%xb%tb37v(i,j),grid%xb%tb37h(i,j))
67          call tb(4,53.0,psfc,ta,gamma,sst,tpw,htpw, &
68             speed,alw,zcld,grid%xb%tb85v(i,j),grid%xb%tb85h(i,j))
69       end do
70    end do
71 
72    if (trace_use) call da_trace_exit("da_transform_xtotb")
73 
74 end subroutine da_transform_xtotb
75 
76