var/da/da_setup_structures/da_qfrmrh.inc

   1 subroutine da_qfrmrh ( p, t, rh, q )
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real, intent(in)  :: p, t, rh
  10    real, intent(out) :: q
  11
  12    real   :: a, b, e, qs
  13
  14    if (trace_use) call da_trace_entry("da_qfrmrh")
  15
  16    a=17.26
  17    b=35.86
  18    if (t <= 263.0) a=21.87
  19    if (t <= 263.0) b= 7.66
  20    e  = 6.11*exp(a*(t-t_triple)/(t-b))
  21    qs = 0.622*e/(p-0.378*e)
  22    q  = qs*rh/100.0
  23
  24    if (trace_use) call da_trace_exit("da_qfrmrh")
  25
  26 end subroutine da_qfrmrh
  27
  28