For ORIG mechanism: Set reaction rates to zero, for ICE and NAT and S…

…TS reactions, when HCl < min-concentration
GEOS-ESM · Apr 25, 2024 · c5a22ea · c5a22ea
1 parent 10555ee
commit c5a22ea
Showing 1 changed file with 11 additions and 11 deletions.
diff --git a/GMI_GridComp/GmiChemistry/StratTrop_Orig/setkin_kcalc.F90 b/GMI_GridComp/GmiChemistry/StratTrop_Orig/setkin_kcalc.F90
@@ -2384,7 +2384,7 @@ FUNCTION sksts_clono2_hcl (tk,ad,pr,sad,clono2,hcl,h2o,ptrop)
           where( hcl > minconc )
             sksts_clono2_hcl = 0.25d0 * gprob_hcl * avgvel * sad / hcl
           elsewhere
-            sksts_clono2_hcl = 0.25d0 * gprob_hcl * avgvel * sad
+            sksts_clono2_hcl = 0.0d0
           end where
 !
           where( sad < 0.0d0 ) sksts_clono2_hcl = 0.0d0
@@ -2540,7 +2540,7 @@ FUNCTION sksts_hocl_hcl (tk,ad,pr,sad,hocl,hcl,h2o,ptrop)
           where( hcl > minconc )
             sksts_hocl_hcl = 0.25d0 * gprob_tot * avgvel * sad / hcl
           elsewhere
-            sksts_hocl_hcl = 0.25d0 * gprob_tot * avgvel * sad
+            sksts_hocl_hcl = 0.0d0
           end where
 !
           where( sad < 0.0d0 ) sksts_hocl_hcl   = 0.0d0
@@ -2654,7 +2654,7 @@ FUNCTION sksts_hobr_hcl (tk,ad,pr,sad,hobr,hcl,h2o,ptrop)
           where( hcl > minconc )
             sksts_hobr_hcl = 0.25d0 * gprob_tot * avgvel * sad / hcl
          elsewhere
-            sksts_hobr_hcl = 0.25d0 * gprob_tot * avgvel * sad
+            sksts_hobr_hcl = 0.0d0
          end where
 !
           where( sad < 0.0d0 ) sksts_hobr_hcl = 0.0d0
@@ -2766,7 +2766,7 @@ FUNCTION sknat_hcl_clono2 (tk,pr,sad,hcl,ptrop)
           sknat_hcl_clono2(:) = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            sknat_hcl_clono2  = sknat_hcl_clono2 / minconc
+            sknat_hcl_clono2  = 0.0d0
           elsewhere
             sknat_hcl_clono2  = sknat_hcl_clono2 / hcl
           end where
@@ -2807,7 +2807,7 @@ FUNCTION sknat_hcl_hocl (tk,pr,sad,hcl,ptrop)
           sknat_hcl_hocl(:) = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            sknat_hcl_hocl  = sknat_hcl_hocl / minconc
+            sknat_hcl_hocl  = 0.0d0
           elsewhere
             sknat_hcl_hocl  = sknat_hcl_hocl / hcl
           end where
@@ -2848,7 +2848,7 @@ FUNCTION sknat_hcl_brono2 (tk,pr,sad,hcl,ptrop)
           sknat_hcl_brono2(:) = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            sknat_hcl_brono2  = sknat_hcl_brono2 / minconc
+            sknat_hcl_brono2  = 0.0d0
           elsewhere
             sknat_hcl_brono2  = sknat_hcl_brono2 / hcl
           end where
@@ -2889,7 +2889,7 @@ FUNCTION sknat_hcl_hobr (tk,pr,sad,hcl,ptrop)
           sknat_hcl_hobr(:) = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            sknat_hcl_hobr  = sknat_hcl_hobr / minconc
+            sknat_hcl_hobr  = 0.0d0
           elsewhere
             sknat_hcl_hobr  = sknat_hcl_hobr / hcl
           end where
@@ -3000,7 +3000,7 @@ FUNCTION skice_hcl_clono2 (tk,pr,sad,hcl,ptrop)
           skice_hcl_clono2(:)  = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            skice_hcl_clono2   = skice_hcl_clono2 / minconc
+            skice_hcl_clono2   = 0.0d0
           elsewhere
             skice_hcl_clono2   = skice_hcl_clono2 / hcl
           end where
@@ -3043,7 +3043,7 @@ FUNCTION skice_hcl_hocl (tk,pr,sad,hcl,ptrop)
           skice_hcl_hocl(:)  = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            skice_hcl_hocl   = skice_hcl_hocl / minconc
+            skice_hcl_hocl   = 0.0d0
           elsewhere
             skice_hcl_hocl   = skice_hcl_hocl / hcl
           end where
@@ -3087,7 +3087,7 @@ FUNCTION skice_hcl_brono2 (tk,pr,sad,hcl,ptrop)
           skice_hcl_brono2(:)  = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            skice_hcl_brono2   = skice_hcl_brono2 / minconc
+            skice_hcl_brono2   = 0.0d0
           elsewhere
             skice_hcl_brono2   = skice_hcl_brono2 / hcl
           end where
@@ -3129,7 +3129,7 @@ FUNCTION skice_hcl_hobr (tk,pr,sad,hcl,ptrop)
           skice_hcl_hobr(:)  = 0.25d0 * gprob * avgvel(:) * sad(:)
 !
           where (hcl < minconc)
-            skice_hcl_hobr   = skice_hcl_hobr / minconc
+            skice_hcl_hobr   = 0.0d0
           elsewhere
             skice_hcl_hobr   = skice_hcl_hobr / hcl
           end where