Changeset 2477 for trunk

Timestamp:

Mar 15, 2021, 9:33:56 AM (4 years ago)

Author:

aslmd

Message:

MESOSCALE workflow: updated .def to latest version: English for run.def, upper atmosphere settings for callphys.def, added callnlte=F to mesoscale.def to allow it for GCM without necessarily doing it for mesoscale

Location:

trunk/MESOSCALE/LMD_MM_MARS/SIMU/MESORUN

Files:

• callphys.def (modified) (3 diffs)
• gcm/run.def (modified) (1 diff)
• mesoscale.def (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SIMU/MESORUN/callphys.def

@@ -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 "clim" 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 
-iaervar = 8
+iaervar = 34
 # Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
 tauvis=0.1
@@ -46 +46 @@
 # call NLTE radiative schemes ?   matters only if callrad=T
 callnlte = .true.
-callnlte = .false.
+#callnlte = .false.
 # NLTE 15um scheme to use.
 # 0-> Old scheme, static oxygen
@@ -158 +158 @@
 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) ; Only used if solvarmod=0
-solarcondate = 1993.4
-#Solar variability as observed for MY? (must bebetween MY23 and MY30)
+#0-> Fixed EUV input  1-> Variability with E10.7 as observed
+solvarmod=1
+# 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
+solvaryear=34
 # value for the UV heating efficiency
 ##(experimental values between 0.19 and 0.23, lower values may

trunk/MESOSCALE/LMD_MM_MARS/SIMU/MESORUN/gcm/run.def

@@ -1 +1 @@
 #
 #----------------------------------------------------------------------- 
-#Parametres de controle du run:                                          
-#------------------------------                                          
+#GCM run control parameters:                                          
+#---------------------------                                     
 
 # planet type
    planet_type=mars
-   
-# Nombre de jours d'integration                                         
-     nday=3
 
-# nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
+# Number of days to run model for                                      
+ nday= 3
+
+# 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=200
-
-# periode de stockage fichier histmoy (en jour)                         
- periodav=60.
-
-# periode de la dissipation (en pas)                                    
+# Apply dissipation every idissip dynamical steps
   idissip=1
 
-# 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            
+# iterate lateral dissipation operator gradiv nitergdiv times 
 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)  
- tetagdiv=1800.
+# time scale (s) for shortest wavelengths for u,v (gradiv)
+ tetagdiv=2500.
 
-# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
- tetagrot=3600.
+# time scale (s) for shortest wavelengths for u,v (nxgradrot)
+ tetagrot=5000.
 
-# temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
- tetatemp=3600.
+# time scale (s) for shortest wavelengths for h (divgrad)
+ tetatemp=5000.
 
-# coefficient pour gamdissip                                            
+# multiplicative constants for dissipation with altitude:
+# coefficient for middle atmosphere (~20-70km)
+dissip_fac_mid = 3
+# coefficient for upper atmosphere (~100km+)
+dissip_fac_up = 30
+
+# 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=10
+# call physics every iphysiq dynamical steps
+  iphysiq=5
 
-# choix d'une grille reguliere                                          
+# Use a regular grid 
   grireg=.true.
 
-# frequence (en pas) de l'ecriture du fichier diagfi               
- ecritphy=40
+# Output in diagfi file every ecritphy dynamical steps 
+ ecritphy=80
 
-# 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= 3
 
-# Sponge:  hauteur de sponge (km)
-  hsponge= 230
-
-# Sponge:  tetasponge (secondes)
+# Sponge layer time scale (s):  tetasponge
   tetasponge = 30000
 

trunk/MESOSCALE/LMD_MM_MARS/SIMU/MESORUN/mesoscale.def

@@ -4 +4 @@
 calllott     = .false.
 TESicealbedo = .false.
+callnlte     = .false.
 ##########################
 # your settings below