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

surendra at geodynamics.org surendra at geodynamics.org
Sat May 29 20:20:23 PDT 2010


Author: surendra
Date: 2010-05-29 20:20:23 -0700 (Sat, 29 May 2010)
New Revision: 16834

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py
Log:
Switched from slip-rate to slip in post-processing for rupture contour.  Also fixed an error in creating off-fault stations

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou	2010-05-29 23:25:42 UTC (rev 16833)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou	2010-05-30 03:20:23 UTC (rev 16834)
@@ -74,7 +74,7 @@
 # ----------------------------------------------------------------------
 # Create nodeset for off-fault stations
 # ----------------------------------------------------------------------
-group "off_fault" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and z_min >= -7.501e+3 and z_coord = 3.000e+3
+group "off_fault" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and z_min >= -7.501e+3 and y_coord = 3.000e+3
 nodeset 25 group off_fault
 nodeset 25 name "off_fault"
 

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-29 23:25:42 UTC (rev 16833)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tabulate_ruptime.py	2010-05-30 03:20:23 UTC (rev 16834)
@@ -52,7 +52,7 @@
     filename = "%s-fault_t%05d.vtk" % (outputRoot,timestamp)
     data = reader.read(filename)
     fields = data['vertex_fields']
-    slipRate = fields['slip_rate'][:,:].squeeze()
+    slipRate = fields['slip'][:,:].squeeze()
 
     # Compute magnitude of slip rate
     slipRateMag = (slipRate[:,0]**2 + slipRate[:,1]**2)**0.5



More information about the CIG-COMMITS mailing list