[cig-commits] r16257 - short/3D/PyLith/trunk/libsrc/friction

[surendra at geodynamics.org]

Author: surendra
Date: 2010-02-12 10:11:16 -0800 (Fri, 12 Feb 2010)
New Revision: 16257

Modified:
   short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc
Log:
Found a bug in SlipWeakening

Modified: short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc	2010-02-12 03:47:46 UTC (rev 16256)
+++ short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc	2010-02-12 18:11:16 UTC (rev 16257)
@@ -258,16 +258,15 @@
   double mu_f = 0.0;
   if (normalTraction < 0.0) {
     // if fault is in compression
-    if (stateVars[s_slipCum] < properties[p_coefD]) {
+    if (stateVars[s_slipCum] < properties[p_d0]) {
 	// if/else linear slip-weakening form of mu_f 
 	mu_f = properties[p_coefS] -
 	  (properties[p_coefS] - properties[p_coefD]) * 
 	  stateVars[s_slipCum] / properties[p_d0];
-	friction = - mu_f * normalTraction;
       } else {
 	mu_f = properties[p_coefD];
-	friction = - mu_f * normalTraction;
       } // if/else
+    friction = - mu_f * normalTraction;
   } // if
 
   PetscLogFlops(5);