/[lmdze]/trunk/phylmd/cv_thermo.f
ViewVC logotype

Diff of /trunk/phylmd/cv_thermo.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 196 by guez, Mon May 23 13:50:39 2016 UTC revision 198 by guez, Tue May 31 16:17:35 2016 UTC
# Line 3  module cv_thermo_m Line 3  module cv_thermo_m
3    ! From LMDZ4/libf/phylmd/cvthermo.h, version 1.1.1.1 2004/05/19 12:53:09    ! From LMDZ4/libf/phylmd/cvthermo.h, version 1.1.1.1 2004/05/19 12:53:09
4    ! Thermodynamical constants for cv_driver    ! Thermodynamical constants for cv_driver
5    
6      use SUPHEC_M, only: rd, rg, rcpd, rcpv, rcw
7    
8    implicit none    implicit none
9    
10    real cpd, cpv, cl, rrv, rrd, lv0, g, rowl, t0    real cpd, cpv, cl, rrv, rrd, rowl, t0
11    real clmcpv, clmcpd, cpdmcp, cpvmcpd, cpvmcl    real, parameter:: clmcpv = rcw - rcpv, clmcpd = rcw - rcpd
12      real, parameter:: cpdmcp = rcpd - rcpv
13      real, parameter:: cpvmcpd = rcpv - rcpd, cpvmcl = rcw - rcpv
14    real eps, epsi, epsim1    real eps, epsi, epsim1
15    real ginv, hrd    real, parameter:: ginv = 1. / rg, hrd = 0.5 * rd
16    real grav  
17      private rd, rg, rcpd, rcpv, rcw
18    
19  contains  contains
20    
# Line 17  contains Line 22  contains
22    
23      ! Set thermodynamical constants for cv_driver      ! Set thermodynamical constants for cv_driver
24    
25      use SUPHEC_M, only: rcpd, rcpv, rcw, rd, rg, rlvtt, rv      use SUPHEC_M, only: rlvtt, rv
26    
27      !-----------------------------------------------------      !-----------------------------------------------------
28    
# Line 26  contains Line 31  contains
31      cl  = RCW      cl  = RCW
32      rrv = RV      rrv = RV
33      rrd = RD      rrd = RD
     lv0 = RLVTT  
     g   = RG  
34      t0  = 273.15      t0  = 273.15
     grav= g  
35    
36      rowl=1000. ! (\`A quelle variable de SUPHEC_M cela correspond-il ?)      rowl = 1000. ! (\`A quelle variable de SUPHEC_M cela correspond-il ?)
37      clmcpv=cl-cpv      eps = rd/rrv
38      clmcpd=cl-cpd      epsi = 1.0/eps
39      cpdmcp=cpd-cpv      epsim1 = epsi - 1.0
     cpvmcpd=cpv-cpd  
     cpvmcl=cl-cpv  
     eps=rrd/rrv  
     epsi=1.0/eps  
     epsim1=epsi-1.0  
     ginv=1.0/grav  
     hrd=0.5*rrd  
40    
41    end SUBROUTINE cv_thermo    end SUBROUTINE cv_thermo
42    
Legend:
Removed from v.196  
changed lines
  Added in v.198
  ViewVC Help
Powered by ViewVC 1.1.21