Changeset 3799 for trunk/LMDZ.MARS

Timestamp:

Jun 5, 2025, 6:10:58 PM (13 days ago)

Author:

yluo

Message:

Mars PCM:
Include H2O2 in the mean molecular weight and specific heat capacity calculations.
YCL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/conc_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4852 +4852 @@
 == 05/06/2025 == JBC
 Handle correctly more variables with different types/shapes/dimensions.
+
+== 05/06/2025 == YCL
+Include H2O2 in the mean molecular weight and specific heat capacity calculations.

trunk/LMDZ.MARS/libf/phymars/conc_mod.F90

@@ -67 +67 @@
  use tracer_mod, only: igcm_co2, igcm_co, igcm_o, igcm_o1d,             &
 &                      igcm_o2, igcm_o3, igcm_h, igcm_h2,               &
-&                      igcm_oh, igcm_ho2, igcm_n2, igcm_ar,             &
+&                      igcm_oh, igcm_ho2, igcm_h2o2, igcm_n2, igcm_ar,  &
 &                      igcm_h2o_vap, igcm_n, igcm_no, igcm_no2, &
 &                      igcm_n2d, igcm_co2plus, igcm_oplus,              &
@@ -185 +185 @@
           cpi(nbq) = 1.065e3  !?
        end if
+       if (igcm_h2o2 /= 0) then
+          nbq = nbq + 1
+          niq(nbq) = igcm_h2o2
+          aki(nbq) = 0.0      !?
+          cpi(nbq) = 1.255e3
+       end if
        if (igcm_n2 /= 0) then
           nbq = nbq + 1