From 87fe79af20e62e3503b21b47a12c3bc452b06fca Mon Sep 17 00:00:00 2001 From: "Noah D. Brenowitz" Date: Wed, 12 May 2021 16:37:24 -0700 Subject: [PATCH] Revert "revert rain evaporation factor change" This reverts commit 8eae74103880ecc723357ecf3d5038ab82c1cc5c. Without this change the rain water becomes far too low. --- FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 b/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 index d8ef9df3d..68831e8ec 100644 --- a/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 +++ b/FV3/gfsphysics/physics/gfdl_cloud_microphys.F90 @@ -1390,7 +1390,8 @@ subroutine revap_racc (ktop, kbot, dt, tz, qv, ql, qr, qi, qs, qg, den, denfac, t2 = tin * tin evap = crevp (1) * t2 * dq * (crevp (2) * sqrt (qden) + crevp (3) * & exp (0.725 * log (qden))) / (crevp (4) * t2 + crevp (5) * qsat * den (k)) - evap = min (qr (k), dt * evap, dqv / (1. + lcpk (k) * dqsdt)) + factor = min (1., re_factor * dqv / qsat) + evap = min (qr (k), dt * evap, factor * dqv / (1. + lcpk (k) * dqsdt)) ! ----------------------------------------------------------------------- ! alternative minimum evap in dry environmental air ! sink = min (qr (k), dim (rh_rain * qsat, qv (k)) / (1. + lcpk (k) * dqsdt))