Fix dimension reading for the MPI version.

Author: timspainNERSC

core/src/ParaGridIO.cpp

@@ -75,7 +75,7 @@ bool ParaGridIO::doOnce()
 
 ParaGridIO::~ParaGridIO() = default;
 
-template <typename N> void ParaGridIO::readDimensions(const N& node)
+template <typename N> void ParaGridIO::readDimensions(const N& node, ModelMetadata& metadata)
 {
     for (auto entry : ModelArray::definedDimensions) {
         auto dimType = entry.first;
@@ -140,11 +140,11 @@ ModelState ParaGridIO::getModelState(const std::string& filePath, ModelMetadata&
         if (ncFile.getGroupCount()) {
             // Read in the legacy file format with groups
             netCDF::NcGroup dataGroup = ncFile.getGroup("data");
-            readDimensions<netCDF::NcGroup>(dataGroup);
+            readDimensions<netCDF::NcGroup>(dataGroup, metadata);
             vars = dataGroup.getVars();
         } else {
             // File format without groups
-            readDimensions<netCDF::NcFile>(ncFile);
+            readDimensions<netCDF::NcFile>(ncFile, metadata);
             vars = ncFile.getVars();
         }
 

core/src/include/ParaGridIO.hpp

@@ -116,7 +116,7 @@ class ParaGridIO : public ParametricGrid::IParaGridIO {
 
     //! A templated function to get the dimensions from a netCDF file, whether
     //! they are stored in the root or a group.
-    template <typename N> void readDimensions(const N& node);
+    template <typename N> void readDimensions(const N& node, ModelMetadata& metadata);
 
     template <typename N>
     static ModelState readForcingTimeGeneric(