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

[brad at geodynamics.org]

Author: brad
Date: 2012-08-30 16:05:53 -0700 (Thu, 30 Aug 2012)
New Revision: 20650

Modified:
   short/3D/PyLith/branches/v1.7-stable/libsrc/pylith/friction/RateStateAgeing.cc
Log:
Prevent zero initial state in RateStateAgeing.

Modified: short/3D/PyLith/branches/v1.7-stable/libsrc/pylith/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-stable/libsrc/pylith/friction/RateStateAgeing.cc	2012-08-30 23:02:08 UTC (rev 20649)
+++ short/3D/PyLith/branches/v1.7-stable/libsrc/pylith/friction/RateStateAgeing.cc	2012-08-30 23:05:53 UTC (rev 20650)
@@ -320,7 +320,7 @@
     const PylithScalar slipRate0 = properties[p_slipRate0];
 
     // Prevent zero value for theta
-    const PylithScalar theta = std::max(stateVars[s_state], 1.0e-10);
+    const PylithScalar theta = std::max(stateVars[s_state], 1.0e-12);
 
     if (slipRate >= slipRateLinear) {
       mu_f = f0 + a*log(slipRate / slipRate0) + b*log(slipRate0*theta/L);