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 |
|
|