1 |
SUBROUTINE cv_thermo |
2 |
|
3 |
! Set thermodynamical constants for convectL |
4 |
|
5 |
use clesphys2, only: iflag_con |
6 |
use SUPHEC_M |
7 |
use cvthermo |
8 |
|
9 |
implicit none |
10 |
|
11 |
!----------------------------------------------------- |
12 |
|
13 |
! original set from convect: |
14 |
if (iflag_con.eq.4) then |
15 |
cpd=1005.7 |
16 |
cpv=1870.0 |
17 |
cl=4190.0 |
18 |
rrv=461.5 |
19 |
rrd=287.04 |
20 |
lv0=2.501E6 |
21 |
g=9.8 |
22 |
t0=273.15 |
23 |
grav=g |
24 |
endif |
25 |
|
26 |
! constants consistent with LMDZ: |
27 |
if (iflag_con.eq.3) then |
28 |
cpd = RCPD |
29 |
cpv = RCPV |
30 |
cl = RCW |
31 |
rrv = RV |
32 |
rrd = RD |
33 |
lv0 = RLVTT |
34 |
g = RG ! not used in convect3 |
35 |
t0 = 273.15 ! convect3 (RTT=273.16) |
36 |
grav= g ! implicitely or explicitely used in convect3 |
37 |
endif |
38 |
|
39 |
rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?) |
40 |
clmcpv=cl-cpv |
41 |
clmcpd=cl-cpd |
42 |
cpdmcp=cpd-cpv |
43 |
cpvmcpd=cpv-cpd |
44 |
cpvmcl=cl-cpv ! for convect3 |
45 |
eps=rrd/rrv |
46 |
epsi=1.0/eps |
47 |
epsim1=epsi-1.0 |
48 |
ginv=1.0/grav |
49 |
hrd=0.5*rrd |
50 |
|
51 |
end SUBROUTINE cv_thermo |