From ed35b2ea37f5579a20e9d816e0baf60163b9526c Mon Sep 17 00:00:00 2001
From: Nathan Arnold <53791641+narnold1@users.noreply.github.com>
Date: Mon, 13 May 2024 18:59:02 -0400
Subject: [PATCH] Ensure tep and qp are always initialized in mpbl_depth

This should fix Issue #930
---
 .../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..5becd6733 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