Sources/phylmd/cv_thermo.f

SUBROUTINE cv_thermo

  ! Set thermodynamical constants for convectL

  use clesphys2, only: iflag_con
  use SUPHEC_M
  use cvthermo

  implicit none

  !-----------------------------------------------------

  ! original set from convect:
  if (iflag_con.eq.4) then
     cpd=1005.7
     cpv=1870.0
     cl=4190.0
     rrv=461.5
     rrd=287.04
     lv0=2.501E6
     g=9.8
     t0=273.15
     grav=g
  endif

  ! constants consistent with LMDZ:
  if (iflag_con.eq.3) then
     cpd = RCPD
     cpv = RCPV
     cl  = RCW
     rrv = RV
     rrd = RD
     lv0 = RLVTT
     g   = RG     ! not used in convect3
     t0  = 273.15 ! convect3 (RTT=273.16)
     grav= g    ! implicitely or explicitely used in convect3
  endif

  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 ! for convect3
  eps=rrd/rrv
  epsi=1.0/eps
  epsim1=epsi-1.0
  ginv=1.0/grav
  hrd=0.5*rrd

end SUBROUTINE cv_thermo
1	SUBROUTINE cv_thermo
2
3	! Set thermodynamical constants for convectL
4
5	use clesphys2, only: iflag_con
6	use SUPHEC_M
7	use cvthermo
8
9	implicit none
10
11	!-----------------------------------------------------
12
13	! original set from convect:
14	if (iflag_con.eq.4) then
15	cpd=1005.7
16	cpv=1870.0
17	cl=4190.0
18	rrv=461.5
19	rrd=287.04
20	lv0=2.501E6
21	g=9.8
22	t0=273.15
23	grav=g
24	endif
25
26	! constants consistent with LMDZ:
27	if (iflag_con.eq.3) then
28	cpd = RCPD
29	cpv = RCPV
30	cl = RCW
31	rrv = RV
32	rrd = RD
33	lv0 = RLVTT
34	g = RG ! not used in convect3
35	t0 = 273.15 ! convect3 (RTT=273.16)
36	grav= g ! implicitely or explicitely used in convect3
37	endif
38
39	rowl=1000.0 !(a quelle variable de SUPHEC_M cela correspond-il?)
40	clmcpv=cl-cpv
41	clmcpd=cl-cpd
42	cpdmcp=cpd-cpv
43	cpvmcpd=cpv-cpd
44	cpvmcl=cl-cpv ! for convect3
45	eps=rrd/rrv
46	epsi=1.0/eps
47	epsim1=epsi-1.0
48	ginv=1.0/grav
49	hrd=0.5*rrd
50
51	end SUBROUTINE cv_thermo