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