Changeset 705 for trunk/LMDZ.MARS/deftank

Timestamp:

Jun 14, 2012, 4:14:11 PM (13 years ago)

Author:

emillour

Message:

Mars GCM:

• Added possibility to run with a varying EUV cycle following real one. The flag solvarmod=1 triggers this behaviour, with companion flag solvaryear=## , where ## is the Mars Year (from 23 to 30). Setting solvarmod=0 reverts to 'old' behaviour, where there is a constant EUV forcing throughout the run, set by the "solarcondate" flag.
• Needs corresponding input data files ("param_v6" subdirectory of "EUV" subdirectory in "datadir").
• Added files jthermcalc_e107.F and param_read_e107.F in "aeronomars", modified files euvheat.F90, hrtherm.F, chemthermos.F90, param_v4.h and param_read.F in "aeronomars" and files inifis.F, physiq.F and callkeys.h in "phymars".

FGG

File:

• trunk/LMDZ.MARS/deftank/callphys.def (modified) (4 diffs)

trunk/LMDZ.MARS/deftank/callphys.def

@@ -45 +45 @@
 # 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 = 1
 # call CO2 NIR absorption ?   matters only if callrad=T
 callnirco2 = .true.
+# NIR NLTE correction ? matters only if callnirco2=T
+nircorr=1
 # call turbulent vertical diffusion ?
 calldifv  = .true.
@@ -59 +66 @@
 # Impose polar cap surface albedos as observed by TES?
 TESicealbedo = .true.
+## Coefficient for Northern cap albedoes
+TESice_Ncoef=1.6
+## Coefficient for Southern cap albedoes
+TESice_Scoef=1.6
+
 
 ## Radiative transfer options :
@@ -104 +116 @@
 # 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.
@@ -112 +122 @@
 # call thermospheric photochemistry ? (only if callthermos=.true.)
 thermochem = .false.
+# call EUV heating ? (only if callthermos=.true.)
+calleuv=.false.
+#Method to include solar variability?
+#0-> Old method   1-> Variability with E10.7 as observed
+solvarmod=1
 # 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       
 
-