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