[cig-commits] r17037 - short/3D/PyLith/branches/v1.5-stable/libsrc/friction

[brad at geodynamics.org]

Author: brad
Date: 2010-07-12 15:41:30 -0700 (Mon, 12 Jul 2010)
New Revision: 17037

Modified:
   short/3D/PyLith/branches/v1.5-stable/libsrc/friction/SlipWeakening.cc
Log:
Adjusted slip-weakening model to include instaneous healing.

Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/friction/SlipWeakening.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/friction/SlipWeakening.cc	2010-07-12 22:41:05 UTC (rev 17036)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/friction/SlipWeakening.cc	2010-07-12 22:41:30 UTC (rev 17037)
@@ -298,11 +298,16 @@
   assert(0 != numStateVars);
   assert(0 != numProperties);
 
-  const double slipPrev = stateVars[s_slipPrev];
- 
-  stateVars[s_slipPrev] = stateVars[s_slipCum];
-  stateVars[s_slipCum] += fabs(slip - slipPrev);
- 
+  if (slipRate != 0.0) {
+    const double slipPrev = stateVars[s_slipPrev];
+
+    stateVars[s_slipPrev] = stateVars[s_slipCum];
+    stateVars[s_slipCum] += fabs(slip - slipPrev);
+  } else {
+    // Sliding has stopped, so reset state variables.
+    stateVars[s_slipPrev] = slip;
+    stateVars[s_slipCum] = 0.0;
+  } // else
 } // _updateStateVars