/[lmdze]/trunk/libf/phylmd/cv_thermo.f90
ViewVC logotype

Contents of /trunk/libf/phylmd/cv_thermo.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 49 - (show annotations)
Wed Aug 24 11:43:14 2011 UTC (12 years, 8 months ago) by guez
File size: 999 byte(s)
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