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

[brad at geodynamics.org]

Author: brad
Date: 2010-04-29 16:52:32 -0700 (Thu, 29 Apr 2010)
New Revision: 16600

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
Log:
Added ability to output slip rate for fault with friction.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-29 04:04:17 UTC (rev 16599)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-29 23:52:32 UTC (rev 16600)
@@ -964,15 +964,17 @@
   const int cohesiveDim = _faultMesh->dimension();
   const int spaceDim = _quadrature->spaceDim();
 
-  const int slipStrLen = strlen("final_slip");
-  const int timeStrLen = strlen("slip_time");
-
   double scale = 0.0;
   int fiberDim = 0;
   if (0 == strcasecmp("slip", name)) {
     const topology::Field<topology::SubMesh>& slip = _fields->get("slip");
     return slip;
 
+  } else if (0 == strcasecmp("slip_rate", name)) {
+    const topology::Field<topology::SubMesh>& slipRate =
+      _fields->get("slip rate");
+    return slipRate;
+
   } else if (cohesiveDim > 0 && 0 == strcasecmp("strike_dir", name)) {
     const ALE::Obj<RealSection>& orientationSection = _fields->get(
       "orientation").section();
@@ -1018,7 +1020,7 @@
     buffer.scale(1.0);
     return buffer;
 
-  } else if (0 == strncasecmp("initial_traction", name, slipStrLen)) {
+  } else if (0 == strcasecmp("initial_traction", name)) {
     assert(0 != _dbInitialTract);
     topology::Field<topology::SubMesh>& buffer =
         _fields->get("buffer (vector)");