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, intent(in)  :: h(kts:kte)
10    real, intent(out) :: dt(kts:kte)
11    real, intent(in)  :: w(kts:kte+1)
12    real, intent(in)  :: terr
13 
14    integer :: k
15 
16    if (trace_use) call da_trace_entry("da_wdt")
17 
18    do k=kte,kts+1,-1
19       if (w(k) >= 0.1) then
20          dt(k)=(h(k)-h(k-1))/w(k)
21       else
22          dt(k)=0.0
23       end if
24    end do
25 
26    if (w(kts) >= 0.1) then
27       dt(kts)=(h(kts)-terr)/w(kts)
28    else
29       dt(kts)=0.0
30    end if
31 
32    if (trace_use) call da_trace_exit("da_wdt")
33 
34 end subroutine da_wdt
35 
36