From f6d0b3929a306cecfd929a3c80159a4d89aa3b08 Mon Sep 17 00:00:00 2001 From: bluefinweiwei Date: Wed, 11 Sep 2024 20:36:42 +0000 Subject: [PATCH] Bugfix Precip diagnostics GFDL MP developer (Linjiong Zhou) provided the code changes. The fix solved the problem of overly produced precip and sensitivity to microphysics timestep modified: module_gfdl_cld_mp.F90 --- .../MP/GFDL_2022_v3/module_gfdl_cld_mp.F90 | 19 +++++++------------ 1 file changed, 7 insertions(+), 12 deletions(-) diff --git a/physics/MP/GFDL_2022_v3/module_gfdl_cld_mp.F90 b/physics/MP/GFDL_2022_v3/module_gfdl_cld_mp.F90 index 41defcd3e..20fb714f5 100644 --- a/physics/MP/GFDL_2022_v3/module_gfdl_cld_mp.F90 +++ b/physics/MP/GFDL_2022_v3/module_gfdl_cld_mp.F90 @@ -1267,7 +1267,7 @@ subroutine mpdrv (hydrostatic, ua, va, wa, delp, pt, qv, ql, qr, qi, qs, qg, & ! unit convert to mm/day ! ----------------------------------------------------------------------- - convt = 86400. * rgrav / dts + convt = 86400. * rgrav / dtm do i = is, ie @@ -1433,7 +1433,7 @@ subroutine mpdrv (hydrostatic, ua, va, wa, delp, pt, qv, ql, qr, qi, qs, qg, & if (fix_negative) & call neg_adj (ks, ke, tz, dp, qvz, qlz, qrz, qiz, qsz, qgz, cond) - condensation (i) = condensation (i) + cond * convt * ntimes + condensation (i) = condensation (i) + cond * convt ! ----------------------------------------------------------------------- ! fast microphysics loop @@ -1946,16 +1946,11 @@ subroutine mp_full (ks, ke, ntimes, tz, qv, ql, qr, qi, qs, qg, dp, dz, u, v, w, snow = snow + s1 * convt graupel = graupel + g1 * convt - !prefluxw = prefluxw + pfw * convt - !prefluxr = prefluxr + pfr * convt - !prefluxi = prefluxi + pfi * convt - !prefluxs = prefluxs + pfs * convt - !prefluxg = prefluxg + pfg * convt - prefluxw = prefluxw + pfw - prefluxr = prefluxr + pfr - prefluxi = prefluxi + pfi - prefluxs = prefluxs + pfs - prefluxg = prefluxg + pfg + prefluxw = prefluxw + pfw * convt + prefluxr = prefluxr + pfr * convt + prefluxi = prefluxi + pfi * convt + prefluxs = prefluxs + pfs * convt + prefluxg = prefluxg + pfg * convt ! ----------------------------------------------------------------------- ! warm rain cloud microphysics