[cig-commits] r22223 - in short/3D/PyLith/trunk: libsrc/pylith/friction unittests/libtests/friction/data

[brad at geodynamics.org]

Author: brad
Date: 2013-06-11 15:05:41 -0700 (Tue, 11 Jun 2013)
New Revision: 22223

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/friction/RateStateAgeing.cc
   short/3D/PyLith/trunk/unittests/libtests/friction/data/RateStateAgeingData.cc
Log:
Fixed sign error in friction derivative for rate-state friction.

Modified: short/3D/PyLith/trunk/libsrc/pylith/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-11 20:29:55 UTC (rev 22222)
+++ short/3D/PyLith/trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-11 22:05:41 UTC (rev 22223)
@@ -367,9 +367,9 @@
     const PylithScalar slipRate0 = properties[p_slipRate0];
 
     if (slipRate >= slipRateLinear) {
-      frictionDeriv = normalTraction * a / (slipRate * _dt);
+      frictionDeriv = -normalTraction * a / (slipRate * _dt);
     } else {
-      frictionDeriv = normalTraction * a / (slipRateLinear * _dt);
+      frictionDeriv = -normalTraction * a / (slipRateLinear * _dt);
     } // else
   } // if    
 

Modified: short/3D/PyLith/trunk/unittests/libtests/friction/data/RateStateAgeingData.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/friction/data/RateStateAgeingData.cc	2013-06-11 20:29:55 UTC (rev 22222)
+++ short/3D/PyLith/trunk/unittests/libtests/friction/data/RateStateAgeingData.cc	2013-06-11 22:05:41 UTC (rev 22223)
@@ -134,8 +134,8 @@
 };
 
 const PylithScalar pylith::friction::RateStateAgeingData::_frictionDeriv[] = {
-  -2.2*0.0125/(0.0011*0.01),
-  -2.3*0.0225/(0.0021*0.01),
+  2.2*0.0125/(0.0011*0.01),
+  2.3*0.0225/(0.0021*0.01),
 };
 
 const PylithScalar pylith::friction::RateStateAgeingData::_slip[] = {