[cig-commits] r16541 - short/3D/PyLith/trunk/libsrc/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-04-12 21:30:01 -0700 (Mon, 12 Apr 2010)
New Revision: 16541

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
Log:
Fixed bug where time step in friction model was not set.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-13 02:27:48 UTC (rev 16540)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-13 04:30:01 UTC (rev 16541)
@@ -341,6 +341,9 @@
   _updateSlipRate(*fields);
   _sensitivitySetup(jacobian);
 
+  // Update time step in friction (can vary).
+  _friction->timeStep(_dt);
+
   const int spaceDim = _quadrature->spaceDim();
 
   // Allocate arrays for vertex values
@@ -626,6 +629,11 @@
   double_array lagrangeTpdtVertex(spaceDim);
   double_array dLagrangeTpdtVertex(spaceDim);
 
+  _updateSlipRate(*fields);
+
+  // Update time step in friction (can vary).
+  _friction->timeStep(_dt);
+
   // Get section information
   double_array slipVertex(spaceDim);
   const ALE::Obj<RealSection>& slipSection = _fields->get("slip").section();