[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