source: trunk/MESOSCALE/LMDZ.MARS.new/in_lmdz_mars_newphys/myGCM/callphys.def.csttau @ 214

Last change on this file since 214 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: 3.9 KB
Line 
1##General options
2##~~~~~~~~~~~~~~~
3#Run with or without tracer transport ?
4tracer=.true.
5
6#Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
7diurnal=.true.
8
9#Seasonal cycle ? if season=False, Ls stays constant, to value set in "start"
10season = .true.   
11
12#write some more output on the screen ?
13lwrite = .false.   
14
15#Save statistics in file "stats.nc" ?
16stats =.false.
17
18#Save EOF profiles in file "profiles" for Climate Database?
19calleofdump = .false.
20
21## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
22## ~~~~~~~~~~~~~
23#  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
24#  =4 Mars Year 24 from TES assimilation (same as =24 for now)
25#  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
26#  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
27#  =26 Mars Year 26 from TES assimilation
28iaervar = 1
29# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
30tauvis=0.5
31# Dust vertical distribution:
32# (=0: old distrib. (Pollack90), =1: top set by "topdustref",
33#  =2: Viking scenario; =3 MGS scenario)
34iddist  = 3
35# Dust top altitude (km). (Matters only if iddist=1)
36topdustref = 55.
37
38## Physical Parameterizations :
39## ~~~~~~~~~~~~~~~~~~~~~~~~~~
40# call radiative transfer ?
41callrad   = .true.
42# call NLTE radiative schemes ?   matters only if callrad=T
43callnlte = .true.
44# call CO2 NIR absorption ?   matters only if callrad=T
45callnirco2 = .true.
46# call turbulent vertical diffusion ?
47calldifv  = .true.
48# call convective adjustment ?
49calladj  = .true.
50# call CO2 condensation ?
51callcond =.true.
52# call thermal conduction in the soil ?
53callsoil = .true.
54# call Lott's gravity wave/subgrid topography scheme ?
55calllott  = .true.
56
57## Radiative transfer options :
58## ~~~~~~~~~~~~~~~~~~~~~~~~~~
59# the rad.transfer is computed every "iradia" physical timestep
60iradia    = 1
61# Output of the exchange coefficient mattrix ? for diagnostic only
62callg2d  = .false.
63# Rayleigh scattering : (should be .false. for now)
64rayleigh  = .false.
65
66## Tracer (dust water, ice and/or chemical species) options (used if tracer=T):
67## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
68# DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
69dustbin    = 0
70# DUST: Radiatively active dust ? (matters if dustbin>0)
71active  = .false.
72# DUST: use mass and number mixing ratios to predict dust size ?
73# (must also have dustbin=1)
74doubleq   = .false.
75# DUST: lifted by GCM surface winds ?
76lifting = .false.
77# DUST: lifted by dust devils ?
78callddevil = .false.
79# DUST: Scavenging by CO2 snowfall ?
80scavenging = .false.
81# DUST/WATERICE: Gravitationnal sedimentation ?
82sedimentation = .false.
83# WATERICE: Radiatively active transported atmospheric water ice ?
84activice   = .false.
85# WATER: Compute water cycle
86water = .false.
87# WATER: current permanent caps at both poles. True IS RECOMMENDED
88#        (with .true., North cap is a source of water and South pole
89#         is a cold trap)
90caps  = .true.
91# PHOTOCHEMISTRY: include chemical species
92photochem  = .false.
93
94## Thermospheric options (relevant if tracer=T) :
95##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
96# call thermosphere ?
97callthermos = .false.
98# WATER: included without cycle (only if water=.false.)
99thermoswater = .false.
100# call thermal conduction ? (only if callthermos=.true.)
101callconduct = .false.
102# call EUV heating ? (only if callthermos=.true.)
103calleuv=.false.
104# call molecular viscosity ? (only if callthermos=.true.)
105callmolvis = .false.
106# call molecular diffusion ? (only if callthermos=.true.)
107callmoldiff = .false.
108# call thermospheric photochemistry ? (only if callthermos=.true.)
109thermochem = .false.
110# date for solar flux calculation: (1985 < date < 2002)
111## (Solar min=1996.4 ave=1993.4 max=1990.6)
112solarcondate = 1993.4
113       
114
115## albedoes of CO2 ice caps from TES ?
116TESicealbedo=.true.
117## Coefficient for Northern cap albedoes
118TESice_Ncoef=1.4
119## Coefficient for Southern cap albedoes
120TESice_Scoef=1.4
121
122
Note: See TracBrowser for help on using the repository browser.