[cig-commits] r19125 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d

[brad at geodynamics.org]

Author: brad
Date: 2011-10-28 16:48:01 -0700 (Fri, 28 Oct 2011)
New Revision: 19125

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/plot_faultth.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg
Log:
Updated results script.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/plot_faultth.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/plot_faultth.py	2011-10-28 22:22:39 UTC (rev 19124)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/plot_faultth.py	2011-10-28 23:48:01 UTC (rev 19125)
@@ -10,7 +10,7 @@
 # ======================================================================
 #
 
-cell = "quad4"
+cell = "tri3"
 dx = 200
 
 inputRoot = "output/%s_%3dm_%s-fault" % (cell,dx,"gradient")
@@ -66,7 +66,7 @@
 h5.close()
 
 
-nrows = 2
+nrows = 3
 ncols = 2
 irow = 1
 icol = 1
@@ -101,6 +101,14 @@
 ax.set_ylabel("Normal Traction (MPa)")
 icol += 1
 
+irow = 3
+icol = 1
+ax = fig.axes(nrows, ncols, irow, icol)
+ax.plot(-slip[:,:,0], -traction[:,:,0])
+ax.set_xlabel("Slip (s)")
+ax.set_ylabel("Shear Traction (MPa)")
+icol += 1
+
 pylab.show()
 
 

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg	2011-10-28 22:22:39 UTC (rev 19124)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg	2011-10-28 23:48:01 UTC (rev 19125)
@@ -26,7 +26,8 @@
 # ----------------------------------------------------------------------
 [pylithapp.mesh_generator]
 reader = pylith.meshio.MeshIOCubit
-#reorder_mesh = True
+reorder_mesh = True
+
 distributor.partitioner = parmetis
 distributor.write_partition = True
 distributor.data_writer.filename = partition.vtk
@@ -120,21 +121,21 @@
 vertex_data_fields=[displacement,velocity]
 output_freq = time_step
 time_step = 0.04999*s
-writer = pylith.meshio.DataWriterHDF5ExtMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
 
 [pylithapp.timedependent.interfaces.fault.output]
 vertex_info_fields = [strike_dir,normal_dir,initial_traction]
 vertex_data_fields = [slip,slip_rate,traction]
 output_freq = time_step
 time_step = 0.04999*s
-writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
 
 [pylithapp.timedependent.materials.elastic.output]
 cell_data_fields = []
 cell_filter = pylith.meshio.CellFilterAvgMesh
 output_freq = time_step
 time_step = 12.0*s
-writer = pylith.meshio.DataWriterHDF5ExtMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
 
 # ----------------------------------------------------------------------
 # PETSc

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py	2011-10-28 22:22:39 UTC (rev 19124)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tabulate_onfaultdata.py	2011-10-28 23:48:01 UTC (rev 19125)
@@ -10,8 +10,8 @@
 # ----------------------------------------------------------------------
 #
 
-cell = "quad4"
-dx = 100
+cell = "tri3"
+dx = 200
 
 inputRoot = "output/%s_%3dm_%s-fault" % (cell,dx,"gradient")
 outdir = "scecfiles/%s_%3dm_%s/" % (cell,dx,"gradient")
@@ -97,7 +97,6 @@
 locName = "st%+04ddp%03d"
 
 lengthScale = 1000.0
-timeScale = 1000.0
 dip = 7.5
 strike = targets[:,1] / lengthScale
 
@@ -111,12 +110,13 @@
     fout.write("# num_timesteps = %8d\n" % nsteps)
     fout.write(locHeader % (strike[iloc], dip))
     fout.write(headerB)
+    # Vertical components hardwired to zero because they don't exist in 2-D.
     data = numpy.transpose((timeStamps, 
                             -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]*0.0,
+                            +slip_rate[:,iloc,1]*0.0,
+                            +traction[:,iloc,1]*0.0/1e+6))
     numpy.savetxt(fout, data, fmt='%14.6e')
     fout.close()

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg	2011-10-28 22:22:39 UTC (rev 19124)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3.cfg	2011-10-28 23:48:01 UTC (rev 19125)
@@ -14,7 +14,7 @@
 [pylithapp.timedependent.materials.elastic]
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
+quadrature.cell.quad_order = 1
 
 
 # ----------------------------------------------------------------------