Changeset 3776

Timestamp:

May 23, 2025, 1:09:29 PM (2 months ago)

Author:

slebonnois

Message:

SL: update of Venus PCM deftank for consistency with current sources

Location:

trunk/LMDZ.VENUS/deftank

Files:

• field_def_physics.xml (modified) (4 diffs)
• gcm-96x96x50-chemistry.def (modified) (2 diffs)
• gcm-96x96x50-notrac.def (modified) (2 diffs)
• gcm.def (modified) (2 diffs)
• iodef.xml (modified) (1 diff)
• physiq-96x96x50-notrac.def (modified) (1 diff)
• physiq-96x96x90-chemistry-IONOSPHERE-IONDIFF.def (modified) (2 diffs)

trunk/LMDZ.VENUS/deftank/field_def_physics.xml

@@ -123 +123 @@
                    long_name="Column density of SO3" 
                    unit="kg/m2" />
-            <field id="col_s2o2"
-                   long_name="Column density of S2O2" 
+            <field id="col_osso_cis"
+                   long_name="Column density of OSSO CIS" 
+                   unit="kg/m2" />
+            <field id="col_osso_trans"
+                   long_name="Column density of OSSO TRANS" 
+                   unit="kg/m2" />
+            <field id="col_s2o2_cyc"
+                   long_name="Column density of S2O2 CYC"
                    unit="kg/m2" />
             <field id="col_ocs"
@@ -284 +290 @@
                    long_name="Kinetic Energy net flux" 
                    unit="W/m2" />
+            <field id="BV2"
+                   long_name="Freq Brunt Väisälä"
+                   unit="Hz^2" />
            <!---- TRACERS ---->
             <field id="age"
@@ -367 +376 @@
                    long_name="SO3" 
                    unit="mol/mol" />
-            <field id="s2o2"
-                   long_name="S2O2" 
+            <field id="osso_cis"
+                   long_name="OSSO CIS" 
+                   unit="mol/mol" />
+            <field id="osso_trans"
+                   long_name="OSSO TRANS" 
+                   unit="mol/mol" />
+            <field id="s2o2_cyc"
+                   long_name="S2O2 CYC"
                    unit="mol/mol" />
             <field id="ocs"
@@ -481 +496 @@
                    long_name="Diff molec N2" 
                    unit="kg/kg" />
-            <field id="BV2"
-                   long_name="Freq Brunt Väisälä"
-                   unit="Hz^2" />
             <field id="iter"
                    long_name="chemical iterations" 

trunk/LMDZ.VENUS/deftank/gcm-96x96x50-chemistry.def

@@ -21 +21 @@
 niterh=2
 ## dissipation time scale (s) for shortest wavelengths for u,v (gradiv)  
-tetagdiv=1.e4
+tetagdiv=3000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
-tetagrot=1.e4
+tetagrot=4000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
-tetatemp=1.e4
+tetatemp=4000.
 ## coefficient for gamdissip                                            
 coefdis=0.
@@ -50 +50 @@
 ## Horizontal dissipation multipliers along the vertical
 dissip_fac_mid=2.
-dissip_fac_up=50.
+dissip_fac_up=30.
 # deltaz et hdelta in km
-dissip_deltaz=30.
-dissip_hdelta=5.
+dissip_deltaz=20.
+dissip_hdelta=10.
 # pupstart in Pa
-dissip_pupstart=1.e4
+dissip_pupstart=1.
 
 ## Sponge layer

trunk/LMDZ.VENUS/deftank/gcm-96x96x50-notrac.def

@@ -21 +21 @@
 niterh=2
 ## dissipation time scale (s) for shortest wavelengths for u,v (gradiv)  
-tetagdiv=1.e4
+tetagdiv=3000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
-tetagrot=1.e4
+tetagrot=4000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
-tetatemp=1.e4
+tetatemp=4000.
 ## coefficient for gamdissip                                            
 coefdis=0.
@@ -50 +50 @@
 ## Horizontal dissipation multipliers along the vertical
 dissip_fac_mid=2.
-dissip_fac_up=50.
+dissip_fac_up=30.
 # deltaz et hdelta in km
-dissip_deltaz=30.
-dissip_hdelta=5.
+dissip_deltaz=20.
+dissip_hdelta=10.
 # pupstart in Pa
-dissip_pupstart=1.e4
+dissip_pupstart=1.
 
 ## Sponge layer

trunk/LMDZ.VENUS/deftank/gcm.def

@@ -21 +21 @@
 niterh=2
 ## dissipation time scale (s) for shortest wavelengths for u,v (gradiv)  
-tetagdiv=1.e4
+tetagdiv=3000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
-tetagrot=1.e4
+tetagrot=4000.
 ## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
-tetatemp=1.e4
+tetatemp=4000.
 ## coefficient for gamdissip                                            
 coefdis=0.
@@ -50 +50 @@
 ## Horizontal dissipation multipliers along the vertical
 dissip_fac_mid=2.
-dissip_fac_up=50.
+dissip_fac_up=30.
 # deltaz et hdelta in km
-dissip_deltaz=30.
-dissip_hdelta=5.
+dissip_deltaz=20.
+dissip_hdelta=10.
 # pupstart in Pa
-dissip_pupstart=1.e4
+dissip_pupstart=1.
 
 ## Sponge layer

trunk/LMDZ.VENUS/deftank/iodef.xml

@@ -17 +17 @@
                     <!-- set info_level between   0 (== none) and
                                                 100 (==lots) -->
-                    <variable id="info_level" type="int">10</variable>
+                    <variable id="info_level" type="int">100</variable>
 
                     <!-- set print_file to true to have XIOS write its

trunk/LMDZ.VENUS/deftank/physiq-96x96x50-notrac.def

@@ -41 +41 @@
 # Reinitialization of tracer abundances 
 reinit_trac=n
+#
+# Parameters for Age of air
+##~~~~~~~~~~~~~~~~~~~~~~~~~
+# Compute age of air (requires that there is an "aoa" tracer; default: n)
+ok_aoa=n
+# Reinitialize age of air tracer (default: y)
+reinit_aoa=n
+# Emission altitude (in model layers) of aoa tracer emission (default: 1)
+lev_aoa=1
 #
 # parameters for chemistry and microphysics

trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry-IONOSPHERE-IONDIFF.def

@@ -14 +14 @@
 nbapp_rad=24000
 ## Number of calls to the chemistry routines (per day)
-nbapp_chem=24000
+nbapp_chem=48000
 ##  Flag for convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(new version JYG), 30 KE(version IPCC AR4), 4 KE vect
 iflag_con=0
@@ -109 +109 @@
 ## include ion diffusion or not | WARNING: it is 2.5-3. time slower with
 #  BE CAREFULL IF (ok_iondiff = y)
-#  ==> nbapp_chem need to be equal to day_step in gcm.def
-ok_iondiff = n
+#  ==> nbapp_chem need to be equal to day_step/iphysiq (gcm.def)
+ok_iondiff = y
 
 #Method to include solar variability?