Merge pull request #583 from xylar/fix-mesh-density-in-converter

Add `meshDensity` in mesh conversion if not already present

Author: xylar

mesh_tools/mesh_conversion_tools_netcdf_c/mpas_mesh_converter.cpp

@@ -476,7 +476,14 @@ int readGridInput(const string inputFilename){/*{{{*/
 
     meshDensity.clear();
     meshDensity.resize(cells.size());
-    ncutil::get_var(inputFilename, "meshDensity", &meshDensity[0]);
+    try {
+        ncutil::get_var(inputFilename, "meshDensity", &meshDensity[0]);
+    } catch (...) {
+        // should be all ones by default
+        for(int i = 0; i < nCells; i++){
+            meshDensity[i] = 1.0;
+        }
+    }
 
     xCellDistance = fabs(xCellRange[1] - xCellRange[0]);
     yCellDistance = fabs(yCellRange[1] - yCellRange[0]);