Merge pull request #34 from jimmielin/hplin/geoschem_14.3_fixes

Fixes for GEOS-Chem v14.3
geoschem · Jun 24, 2024 · 2c0f0bd · 2c0f0bd
2 parents 4126267 + e53afbe
commit 2c0f0bd
Show file tree

Hide file tree

Showing 3 changed files with 266 additions and 192 deletions.
diff --git a/src/chemistry/geoschem/chem_mods.F90 b/src/chemistry/geoschem/chem_mods.F90
@@ -62,7 +62,7 @@ module chem_mods
                             rxntot = 212, & ! number of total reactions
                             gascnt = 172, & ! number of gas phase reactions
                             nabscol = 2, & ! number of absorbing column densities
-                            gas_pcnst = 271, & ! number of "gas phase" species (=solsym length)
+                            gas_pcnst = 357, & ! number of "gas phase" species (=solsym length)
                                                ! Includes GC advected species, MAM aerosols,
                                                ! and CO2, and any non-advected species added
                                                ! to solsym within mo_sim_dat.F90.
@@ -81,7 +81,7 @@ module chem_mods
                             clsze = 1, & ! loop length for implicit chemistry
                             rxt_tag_cnt = 0, & ! number of tagged reactions (unused in GEOS-Chem)
                             enthalpy_cnt = 0, &
-                            nslvd = 88  ! number of short-lived (non-advected) species
+                            nslvd = 88  ! number of short-lived (non-advected) species minus CO2
       integer :: clscnt(5) = 0
       integer :: cls_rxt_cnt(4,5) = 0
       integer :: clsmap(gas_pcnst,5) = 0

diff --git a/src/chemistry/geoschem/chemistry.F90 b/src/chemistry/geoschem/chemistry.F90
@@ -1907,7 +1907,7 @@ subroutine chem_timestep_tend( state, ptend, cam_in, cam_out, dT, pbuf, fh2o )
     use Toms_Mod,            only : Compute_Overhead_O3
     use UCX_Mod,             only : Set_H2O_Trac
     use Unitconv_Mod,        only : Convert_Spc_Units, UNIT_STR
-    use Unitconv_Mod,        only : KG_SPECIES_PER_KG_DRY_AIR
+    use Unitconv_Mod,        only : KG_SPECIES_PER_KG_DRY_AIR, KG_SPECIES_PER_M2
     use Wetscav_Mod,         only : Setup_Wetscav
 
     REAL(r8),            INTENT(IN)    :: dT          ! Time step
@@ -3744,7 +3744,7 @@ subroutine chem_timestep_tend( state, ptend, cam_in, cam_out, dT, pbuf, fh2o )
                             State_Chm  = State_Chm(LCHNK),            &
                             State_Grid = State_Grid(LCHNK),           &
                             State_Met  = State_Met(LCHNK),            &
-                            OutUnit    = KG_SPECIES_PER_KG_DRY_AIR,   &
+                            OutUnit    = KG_SPECIES_PER_M2,           &
                             OrigUnit   = OrigUnit,                    &
                             RC         = RC                          )