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