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

[brad at geodynamics.org]

Author: brad
Date: 2010-04-21 17:07:59 -0700 (Wed, 21 Apr 2010)
New Revision: 16573

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
Log:
Do not allow fault interpenetration with friction.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-21 20:08:52 UTC (rev 16572)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-04-22 00:07:59 UTC (rev 16573)
@@ -518,8 +518,11 @@
         dSlipVertex[iDim] += 
 	  orientationVertex[iDim*spaceDim+kDim] * dispRelVertex[kDim];
 
-    // Set fault opening to zero if fault is under compression.
+    // Do not allow fault interpenetration and set fault opening to
+    // zero if fault is under compression.
     const int indexN = spaceDim - 1;
+    if (dSlipVertex[indexN] < 0.0)
+      dSlipVertex[indexN] = 0.0;
     const double lagrangeTpdtNormal = lagrangeTVertex[indexN] + 
       lagrangeTIncrVertex[indexN] + dLagrangeTpdtVertex[indexN];
     if (lagrangeTpdtNormal < 0.0)