da_wdt.inc

1 subroutine da_wdt(H,W,TERR,DT)
2 
3    !----------------------------------------------------------------------
4    ! Purpose: Calculate DT
5    !----------------------------------------------------------------------
6 
7    implicit none
8 
9    real, dimension(kts:kte),intent(in)   :: h
10    real, dimension(kts:kte),intent(out)  :: dt
11    real, dimension(kts:kte+1),intent(in) :: w
12    real,intent(in)                       :: terr
13 
14    integer                     :: k
15 
16    do k=kte,kts+1,-1
17       if (w(k) >= 0.1) then
18          dt(k)=(h(k)-h(k-1))/w(k)
19       else
20          dt(k)=0.
21       end if
22    end do
23 
24    if (w(kts) >= 0.1) then
25       dt(kts)=(h(kts)-terr)/w(kts)
26    else
27       dt(kts)=0.0
28    end if
29 
30 end subroutine da_wdt
31 
32