source: trunk/MESOSCALE/LMD_LES_MARS/modif_mars/module_model_constants.F @ 134

Last change on this file since 134 was 54, checked in by aslmd, 14 years ago

LMD_LES_MARS: Tests with new physics in unified frame
M 53 mesoscale/LMD_LES_MARS/modif_mars/module_model_constants.F
M 53 mesoscale/LMD_LES_MARS/modif_mars/module_first_rk_step_part1.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
Modifications made to cope with new soil scheme

M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F
Commented aeronomars in old physics so that the model can compile with lots of grid points

LMD_MM_MARS: Implement tracers with the new physics. Test with JBM runs with 5 tracers

including radiatively active dust and water ice

TODO: propagate reading "traceur.def" and few initialisations from

testphys1d.F into module_lmd_driver.F

M 53 mars/libf/phymars/dimradmars.h
M 53 mars/libf/phymars/callradite.F
Change it back to both dust and ice radiatively active

M 53 mesoscale/LMD_MM_MARS/makemeso
Allow to "fresh start" with the inclusion of "clean -a" in the script

M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
M 53 mesoscale/LMD_MM_MARS/SIMU/runmeso
Added the option mars==11 with 5 tracers: CO2 H2Ovap H2Oice DUST DUSTN
-- Note that the order matters for H2Ovap and H2Oice

LMD_MM_MARS: Create a folder to run LMD GCM with exact same physics as mesoscale

A 0 mesoscale/LMDZ.MARS/libf/phymars/physiq.F
All files in mesoscale/LMDZ.MARS/libf are links except this one
--> Because specific WRITEDIAGFI commands are set to output what is needed PREP_MARS
--> Something must be done so that qCO2, qdust, qdustN can be propagated too !!!

A 0 mesoscale/LMDZ.MARS/libf_gcm
A 0 mesoscale/LMDZ.MARS/in_lmdz_mars_newphys
Two important links that must not be broken

A 0 mesoscale/LMDZ.MARS/myGCM
This folder is for runs. Only links here -- except temp files which are not synchronized.

M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/compile
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/run.def
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/traceur.def
M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def
A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def.csttau
Those files are needed to run the GCM in order to prepare initial & boundary

conditions for the mesoscale... but basically it is the exact same GCM

--> compile is a convenient script

A 0 mesoscale/LMDZ.MARS/makegcm
A 0 mesoscale/LMDZ.MARS/makegcm_g95
A 0 mesoscale/LMDZ.MARS/create_make_gcm
This is supposed to be merged one day with files in trunk/mars/

A 0 mesoscale/LMDZ.MARS/myGCM/DEFS_JB
Necessary files to restart a GCM run corresponding to new water cycle in JB's thesis
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          = 191.  ! gas constant m2 s-2 K-1
22   REAL    , PARAMETER :: cp           = 744.5   ! r= 8.314511E+0 *1000.E+0/mugaz
23#else
24   REAL    , PARAMETER :: r_d          = 191.
25   REAL    , PARAMETER :: cp           = 744.5
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          = 3.9*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   real    , parameter :: mwdry        = 43.49 ! molecular weight of dry air (g/mole)
41
42   REAL    , PARAMETER :: p1000mb      = 610.
43   REAL    , PARAMETER :: t0           = 220. ! earth : 300 ...  remplacer par 220 ?
44   REAL    , PARAMETER :: p0           = 610.
45   REAL    , PARAMETER :: cpovcv       = cp/(cp-r_d)
46   REAL    , PARAMETER :: cvovcp       = 1./cpovcv
47   REAL    , PARAMETER :: rvovrd       = r_v/r_d
48
49   REAL    , PARAMETER :: reradius     = 1./3397200.    ! Rayon de mars (m)  ~3397200 m
50
51!-------------------------------
52!****MARS
53! here some LMD stuff
54! >> TODO: possible modification with startfi.nc values ?
55   REAL    , PARAMETER :: wdaysec       = 88800.     ! duree du sol (s)  ~88775 s
56   REAL    , PARAMETER :: wmugaz        = 43.49      ! Masse molaire de l'atm (g.mol-1) ~43.49
57   REAL    , PARAMETER :: womeg         = 7.0721E-5  ! omega (rad.s-1)
58   REAL    , PARAMETER :: wyear_day     = 669        ! Duree de l'annee (sols) ~668.6
59   REAL    , PARAMETER :: wperiheli     = 206.66     ! Dist.min. soleil-mars (Mkm) ~206.66     
60   REAL    , PARAMETER :: waphelie      = 249.22     ! Dist.max. soleil-mars (Mkm) ~249.22
61   REAL    , PARAMETER :: wperi_day     = 485.       ! Date du perihelie (sols depuis printemps)
62   REAL    , PARAMETER :: wobliquit     = 25.2       ! Obliquite de la planete (deg) ~25.2
63   REAL    , PARAMETER :: wz0           = 1.e-2      ! surface roughness (m) ~0.01
64   REAL    , PARAMETER :: wlmixmin      = 30         ! longueur de melange ~100
65   REAL    , PARAMETER :: wemin_turb    = 1.e-6      ! energie minimale ~1.e-8
66   REAL    , PARAMETER :: wemissiv      = 0.95       ! Emissivite du sol martien ~.95
67   REAL    , PARAMETER :: wemissiceN    = 0.95       ! Emissivite calotte nord
68   REAL    , PARAMETER :: wemissiceS    = 0.95       ! Emissivite calotte sud
69   REAL    , PARAMETER :: walbediceN    = 0.65        ! Albedo calotte nord !0.5
70   REAL    , PARAMETER :: walbediceS    = 0.65        ! Albedo calotte sud  !0.5
71   REAL    , PARAMETER :: wiceradiusN   = 100.e-6    ! mean scat radius of CO2 snow (north)
72   REAL    , PARAMETER :: wiceradiusS   = 100.e-6    ! mean scat radius of CO2 snow (south)
73   REAL    , PARAMETER :: wdtemisiceN   = 0.4  ! time scale for snow metamorphism (north) !2 
74   REAL    , PARAMETER :: wdtemisiceS   = 0.4  ! time scale for snow metamorphism (south) !2
75#ifdef NEWPHYS
76   REAL    , PARAMETER :: wvolcapa      = 1.e6       ! volumetric capacity of soil (new soil model)
77#endif
78!
79!-------------------------------
80
81   REAL    , PARAMETER :: asselin      = .025
82!   REAL    , PARAMETER :: asselin      = .0
83   REAL    , PARAMETER :: cb           = 25.
84
85   REAL    , PARAMETER :: XLV0         = 3.15E6
86   REAL    , PARAMETER :: XLV1         = 2370.
87   REAL    , PARAMETER :: XLS0         = 2.905E6
88   REAL    , PARAMETER :: XLS1         = 259.532
89
90   REAL    , PARAMETER :: XLS          = 2.85E6
91   REAL    , PARAMETER :: XLV          = 2.5E6
92   REAL    , PARAMETER :: XLF          = 3.50E5
93
94   REAL    , PARAMETER :: rhowater     = 1000.
95   REAL    , PARAMETER :: rhosnow      = 100.
96   REAL    , PARAMETER :: rhoair0      = 0.02   !Earth Surface density: 1.217 kg/m3
97                                                !Mars Surface density: 0.020 kg/m3
98   REAL    , PARAMETER :: DEGRAD       = 3.1415926/180.
99   REAL    , PARAMETER :: DPD          = 360./669. ! longitude solaire d'un jour?
100
101   REAL    , PARAMETER ::  SVP1=0.6112
102   REAL    , PARAMETER ::  SVP2=17.67
103   REAL    , PARAMETER ::  SVP3=29.65
104   REAL    , PARAMETER ::  SVPT0=273.15
105   REAL    , PARAMETER ::  EP_1=R_v/R_d-1.
106   REAL    , PARAMETER ::  EP_2=R_d/R_v
107   REAL    , PARAMETER ::  KARMAN=0.4  ! von karman constant
108   REAL    , PARAMETER ::  EOMEG=7.0721E-5 ! angular rotation rate (rad.s-1)
109   REAL    , PARAMETER ::  STBOLT=5.67051E-8 ! stefan-boltzmann constant
110
111
112
113
114
115   REAL    , PARAMETER ::  prandtl = 1./3.0
116                                         ! constants for w-damping option
117   REAL    , PARAMETER ::  w_alpha = 0.3 ! strength m/s/s
118   REAL    , PARAMETER ::  w_beta  = 1.0 ! activation cfl number
119
120       REAL , PARAMETER ::  pq0=379.90516
121       REAL , PARAMETER ::  epsq2=0.2
122!      REAL , PARAMETER ::  epsq2=0.02
123       REAL , PARAMETER ::  a2=17.2693882
124       REAL , PARAMETER ::  a3=273.16
125       REAL , PARAMETER ::  a4=35.86
126       REAL , PARAMETER ::  epsq=1.e-12
127       REAL , PARAMETER ::  p608=rvovrd-1.
128!#if ( NMM_CORE == 1 )
129       REAL , PARAMETER ::  climit=1.e-20
130       REAL , PARAMETER ::  cm1=2937.4
131       REAL , PARAMETER ::  cm2=4.9283
132       REAL , PARAMETER ::  cm3=23.5518
133!       REAL , PARAMETER ::  defc=8.0
134!       REAL , PARAMETER ::  defm=32.0
135       REAL , PARAMETER ::  defc=0.0
136       REAL , PARAMETER ::  defm=99999.0
137       REAL , PARAMETER ::  epsfc=1./1.05
138       REAL , PARAMETER ::  epswet=0.0
139       REAL , PARAMETER ::  fcdif=1./3.
140       REAL , PARAMETER ::  fcm=0.00003
141       REAL , PARAMETER ::  gma=-r_d*(1.-rcp)*0.5
142       REAL , PARAMETER ::  p400=40000.0
143       REAL , PARAMETER ::  phitp=15000.0
144       REAL , PARAMETER ::  pi2=2.*3.1415926
145       REAL , PARAMETER ::  plbtm=105000.0
146       REAL , PARAMETER ::  plomd=64200.0
147       REAL , PARAMETER ::  pmdhi=35000.0
148       REAL , PARAMETER ::  q2ini=0.50
149       REAL , PARAMETER ::  rfcp=0.25/cp
150       REAL , PARAMETER ::  rhcrit_land=0.75
151       REAL , PARAMETER ::  rhcrit_sea=0.80
152       REAL , PARAMETER ::  rlag=14.8125
153       REAL , PARAMETER ::  rlx=0.90
154       REAL , PARAMETER ::  scq2=50.0
155       REAL , PARAMETER ::  slopht=0.001
156       REAL , PARAMETER ::  tlc=2.*0.703972477
157       REAL , PARAMETER ::  wa=0.15
158       REAL , PARAMETER ::  wght=0.35
159       REAL , PARAMETER ::  wpc=0.075
160       REAL , PARAMETER ::  z0land=0.10
161       REAL , PARAMETER ::  z0max=0.008
162       REAL , PARAMETER ::  z0sea=0.001
163!#endif
164
165
166   !  Earth
167
168   !  The value for P2SI *must* be set to 1.0 for Earth
169   !  Although, now we may not need this declaration here (see above)
170   !REAL    , PARAMETER :: P2SI         = 1.0
171
172   !  Orbital constants:
173
174   INTEGER , PARAMETER :: PLANET_YEAR = 365
175   REAL , PARAMETER :: OBLIQUITY = 23.5
176   REAL , PARAMETER :: ECCENTRICITY = 0.014
177   REAL , PARAMETER :: SEMIMAJORAXIS = 1.0 ! In AU
178   ! Don't know the following values, so we'll fake them for now
179   REAL , PARAMETER :: zero_date = 0.0   ! Time of perihelion passage
180   !  Fraction into the year (from perhelion) of the
181   !  occurrence of the Northern Spring Equinox
182   REAL , PARAMETER :: EQUINOX_FRACTION= 0.0
183
184 CONTAINS
185   SUBROUTINE init_module_model_constants
186   END SUBROUTINE init_module_model_constants
187 END MODULE module_model_constants
Note: See TracBrowser for help on using the repository browser.