Correct the dimension ordering, which was wrong.

core/src/ParaGridIO.cpp

@@ -102,7 +102,7 @@ ModelState ParaGridIO::getModelState(const std::string& filePath)
         std::vector<netCDF::NcDim> varDims = var.getDims();
         std::string dimKey = "";
         for (netCDF::NcDim& dim : varDims) {
-            dimKey = dim.getName() + dimKey;
+            dimKey += dim.getName();
         }
         if (!dimensionKeys.count(dimKey)) {
             throw std::out_of_range(
@@ -215,12 +215,6 @@ void ParaGridIO::dumpModelState(
         dimMap.at(entry.second).push_back(ncFromMAMap.at(entry.first));
     }
 
-    // Reverse the order of the dimensions to translate between column-major ModelArray and
-    // row-major netCDF
-    for (auto& [type, v] : dimMap) {
-        std::reverse(v.begin(), v.end());
-    }
-
     std::set<std::string> restartFields = { hiceName, ciceName, hsnowName, ticeName, sstName,
         sssName, maskName, coordsName }; // TODO and others
     // Loop through either the above list (isRestart) or all provided fields(!isRestart)