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

brad at geodynamics.org brad at geodynamics.org
Fri Mar 1 14:08:39 PST 2013


Author: brad
Date: 2013-03-01 14:08:39 -0800 (Fri, 01 Mar 2013)
New Revision: 21429

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/tabulate_offfaultdata.py
Log:
Set targets for off-fault data.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/tabulate_offfaultdata.py	2013-03-01 19:22:27 UTC (rev 21428)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv24/tabulate_offfaultdata.py	2013-03-01 22:08:39 UTC (rev 21429)
@@ -8,7 +8,7 @@
 # ----------------------------------------------------------------------
 #
 
-sim = "tpv23"
+sim = "tpv24"
 cell = "tet4"
 dx = 200
 
@@ -21,13 +21,16 @@
 import time
 
 # ----------------------------------------------------------------------
-targets = numpy.array([[-3.0e+3,  -5.0e+3,      0.0],
-                       [-3.0e+3,  +5.0e+3,      0.0],
-                       [-3.0e+3, +15.0e+3,      0.0],
-                       [+3.0e+3,  -5.0e+3,      0.0],
-                       [+3.0e+3,  +5.0e+3,      0.0],
-                       [+3.0e+3, +15.0e+3,      0.0],
+targets = numpy.array([[-3.0e+3,  -2.0e+3,      0.0],
+                       [+3.0e+3,  -2.0e+3,      0.0],
+                       [-3.0e+3,  +2.0e+3,      0.0],
+                       [+0.6e+3,  +2.0e+3,      0.0],
+                       [+4.2e+3,  +2.0e+3,      0.0],
+                       [-3.0e+3,  +8.0e+3,      0.0],
+                       [+2.3e+3,  +8.0e+3,      0.0],
+                       [+7.6e+3,  +8.0e+3,      0.0],
     ]) 
+targets[:,1] += 2.0e+3
 
 
 h5 = h5py.File("%s.h5" % (inputRoot), 'r', driver="sec2")
@@ -86,6 +89,7 @@
 dip = 0.0
 body = -targets[:,0] / lengthScale
 strike = (targets[:,1]) / lengthScale
+strike -= 2.0
 
 for iloc in xrange(ntargets):
     pt = locName % (round(10*body[iloc]), 



More information about the CIG-COMMITS mailing list