diff --git a/src/DataIO_NetCDF.cpp b/src/DataIO_NetCDF.cpp
index d20e982812..9cad754624 100644
--- a/src/DataIO_NetCDF.cpp
+++ b/src/DataIO_NetCDF.cpp
@@ -1428,8 +1428,8 @@ int DataIO_NetCDF::writeData_vector_scalar(DataSet const* ds) {
 }
 
 /** Write 1D DataSets that share an index dimension. */
-int DataIO_NetCDF::writeData_1D(DataSet const* ds, Dimension const& dim, SetArray const& sets) {
-  if (ds->Type() == DataSet::PH)
+int DataIO_NetCDF::writeData_1D(DataSet const* firstSet, Dimension const& dim, SetArray const& sets) {
+  if (firstSet->Type() == DataSet::PH)
     mprintf("Warning: Currently only State information saved for pH sets.\n");
   if (debug_ > 0) {
     mprintf("DEBUG: Sets for dimension '%s' %f %f:", dim.label(), dim.Min(), dim.Step());
@@ -1439,7 +1439,7 @@ int DataIO_NetCDF::writeData_1D(DataSet const* ds, Dimension const& dim, SetArra
   }
   // Define the dimension. Ensure name is unique by appending an index.
   if (EnterDefineMode(ncid_)) return 1;
-  int dimIdx = defineDim( "length", ds->Size(), ds->Meta().Legend() );
+  int dimIdx = defineDim( "length", firstSet->Size(), firstSet->Meta().Legend() );
   if (dimIdx < 0) return 1;
   NcDim lengthDim = Dimensions_[dimIdx];
   
@@ -1466,6 +1466,8 @@ int DataIO_NetCDF::writeData_1D(DataSet const* ds, Dimension const& dim, SetArra
         return 1;
     }
     // Add variable(s)
+    if (debug_ > 0)
+      mprintf("DEBUG: Defining variable '%s' for set '%s'\n", it->DS()->Meta().PrintName().c_str(), it->DS()->legend());
     variables.push_back( VarArray() );
     VarArray& set_vars = variables.back();
     if (it->DS()->Type() == DataSet::STRING) {
@@ -1529,22 +1531,23 @@ int DataIO_NetCDF::writeData_1D(DataSet const* ds, Dimension const& dim, SetArra
       set_vars.push_back( mat3x3Var );
     } else if (it->DS()->Type() == DataSet::XYMESH) {
       // ----- XY mesh set -------------
+      DataSet_Mesh const& meshSet = static_cast<DataSet_Mesh const&>( *(it->DS()) );
       // Define the Y variable
-      NcVar yVar = defineVar(Dimensions_[dimIdx].DID(), NC_DOUBLE, ds->Meta().PrintName(), "Y");
+      NcVar yVar = defineVar(Dimensions_[dimIdx].DID(), NC_DOUBLE, meshSet.Meta().PrintName(), "Y");
       if (yVar.Empty()) {
-        mprinterr("Error: Could not define Y variable for set '%s'\n", ds->legend());
+        mprinterr("Error: Could not define Y variable for set '%s'\n", meshSet.legend());
         return 1;
       }
       set_vars.push_back( yVar );
       // Define the X variable.
-      NcVar xVar = defineVar(Dimensions_[dimIdx].DID(), NC_DOUBLE, ds->Meta().PrintName(), "X");
+      NcVar xVar = defineVar(Dimensions_[dimIdx].DID(), NC_DOUBLE, meshSet.Meta().PrintName(), "X");
       if (xVar.Empty()) {
-        mprinterr("Error: Could not define Y variable for '%s'\n", ds->legend());
+        mprinterr("Error: Could not define X variable for '%s'\n", meshSet.legend());
         return 1;
       }
       set_vars.push_back( xVar );
       // Add DataSet info to variable
-      if (AddDataSetInfo(ds, std::vector<int>(1, xVar.VID()), ncid_, yVar.VID())) return 1;
+      if (AddDataSetInfo(it->DS(), std::vector<int>(1, xVar.VID()), ncid_, yVar.VID())) return 1;
     } else {
       // ----- All other 1D sets -------
       set_vars.push_back( defineVar(lengthDim.DID(), dtype, it->DS()->Meta().PrintName(), "Y") );
diff --git a/src/NC_Routines.cpp b/src/NC_Routines.cpp
index b59fd39142..3edc6d7dd7 100644
--- a/src/NC_Routines.cpp
+++ b/src/NC_Routines.cpp
@@ -69,7 +69,7 @@ NC::FormatType NC::GetFormatType(std::string const& fnameIn) {
 // NC::CheckErr()
 bool NC::CheckErr(int ncerr) {
   if ( ncerr != NC_NOERR ) {
-    mprintf("%s\n", nc_strerror(ncerr));
+    mprinterr("Error: %s\n", nc_strerror(ncerr));
     return true;
   }
   return false;
diff --git a/src/Version.h b/src/Version.h
index 503351e627..24ea04b037 100644
--- a/src/Version.h
+++ b/src/Version.h
@@ -12,7 +12,7 @@
  * Whenever a number that precedes <revision> is incremented, all subsequent
  * numbers should be reset to 0.
  */
-#define CPPTRAJ_INTERNAL_VERSION "V6.29.11"
+#define CPPTRAJ_INTERNAL_VERSION "V6.29.12"
 /// PYTRAJ relies on this
 #define CPPTRAJ_VERSION_STRING CPPTRAJ_INTERNAL_VERSION
 #endif