wrfv2_fire/var/da/da_setup_structures/da_lcl.inc

   1 subroutine da_lcl(p0, z0, t0, q0, plcl, zlcl, tlcl, qlcl)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real, intent(in)   :: p0, z0, t0, q0
  10    real, intent(out)  :: plcl, zlcl, tlcl, qlcl
  11
  12    real   :: dp, qs, eps
  13
  14    if (trace_use) call da_trace_entry("da_lcl")
  15
  16    dp=5.0
  17    plcl=300.0
  18
  19    do
  20       tlcl=t0*((plcl/p0)**0.286)
  21
  22       call da_qfrmrh(plcl, tlcl, 100.0, qs)
  23
  24       eps=qs-q0
  25
  26       if (eps >= 0.0) then
  27          zlcl=(1004.0/gravity)*(t0-tlcl)+z0
  28          qlcl=qs
  29          return
  30       else
  31          plcl=plcl+dp
  32
  33          if (plcl >= p0) then
  34             zlcl=z0
  35             qlcl=q0
  36             plcl=p0
  37             return
  38         end if
  39       end if
  40    end do
  41
  42    if (trace_use) call da_trace_exit("da_lcl")
  43
  44 end subroutine da_lcl
  45
  46