[cig-commits] r18473 - short/3D/PyLith/trunk/libsrc/pylith/meshio

[brad at geodynamics.org]

Author: brad
Date: 2011-05-26 18:16:02 -0700 (Thu, 26 May 2011)
New Revision: 18473

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc
Log:
Fixed inconsistency in dimensions of time array b/t DataWriterHDF5 and DataWriterHDF5Ext.

Modified: short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc	2011-05-27 00:01:39 UTC (rev 18472)
+++ short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc	2011-05-27 01:16:02 UTC (rev 18473)
@@ -572,27 +572,31 @@
   assert(_h5->hasGroup(group));
   std::string datasetFullName = std::string(group) + "/time";
 
-  const int ndims = 1;
+  const int ndims = 3;
 
   // Each time stamp has a size of 1.
-  hsize_t dimsChunk[1];
+  hsize_t dimsChunk[3]; // Use 3 dims for compatibility with PETSc viewer
   dimsChunk[0] = 1;
+  dimsChunk[1] = 1;
+  dimsChunk[2] = 1;
 
-  
-
   if (!_h5->hasDataset(datasetFullName.c_str())) {
     // Create dataset
     // Dataset has unknown size.
-    hsize_t dims[1];
+    hsize_t dims[3];
     dims[0] = H5S_UNLIMITED;
+    dims[1] = 1;
+    dims[2] = 1;
     _h5->createDataset(group, "time", dims, dimsChunk, ndims, 
 		       H5T_NATIVE_DOUBLE);
   } // if
   
   // Write time stamp as chunk to HDF5 file.
   // Current dimensions of dataset.
-  hsize_t dims[1];
+  hsize_t dims[3];
   dims[0] = _tstampIndex+1;
+  dims[1] = 1;
+  dims[2] = 1;
   _h5->writeDatasetChunk(group, "time", &t, dims, dimsChunk, ndims, 
 			 _tstampIndex, H5T_NATIVE_DOUBLE);