Changeset 1949 for trunk/LMDZ.MARS/deftank

Timestamp:

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

Author:

emillour

Message:

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

Location:

trunk/LMDZ.MARS/deftank

Files:

• callphys.def.MCD5 (modified) (2 diffs)
• callphys.def.notracers (moved) (moved from trunk/LMDZ.MARS/deftank/callphys.def)
• callphys.def.watercycle (moved) (moved from trunk/LMDZ.MARS/deftank/callphys.def.outliers) (6 diffs)
• run.def.32x24x25 (modified) (1 diff)
• run.def.32x24x32 (modified) (1 diff)
• run.def.64x48x25 (deleted)
• run.def.64x48x32 (modified) (1 diff)
• traceur.def.scavenging (deleted)
• traceur.def.watercycle (modified) (1 diff)

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

@@ -5 +5 @@
 
 #Directory where external input files are:
-datadir=/planeto/emlmd/run_mcd5.v0_new/datadir
+datadir=/users/lmdz/WWW/planets/mars/datadir
 
 #Diurnal cycle ?  if diurnal=False, diurnal averaged solar heating
@@ -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 

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       
 
+

trunk/LMDZ.MARS/deftank/run.def.32x24x25

@@ -3 +3 @@
 #GCM run control parameters:                                          
 #---------------------------                                     
+
+# planet type
+   planet_type=mars
 
 # Number of days to run model for                                      

trunk/LMDZ.MARS/deftank/run.def.32x24x32

@@ -3 +3 @@
 #GCM run control parameters:                                          
 #---------------------------                                     
+
+# planet type
+   planet_type=mars
 
 # Number of days to run model for                                      

trunk/LMDZ.MARS/deftank/run.def.64x48x32

@@ -1 +1 @@
 #
 #----------------------------------------------------------------------- 
-#Parametres de controle du run:                                          
-#------------------------------                                          
+#GCM run control parameters:                                          
+#---------------------------                                          
 
-# Nombre de jours d'integration                                         
+# planet type
+   planet_type=mars
+
+# Number of days to run model for                                      
      nday=9999
 
-# nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+# Number of dynamical steps per day (must be a multiple of iperiod)
  day_step = 960
 
-# periode pour le pas Matsuno (en pas)                                  
+# Apply a Matsuno step every iperiod dynamical step
   iperiod=5
 
-# periode de sortie des variables de controle (en pas)                  
+# Control output information in the dynamics every iconser dynamical steps
   iconser=120
 
-# periode d'ecriture du fichier histoire (en jour)                      
-    iecri=100
-
-# periode de stockage fichier histmoy (en jour)                         
- periodav=60.
-
-# periode de la dissipation (en pas)                                    
+# Apply dissipation every idissip dynamical steps
   idissip=5
 
-# choix de l'operateur de dissipation (star ou  non star )              
+# dissipation operator to use (star or non-star)
  lstardis=.true.
 
-# avec ou sans coordonnee hybrides
+# use hybrid vertical coordinate (else will use sigma levels)
  hybrid=.true.
 
-# nombre d'iterations de l'operateur de dissipation   gradiv            
+# use hybrid vertical coordinate (else will use sigma levels)
 nitergdiv=1
 
-# nombre d'iterations de l'operateur de dissipation  nxgradrot          
+# iterate lateral dissipation operator nxgradrot nitergrot times         
 nitergrot=2
 
-# nombre d'iterations de l'operateur de dissipation  divgrad            
+# iterate lateral dissipation operator divgrad niterh times          
    niterh=2
 
-# temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
+# time scale (s) for shortest wavelengths for u,v (gradiv)
  tetagdiv= 4000.
 
-# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+# time scale (s) for shortest wavelengths for u,v (nxgradrot)
  tetagrot=5000.
 
-# temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
+# time scale (s) for shortest wavelengths for h (divgrad)
  tetatemp=5000.
 
-# coefficient pour gamdissip                                            
+# coefficient for gamdissip                                            
   coefdis=0.
 
-# choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+# time marching scheme (Matsuno if purmats is true, else Matsuno-Leapfrog) 
   purmats=.false.
 
-# avec ou sans physique                                                 
+# run with (true) or without (false) physics
    physic=.true.
 
-# periode de la physique (en pas)                                       
-  iphysiq=20
+# call physics every iphysiq dynamical steps
+  iphysiq=10
 
-# choix d'une grille reguliere                                          
+# Use a regular grid 
   grireg=.true.
 
-# frequence (en pas) de l'ecriture du fichier diagfi               
+# Output in diagfi file every ecritphy dynamical steps 
  ecritphy=240
 
-# longitude en degres du centre du zoom                                 
+# longitude (degrees) of zoom center
    clon=63.
 
-# latitude en degres du centre du zoom                                  
+# latitude (degrees) of zoom center
    clat=0.
 
-# facteur de grossissement du zoom,selon longitude                      
+# enhancement factor of zoom, along longitudes
   grossismx=1.
 
-# facteur de grossissement du zoom ,selon latitude                      
+# enhancement factor of zoom, along latitudes                      
  grossismy=1.
 
-#  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+#  Use an hyperbolic function f(y) if .true., else use a sine      
   fxyhypb=.false.
 
-# extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+# extention along longitudes of zoom region (fraction of global domain)
    dzoomx= 0.
 
-# extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+# extention along latitudes of zoom region (fraction of global domain)
    dzoomy=0.
 
-#  raideur du zoom en  X
+# zoom stiffness along longitudes
     taux=2.
 
-#  raideur du zoom en  Y
+# zoom stiffness along latitudes
     tauy=2.
 
-#  Fonction  f(y) avec y = Sin(latit.) si = .TRUE. ,  Sinon  y = latit.
+#  Fonction  f(y) as y = Sin(latitude) if = .true. ,  else  y = latitude
   ysinus= .false.
 
-# Avec sponge layer
+# Use a sponge layer
   callsponge  = .true.
  
+# Sponge layer extends over topmost nsponge layers
+  nsponge = 3
+
 # Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
-  mode_sponge= 2
+  mode_sponge= 3
 
-# Sponge:  tetasponge (secondes)
+# Sponge layer time scale (s):  tetasponge
   tetasponge = 30000
 

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

@@ -1 @@
-2
+7
+co2
+dust_number
+dust_mass
+ccn_number
+ccn_mass
 h2o_ice
 h2o_vap