[cig-commits] r21435 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22
brad at geodynamics.org
brad at geodynamics.org
Fri Mar 1 14:57:32 PST 2013
Author: brad
Date: 2013-03-01 14:57:32 -0800 (Fri, 01 Mar 2013)
New Revision: 21435
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_ruptime.py
Log:
Updated tabulation scripts for separate dirs for each sim.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_offfaultdata.py 2013-03-01 22:56:49 UTC (rev 21434)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_offfaultdata.py 2013-03-01 22:57:32 UTC (rev 21435)
@@ -8,7 +8,7 @@
# ----------------------------------------------------------------------
#
-sim = "tpv22"
+sim = "tpv23"
cell = "tet4"
dx = 200
@@ -91,7 +91,7 @@
pt = locName % (round(10*body[iloc]),
round(10*strike[iloc]),
round(10*dip))
- filename = "%sbody%s.dat" % (outputRoot, pt)
+ filename = "%s/body%s.dat" % (outputRoot, pt)
fout = open(filename, 'w')
fout.write(headerA)
fout.write("# time_step = %14.6E\n" % dt)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_onfaultdata.py 2013-03-01 22:56:49 UTC (rev 21434)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_onfaultdata.py 2013-03-01 22:57:32 UTC (rev 21435)
@@ -8,12 +8,12 @@
# ----------------------------------------------------------------------
#
-sim = "tpv23"
+sim = "tpv22"
cell = "tet4"
dx = 200
inputRoot = "output/%s_%s_%03dm" % (sim, cell,dx)
-outputRoot = "scecfiles/%s_%s_%03dm" % (sim, cell,dx)
+outputRoot = "scecfiles/%s_%s_%03dm/" % (sim, cell,dx)
# ----------------------------------------------------------------------
import h5py
@@ -85,7 +85,7 @@
for iloc in xrange(ntargets):
pt = locName % (round(10*strike[iloc]),
round(10*dip[iloc]))
- filename = "%s-%s-%s.dat" % (outputRoot, fault, pt)
+ filename = "%s/%s-%s.dat" % (outputRoot, fault, pt)
fout = open(filename, 'w');
fout.write(headerA)
fout.write("# time_step = %14.6E\n" % dt)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_ruptime.py 2013-03-01 22:56:49 UTC (rev 21434)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv22/tabulate_ruptime.py 2013-03-01 22:57:32 UTC (rev 21435)
@@ -8,12 +8,12 @@
# ----------------------------------------------------------------------
#
-sim = "tpv22"
+sim = "tpv23"
cell = "tet4"
dx = 200
inputRoot = "output/%s_%s_%03dm" % (sim, cell, dx)
-outputRoot = "scecfiles/%s_%s_%03dm" % (sim, cell, dx)
+outputRoot = "scecfiles/%s_%s_%03dm/" % (sim, cell, dx)
# ----------------------------------------------------------------------
import h5py
@@ -79,7 +79,7 @@
distDip = -vertices[:,2]
distStrike = vertices[:,1]
- filename = "%s_ruptime_%s.dat" % (outputRoot, fault)
+ filename = "%s/ruptime_%s.dat" % (outputRoot, fault)
fout = open(filename, "w")
fout.write(headerA)
data = numpy.transpose((distStrike, distDip, rupTime))
More information about the CIG-COMMITS
mailing list