From e93eb34614d4e3a2671fa10e146882f01e13d803 Mon Sep 17 00:00:00 2001 From: Atanas Trayanov Date: Tue, 14 May 2024 09:28:18 -0400 Subject: [PATCH] Fixed uninitialed tep and qp when wstar <= 0.001 --- .../GEOSturbulence_GridComp/LockEntrain.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSturbulence_GridComp/LockEntrain.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSturbulence_GridComp/LockEntrain.F90 index 5d87dd086..1f079aaa0 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSturbulence_GridComp/LockEntrain.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSturbulence_GridComp/LockEntrain.F90 @@ -1265,7 +1265,8 @@ subroutine mpbl_depth(i,j,icol,jcol,nlev,tpfac, entrate, pceff, vscale, pertopt, qp = q(i,j,nlev) + 2.*evap(i,j)/(zrho*wstar) ! print *,'tpert=',2.*sh(i,j)/(zrho*wstar*MAPL_CP) else - + tep = t(i,j,nlev) + 0.4 + qp = q(i,j,nlev) end if else ! tpfac scales up bstar by inv. ratio of ! heat-bubble area to stagnant area