17 |
|
|
18 |
contains |
contains |
19 |
|
|
20 |
elemental REAL function FOEEW(T, DEL) |
elemental REAL function FOEEW(T, ICE) |
21 |
|
|
22 |
use yoethf_m, only: R3LES, R3IES, R4LES, R4IES |
use yoethf_m, only: R3LES, R3IES, R4LES, R4IES |
23 |
use SUPHEC_M, only: rtt |
use SUPHEC_M, only: rtt |
24 |
|
|
25 |
REAL, intent(in):: T |
REAL, intent(in):: T |
26 |
logical, intent(in):: DEL ! ice, else liquid |
logical, intent(in):: ICE ! else liquid |
27 |
|
|
28 |
!----------------------- |
!----------------------- |
29 |
|
|
30 |
FOEEW = exp(merge(R3IES / (T - R4IES), R3lES / (T - R4lES), del) & |
FOEEW = exp(merge(R3IES / (T - R4IES), R3lES / (T - R4lES), ice) & |
31 |
* (T - RTT)) |
* (T - RTT)) |
32 |
|
|
33 |
end function FOEEW |
end function FOEEW |
34 |
|
|
35 |
!****************************************** |
!****************************************** |
36 |
|
|
37 |
REAL function FOEDE(T, DEL, P5ARG, QS, PCOARG) |
REAL function FOEDE(T, ICE, P5ARG, QS, PCOARG) |
38 |
|
|
39 |
use yoethf_m, only: R4LES, R4IES |
use yoethf_m, only: R4LES, R4IES |
40 |
|
|
41 |
REAL, intent(in):: T |
REAL, intent(in):: T |
42 |
logical, intent(in):: DEL ! ice, else liquid |
logical, intent(in):: ICE ! else liquid |
43 |
real, intent(in):: P5ARG, QS, PCOARG |
real, intent(in):: P5ARG, QS, PCOARG |
44 |
|
|
45 |
!----------------------- |
!----------------------- |
46 |
|
|
47 |
FOEDE = QS * PCOARG * P5ARG / (T - merge(R4IES, R4lES, del))**2 |
FOEDE = QS * PCOARG * P5ARG / (T - merge(R4IES, R4lES, ice))**2 |
48 |
|
|
49 |
end function FOEDE |
end function FOEDE |
50 |
|
|
100 |
|
|
101 |
!----------------------- |
!----------------------- |
102 |
|
|
103 |
dqsatl = RLVTT / RCPD * qs * LOG(10.) & |
dqsatl = RLVTT / RCPD * qs * LOG(10.) * (2948.964 / t**2 - 5.028 & |
104 |
* (2948.964/t**2-5.028/LOG(10.)/t & |
/ LOG(10.) / t + 25.21935 * 2999.924 / t**2 * EXP(- 2999.924 / t) & |
105 |
+25.21935*2999.924/t**2*EXP(-2999.924/t) & |
+ 29810.16 * 0.0699382 * EXP(- 0.0699382 * t)) |
|
+29810.16*0.0699382*EXP(-0.0699382*t)) |
|
106 |
|
|
107 |
end function dqsatl |
end function dqsatl |
108 |
|
|