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

[brad at geodynamics.org]

Author: brad
Date: 2013-02-28 16:15:32 -0800 (Thu, 28 Feb 2013)
New Revision: 21420

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

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py	2013-02-28 23:51:31 UTC (rev 21419)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/create_friction.py	2013-03-01 00:15:32 UTC (rev 21420)
@@ -53,6 +53,12 @@
 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
 
+# Fault edges
+mask = z < -14.999e+3
+coefStatic[mask] = 1.0e+6
+mask = numpy.abs(y) < 13.999e+3
+coefStatic[mask] = 1.0e+6
+
 cs = CSCart()
 cs._configure()
 cs.initialize()
@@ -100,7 +106,7 @@
 for i in xrange(nptsz):
     istart = i
     istop = i+nptsx*nptsz
-    points[istart:istop:nptsz,1] = x
+    points[istart:istop:nptsz,0] = x
 for i in xrange(nptsx):
     istart = i*nptsz
     istop = (i+1)*nptsz
@@ -113,6 +119,14 @@
 mask = points[:,2] > -4000.0
 cohesion = mask*(0.30+0.000675*(+4.0e+3+points[:,2])) + ~mask*0.30
 
+# Fault edges
+mask = z < -14.999e+3
+coefStatic[mask] = 1.0e+6
+mask = numpy.abs(x) > 5.999e+3
+coefStatic[mask] = 1.0e+6
+mask = numpy.abs(x) < 99.0
+coefStatic[mask] = 1.0e+6
+
 writer = SimpleIOAscii()
 writer.inventory.filename = "empty"
 writer._configure()

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/pylithapp.cfg	2013-02-28 23:51:31 UTC (rev 21419)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/pylithapp.cfg	2013-03-01 00:15:32 UTC (rev 21420)
@@ -112,7 +112,7 @@
 
 [pylithapp.timedependent.interfaces.fault_main]
 id = 100
-label = fault_main
+label = main_fault
 
 friction = pylith.friction.SlipWeakeningTimeStable
 friction.label = Slip Weakening on main fault
@@ -130,7 +130,7 @@
 
 [pylithapp.timedependent.interfaces.fault_branch]
 id = 101
-label = fault_branch
+label = branch_fault
 
 friction = pylith.friction.SlipWeakening
 friction.label = Slip Weakening on branch fault