Add bug fix to make sure that arrays on level edges are not going out…

… of bounds (issue GEOS-ESM/GMAO_Shared#99)
mathomp4 · Feb 16, 2021 · 5a5247b · 5a5247b
1 parent f4d8957
commit 5a5247b
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/Lightning_mod.F90 b/Lightning_mod.F90
@@ -913,8 +913,8 @@ subroutine hemcoFlashrate (cellArea, lwi, lonslocal, latslocal, airTemp, ple, ge
       ALLOCATE(ltop2d (1:IM, 1:JM),   STAT=RC)
 
       ! edge layer fields
-      pleHemco(:,:,1:KM) = ple(:,:,KM:1:-1)
-      cnvMfcHemco(:,:,1:KM) = cnvMfc(:,:,KM:1:-1)
+      pleHemco(:,:,1:KM) = ple(:,:,LM:0:-1)
+      cnvMfcHemco(:,:,1:KM) = cnvMfc(:,:,LM:0:-1)
 
       ! fields at grid box center
       airTempHemco(:,:,1:LM) = airTemp(:,:,LM:1:-1)
@@ -925,7 +925,7 @@ subroutine hemcoFlashrate (cellArea, lwi, lonslocal, latslocal, airTemp, ple, ge
 
 
       do levCount=1,LM
-         gridBoxHeightHemco(:,:,levCount) = geoPotHeight(:,:,levCount) - geoPotHeight(:,:,levCount+1)
+         gridBoxHeightHemco(:,:,levCount) = geoPotHeight(:,:,levCount-1) - geoPotHeight(:,:,levCount)
       enddo
       gridBoxHeightHemco2(:,:,1:LM) = gridBoxHeightHemco(:,:,LM:1:-1)