2 |
|
|
3 |
! Set thermodynamical constants for convectL |
! Set thermodynamical constants for convectL |
4 |
|
|
|
use clesphys2, only: iflag_con |
|
5 |
use SUPHEC_M |
use SUPHEC_M |
6 |
use cvthermo |
use cvthermo |
7 |
|
|
9 |
|
|
10 |
!----------------------------------------------------- |
!----------------------------------------------------- |
11 |
|
|
12 |
! original set from convect: |
cpd = RCPD |
13 |
if (iflag_con.eq.4) then |
cpv = RCPV |
14 |
cpd=1005.7 |
cl = RCW |
15 |
cpv=1870.0 |
rrv = RV |
16 |
cl=4190.0 |
rrd = RD |
17 |
rrv=461.5 |
lv0 = RLVTT |
18 |
rrd=287.04 |
g = RG |
19 |
lv0=2.501E6 |
t0 = 273.15 |
20 |
g=9.8 |
grav= g |
|
t0=273.15 |
|
|
grav=g |
|
|
endif |
|
|
|
|
|
! constants consistent with LMDZ: |
|
|
if (iflag_con.eq.3) then |
|
|
cpd = RCPD |
|
|
cpv = RCPV |
|
|
cl = RCW |
|
|
rrv = RV |
|
|
rrd = RD |
|
|
lv0 = RLVTT |
|
|
g = RG ! not used in convect3 |
|
|
t0 = 273.15 ! convect3 (RTT=273.16) |
|
|
grav= g ! implicitely or explicitely used in convect3 |
|
|
endif |
|
21 |
|
|
22 |
rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?) |
rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?) |
23 |
clmcpv=cl-cpv |
clmcpv=cl-cpv |
24 |
clmcpd=cl-cpd |
clmcpd=cl-cpd |
25 |
cpdmcp=cpd-cpv |
cpdmcp=cpd-cpv |
26 |
cpvmcpd=cpv-cpd |
cpvmcpd=cpv-cpd |
27 |
cpvmcl=cl-cpv ! for convect3 |
cpvmcl=cl-cpv |
28 |
eps=rrd/rrv |
eps=rrd/rrv |
29 |
epsi=1.0/eps |
epsi=1.0/eps |
30 |
epsim1=epsi-1.0 |
epsim1=epsi-1.0 |