[cig-commits] r22319 - short/3D/PyLith/trunk/libsrc/pylith/faults

[knepley at geodynamics.org]

Author: knepley
Date: 2013-06-15 23:55:14 -0700 (Sat, 15 Jun 2013)
New Revision: 22319

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
Log:
Now restrict fault mesh with material id

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-06-16 06:27:22 UTC (rev 22318)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/CohesiveTopology.cc	2013-06-16 06:55:14 UTC (rev 22319)
@@ -674,6 +674,7 @@
   PYLITH_METHOD_BEGIN;
 
   assert(faultMesh);
+  const char    *labelname = "material-id";
   PetscErrorCode err;
 
   faultMesh->coordsys(mesh.coordsys());
@@ -681,7 +682,8 @@
   PetscDM dmMesh = mesh.dmMesh();assert(dmMesh);
   PetscDM dmFaultMesh;
 
-  err = DMPlexCreateCohesiveSubmesh(dmMesh, constraintCell ? PETSC_TRUE : PETSC_FALSE, &dmFaultMesh);PYLITH_CHECK_ERROR(err);
+
+  err = DMPlexCreateCohesiveSubmesh(dmMesh, constraintCell ? PETSC_TRUE : PETSC_FALSE, labelname, materialId, &dmFaultMesh);PYLITH_CHECK_ERROR(err);
   std::string meshLabel = "fault_" + std::string(label);
 
   PetscReal lengthScale = 1.0;