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