From 5010ceda52091d14ebb9c38b8b023fafdb8d7a98 Mon Sep 17 00:00:00 2001 From: Atanas Trayanov Date: Wed, 11 Feb 2026 15:23:00 -0500 Subject: [PATCH] Initialized iTable and rTable after allocation. Increased the initial size of hash to speed-up the code --- .../Utils/Raster/makebcs/CombineRasters.F90 | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Utils/Raster/makebcs/CombineRasters.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Utils/Raster/makebcs/CombineRasters.F90 index 001bce06c..8d2bc31b6 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Utils/Raster/makebcs/CombineRasters.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSsurface_GridComp/Utils/Raster/makebcs/CombineRasters.F90 @@ -8,6 +8,7 @@ program mkOverlaySimple use MAPL_HashMod use MAPL_ExceptionHandling use MAPL_Constants + use, intrinsic :: iso_fortran_env, only: REAL64 ! Overlay atmosphere, land, and ocean rasters, creating a .idx file. ! The ocean raster should be defined everywhere, or at least, everywhere @@ -210,8 +211,11 @@ program mkOverlaySimple allocate(iTable(0:nvars,maxtiles),stat=status) VERIFY_(STATUS) + iTable = 0 + allocate(rTable(1:rvars,maxtiles),stat=status) VERIFY_(STATUS) + rTable = 0.0_REAL64 allocate(rst1(nx,ny), stat=status) VERIFY_(STATUS) @@ -279,7 +283,7 @@ program mkOverlaySimple end if ip = 0 - Hash = MAPL_HashCreate(8*1024) + Hash = MAPL_HashCreate(512*1024) if(Verb) write (6, '(A)', advance='NO') ' Started Overlay'