Ensure tep and qp are always initialized in mpbl_depth

This should fix Issue #930
GEOS-ESM · May 13, 2024 · ed35b2e · ed35b2e
1 parent db55c30
commit ed35b2e
Showing 1 changed file with 2 additions and 1 deletion.
diff --git 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