From 45ed3770b3a390a5285499fd727c46cee4c1e53a Mon Sep 17 00:00:00 2001 From: "Chunxi.Zhang-NOAA" Date: Wed, 25 Jan 2023 21:54:32 +0000 Subject: [PATCH] Bug fix for cloud effective radius for convective clouds (HR1) --- physics/radiation_clouds.f | 14 ++++++++++---- 1 file changed, 10 insertions(+), 4 deletions(-) diff --git a/physics/radiation_clouds.f b/physics/radiation_clouds.f index bf255ce00..7255f1578 100644 --- a/physics/radiation_clouds.f +++ b/physics/radiation_clouds.f @@ -2212,10 +2212,16 @@ subroutine progcld_thompson_wsm6 & !> The total condensate includes convective condensate. do k = 1, NLAY-1 do i = 1, IX - cwp(i,k) = max(0.0, (clw(i,k,ntcw)+cnvw(i,k)* - & (1.-tem2d(i,k))) * gfac * delp(i,k)) - cip(i,k) = max(0.0, (clw(i,k,ntiw) + cnvw(i,k)* - & tem2d(i,k)) *gfac * delp(i,k)) + tem1 = cnvw(i,k)*(1.-tem2d(i,k)) + cwp(i,k) = max(0.0, (clw(i,k,ntcw)+tem1) * + & gfac * delp(i,k)) + if(tem1 > 1.e-12 .and. clw(i,k,ntcw) < 1.e-12) + & rew(i,k)=reliq_def + tem2 = cnvw(i,k)*tem2d(i,k) + cip(i,k) = max(0.0, (clw(i,k,ntiw) + tem2 ) + & *gfac * delp(i,k)) + if(tem2 > 1.e-12 .and. clw(i,k,ntiw) < 1.e-12) + & rei(i,k)=reice_def crp(i,k) = max(0.0, clw(i,k,ntrw) * gfac * delp(i,k)) csp(i,k) = max(0.0, clw(i,k,ntsw) * gfac * delp(i,k)) enddo