Debug grid extraction for 3D and 4D

GeoRegionsEcosystem · Jun 22, 2023 · 037e184 · 037e184 · natgeo-wong · Jun 22, 2023
1 parent dacf22a
commit 037e184
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "GeoRegions"
 uuid = "b001f823-fa75-4bff-bf55-6610c8f3688a"
 authors = ["natgeo-wong <[email protected]>"]
-version = "5.2.2"
+version = "5.2.3"
 
 [deps]
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"

diff --git a/src/extract/grid.jl b/src/extract/grid.jl
@@ -143,7 +143,7 @@ function extractGrid(
     ilat  = ggrd.ilat; nlat = length(ggrd.ilat)
                        n3D  = size(odata,3)
     mask  = ggrd.mask
-	ndata = zeros(nlon,nlat)
+	ndata = zeros(nlon,nlat,n3D)
 	for i3D = 1 : n3D, glat in 1 : nlat, glon in 1 : nlon
 		ndata[glon,glat,i3D] = odata[ilon[glon],ilat[glat],i3D] * mask[glon,glat]
 	end
@@ -179,7 +179,7 @@ function extractGrid(
 	ilat  = ggrd.ilat; nlat = length(ggrd.ilat)
                        n3D  = size(odata,3)
                        n4D  = size(odata,4)
-	ndata = zeros(nlon,nlat,n3D)
+	ndata = zeros(nlon,nlat,n3D,n4D)
 	for i4D = 1 : n4D, i3D = 1 : n3D, glat in 1 : nlat, glon in 1 : nlon
 		ndata[glon,glat,i3D,i4D] = odata[ilon[glon],ilat[glat],i3D,i4D]
 	end
@@ -216,7 +216,7 @@ function extractGrid(
                        n3D  = size(odata,3)
                        n4D  = size(odata,4)
     mask  = ggrd.mask
-	ndata = zeros(nlon,nlat)
+	ndata = zeros(nlon,nlat,n3D,n4D)
 	for i4D = 1 : n4D, i3D = 1 : n3D, glat in 1 : nlat, glon in 1 : nlon
 		ndata[glon,glat,i3D,i4D] = odata[ilon[glon],ilat[glat],i3D,i4D] * mask[glon,glat]
 	end