var/da/da_physics/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