In routine read_grid_file, initialize each array after allocation.

Fixes ufs-community#105.
GeorgeGayno-NOAA · Aug 6, 2021 · ccc6692 · ccc6692
1 parent dbe2ba5
commit ccc6692
Showing 1 changed file with 12 additions and 2 deletions.
diff --git a/sorc/grid_tools.fd/filter_topo.fd/filter_topo.F90 b/sorc/grid_tools.fd/filter_topo.fd/filter_topo.F90
@@ -705,8 +705,10 @@ subroutine read_grid_file(regional)
           jsd=js-ng; jed=je+ng
 
           allocate(tmpvar(ni+1,nj+1))
-          allocate(geolon_c(isd:ied+1,jsd:jed+1,6))
-          allocate(geolat_c(isd:ied+1,jsd:jed+1,6))
+          allocate(geolon_c(isd:ied+1,jsd:jed+1,ntiles))
+          allocate(geolat_c(isd:ied+1,jsd:jed+1,ntiles))
+          geolon_c(:,:,:) = 0.
+          geolat_c(:,:,:) = 0.
        else if ( t == 7 ) then ! nested grid
           if(mod(ni,2) .NE. 0 .or. mod(nj,2) .NE. 0) &
                call handle_err(-1, "read_grid_file: ni and nj must be even")
@@ -798,6 +800,8 @@ subroutine read_grid_file(regional)
     !--- compute dx, dy
     allocate(dx(isd:ied,jsd:jed+1,ntiles))
     allocate(dy(isd:ied+1,jsd:jed,ntiles))
+    dx(:,:,:) = -1.e25
+    dy(:,:,:) = -1.e25
     do t = 1, ntiles
        do j = js, je+1 ; do i = is, ie
           g1(1) = geolon_c(i  ,j,t)
@@ -846,6 +850,8 @@ subroutine read_grid_file(regional)
     !--- compute dxa and dya -----
     allocate(dxa(isd:ied,jsd:jed,ntiles))
     allocate(dya(isd:ied,jsd:jed,ntiles))
+    dxa(:,:,:) = -1.e25
+    dya(:,:,:) = -1.e25
     do t = 1, ntiles
        do j=js,je ; do i=is,ie
           g1(1) = geolon_c(i,j,t); g1(2) = geolat_c(i,j,t)
@@ -875,6 +881,8 @@ subroutine read_grid_file(regional)
     !--- compute dxc and dyc
     allocate(dxc(isd:ied+1,jsd:jed,ntiles))
     allocate(dyc(isd:ied,jsd:jed+1,ntiles))
+    dxc(:,:,:) = -1.e25
+    dyc(:,:,:) = -1.e25
     do t = 1, ntiles
        do j=jsd,jed
           do i=isd+1,ied
@@ -901,6 +909,7 @@ subroutine read_grid_file(regional)
 
     !--- compute area
     allocate(area(isd:ied,jsd:jed,ntiles))
+    area(:,:,:) = 0.
     do t = 1, ntiles
        do j=js,je
           do i=is,ie
@@ -925,6 +934,7 @@ subroutine read_grid_file(regional)
     allocate(sin_sg(4,isd:ied,jsd:jed,ntiles))
     allocate(cos_sg(4,isd:ied,jsd:jed))
     allocate(grid3(3, npx, npy))
+    grid3(:,:,:) = -1.e25
     cos_sg(:,:,:) =  big_number
     sin_sg(:,:,:,:) = tiny_number