[cig-commits] r22207 - in short/3D/PyLith/trunk/examples/3d/hex8: . spatialdb

[brad at geodynamics.org]

Author: brad
Date: 2013-06-10 17:54:50 -0700 (Mon, 10 Jun 2013)
New Revision: 22207

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/greensfns.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/impulse_amplitude.spatialdb
   short/3D/PyLith/trunk/examples/3d/hex8/step21.cfg
Log:
Minor cleanup of step21 (Green's functions). Simplify spatialdata file and compute two shear components of slip. Added output of impulse amplitude in fault info file.

Modified: short/3D/PyLith/trunk/examples/3d/hex8/greensfns.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/greensfns.cfg	2013-06-11 00:53:25 UTC (rev 22206)
+++ short/3D/PyLith/trunk/examples/3d/hex8/greensfns.cfg	2013-06-11 00:54:50 UTC (rev 22207)
@@ -40,9 +40,9 @@
 
 # Set the parameters for the fault interface condition.
 [greensfns.interfaces.fault]
-# Generate impulses for lateral slip only, no fault opening.
-# Fault DOF 0 corresponds to left-lateral slip.
-impulse_dof = [0]
+# Generate impulses for lateral slip and reverse slip, no fault opening.
+# Fault DOF 0 and 1 correspond to left-lateral and reverse slip.
+impulse_dof = [0, 1]
 
 # Set the amplitude of the slip impulses (amplitude is nonzero on only
 # a subset of the fault)

Modified: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/impulse_amplitude.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/impulse_amplitude.spatialdb	2013-06-11 00:53:25 UTC (rev 22206)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/impulse_amplitude.spatialdb	2013-06-11 00:54:50 UTC (rev 22207)
@@ -1,14 +1,14 @@
 // -*- C++ -*- (syntax highlighting)
 //
 // This spatial database specifies the distribution of slip on the
-// fault surface for the forward problem.
+// fault surface for the Green's function impulses.
 //
 #SPATIAL.ascii 1
 SimpleDB {
   num-values = 1
   value-names =  slip
   value-units =  m
-  num-locs = 35
+  num-locs = 10
   data-dim = 2 // Data is specified on a plane.
   space-dim = 3
   cs-data = cartesian {
@@ -21,38 +21,15 @@
 // (2) y coordinate (km)
 // (3) z coordinate (km)
 // (4) slip (m)
+0.0  -2.1   0.0    0.0
 0.0  -2.0   0.0   -1.0
-0.0  -2.0  -1.0   -1.0
-0.0  -2.0  -2.0   -1.0
+0.0  +2.0   0.0   -1.0
+0.0  +2.1   0.0    0.0
+
+0.0  -2.1  -3.0    0.0
 0.0  -2.0  -3.0   -1.0
-0.0  -1.0   0.0   -1.0
-0.0  -1.0  -1.0   -1.0
-0.0  -1.0  -2.0   -1.0
-0.0  -1.0  -3.0   -1.0
-0.0   0.0   0.0   -1.0
-0.0   0.0  -1.0   -1.0
-0.0   0.0  -2.0   -1.0
-0.0   0.0  -3.0   -1.0
-0.0   1.0   0.0   -1.0
-0.0   1.0  -1.0   -1.0
-0.0   1.0  -2.0   -1.0
-0.0   1.0  -3.0   -1.0
-0.0   2.0   0.0   -1.0
-0.0   2.0  -1.0   -1.0
-0.0   2.0  -2.0   -1.0
-0.0   2.0  -3.0   -1.0
-0.0  -3.0   0.0    0.0
-0.0  -3.0  -1.0    0.0
-0.0  -3.0  -2.0    0.0
-0.0  -3.0  -3.0    0.0
-0.0   3.0   0.0    0.0
-0.0   3.0  -1.0    0.0
-0.0   3.0  -2.0    0.0
-0.0   3.0  -3.0    0.0
-0.0  -3.0  -4.0    0.0
-0.0  -2.0  -4.0    0.0
-0.0  -1.0  -4.0    0.0
-0.0   0.0  -4.0    0.0
-0.0   1.0  -4.0    0.0
-0.0   2.0  -4.0    0.0
-0.0   3.0  -4.0    0.0
+0.0  +2.0  -3.0   -1.0
+0.0  +2.1  -3.0    0.0
+
+0.0  -2.0  -3.1    0.0
+0.0  +2.0  -3.1    0.0

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step21.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step21.cfg	2013-06-11 00:53:25 UTC (rev 22206)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step21.cfg	2013-06-11 00:54:50 UTC (rev 22207)
@@ -119,7 +119,7 @@
 [pylithapp.problem.interfaces.fault.output]
 writer = pylith.meshio.DataWriterHDF5
 writer.filename = output/step21-fault.h5
-vertex_info_fields = [area]
+vertex_info_fields = [area, impulse_amplitude]
 
 # We will not be writing any cell_info or cell_data fields for material
 # output.