/[lmdze]/trunk/Sources/phylmd/FCTTRE.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/FCTTRE.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/FCTTRE.f revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC trunk/Sources/phylmd/FCTTRE.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC
# Line 17  module FCTTRE Line 17  module FCTTRE
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    
# Line 100  contains Line 100  contains
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    
Legend:
Removed from v.103  
changed lines
  Added in v.134
  ViewVC Help
Powered by ViewVC 1.1.21