Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

module_model_constants.F_gcm @ 3547

Last change on this file since 3547 was 31, checked in by aslmd, 14 years ago
LMD_MM_MARS: teste compilation, scripts additionnels et initialisation nouvelle physique
File size: 7.9 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format