60 real,
parameter :: rd = 287.0
61 real,
parameter :: rv = 461.0
62 real,
parameter :: rddrv = rd/rv
67 if ( t.ge.253.15 )
then
75 if (p-(1.0-rddrv)*fespt == 0.)
then
78 f_qvsat = rddrv * fespt / (p-(1.0-rddrv)*fespt)
102 real,
parameter :: satfwa = 1.0007
103 real,
parameter :: satfwb = 3.46e-8
105 real,
parameter :: satewa = 611.21
106 real,
parameter :: satewb = 17.502
107 real,
parameter :: satewc = 32.18
109 real,
parameter :: satfia = 1.0003
110 real,
parameter :: satfib = 4.18e-8
112 real,
parameter :: sateia = 611.15
113 real,
parameter :: sateib = 22.452
114 real,
parameter :: sateic = 0.6
116 f = satfwa + satfwb * p
117 f_esl = f * satewa * exp( satewb*(t-273.15)/(t-satewc) )
139 real,
parameter :: satfwa = 1.0007
140 real,
parameter :: satfwb = 3.46e-8
142 real,
parameter :: satewa = 611.21
143 real,
parameter :: satewb = 17.502
144 real,
parameter :: satewc = 32.18
146 real,
parameter :: satfia = 1.0003
147 real,
parameter :: satfib = 4.18e-8
149 real,
parameter :: sateia = 611.15
150 real,
parameter :: sateib = 22.452
151 real,
parameter :: sateic = 0.6
153 f = satfia + satfib * p
154 f_esi = f * sateia * exp( sateib*(t-273.15)/(t-sateic) )
real function f_esi(p, t)
real function f_esl(p, t)
real function f_qvsat(p, t)