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

[knepley at geodynamics.org]

Author: knepley
Date: 2013-02-13 04:46:07 -0800 (Wed, 13 Feb 2013)
New Revision: 21361

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc
Log:
Fix (for now) the determination of dof for HSFExt output

Modified: short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc	2013-02-13 00:31:10 UTC (rev 21360)
+++ short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterHDF5Ext.cc	2013-02-13 12:46:07 UTC (rev 21361)
@@ -464,11 +464,15 @@
     err = VecView(vector, binaryViewer);CHECK_PETSC_ERROR(err);
     ++_datasets[field.label()].numTimeSteps;
 
-    PetscSection section = field.petscSection();
+    DM           dm;
+    PetscSection section;// = field.petscSection();
     PetscInt     dof     = 0, n, numLocalVertices = 0, numVertices, vStart;
     IS           globalVertexNumbers;
 
+    err = VecGetDM(vector, &dm);CHECK_PETSC_ERROR(err);
+    err = DMGetDefaultSection(dm, &section);CHECK_PETSC_ERROR(err);
     assert(section);
+
     err = DMPlexGetDepthStratum(dmMesh, 0, &vStart, PETSC_NULL);CHECK_PETSC_ERROR(err);
     err = DMPlexGetVertexNumbering(dmMesh, &globalVertexNumbers);CHECK_PETSC_ERROR(err);
     err = ISGetSize(globalVertexNumbers, &n);CHECK_PETSC_ERROR(err);