From c5a22ead35a0a5026bbbbb8578ab593372bf6b70 Mon Sep 17 00:00:00 2001
From: Mike Manyin <michael.manyin@nasa.gov>
Date: Thu, 25 Apr 2024 14:58:47 -0400
Subject: [PATCH] For ORIG mechanism: Set reaction rates to zero, for ICE and
 NAT and STS reactions, when HCl < min-concentration

---
 .../StratTrop_Orig/setkin_kcalc.F90           | 22 +++++++++----------
 1 file changed, 11 insertions(+), 11 deletions(-)

diff --git a/GMI_GridComp/GmiChemistry/StratTrop_Orig/setkin_kcalc.F90 b/GMI_GridComp/GmiChemistry/StratTrop_Orig/setkin_kcalc.F90
index 701cad0..0a98b07 100644
--- 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