Parent Directory | Revision Log
LMDZE now uses library Jumble. Removed all calls to "flinget". Replaced calls to "flinget", "flininfo", "flinopen_nozoom" by calls to NetCDF95 and Jumble. Split file "cv_driver.f" into "cv_driver.f90", "cv_flag.f90" and "cv_thermo.f90". Bug fix: "QANCIEN" was read twice in "phyeytat0". In "physiq", initialization of "d_t", "d_u", "d_v" was useless.
1 | guez | 49 | SUBROUTINE cv_thermo(iflag_con) |
2 | |||
3 | ! Set thermodynamical constants for convectL | ||
4 | |||
5 | use SUPHEC_M | ||
6 | use cvthermo | ||
7 | |||
8 | implicit none | ||
9 | |||
10 | integer, intent(in):: iflag_con | ||
11 | |||
12 | !----------------------------------------------------- | ||
13 | |||
14 | ! original set from convect: | ||
15 | if (iflag_con.eq.4) then | ||
16 | cpd=1005.7 | ||
17 | cpv=1870.0 | ||
18 | cl=4190.0 | ||
19 | rrv=461.5 | ||
20 | rrd=287.04 | ||
21 | lv0=2.501E6 | ||
22 | g=9.8 | ||
23 | t0=273.15 | ||
24 | grav=g | ||
25 | endif | ||
26 | |||
27 | ! constants consistent with LMDZ: | ||
28 | if (iflag_con.eq.3) then | ||
29 | cpd = RCPD | ||
30 | cpv = RCPV | ||
31 | cl = RCW | ||
32 | rrv = RV | ||
33 | rrd = RD | ||
34 | lv0 = RLVTT | ||
35 | g = RG ! not used in convect3 | ||
36 | t0 = 273.15 ! convect3 (RTT=273.16) | ||
37 | grav= g ! implicitely or explicitely used in convect3 | ||
38 | endif | ||
39 | |||
40 | rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?) | ||
41 | clmcpv=cl-cpv | ||
42 | clmcpd=cl-cpd | ||
43 | cpdmcp=cpd-cpv | ||
44 | cpvmcpd=cpv-cpd | ||
45 | cpvmcl=cl-cpv ! for convect3 | ||
46 | eps=rrd/rrv | ||
47 | epsi=1.0/eps | ||
48 | epsim1=epsi-1.0 | ||
49 | ginv=1.0/grav | ||
50 | hrd=0.5*rrd | ||
51 | |||
52 | end SUBROUTINE cv_thermo |
ViewVC Help | |
Powered by ViewVC 1.1.21 |