[cig-commits] r21381 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22

[brad at geodynamics.org]

Author: brad
Date: 2013-02-20 15:10:36 -0800 (Wed, 20 Feb 2013)
New Revision: 21381

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py
Log:
Fixed some errors in setting up friction.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py	2013-02-20 22:53:16 UTC (rev 21380)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/create_friction.py	2013-02-20 23:10:36 UTC (rev 21381)
@@ -46,11 +46,11 @@
 d0 = 0.30*numpy.ones( (nptsy*nptsz), dtype=numpy.float64)
 
 mask = points[:,2] > -5000.0
-cohesion = mask*0.0014*(-50000-points[:,2]) + ~mask*0.0
+cohesion = mask*0.0014*(+5.0e+3+points[:,2]) + ~mask*0.0
 
 rcrit = 3000.0
 mask = r <= rcrit
-weakTime = mask*(r/(0.7*vs) + 0.081*rcrit/(0.7*vs)*(1.0/(1.0-(r/rcrit)**2-1))) + ~mask*1.0e+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
 
 cs = CSCart()
 cs._configure()