Get global X and Y dimensions from partition metadata file

They are needed to write to a correct location in parallel mode.
nextsimhub · Nov 9, 2023 · 2ff91ea · 2ff91ea
1 parent 791c886
commit 2ff91ea
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/core/src/RectGridIO.cpp b/core/src/RectGridIO.cpp
@@ -161,6 +161,8 @@ void RectGridIO::readPartitionData(const std::string& partitionFile, ModelMetada
             + std::to_string(nBoxes) + "\n";
         throw std::runtime_error(errorMsg);
     }
+    metadata.globalExtentX = ncFile.getDim("globalX").getSize();
+    metadata.globalExtentY = ncFile.getDim("globalY").getSize();
     netCDF::NcGroup bboxGroup(ncFile.getGroup(bboxName));
     std::vector<size_t> index(1, metadata.mpiMyRank);
     bboxGroup.getVar("global_x").getVar(index, &metadata.localCornerX);

diff --git a/core/src/include/ModelMetadata.hpp b/core/src/include/ModelMetadata.hpp
@@ -61,7 +61,7 @@ class ModelMetadata {
     MPI_Comm mpiComm;
     int mpiSize = 0;
     int mpiMyRank = -1;
-    int localCornerX, localCornerY, localExtentX, localExtentY;
+    int localCornerX, localCornerY, localExtentX, localExtentY, globalExtentX, globalExtentY;
 #endif
 
 private: