Changeset 4141 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 18, 2026, 10:58:43 PM (8 days ago)

Author:

yluo

Message:

Mars PCM:
Enabled output of vmr_<tracer> for all chemical tracers, including gas-phase species, charged species, and H2O ice.
YCL

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (1 diff)

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

@@ -3138 +3138 @@
                     call wstats(ngrid,"vmr_"//trim(noms(iq)),
      $                        "Volume mixing ratio","mol/mol",3,vmr)
-                    if ((noms(iq).eq."o")
-     $             .or. (noms(iq).eq."co2")
+                    if ((noms(iq).eq."co2")
+     $             .or. (noms(iq).eq."co")
+     $             .or. (noms(iq).eq."o")
+     $             .or. (noms(iq).eq."o1d")
+     $             .or. (noms(iq).eq."o2")
      $             .or. (noms(iq).eq."o3")
+     $             .or. (noms(iq).eq."h")
+     $             .or. (noms(iq).eq."h2")
+     $             .or. (noms(iq).eq."oh")
+     $             .or. (noms(iq).eq."ho2")
+     $             .or. (noms(iq).eq."h2o2")
+     $             .or. (noms(iq).eq."n2")
      $             .or. (noms(iq).eq."ar")
-     $             .or. (noms(iq).eq."o2")
-     $             .or. (noms(iq).eq."h2o_vap") ) then
+     $             .or. (noms(iq).eq."h2o_ice")
+     $             .or. (noms(iq).eq."h2o_vap")
+     $             .or. (noms(iq).eq."n")
+     $             .or. (noms(iq).eq."no")
+     $             .or. (noms(iq).eq."no2")
+     $             .or. (noms(iq).eq."n2d")
+     $             .or. (noms(iq).eq."he")
+     $             .or. (noms(iq).eq."co2plus")
+     $             .or. (noms(iq).eq."oplus")
+     $             .or. (noms(iq).eq."o2plus")
+     $             .or. (noms(iq).eq."coplus")
+     $             .or. (noms(iq).eq."cplus")
+     $             .or. (noms(iq).eq."nplus")
+     $             .or. (noms(iq).eq."noplus")
+     $             .or. (noms(iq).eq."n2plus")
+     $             .or. (noms(iq).eq."hplus")
+     $             .or. (noms(iq).eq."hco2plus")
+     $             .or. (noms(iq).eq."hcoplus")
+     $             .or. (noms(iq).eq."h2oplus")
+     $             .or. (noms(iq).eq."h3oplus")
+     $             .or. (noms(iq).eq."ohplus")
+     $             .or. (noms(iq).eq."elec") ) then
                       call write_output("vmr_"//trim(noms(iq)),
      $                         "Volume mixing ratio","mol/mol",vmr(:,:))