[cig-commits] r19092 - short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/friction

[brad at geodynamics.org]

Author: brad
Date: 2011-10-18 13:07:12 -0700 (Tue, 18 Oct 2011)
New Revision: 19092

Modified:
   short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/friction/RateStateAgeing.cc
Log:
Removed hardwired minimum slip rate. Use data member.

Modified: short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/friction/RateStateAgeing.cc	2011-10-18 20:06:33 UTC (rev 19091)
+++ short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/friction/RateStateAgeing.cc	2011-10-18 20:07:12 UTC (rev 19092)
@@ -316,7 +316,7 @@
 
     // Since regulatized friction -> 0 as slipRate -> 0, limit slip
     // rate to some minimum value
-    const double slipRateEff = std::max(1.0e-14, slipRate);
+    const double slipRateEff = std::max(_minSlipRate, slipRate);
 
     const double slipRate0 = properties[p_slipRate0];
     const double a = properties[p_a];
@@ -374,7 +374,7 @@
 
   // Since regulatized friction -> 0 as slipRate -> 0, limit slip
   // rate to some minimum value
-  const double slipRateEff = std::max(1.0e-14, slipRate);
+  const double slipRateEff = std::max(_minSlipRate, slipRate);
 
   const double dt = _dt;
   const double thetaTVertex = stateVars[s_state];