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

[knepley at geodynamics.org]

Author: knepley
Date: 2008-02-14 17:13:33 -0800 (Thu, 14 Feb 2008)
New Revision: 9355

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
Log:
Fix for boundary and fault section output


Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-02-14 23:44:22 UTC (rev 9354)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-02-15 01:13:33 UTC (rev 9355)
@@ -192,13 +192,15 @@
   assert(0 != name);
 
   try {
+    // Correctly handle boundary and fault meshes
+    //const int depth = mesh->depth();
+    const int depth = mesh->getDimension();
     const std::string labelName = 
       (mesh->hasLabel("censored depth")) ? "censored depth" : "depth";
     const ALE::Obj<Mesh::numbering_type>& numbering = 
-      mesh->getFactory()->getNumbering(mesh, labelName, mesh->depth());
+      mesh->getFactory()->getNumbering(mesh, labelName, depth);
     const int fiberDim = 
-      field->getFiberDimension(*mesh->getLabelStratum(labelName, 
-						      mesh->depth())->begin());
+      field->getFiberDimension(*mesh->getLabelStratum(labelName, depth)->begin());
     const int enforceDim = (fieldType != VECTOR_FIELD) ? fiberDim : 3;
 
     if (!_wroteCellHeader) {