[cig-commits] r16809 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205

[surendra at geodynamics.org]

Author: surendra
Date: 2010-05-27 11:59:55 -0700 (Thu, 27 May 2010)
New Revision: 16809

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
Log:
Fixed some minor errors in post-processing scripts for tpv205

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg	2010-05-27 16:53:44 UTC (rev 16808)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/pylithapp.cfg	2010-05-27 18:59:55 UTC (rev 16809)
@@ -68,7 +68,7 @@
 dt = 0.01*s
 
 
-[pylithapp.timedependent.implicit]
+[pylithapp.timedependent.formulation]
 output = [domain,subdomain]
 output.subdomain = pylith.meshio.OutputSolnSubset
 

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py	2010-05-27 16:53:44 UTC (rev 16808)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_onfaultdata.py	2010-05-27 18:59:55 UTC (rev 16809)
@@ -10,7 +10,7 @@
 # ----------------------------------------------------------------------
 #
 
-cell = "tet4"
+cell = "hex8"
 dx = 200
 dt = 0.05
 
@@ -145,8 +145,8 @@
                             -slip[:,iloc,0],
                             -slip_rate[:,iloc,0],
                             -traction[:,iloc,0]/1e+6,
-                            -slip[:,iloc,1],
-                            -slip_rate[:,iloc,1],
-                            -traction[:,iloc,1]/1e+6))
+                            +slip[:,iloc,1],
+                            +slip_rate[:,iloc,1],
+                            +traction[:,iloc,1]/1e+6))
     numpy.savetxt(fout, data, fmt='%14.6e')
     fout.close()

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py	2010-05-27 16:53:44 UTC (rev 16808)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py	2010-05-27 18:59:55 UTC (rev 16809)
@@ -38,7 +38,7 @@
 nsteps = timestamps.shape[0]
 
 # Set default rupture time to a large value (1.0e+30)
-rupTime = 1.0e+30 * numpy.ones( (npts,), dtype=numpy.float64)
+rupTime = 1.0e+9 * numpy.ones( (npts,), dtype=numpy.float64)
 
 # Create buffer for current rupture time
 tmpTime = numpy.zeros( (npts,), dtype=numpy.float64)
@@ -87,8 +87,8 @@
 
 lengthScale = 1000.0
 timeScale = 1000.0
-distDip = -vertices[:,2] / lengthScale
-distStrike = vertices[:,0] / lengthScale
+distDip = -vertices[:,2]
+distStrike = vertices[:,0]
 time =  timestamps / timeScale
 print "time", time