--- trunk/Sources/phylmd/cv_thermo.f 2016/03/16 11:11:27 182 +++ trunk/Sources/phylmd/cv_thermo.f 2016/06/06 17:42:15 201 @@ -1,34 +1,19 @@ -SUBROUTINE cv_thermo +module cv_thermo_m - ! Set thermodynamical constants for convectL + ! From LMDZ4/libf/phylmd/cvthermo.h, version 1.1.1.1 2004/05/19 12:53:09 + ! Thermodynamical constants for cv_driver - use SUPHEC_M - use cvthermo + use SUPHEC_M, only: rd, rg, rcpd, rcpv, rcw, rv implicit none - !----------------------------------------------------- + real, parameter:: rowl = 1000., t0 = 273.15 + real, parameter:: clmcpv = rcw - rcpv, clmcpd = rcw - rcpd + real, parameter:: cpdmcp = rcpd - rcpv + real, parameter:: cpvmcpd = rcpv - rcpd, cpvmcl = rcw - rcpv + real, parameter:: eps = rd / rv, epsi = 1. / eps, epsim1 = epsi - 1. + real, parameter:: ginv = 1. / rg, hrd = 0.5 * rd - cpd = RCPD - cpv = RCPV - cl = RCW - rrv = RV - rrd = RD - lv0 = RLVTT - g = RG - t0 = 273.15 - grav= g - - rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?) - clmcpv=cl-cpv - clmcpd=cl-cpd - cpdmcp=cpd-cpv - 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 + private rd, rg, rcpd, rcpv, rcw, rv -end SUBROUTINE cv_thermo +end module cv_thermo_m