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

[brad at geodynamics.org]

Author: brad
Date: 2010-02-21 14:02:11 -0800 (Sun, 21 Feb 2010)
New Revision: 16301

Modified:
   short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc
Log:
Fixed error.

Modified: short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc	2010-02-21 20:43:07 UTC (rev 16300)
+++ short/3D/PyLith/trunk/libsrc/friction/RateStateAgeing.cc	2010-02-21 22:02:11 UTC (rev 16301)
@@ -314,7 +314,6 @@
 
   const double dt = _dt;
   const double thetaTVertex = stateVars[s_state];
-  const double thetaTpdtVertex;
   const double L = properties[p_L];
   const double expTerm = exp(-slipRate * dt / L);
   const double vDtL = slipRate * dt / L;
@@ -324,6 +323,7 @@
   // Above ODE is integrated from t->t+dt keeping slipRate constant gives
   // thetaTpdt = thetaT * exp(- slipRate * theta / L)
   //             + L / slipRate * (1 -  exp(- slipRate * theta / L))
+  double thetaTpdtVertex = 0.0;
   if (vDtL < 0.00001)
     // As (slipRate * dt / L) --> 0, exp(-slipRate * dt / L) --> 1
     // So using first three term in the Taylor series expansion of