[cig-commits] r22224 - in short/3D/PyLith/branches/v1.7-trunk: libsrc/pylith/friction unittests/libtests/friction/data

[brad at geodynamics.org]

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

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

Modified: short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-11 22:05:41 UTC (rev 22223)
+++ short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/friction/RateStateAgeing.cc	2013-06-11 22:05:50 UTC (rev 22224)
@@ -371,9 +371,9 @@
     // We approximate the derivative wrt slip by taking the
     // derivative wrt slip rate and multiplying by dt. 
     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/branches/v1.7-trunk/unittests/libtests/friction/data/RateStateAgeingData.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/friction/data/RateStateAgeingData.cc	2013-06-11 22:05:41 UTC (rev 22223)
+++ short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/friction/data/RateStateAgeingData.cc	2013-06-11 22:05:50 UTC (rev 22224)
@@ -134,8 +134,8 @@
 };
 
 const PylithScalar pylith::friction::RateStateAgeingData::_frictionDeriv[] = {
-  -2.2e+5*0.0125/(0.0011*0.01),
-  -2.3e+5*0.0225/(0.0021*0.01),
+  2.2e+5*0.0125/(0.0011*0.01),
+  2.3e+5*0.0225/(0.0021*0.01),
 };
 
 const PylithScalar pylith::friction::RateStateAgeingData::_slip[] = {