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    A=17.26
15    B=35.86
16    if (T <= 263.) A=21.87
17    if (T <= 263.) B= 7.66
18    E=6.11*EXP(A*(T-t_triple)/(T-B))
19    QS=0.622*E/(P-0.378*E)
20    Q=QS*RH/100.
21 
22 end subroutine da_qfrmrh
23 
24