From 4f8f61501ff9ec315557718e851a4c81d6a373ea Mon Sep 17 00:00:00 2001 From: JungHoon Shin Date: Mon, 29 Jul 2024 15:58:55 +0000 Subject: [PATCH] JungHoon: Increasing the absolute value of entrainment/detrainment coefficients. Simply changed 1.75e-3 to 2.50e-3 in line 2303 (entrainment part) and -1.75e-3 to -2.50e-3 in line 2179 (detrainment part) of the code --- physics/CONV/nTiedtke/cu_ntiedtke.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/physics/CONV/nTiedtke/cu_ntiedtke.F90 b/physics/CONV/nTiedtke/cu_ntiedtke.F90 index 06246457c..e2f4730c0 100644 --- a/physics/CONV/nTiedtke/cu_ntiedtke.F90 +++ b/physics/CONV/nTiedtke/cu_ntiedtke.F90 @@ -2176,7 +2176,7 @@ subroutine cuascn & jl = jlx(jll) zdmfde(jl) = min(zdmfde(jl),0.75*pmfu(jl,jk+1)) if ( jk == kcbot(jl) ) then - zoentr(jl) = -1.75e-3*(min(1.,pqen(jl,jk)/pqsen(jl,jk)) - & + zoentr(jl) = -2.50e-3*(min(1.,pqen(jl,jk)/pqsen(jl,jk)) - & 1.)*(pgeoh(jl,jk)-pgeoh(jl,jk+1))*zrg zoentr(jl) = min(0.4,zoentr(jl))*pmfu(jl,jk+1) end if @@ -2300,7 +2300,7 @@ subroutine cuascn & end if if ( zbuo(jl,jk) > -0.2 ) then ikb = kcbot(jl) - zoentr(jl) = 1.75e-3*(0.3-(min(1.,pqen(jl,jk-1) / & + zoentr(jl) = 2.50e-3*(0.3-(min(1.,pqen(jl,jk-1) / & pqsen(jl,jk-1))-1.))*(pgeoh(jl,jk-1)-pgeoh(jl,jk)) * & zrg*min(1.,pqsen(jl,jk)/pqsen(jl,ikb))**3 zoentr(jl) = min(0.4,zoentr(jl))*pmfu(jl,jk)