Debug grid extraction for 3D and 4D

GeoRegionsEcosystem · Jun 22, 2023 · 7c3ae0f · 7c3ae0f · natgeo-wong · Jun 22, 2023
1 parent 99e2679
commit 7c3ae0f
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.1.2"
+version = "5.1.3"
 
 [deps]
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"

diff --git a/src/ExtractGrid.jl b/src/ExtractGrid.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