[cig-commits] r21422 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24

[brad at geodynamics.org]

Author: brad
Date: 2013-02-28 16:33:05 -0800 (Thu, 28 Feb 2013)
New Revision: 21422

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py
Log:
Small fix.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py	2013-03-01 00:19:53 UTC (rev 21421)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py	2013-03-01 00:33:05 UTC (rev 21422)
@@ -54,9 +54,9 @@
 weakTime = mask*(r/(0.7*vs) + 0.081*rcrit/(0.7*vs)*(1.0/(1.0-(r/rcrit)**2)-1)) + ~mask*1.0e+9
 
 # Fault edges
-mask = z < -14.999e+3
+mask = points[:,2] < -14.999e+3
 coefStatic[mask] = 1.0e+6
-mask = numpy.abs(y) < 13.999e+3
+mask = numpy.abs(points[:,1]) > 13.999e+3
 coefStatic[mask] = 1.0e+6
 
 cs = CSCart()
@@ -120,11 +120,11 @@
 cohesion = mask*(0.30+0.000675*(+4.0e+3+points[:,2])) + ~mask*0.30
 
 # Fault edges
-mask = z < -14.999e+3
+mask = points[:,2] < -14.999e+3
 coefStatic[mask] = 1.0e+6
-mask = numpy.abs(x) > 5.999e+3
+mask = numpy.abs(points[:,0]) > 5.999e+3
 coefStatic[mask] = 1.0e+6
-mask = numpy.abs(x) < 99.0
+mask = numpy.abs(points[:,0]) < 99.0
 coefStatic[mask] = 1.0e+6
 
 writer = SimpleIOAscii()