Changeset 3058

Timestamp:

Sep 28, 2023, 10:48:50 AM (14 months ago)

Author:

emillour

Message:

Mars PCM:
Update reference deftank files:

• remove obsolete callphys.def.notracers
• rename z2sig.def to z2sig.def.MCD6
• update callphys.def.GCM6 to include non-oro GW parametrization
• update callphys.def.MCD6 to disable Frost metamorphism option

EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• deftank/README (modified) (1 diff)
• deftank/callphys.def.GCM6 (modified) (3 diffs)
• deftank/callphys.def.MCD6 (modified) (1 diff)
• deftank/callphys.def.notracers (deleted)
• deftank/z2sig.def.MCD6 (moved) (moved from trunk/LMDZ.MARS/deftank/z2sig.def)
• libf/phymars/vdifc_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4224 +4224 @@
 Fix ticket #141 : when doing a log interpolation in zrecast (e.g. for density), negative and null values are changed into a very small positive value (below the numerical precision) so that the log interpolation yields a null value.
 + add '_FillValue' attribute (redundant with 'missing_value') checks and outputs in zrecast to comply with newest NCO versions (part of ticket #152)
+
+== 28/09/2023 == EM
+Update reference deftank files:
+- remove obsolete callphys.def.notracers
+- rename z2sig.def to z2sig.def.MCD6
+- update callphys.def.GCM6 to include non-oro GW parametrization
+- update callphys.def.MCD6 to disable Frost metamorphism option
+
+

trunk/LMDZ.MARS/deftank/README

@@ -15 +15 @@
 - run.def.lonxlatxalt illustrate cases for different spatial grids
   (typically the GCM is run at 64x48 horizontal resolution)
-- z2sig.def is the current standard vertical grid definition (use 54 layers
+- z2sig.def.MCD6 is the current standard vertical grid definition (use 54 layers
   if running up to ~120km, e.g. water & dust cycle simulations, 73 layers
-  if including chemistry and thermosphere up to ~300km)
+  if including chemistry and thermosphere up to ~250km)
 - *.GCM5 files correspond to GCM5 version using z2sig.def.MCD5
   up to ~120km (32 levels) or ~300km (49 levels)
 - *.MCD5 files correspond to what was used to generate MCDv5.x simulations
 - *.GCM6 files correspond to GCM6 version with improved dust and water cycles
-  up to ~120km (54 levels)
+  up to ~120km (54 levels) using z2sig.def.MCD6
 - *.MCD6 files correspond to what was used to generate MCDv6.1 simulations
+  up to ~250km (73 levels) using z2sig.def.MCD6
 
 >> XML FILES

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

@@ -2 +2 @@
 ##~~~~~~~~~~~~~~~
 #Directory where external input files are:
-datadir=/users/lmdz/WWW/planets/mars/datadir
+datadir=datadir
 
 #Run with or without tracer transport ?
@@ -38 +38 @@
 #  ...
 #  =35 Mars Year 35 from observers assimilation
-iaervar = 35
+iaervar = 8
 # Dust opacity at 610 Pa (when constant, i.e. for the iaervar=1 case)
 tauvis=0.2
@@ -82 +82 @@
 # call thermal conduction in the soil ?
 callsoil = .true.
+
+##Gravity Waves
 # call Lott's gravity wave/subgrid topography scheme ?
 calllott  = .true.
+# call Lott's non-orographic gravity waves scheme?
+calllott_nonoro=.true.
+#Eliassen-Palm FLux(only if calllott_nonoro=.true.)
+nonoro_gwd_epflux_max=5.E-4
+# saturation parameter non-orographic gravity waves(only if calllott_nonoro=.true.)
+nonoro_gwd_sat=1.5
+# value of the dissaption coefficiet(only if calllott_nonoro=.true.)
+nonoro_gwd_rdiss=0.15
+# value of the max wave number
+nonoro_gwd_kmax=1.E-4
+# value of the min wave number
+nonoro_gwd_kmin=7.E-6
+# value to control the launch altitude
+nonoro_gwd_xlaunch=0.6
+
 # Impose polar cap surface albedos as observed by TES?
 TESicealbedo = .true.

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

@@ -189 +189 @@
 frost_albedo_threshold = 0.005
 # Frost metamorphism into old ice
-refill_watercap = .true.
+refill_watercap = .false.
 frost_metam_threshold = 0.05
 # WATER: Water ice thermal inertia (old ice and fresh frost)

trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

@@ -831 +831 @@
            else
 #endif
-            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zcdh,
+            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zicdh,
      &                    pqsurf_tmp(:,igcm_co2),pdqsdif_tmp)
 #ifndef MESOSCALE