Changeset 4080 for trunk/LMDZ.VENUS/deftank

Timestamp:

Feb 20, 2026, 11:39:45 AM (3 weeks ago)

Author:

streel

Message:

!!! Important : The microphysical scheme is still a work in progress, and I was trying to fix a bug prior this commit. Ehouarn is planning to review the tracers in the model and to facilitate his work I am now commiting an unfinished version of the microphysical scheme.
For users, the Aurélien Stolzenbach parametrization remains untouched with this commit, meaning that this commit will act as a transparent one for your current simulations. Thanks for your understanding. !!!

Patch note :

First deployment of the microphysical scheme.

Modifications of the Venus Physics

• Removal of the old files / Addition of the new ones linked to the microphysical scheme in the cloud venus folder.
• Modification of cles_phys_mod and conf_phys : Additions of the new keywords used by the microphysical scheme
• Modification of chemparam_mod / phys_state_var_mod : Adding the new tracers and parameters used by the microphysical scheme.
• Modification of phytrac_chimie : Adding the call to the microphysical scheme and computation of the tendencies.
• Modification of physiq_mod : New sedimentation of the microphysical tracers replacing the older ones. Computation of the final mixing ratios of the microphysical tracers. Outputs of the new variables.

New definitions files for the simulations

• New tracer_mic.def, traceur_ion_mic.def
• Update of the differents physiq.def (with a new microphysical sub part)
• Update of the xml files, with the new possible outputs

Miscellaneous

• Fixed a typo in a comment in photolysis_mod

All the information related to the microphysical scheme will be sum-up and updated on the wiki page : ​https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Venus_PCM_-_Microphysical_scheme

NS

Location:

trunk/LMDZ.VENUS/deftank

Files:

• field_def_physics.xml (modified) (1 diff)
• file_def_physics-chemcld.xml (added)
• physiq-96x96x50-chemistry.def (modified) (1 diff)
• physiq-96x96x78-chemistry.def (modified) (1 diff)
• physiq-96x96x90-chemistry-IONOSPHERE-IONDIFF.def (modified) (1 diff)
• physiq-96x96x90-chemistry.def (modified) (1 diff)
• physiq.def (modified) (1 diff)
• traceur-chemcld.def (modified) (2 diffs)

trunk/LMDZ.VENUS/deftank/field_def_physics.xml

@@ -510 +510 @@
             <field id="o2_emis"
                    long_name="Volume emission rate O2 emission"
-                   unit="photons.cm-3.s-1" />
+                    unit="photons.cm-3.s-1" />
+            <field id="h2otot"
+                   long_name="Total vmr of H2O"
+                    unit="mol/mol" />
+            <field id="h2so4tot"
+                   long_name="Total vmr of H2SO4"
+                    unit="mol/mol" />
 
-       </field_group>
+           <!---- MICROPHYSICS TRACERS ---->
+            <field id="m0_aer"
+                   long_name="Order 0 moment of the aerosols"
+                    unit="m-3" />
+            <field id="m0_mode1drop"
+                   long_name="Order 0 moment of mode 1 droplets"
+                    unit="m-3" />
+            <field id="m0_mode1ccn"
+                   long_name="Order 0 moment of mode 1 ccn"
+                    unit="m-3" />
+            <field id="m0_mode2drop"
+                   long_name="Order 0 moment of mode 2 droplets"
+                    unit="m-3" />
+            <field id="m0_mode2ccn"
+                   long_name="Order 0 moment of mode2 ccn"
+                    unit="m-3" />
+            <field id="m3_aer"
+                   long_name="Order 3 moment of the aerosols"
+                    unit="m-3/m-3" />
+            <field id="m3_mode1sa"
+                   long_name="Order 3 moment of mode 1 Sulfuric Acid"
+                    unit="m-3/m-3" />
+            <field id="m3_mode1wv"
+                   long_name="Order 3 moment of mode 1 Water Vapor"
+                    unit="m-3/m-3" />
+            <field id="m3_mode1ccn"
+                   long_name="Order 3 moment of mode 1 ccn"
+                    unit="m-3/m-3" />
+            <field id="m3_mode2sa"
+                   long_name="Order 3 moment of mode 2 Sulfuric Acid"
+                    unit="m-3/m-3" />
+            <field id="m3_mode2wv"
+                   long_name="Order 3 moment of mode 2 Water Vapor"
+                    unit="m-3/m-3" />
+            <field id="m3_mode2ccn"
+                   long_name="Order 3 moment of mode 2 ccn"
+                    unit="m-3/m-3" />
+            <field id="r1"
+                   long_name="Radius of mode 1 droplets"
+                    unit="m" />
+            <field id="r2"
+                   long_name="Radius of mode 2 droplets"
+                    unit="m" />
+            <field id="RHtab"
+                   long_name="Relative humidity of H2SO4"
+                    unit="" />
+            <field id="rhotab"
+                   long_name="Density of droplets"
+                    unit="kg/m3" />
+            <field id="WSAVtab"
+                   long_name="mass fraction of SA of droplets"
+                    unit="" />
+            <field id="dM0m1_hom"
+                   long_name="Tendencie of the numb of drop by hom nuc"
+                    unit="/m3" />
+            <field id="dM3m1_hom_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by hom nuc"
+                    unit="m3/m3" />
+            <field id="dM3m1_hom_h2o"
+                   long_name="Tendencie of the volume of h2o by hom nuc"
+                    unit="m3/m3" />
+            <field id="dM0m1_het"
+                   long_name="Tendencie of the numb of drop by het nuc"
+                    unit="/m3" />
+            <field id="dM0m1_het_ccn"
+                   long_name="Tendencie of the numb of ccn by het nuc"
+                    unit="/m3" />
+            <field id="dM3m1_het_ccn"
+                   long_name="Tendencie of the volume of ccn by het nuc"
+                    unit="m3/m3" />
+            <field id="dM3m1_het_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by het nuc"
+                    unit="m3/m3" />
+            <field id="dM3m1_het_h2o"
+                   long_name="Tendencie of the volume of h2o by het nuc"
+                    unit="m3/m3" />
+            <field id="dM0_het_aer"
+                   long_name="Tendencie of the numb of aer by het nuc"
+                    unit="/m3" />
+            <field id="dM3_het_aer"
+                   long_name="Tendencie of the volume of aer by het nuc"
+                    unit="/m3" />
+            <field id="dM3m1_cond_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by cond"
+                    unit="/m3" />
+            <field id="dM3m1_cond_h2o"
+                   long_name="Tendencie of the volume of h2o by cond"
+                    unit="/m3" />
+            <field id="dM3m2_cond_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by cond"
+                    unit="/m3" />
+            <field id="dM3m2_cond_h2o"
+                   long_name="Tendencie of the volume of h2o by cond"
+                    unit="/m3" />
+            <field id="dM0m1_coag"
+                   long_name="Tendencie of the numb of drop by coag"
+                    unit="/m3" />
+            <field id="dM0m1_coag_ccn"
+                   long_name="Tendencie of the numb of ccn by coag"
+                    unit="/m3" />
+            <field id="dM3m1_coag_ccn"
+                   long_name="Tendencie of the volume of ccn by coag"
+                    unit="m3/m3" />
+            <field id="dM3m1_coag_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by coag"
+                    unit="m3/m3" />
+            <field id="dM3m1_coag_h2o"
+                   long_name="Tendencie of the volume of h2o by coag"
+                    unit="m3/m3" />
+            <field id="dM0m2_coag"
+                   long_name="Tendencie of the numb of drop by coag"
+                    unit="/m3" />
+            <field id="dM0m2_coag_ccn"
+                   long_name="Tendencie of the numb of ccn by coag"
+                    unit="/m3" />
+            <field id="dM3m2_coag_ccn"
+                   long_name="Tendencie of the volume of ccn by coag"
+                    unit="m3/m3" />
+            <field id="dM3m2_coag_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by coag"
+                    unit="m3/m3" />
+            <field id="dM3m2_coag_h2o"
+                   long_name="Tendencie of the volume of h2o by coag"
+                    unit="m3/m3" />
+            <field id="dM0m1_merge"
+                   long_name="Tendencie of the numb of drop by mode-merge"
+                    unit="/m3" />
+            <field id="dM0m1_merge_ccn"
+                   long_name="Tendencie of the numb of ccn by mode-merge"
+                    unit="/m3" />
+            <field id="dM3m1_merge_ccn"
+                   long_name="Tendencie of the volume of ccn by mode-merge"
+                    unit="m3/m3" />
+            <field id="dM3m1_merge_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by mode-merge"
+                    unit="m3/m3" />
+            <field id="dM3m1_merge_h2o"
+                   long_name="Tendencie of the volume of h2o by mode-merge"
+                    unit="m3/m3" />
+            <field id="dM0m2_merge"
+                   long_name="Tendencie of the numb of drop by mode-merge"
+                    unit="/m3" />
+            <field id="dM0m2_merge_ccn"
+                   long_name="Tendencie of the numb of ccn by mode-merge"
+                    unit="/m3" />
+            <field id="dM3m2_merge_ccn"
+                   long_name="Tendencie of the volume of ccn by mode-merge"
+                    unit="m3/m3" />
+            <field id="dM3m2_merge_h2so4"
+                   long_name="Tendencie of the volume of h2so4 by mode-merge"
+                    unit="m3/m3" />
+            <field id="dM3m2_merge_h2o"
+                   long_name="Tendencie of the volume of h2o by mode-merge"
+                    unit="m3/m3" />
+            <field id="satm1"
+                   long_name="sat ratio of mode 1"
+                    unit="" />
+            <field id="satm2"
+                   long_name="sat ratio of mode 2"
+                    unit="" />
+        </field_group>
 
     </field_definition>

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

@@ -60 +60 @@
 ok_sedim=y
 nb_mode=3
+#
+# parameters for microphysics (cl_scheme = 2)
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# Verboose debug mode
+debug_cloud = n
+#Keys to calls the different microphysics processes or not
+mode_merging = y
+hom_nuclea = y
+het_nuclea = y
+cond_evap = y
+coag = y
+#
+# Regarding the mode 1 and 2 size distributions
+sigma_mode1 = 1.56
+sigma_mode2 = 1.29
+#
+# Aerosol production (precursors off the droplets)
+prod_aer = y
+#Characteristics of the size distribution of the aer
+rad_aer = 0.125e-6
+sigma_aer = 1.56
+rho_aer = 1800
+# Vertical distribution and production rate
+alt_aer = 62000
+sigmaz_aer = 2000
+tx_prod_aer = 1.D-10
 #
 # parameters for the boundary layer

trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def

@@ -60 +60 @@
 ok_sedim=y
 nb_mode=3
+#
+# parameters for microphysics (cl_scheme = 2)
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# Verboose debug mode
+debug_cloud = n
+#Keys to calls the different microphysics processes or not
+mode_merging = y
+hom_nuclea = y
+het_nuclea = y
+cond_evap = y
+coag = y
+#
+# Regarding the mode 1 and 2 size distributions
+sigma_mode1 = 1.56
+sigma_mode2 = 1.29
+#
+# Aerosol production (precursors off the droplets)
+prod_aer = y
+#Characteristics of the size distribution of the aer
+rad_aer = 0.125e-6
+sigma_aer = 1.56
+rho_aer = 1800
+# Vertical distribution and production rate
+alt_aer = 62000
+sigmaz_aer = 2000
+tx_prod_aer = 1.D-10
 #
 # parameters for the boundary layer

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

@@ -60 +60 @@
 ok_sedim=y
 nb_mode=3
+# Verboose microphysical model
+debug_cloud=n
 #
 # parameters for the boundary layer

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

@@ -60 +60 @@
 ok_sedim=y
 nb_mode=3
+#
+# parameters for microphysics (cl_scheme = 2)
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# Verboose debug mode
+debug_cloud = n
+#Keys to calls the different microphysics processes or not
+mode_merging = y
+hom_nuclea = y
+het_nuclea = y
+cond_evap = y
+coag = y
+#
+# Regarding the mode 1 and 2 size distributions
+sigma_mode1 = 1.56
+sigma_mode2 = 1.29
+#
+# Aerosol production (precursors off the droplets)
+prod_aer = y
+#Characteristics of the size distribution of the aer
+rad_aer = 0.125e-6
+sigma_aer = 1.56
+rho_aer = 1800
+# Vertical distribution and production rate
+alt_aer = 62000
+sigmaz_aer = 2000
+tx_prod_aer = 1.D-10
 #
 # parameters for the boundary layer

trunk/LMDZ.VENUS/deftank/physiq.def

@@ -65 +65 @@
 ok_sedim=n
 nb_mode=3
+#
+# parameters for microphysics (cl_scheme = 2)
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# Verboose debug mode
+debug_cloud = n
+#Keys to calls the different microphysics processes or not
+mode_merging = y
+hom_nuclea = y
+het_nuclea = y
+cond_evap = y
+coag = y
+#
+# Regarding the mode 1 and 2 size distributions
+sigma_mode1 = 1.56
+sigma_mode2 = 1.29
+#
+# Aerosol production (precursors off the droplets)
+prod_aer = y
+#Characteristics of the size distribution of the aer
+rad_aer = 0.125e-6
+sigma_aer = 1.56
+rho_aer = 1800
+# Vertical distribution and production rate
+alt_aer = 62000
+sigmaz_aer = 2000
+tx_prod_aer = 1.D-10
 #
 # parameters for the boundary layer

trunk/LMDZ.VENUS/deftank/traceur-chemcld.def

@@ -1 @@
-48
+55
 10 10 co2
 10 10 co
@@ -36 +36 @@
 10 10 oscl
 10 10 n2
-10 10 M0_aer
-10 10 M3_aer
-10 10 M0_m1drop
-10 10 M0_m1ccn
-10 10 M3_m1sa
-10 10 M3_m1w
-10 10 M3_m1ccn
-10 10 M0_m2drop
-10 10 M0_m2ccn
-10 10 M3_m2sa
-10 10 M3_m2w
-10 10 M3_m2ccn
+10 10 he
+10 10 no2
+10 10 no
+10 10 n
+10 10 n2d
+10 10 h2oliq
+10 10 h2so4liq
+10 10 m0_aer
+10 10 m3_aer
+10 10 m0_mode1drop
+10 10 m0_mode1ccn
+10 10 m3_mode1sa
+10 10 m3_mode1wv
+10 10 m3_mode1ccn
+10 10 m0_mode2drop
+10 10 m0_mode2ccn
+10 10 m3_mode2sa
+10 10 m3_mode2wv
+10 10 m3_mode2ccn