Changeset 1042 for trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0.nothermals

Timestamp:

Sep 19, 2013, 2:15:35 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE. def files closer to callphys.def.MCD5

File:

• trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0.nothermals (copied) (copied from trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0) (6 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0.nothermals

@@ -3 +3 @@
 #Run with or without tracer transport ?
 tracer=.false.
+
+#Directory where external input files are:
+datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
+#datadir = /san/home/spiga/DATAFILE
+#datadir = /d5/aslmd/datafile
+#datadir=/users/forget/WWW/datagcm/datafile
 
 #Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
@@ -23 +29 @@
 #  =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario,
 #  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
-#  =24 Mars Year 24 from TES assimilation (ie: MCD reference case)
-#  =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm)
-#  =26 Mars Year 26 from TES assimilation
-iaervar = 26 
-# Dust opacity at 700 Pa (when constant, i.e. for the iaervar=1 case)
-tauvis=0.2
+#  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
+#  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
+#  =30 Mars Year 30 
+iaervar = 24
+iaervar = 26
+# Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
+tauvis=0.1
 # Dust vertical distribution: 
 # (=0: old distrib. (Pollack90), =1: top set by "topdustref",
@@ -41 +48 @@
 callrad   = .true.
 # call NLTE radiative schemes ?   matters only if callrad=T
+callnlte = .true.
 callnlte = .false.
+# NLTE 15um scheme to use.
+# 0-> Old scheme, static oxygen
+# 1-> Old scheme, dynamic oxygen
+# 2-> New scheme
+nltemodel = 2 
+nltemodel = 0
 # call CO2 NIR absorption ?   matters only if callrad=T
 callnirco2 = .true.
+# NIR NLTE correction ? matters only if callnirco2=T
+nircorr=1
+nircorr=0
 # call turbulent vertical diffusion ?
 calldifv  = .true.
 # call convective adjustment ?
 calladj  = .true.
+# Thermals
+calltherm = .true.
+callrichsl = .true.
+calltherm = .false.
+callrichsl = .false.
 # call CO2 condensation ?
 callcond =.true.
@@ -56 +78 @@
 # Impose polar cap surface albedos as observed by TES?
 TESicealbedo = .true.
+TESicealbedo = .false.
+## Coefficient for Northern cap albedoes
+TESice_Ncoef=1.6
+## Coefficient for Southern cap albedoes
+TESice_Scoef=1.6
+
 
 ## Radiative transfer options :
@@ -69 +97 @@
 ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # DUST: Transported dust ? (if >0, use 'dustbin' dust bins)
+dustbin    = 2
 dustbin    = 0
 # DUST: Radiatively active dust ? (matters if dustbin>0)
+active  = .true.
 active  = .false.
 # DUST: use mass and number mixing ratios to predict dust size ?
+doubleq   = .true.
 doubleq   = .false. 
 # DUST: lifted by GCM surface winds ?
-lifting = .false. 
+lifting = .true.
+lifting = .false.
 # DUST: lifted by dust devils ?
 callddevil = .false.
-# DUST: Scavenging by CO2 snowfall ?
+# DUST: Scavenging by H2O/CO2 snowfall ?
+scavenging = .true.
 scavenging = .false.
 # DUST/WATERICE: Gravitationnal sedimentation ?
+sedimentation = .true.
 sedimentation = .false.
 # WATERICE: Radiatively active transported atmospheric water ice ?
+activice   = .true.
 activice   = .false.
 # WATER: Compute water cycle
+water = .true.
 water = .false.
+# WATER: Microphysical scheme for water-ice clouds?
+microphys = .true.
+microphys = .false.
+# WATER: parameter contact
+mteta = 0.95
+# WATER: Effective variance for sedimentation for the log-normal 
+#        distribution of ice particles ?
+nuice_sed=0.1
 # WATER: current permanent caps at both poles. True IS RECOMMENDED
 #        (with .true., North cap is a source of water and South pole
 #         is a cold trap)
+caps  = .true.
 caps  = .false.
+# WATER: Water ice albedo ?
+albedo_h2o_ice = 0.35
+# WATER: Water ice thermal inertia
+inert_h2o_ice = 800
+# WATER: Frost thickness threshold for albedo
+frost_albedo_threshold = 0.005
 # PHOTOCHEMISTRY: include chemical species 
+photochem  = .true.
 photochem  = .false.
 
@@ -96 +148 @@
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # call thermosphere ?
+callthermos = .true.
 callthermos = .false.
 # WATER: included without cycle (only if water=.false.)
 thermoswater = .false.
 # call thermal conduction ? (only if callthermos=.true.)
+callconduct = .true.
 callconduct = .false.
-# call EUV heating ? (only if callthermos=.true.)
-calleuv=.false.
 # call molecular viscosity ? (only if callthermos=.true.)
+callmolvis = .true.
 callmolvis = .false.
 # call molecular diffusion ? (only if callthermos=.true.)
+callmoldiff = .true.
 callmoldiff = .false.
 # call thermospheric photochemistry ? (only if callthermos=.true.)
+thermochem = .true.
 thermochem = .false.
+# call EUV heating ? (only if callthermos=.true.)
+calleuv=.true.
+calleuv=.false.
+#Method to include solar variability?
+#0-> Old method (fixed EUV input)  1-> Variability with E10.7 as observed
+solvarmod=0
 # date for solar flux calculation: (1985 < date < 2002)
-## (Solar min=1996.4 ave=1993.4 max=1990.6)
+## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
 solarcondate = 1993.4
-       
+#Solar variability as observed for MY? (must bebetween MY23 and MY30)
+# (only matters if solvarmod=1)
+solvaryear=24
+# value for the UV heating efficiency
+##(experimental values between 0.19 and 0.23, lower values may
+## be used to compensate for low 15 um cooling)
+euveff = 0.21       
 
 
+