Changeset 1949 for trunk/LMDZ.MARS/deftank/callphys.def.watercycle

Timestamp:

Jun 19, 2018, 12:06:14 PM (7 years ago)

Author:

emillour

Message:

Mars GCM:
Some cleanup in the reference def files.
EM

File:

• trunk/LMDZ.MARS/deftank/callphys.def.watercycle (moved) (moved from trunk/LMDZ.MARS/deftank/callphys.def.outliers) (6 diffs)

trunk/LMDZ.MARS/deftank/callphys.def.watercycle

@@ -20 +20 @@
 
 #Directory where external input files are:
-datadir=/users/forget/WWW/datagcm/datafile/
+datadir=/users/lmdz/WWW/planets/mars/datadir
 
 ## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F)
@@ -27 +27 @@
 #  =4 Mars Year 24 from TES assimilation (old version of MY24; dust_tes.nc file)
 #  =6 "cold" (low dust) scenario ; =7 "warm" (high dust) scenario
+#  =8 "climatology" (our best guess of a typical Mars year) scenario
 #  =24 Mars Year 24  ; =25 Mars Year 25 (year with a global dust storm) ; ...
 #  =30 Mars Year 30 
@@ -45 +46 @@
 # call NLTE radiative schemes ?   matters only if callrad=T
 callnlte = .true.
+# NLTE 15um scheme to use.
+# 0-> Old scheme, static oxygen
+# 1-> Old scheme, dynamic oxygen
+# 2-> New scheme
+nltemodel = 0
 # call CO2 NIR absorption ?   matters only if callrad=T
 callnirco2 = .true.
+# NIR NLTE correction ? matters only if callnirco2=T
+nircorr=0
 # call turbulent vertical diffusion ?
 calldifv  = .true.
 # call convective adjustment ?
 calladj  = .true.
+# Thermals
+calltherm = .true.
+callrichsl = .true.
 # call CO2 condensation ?
 callcond =.true.
@@ -60 +71 @@
 TESicealbedo = .true.
 ## Coefficient for Northern cap albedoes
-TESice_Ncoef=1.4
+TESice_Ncoef=1.6
 ## Coefficient for Southern cap albedoes
-TESice_Scoef=1.4
+TESice_Scoef=1.6
 
 
@@ -130 +141 @@
 # call thermal conduction ? (only if callthermos=.true.)
 callconduct = .false.
-# call EUV heating ? (only if callthermos=.true.)
-calleuv=.false.
 # call molecular viscosity ? (only if callthermos=.true.)
 callmolvis = .false.
@@ -138 +147 @@
 # call thermospheric photochemistry ? (only if callthermos=.true.)
 thermochem = .false.
-# date for solar flux calculation: (1985 < date < 2002)
-## (Solar min=1996.4 ave=1993.4 max=1990.6)
-solarcondate = 1993.4
-       
-## Thermals
-## ~~~~~~~~
-calltherm = .true.
-callrichsl = .true.
+# call EUV heating ? (only if callthermos=.true.)
+calleuv=.false.
+#Method to include solar variability?
+#0-> Fixed EUV input  1-> Variability with E10.7 as observed
+solvarmod=0
+# fixed E10.7 value (for solvarmod=0)
+# (min=80 , ave=140, max=320)
+fixed_euv_value=140
+#Solar variability as observed for MY? (must be between MY23 and MY32)
+# (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       
 
+