Change saturation

ai2cm · Mar 5, 2021 · 7b39eb3 · 7b39eb3
1 parent fe24ba8
commit 7b39eb3
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/FV3/atmos_cubed_sphere/model/fv_cmp.F90 b/FV3/atmos_cubed_sphere/model/fv_cmp.F90
@@ -406,7 +406,7 @@ subroutine fv_sat_adj (mdt, zvir, is, ie, js, je, ng, hydrostatic, consv_te, te0
                 ! factor = - min (1., fac_l2v * sqrt (max (0., ql (i, j)) / 1.e-5) * 10. * (1. - qv (i, j) / wqsat (i)))
                 ! factor = - fac_l2v
                 ! factor = - 1
-                factor = - min (1., fac_l2v * 10. * (1. - qv (i, j) / wqsat (i))) ! the rh dependent factor = 1 at 90%
+                factor = - fac_l2v * min (1., 10. * (1. - qv (i, j) / wqsat (i))) ! the rh dependent factor = 1 at 90%
                 src (i) = - min (ql (i, j), factor * dq0)
             endif
             qv (i, j) = qv (i, j) - src (i)
@@ -449,7 +449,7 @@ subroutine fv_sat_adj (mdt, zvir, is, ie, js, je, ng, hydrostatic, consv_te, te0
                     ! factor = - min (1., fac_l2v * sqrt (max (0., ql (i, j)) / 1.e-5) * 10. * (1. - qv (i, j) / wqsat (i))) ! the rh dependent factor = 1 at 90%
                     ! factor = - fac_l2v
                     ! factor = - 1
-                    factor = - min (1., fac_l2v * 10. * (1. - qv (i, j) / wqsat (i))) ! the rh dependent factor = 1 at 90%
+                    factor = - fac_l2v * min (1., 10. * (1. - qv (i, j) / wqsat (i))) ! the rh dependent factor = 1 at 90%
                     src (i) = - min (ql (i, j), factor * dq0)
                 endif
                 adj_fac = 1.

diff --git a/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 b/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90
@@ -2083,7 +2083,7 @@ subroutine subgrid_z_proc (ktop, kbot, p1, den, denfac, dts, rh_adj, tz, qv, &
             ! factor = min (1., fac_l2v * sqrt (max (0., ql (k)) / 1.e-5) * 10. * dq0 / qsw)
             ! factor = fac_l2v
             ! factor = 1
-            factor = min (1., fac_l2v * (10. * dq0 / qsw)) ! the rh dependent factor = 1 at 90%
+            factor = fac_l2v * min (1., 10. * dq0 / qsw) ! the rh dependent factor = 1 at 90%
             evap = min (ql (k), factor * dq0 / (1. + tcp3 (k) * dwsdt))
         else ! condensate all excess vapor into cloud water
             ! -----------------------------------------------------------------------