source: lmdz_wrf/trunk/WRFV3/share/module_model_constants.F @ 2661

Last change on this file since 2661 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 5.6 KB
Line 
1!WRF:MODEL_LAYER:CONSTANTS
2!
3
4 MODULE module_model_constants
5
6   !  2. Following are constants for use in defining real number bounds.
7
8   !  A really small number.
9
10   REAL    , PARAMETER :: epsilon         = 1.E-15
11
12   !  4. Following is information related to the physical constants.
13
14   !  These are the physical constants used within the model.
15
16! JM NOTE -- can we name this grav instead?
17   REAL    , PARAMETER :: g = 9.81  ! acceleration due to gravity (m {s}^-2)
18
19#if ( NMM_CORE == 1 )
20   REAL    , PARAMETER :: r_d          = 287.04
21   REAL    , PARAMETER :: cp           = 1004.6
22#else
23   REAL    , PARAMETER :: r_d          = 287.
24   REAL    , PARAMETER :: cp           = 7.*r_d/2.
25#endif
26
27   REAL    , PARAMETER :: r_v          = 461.6
28   REAL    , PARAMETER :: cv           = cp-r_d
29   REAL    , PARAMETER :: cpv          = 4.*r_v
30   REAL    , PARAMETER :: cvv          = cpv-r_v
31   REAL    , PARAMETER :: cvpm         = -cv/cp
32   REAL    , PARAMETER :: cliq         = 4190.
33   REAL    , PARAMETER :: cice         = 2106.
34   REAL    , PARAMETER :: psat         = 610.78
35   REAL    , PARAMETER :: rcv          = r_d/cv
36   REAL    , PARAMETER :: rcp          = r_d/cp
37   REAL    , PARAMETER :: rovg         = r_d/g
38   REAL    , PARAMETER :: c2           = cp * rcv
39   real    , parameter :: mwdry        = 28.966 ! molecular weight of dry air (g/mole)
40
41   REAL    , PARAMETER :: p1000mb      = 100000.
42   REAL    , PARAMETER :: t0           = 300.
43   REAL    , PARAMETER :: p0           = p1000mb
44   REAL    , PARAMETER :: cpovcv       = cp/(cp-r_d)
45   REAL    , PARAMETER :: cvovcp       = 1./cpovcv
46   REAL    , PARAMETER :: rvovrd       = r_v/r_d
47
48   REAL    , PARAMETER :: reradius     = 1./6370.0e03
49
50   REAL    , PARAMETER :: asselin      = .025
51!   REAL    , PARAMETER :: asselin      = .0
52   REAL    , PARAMETER :: cb           = 25.
53
54   REAL    , PARAMETER :: XLV0         = 3.15E6
55   REAL    , PARAMETER :: XLV1         = 2370.
56   REAL    , PARAMETER :: XLS0         = 2.905E6
57   REAL    , PARAMETER :: XLS1         = 259.532
58
59   REAL    , PARAMETER :: XLS          = 2.85E6
60   REAL    , PARAMETER :: XLV          = 2.5E6
61   REAL    , PARAMETER :: XLF          = 3.50E5
62
63   REAL    , PARAMETER :: rhowater     = 1000.
64   REAL    , PARAMETER :: rhosnow      = 100.
65   REAL    , PARAMETER :: rhoair0      = 1.28
66!
67   REAL    , PARAMETER :: n_ccn0       = 1.0E8
68!
69   REAL    , PARAMETER :: DEGRAD       = 3.1415926/180.
70   REAL    , PARAMETER :: DPD          = 360./365.
71
72   REAL    , PARAMETER ::  SVP1=0.6112
73   REAL    , PARAMETER ::  SVP2=17.67
74   REAL    , PARAMETER ::  SVP3=29.65
75   REAL    , PARAMETER ::  SVPT0=273.15
76   REAL    , PARAMETER ::  EP_1=R_v/R_d-1.
77   REAL    , PARAMETER ::  EP_2=R_d/R_v
78   REAL    , PARAMETER ::  KARMAN=0.4
79   REAL    , PARAMETER ::  EOMEG=7.2921E-5
80   REAL    , PARAMETER ::  STBOLT=5.67051E-8
81
82   REAL    , PARAMETER ::  prandtl = 1./3.0
83                                         ! constants for w-damping option
84   REAL    , PARAMETER ::  w_alpha = 0.3 ! strength m/s/s
85   REAL    , PARAMETER ::  w_beta  = 1.0 ! activation cfl number
86
87       REAL , PARAMETER ::  pq0=379.90516
88       REAL , PARAMETER ::  epsq2=0.2
89       REAL , PARAMETER ::  a2=17.2693882
90       REAL , PARAMETER ::  a3=273.16
91       REAL , PARAMETER ::  a4=35.86
92       REAL , PARAMETER ::  epsq=1.e-12
93       REAL , PARAMETER ::  p608=rvovrd-1.
94!#if ( NMM_CORE == 1 )
95       REAL , PARAMETER ::  climit=1.e-20
96       REAL , PARAMETER ::  cm1=2937.4
97       REAL , PARAMETER ::  cm2=4.9283
98       REAL , PARAMETER ::  cm3=23.5518
99!       REAL , PARAMETER ::  defc=8.0
100!       REAL , PARAMETER ::  defm=32.0
101       REAL , PARAMETER ::  defc=0.0
102       REAL , PARAMETER ::  defm=99999.0
103       REAL , PARAMETER ::  epsfc=1./1.05
104       REAL , PARAMETER ::  epswet=0.0
105       REAL , PARAMETER ::  fcdif=1./3.
106#ifdef HWRF
107       REAL , PARAMETER ::  fcm=0.0
108#else
109       REAL , PARAMETER ::  fcm=0.00003
110#endif
111       REAL , PARAMETER ::  gma=-r_d*(1.-rcp)*0.5
112       REAL , PARAMETER ::  p400=40000.0
113       REAL , PARAMETER ::  phitp=15000.0
114       REAL , PARAMETER ::  pi2=2.*3.1415926
115       REAL , PARAMETER ::  plbtm=105000.0
116       REAL , PARAMETER ::  plomd=64200.0
117       REAL , PARAMETER ::  pmdhi=35000.0
118       REAL , PARAMETER ::  q2ini=0.50
119       REAL , PARAMETER ::  rfcp=0.25/cp
120       REAL , PARAMETER ::  rhcrit_land=0.75
121       REAL , PARAMETER ::  rhcrit_sea=0.80
122       REAL , PARAMETER ::  rlag=14.8125
123       REAL , PARAMETER ::  rlx=0.90
124       REAL , PARAMETER ::  scq2=50.0
125       REAL , PARAMETER ::  slopht=0.001
126       REAL , PARAMETER ::  tlc=2.*0.703972477
127       REAL , PARAMETER ::  wa=0.15
128       REAL , PARAMETER ::  wght=0.35
129       REAL , PARAMETER ::  wpc=0.075
130       REAL , PARAMETER ::  z0land=0.10
131#ifdef HWRF
132       REAL , PARAMETER ::  z0max=0.01
133#else
134       REAL , PARAMETER ::  z0max=0.008
135#endif
136       REAL , PARAMETER ::  z0sea=0.001
137!#endif
138
139
140   !  Earth
141
142   !  The value for P2SI *must* be set to 1.0 for Earth
143   !  Although, now we may not need this declaration here (see above)
144   !REAL    , PARAMETER :: P2SI         = 1.0
145
146   !  Orbital constants:
147
148   INTEGER , PARAMETER :: PLANET_YEAR = 365
149   REAL , PARAMETER :: OBLIQUITY = 23.5
150   REAL , PARAMETER :: ECCENTRICITY = 0.014
151   REAL , PARAMETER :: SEMIMAJORAXIS = 1.0 ! In AU
152   ! Don't know the following values, so we'll fake them for now
153   REAL , PARAMETER :: zero_date = 0.0   ! Time of perihelion passage
154   !  Fraction into the year (from perhelion) of the
155   !  occurrence of the Northern Spring Equinox
156   REAL , PARAMETER :: EQUINOX_FRACTION= 0.0
157
158 CONTAINS
159   SUBROUTINE init_module_model_constants
160   END SUBROUTINE init_module_model_constants
161 END MODULE module_model_constants
Note: See TracBrowser for help on using the repository browser.