[cig-commits] r18839 - short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2011-08-19 16:28:05 -0700 (Fri, 19 Aug 2011)
New Revision: 18839

Modified:
   short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc
Log:
Fixed small bug (typo in previous attempted fix).

Modified: short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc	2011-08-19 16:47:18 UTC (rev 18838)
+++ short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/FaultCohesiveDyn.cc	2011-08-19 23:28:05 UTC (rev 18839)
@@ -558,11 +558,10 @@
     // Do not allow fault interpenetration and set fault opening to
     // zero if fault is under compression.
     const int indexN = spaceDim - 1;
-    if (slipVertex[indexN] + dSlipVertex[indexN] < 0.0)
-      dSlipVertex[indexN] = -dSlipVertex[indexN];
     const double lagrangeTpdtNormal = lagrangeTVertex[indexN] + 
       lagrangeTIncrVertex[indexN] + dLagrangeTpdtVertex[indexN];
-    if (lagrangeTpdtNormal < 0.0)
+    if (lagrangeTpdtNormal < 0.0 || 
+	slipVertex[indexN] + dSlipVertex[indexN] < 0.0)
       dSlipVertex[indexN] = -slipVertex[indexN];
 
     // Set change in slip.