1 | !WRF:MODEL_LAYER:CONSTANTS |
---|
2 | ! |
---|
3 | !PLANET MARS |
---|
4 | |
---|
5 | MODULE module_model_constants |
---|
6 | |
---|
7 | ! 2. Following are constants for use in defining real number bounds. |
---|
8 | |
---|
9 | ! A really small number. |
---|
10 | |
---|
11 | REAL , PARAMETER :: epsilon = 1.E-15 |
---|
12 | |
---|
13 | ! 4. Following is information related to the physical constants. |
---|
14 | |
---|
15 | ! These are the physical constants used within the model. |
---|
16 | |
---|
17 | ! JM NOTE -- can we name this grav instead? |
---|
18 | REAL , PARAMETER :: g = 3.72 ! acceleration due to gravity (m {s}^-2) |
---|
19 | |
---|
20 | #if ( NMM_CORE == 1 ) |
---|
21 | REAL , PARAMETER :: r_d = 192. ! gas constant m2 s-2 K-1 |
---|
22 | REAL , PARAMETER :: cp = 844.6 ! r= 8.314511E+0 *1000.E+0/mugaz |
---|
23 | #else |
---|
24 | REAL , PARAMETER :: r_d = 192. |
---|
25 | REAL , PARAMETER :: cp = 844.6 |
---|
26 | #endif |
---|
27 | |
---|
28 | REAL , PARAMETER :: r_v = 461.6 ! gas constant for water vapor |
---|
29 | REAL , PARAMETER :: cv = cp-r_d |
---|
30 | REAL , PARAMETER :: cpv = 4.4*r_v ! cp / R = 4.4 ?? ! rcp = r/cp ~0.256793 |
---|
31 | REAL , PARAMETER :: cvv = cpv-r_v |
---|
32 | REAL , PARAMETER :: cvpm = -cv/cp |
---|
33 | REAL , PARAMETER :: cliq = 4190. |
---|
34 | REAL , PARAMETER :: cice = 2106. |
---|
35 | REAL , PARAMETER :: psat = 610.78 |
---|
36 | REAL , PARAMETER :: rcv = r_d/cv |
---|
37 | REAL , PARAMETER :: rcp = r_d/cp |
---|
38 | REAL , PARAMETER :: rovg = r_d/g |
---|
39 | REAL , PARAMETER :: c2 = cp * rcv |
---|
40 | |
---|
41 | REAL , PARAMETER :: p1000mb = 610. |
---|
42 | REAL , PARAMETER :: t0 = 220. ! earth : 300 ... à remplacer par 220 ? |
---|
43 | REAL , PARAMETER :: p0 = 610. |
---|
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./3397200. ! Rayon de mars (m) ~3397200 m |
---|
49 | |
---|
50 | !------------------------------- |
---|
51 | !****MARS |
---|
52 | ! here some LMD stuff |
---|
53 | ! >> TODO: possible modification with startfi.nc values ? |
---|
54 | REAL , PARAMETER :: wdaysec = 88800. ! duree du sol (s) ~88775 s |
---|
55 | REAL , PARAMETER :: wmugaz = 43.49 ! Masse molaire de l'atm (g.mol-1) ~43.49 |
---|
56 | REAL , PARAMETER :: womeg = 7.0721E-5 ! omega (rad.s-1) |
---|
57 | REAL , PARAMETER :: wyear_day = 669 ! Duree de l'annee (sols) ~668.6 |
---|
58 | REAL , PARAMETER :: wperiheli = 206.66 ! Dist.min. soleil-mars (Mkm) ~206.66 |
---|
59 | REAL , PARAMETER :: waphelie = 249.22 ! Dist.max. soleil-mars (Mkm) ~249.22 |
---|
60 | REAL , PARAMETER :: wperi_day = 485. ! Date du perihelie (sols depuis printemps) |
---|
61 | REAL , PARAMETER :: wobliquit = 25.2 ! Obliquite de la planete (deg) ~25.2 |
---|
62 | REAL , PARAMETER :: wz0 = 1.e-2 ! surface roughness (m) ~0.01 |
---|
63 | REAL , PARAMETER :: wlmixmin = 30 ! longueur de melange ~100 |
---|
64 | REAL , PARAMETER :: wemin_turb = 1.e-6 ! energie minimale ~1.e-8 |
---|
65 | REAL , PARAMETER :: wemissiv = 0.95 ! Emissivite du sol martien ~.95 |
---|
66 | REAL , PARAMETER :: wemissiceN = 0.95 ! Emissivite calotte nord |
---|
67 | REAL , PARAMETER :: wemissiceS = 0.95 ! Emissivite calotte sud |
---|
68 | REAL , PARAMETER :: walbediceN = 0.65 ! Albedo calotte nord !0.5 |
---|
69 | REAL , PARAMETER :: walbediceS = 0.65 ! Albedo calotte sud !0.5 |
---|
70 | REAL , PARAMETER :: wiceradiusN = 100.e-6 ! mean scat radius of CO2 snow (north) |
---|
71 | REAL , PARAMETER :: wiceradiusS = 100.e-6 ! mean scat radius of CO2 snow (south) |
---|
72 | REAL , PARAMETER :: wdtemisiceN = 0.4 ! time scale for snow metamorphism (north) !2 |
---|
73 | REAL , PARAMETER :: wdtemisiceS = 0.4 ! time scale for snow metamorphism (south) !2 |
---|
74 | #ifdef NEWPHYS |
---|
75 | REAL , PARAMETER :: wvolcapa = 1.e6 ! volumetric capacity of soil (new soil model) |
---|
76 | #endif |
---|
77 | ! |
---|
78 | !------------------------------- |
---|
79 | |
---|
80 | ! REAL , PARAMETER :: asselin = .025 |
---|
81 | REAL , PARAMETER :: asselin = .0 |
---|
82 | REAL , PARAMETER :: cb = 25. |
---|
83 | |
---|
84 | REAL , PARAMETER :: XLV0 = 3.15E6 |
---|
85 | REAL , PARAMETER :: XLV1 = 2370. |
---|
86 | REAL , PARAMETER :: XLS0 = 2.905E6 |
---|
87 | REAL , PARAMETER :: XLS1 = 259.532 |
---|
88 | |
---|
89 | REAL , PARAMETER :: XLS = 2.85E6 |
---|
90 | REAL , PARAMETER :: XLV = 2.5E6 |
---|
91 | REAL , PARAMETER :: XLF = 3.50E5 |
---|
92 | |
---|
93 | REAL , PARAMETER :: rhowater = 1000. |
---|
94 | REAL , PARAMETER :: rhosnow = 100. |
---|
95 | REAL , PARAMETER :: rhoair0 = 0.02 !Earth Surface density: 1.217 kg/m3 |
---|
96 | !Mars Surface density: 0.020 kg/m3 |
---|
97 | REAL , PARAMETER :: DEGRAD = 3.1415926/180. |
---|
98 | REAL , PARAMETER :: DPD = 360./669. ! longitude solaire d'un jour? |
---|
99 | |
---|
100 | REAL , PARAMETER :: SVP1=0.6112 |
---|
101 | REAL , PARAMETER :: SVP2=17.67 |
---|
102 | REAL , PARAMETER :: SVP3=29.65 |
---|
103 | REAL , PARAMETER :: SVPT0=273.15 |
---|
104 | REAL , PARAMETER :: EP_1=R_v/R_d-1. |
---|
105 | REAL , PARAMETER :: EP_2=R_d/R_v |
---|
106 | REAL , PARAMETER :: KARMAN=0.4 ! von karman constant |
---|
107 | REAL , PARAMETER :: EOMEG=7.0721E-5 ! angular rotation rate (rad.s-1) |
---|
108 | REAL , PARAMETER :: STBOLT=5.67051E-8 ! stefan-boltzmann constant |
---|
109 | |
---|
110 | ! proportionality constants for eddy viscosity coefficient calc |
---|
111 | REAL , PARAMETER :: c_s = .25 ! turbulence parameterization constant, for smagorinsky |
---|
112 | ! REAL , PARAMETER :: c_s = .1 ! turbulence parameterization constant, for smagorinsky |
---|
113 | REAL , PARAMETER :: c_k = .15 ! turbulence parameterization constant, for TKE |
---|
114 | REAL , PARAMETER :: prandtl = 1./3.0 |
---|
115 | ! constants for w-damping option |
---|
116 | ! REAL , PARAMETER :: w_alpha = 0.3 ! strength m/s/s |
---|
117 | REAL , PARAMETER :: w_alpha = 2.0 ! strength m/s/s |
---|
118 | ! REAL , PARAMETER :: w_beta = 1.0 ! activation cfl number |
---|
119 | REAL , PARAMETER :: w_beta = 0.8 ! activation cfl number |
---|
120 | |
---|
121 | REAL , PARAMETER :: pq0=379.90516 |
---|
122 | REAL , PARAMETER :: epsq2=0.2 |
---|
123 | !try both ? |
---|
124 | ! REAL , PARAMETER :: epsq2=0.02 |
---|
125 | REAL , PARAMETER :: a2=17.2693882 |
---|
126 | REAL , PARAMETER :: a3=273.16 |
---|
127 | REAL , PARAMETER :: a4=35.86 |
---|
128 | REAL , PARAMETER :: epsq=1.e-12 |
---|
129 | REAL , PARAMETER :: p608=rvovrd-1. |
---|
130 | !#if ( NMM_CORE == 1 ) |
---|
131 | REAL , PARAMETER :: climit=1.e-20 |
---|
132 | REAL , PARAMETER :: cm1=2937.4 |
---|
133 | REAL , PARAMETER :: cm2=4.9283 |
---|
134 | REAL , PARAMETER :: cm3=23.5518 |
---|
135 | ! REAL , PARAMETER :: defc=8.0 |
---|
136 | ! REAL , PARAMETER :: defm=32.0 |
---|
137 | REAL , PARAMETER :: defc=0.0 |
---|
138 | REAL , PARAMETER :: defm=99999.0 |
---|
139 | REAL , PARAMETER :: epsfc=1./1.05 |
---|
140 | REAL , PARAMETER :: epswet=0.0 |
---|
141 | REAL , PARAMETER :: fcdif=1./3. |
---|
142 | ! REAL , PARAMETER :: fcm=0.003 |
---|
143 | REAL , PARAMETER :: fcm=0.00003 |
---|
144 | ! REAL , PARAMETER :: fcm=0.0 |
---|
145 | REAL , PARAMETER :: gma=-r_d*(1.-rcp)*0.5 |
---|
146 | REAL , PARAMETER :: p400=40000.0 |
---|
147 | REAL , PARAMETER :: phitp=15000.0 |
---|
148 | REAL , PARAMETER :: pi2=2.*3.1415926 |
---|
149 | REAL , PARAMETER :: plbtm=105000.0 |
---|
150 | REAL , PARAMETER :: plomd=64200.0 |
---|
151 | REAL , PARAMETER :: pmdhi=35000.0 |
---|
152 | REAL , PARAMETER :: q2ini=0.50 |
---|
153 | REAL , PARAMETER :: rfcp=0.25/cp |
---|
154 | REAL , PARAMETER :: rhcrit_land=0.75 |
---|
155 | REAL , PARAMETER :: rhcrit_sea=0.80 |
---|
156 | REAL , PARAMETER :: rlag=14.8125 |
---|
157 | REAL , PARAMETER :: rlx=0.90 |
---|
158 | REAL , PARAMETER :: scq2=50.0 |
---|
159 | REAL , PARAMETER :: slopht=0.001 |
---|
160 | REAL , PARAMETER :: tlc=2.*0.703972477 |
---|
161 | REAL , PARAMETER :: wa=0.15 |
---|
162 | REAL , PARAMETER :: wght=0.35 |
---|
163 | REAL , PARAMETER :: wpc=0.075 |
---|
164 | REAL , PARAMETER :: z0land=0.10 |
---|
165 | ! REAL , PARAMETER :: z0max=0.01 |
---|
166 | REAL , PARAMETER :: z0max=0.008 |
---|
167 | REAL , PARAMETER :: z0sea=0.001 |
---|
168 | !#endif |
---|
169 | |
---|
170 | |
---|
171 | CONTAINS |
---|
172 | SUBROUTINE init_module_model_constants |
---|
173 | END SUBROUTINE init_module_model_constants |
---|
174 | END MODULE module_model_constants |
---|