[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