[cig-commits] r19830 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210
brad at geodynamics.org
brad at geodynamics.org
Tue Mar 20 16:40:21 PDT 2012
Author: brad
Date: 2012-03-20 16:40:21 -0700 (Tue, 20 Mar 2012)
New Revision: 19830
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tet4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tractions.spatialdb
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb
Log:
Small correction to nucleation stresses in tpv10. Added tpv11.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py 2012-03-20 23:40:21 UTC (rev 19830)
@@ -0,0 +1,55 @@
+#!/usr/bin/env python
+#
+# ======================================================================
+#
+# Brad T. Aagaard
+# U.S. Geological Survey
+#
+# ======================================================================
+#
+
+sim = "tpv11"
+cell = "tet4"
+dx = 200
+
+inputRoot = "output/%s_%s_%03dm-fault_info" % (sim,cell,dx)
+
+# ----------------------------------------------------------------------
+import tables
+import numpy
+import pylab
+
+# ----------------------------------------------------------------------
+h5 = tables.openFile("%s.h5" % inputRoot, 'r')
+vertices = h5.root.geometry.vertices[:]
+traction = h5.root.vertex_fields.initial_traction[:].squeeze()
+h5.close()
+
+from math import pi, sin
+dipDist = vertices[:,2]/sin(60.0*pi/180.0)
+
+# Expected tractions
+if sim == "tpv10":
+ mask = numpy.bitwise_and(dipDist <= -10.5e+3, dipDist >= -13.5e+3)
+ shearE = mask*(0.2e+6 + ((0.760+0.0057)*7378*-dipDist)) + \
+ ~mask*0.55*7378*-dipDist
+ normalE = 7378*dipDist
+
+elif sim == "tpv11": # No 2-D case
+ mask = numpy.bitwise_and(dipDist <= -10.5e+3, dipDist >= -13.5e+3)
+ shearE = mask*(0.2e+6 + ((0.570+0.0057)*7378*-dipDist)) + \
+ ~mask*0.55*7378*-dipDist
+ normalE = 7378*dipDist
+
+elif sim == "tpv12":
+ mask = dipDist > -13.8e+3
+ normalE = mask*7390.01*dipDist + ~mask*14427.98*dipDist
+ shearE = mask*-0.549847*normalE + ~mask*0.0
+
+pylab.plot(shearE, dipDist, 'k*',
+ traction[:,1], dipDist, 'rx',
+ normalE, dipDist, 'k*',
+ traction[:,2], dipDist, 'rx')
+pylab.show()
+
+# End of file
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py
___________________________________________________________________
Name: svn:executable
+ *
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg 2012-03-20 23:04:38 UTC (rev 19829)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg 2012-03-20 23:40:21 UTC (rev 19830)
@@ -50,7 +50,7 @@
[pylithapp.timedependent.formulation.time_step]
#total_time = 15.01*s
-total_time = 0.01*s
+total_time = 0.0*s
dt = 0.01*s
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb 2012-03-20 23:04:38 UTC (rev 19829)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb 2012-03-20 23:40:21 UTC (rev 19830)
@@ -54,10 +54,10 @@
-5.250 -6.000 0.0 0.0 42.60795000 -77.4690000
-5.250 -3.000 0.0 0.0 42.60795000 -77.4690000
-5.250 -1.501 0.0 0.0 42.60795000 -77.4690000
--5.250 -1.500 0.0 0.0 51.32321250 -77.4690000
--5.250 -1.499 0.0 0.0 60.03847500 -77.4690000
--5.250 +1.499 0.0 0.0 60.03847500 -77.4690000
--5.250 +1.500 0.0 0.0 51.32321250 -77.4690000
+-5.250 -1.500 0.0 0.0 46.83546582 -77.4690000
+-5.250 -1.499 0.0 0.0 51.06298165 -77.4690000
+-5.250 +1.499 0.0 0.0 51.06298165 -77.4690000
+-5.250 +1.500 0.0 0.0 46.83546582 -77.4690000
-5.250 +1.501 0.0 0.0 42.60795000 -77.4690000
-5.250 +3.000 0.0 0.0 42.60795000 -77.4690000
-5.250 +6.000 0.0 0.0 42.60795000 -77.4690000
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11.cfg (from rev 19829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11.cfg 2012-03-20 23:40:21 UTC (rev 19830)
@@ -0,0 +1,23 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
+db_initial_tractions.iohandler.filename = tpv11_tractions.spatialdb
+db_initial_tractions.query_type = linear
+
+friction.db_properties.iohandler.filename = tpv11_friction.spatialdb
+friction.db_properties.query_type = nearest
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,dip_dir,normal_dir,initial_traction]
+
+
+# End of file
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_friction.spatialdb (from rev 19829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_friction.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_friction.spatialdb 2012-03-20 23:40:21 UTC (rev 19830)
@@ -0,0 +1,22 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
+ value-units = none none m MPa
+ num-locs = 8
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
+-7.499 -15.000 0.0 10000.0 0.448 0.5 0.2
+-7.499 -14.999 0.0 0.570 0.448 0.5 0.2
+-7.499 +14.999 0.0 0.570 0.448 0.5 0.2
+-7.499 +15.000 0.0 10000.0 0.448 0.5 0.2
+
+-7.500 -15.000 0.0 10000.0 0.448 0.5 0.2
+-7.500 -14.999 0.0 10000.0 0.448 0.5 0.2
+-7.500 +14.999 0.0 10000.0 0.448 0.5 0.2
+-7.500 +15.000 0.0 10000.0 0.448 0.5 0.2
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tet4_200m.cfg (from rev 19829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tet4_200m.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tet4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tet4_200m.cfg 2012-03-20 23:40:21 UTC (rev 19830)
@@ -0,0 +1,28 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.6*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv11_tet4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv11_tet4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv11_tet4_200m-elastic.h5
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tractions.spatialdb (from rev 19829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv10_tractions.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv11_tractions.spatialdb 2012-03-20 23:40:21 UTC (rev 19830)
@@ -0,0 +1,129 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear-leftlateral traction-shear-updip traction-normal
+ value-units = MPa MPa MPa
+ num-locs = 108
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
++0.000 -15.000 0.0 0.0 0.0 0.0
++0.000 -6.000 0.0 0.0 0.0 0.0
++0.000 -3.000 0.0 0.0 0.0 0.0
++0.000 -1.501 0.0 0.0 0.0 0.0
++0.000 -1.500 0.0 0.0 0.0 0.0
++0.000 -1.499 0.0 0.0 0.0 0.0
++0.000 +1.499 0.0 0.0 0.0 0.0
++0.000 +1.500 0.0 0.0 0.0 0.0
++0.000 +1.501 0.0 0.0 0.0 0.0
++0.000 +3.000 0.0 0.0 0.0 0.0
++0.000 +6.000 0.0 0.0 0.0 0.0
++0.000 +15.000 0.0 0.0 0.0 0.0
+
+-3.000 -15.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -6.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -3.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.501 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.500 0.0 0.0 24.34740000 -44.2680000
+-3.000 -1.499 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.499 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.500 0.0 0.0 24.34740000 -44.2680000
+-3.000 +1.501 0.0 0.0 24.34740000 -44.2680000
+-3.000 +3.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 +6.000 0.0 0.0 24.34740000 -44.2680000
+-3.000 +15.000 0.0 0.0 24.34740000 -44.2680000
+
+-5.249 -15.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -6.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -3.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.501 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.500 0.0 0.0 42.59983420 -77.4542440
+-5.249 -1.499 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.499 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.500 0.0 0.0 42.59983420 -77.4542440
+-5.249 +1.501 0.0 0.0 42.59983420 -77.4542440
+-5.249 +3.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 +6.000 0.0 0.0 42.59983420 -77.4542440
+-5.249 +15.000 0.0 0.0 42.59983420 -77.4542440
+
+-5.250 -15.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -6.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -3.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 -1.501 0.0 0.0 42.60795000 -77.4690000
+-5.250 -1.500 0.0 0.0 43.15568832 -77.4690000
+-5.250 -1.499 0.0 0.0 43.70342665 -77.4690000
+-5.250 +1.499 0.0 0.0 43.70342665 -77.4690000
+-5.250 +1.500 0.0 0.0 43.15568832 -77.4690000
+-5.250 +1.501 0.0 0.0 42.60795000 -77.4690000
+-5.250 +3.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 +6.000 0.0 0.0 42.60795000 -77.4690000
+-5.250 +15.000 0.0 0.0 42.60795000 -77.4690000
+
+-5.251 -15.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -6.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -3.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 -1.501 0.0 0.0 42.61606580 -77.4837560
+-5.251 -1.500 0.0 0.0 43.71173206 -77.4837560
+-5.251 -1.499 0.0 0.0 44.80739833 -77.4837560
+-5.251 +1.499 0.0 0.0 44.80739833 -77.4837560
+-5.251 +1.500 0.0 0.0 43.71173206 -77.4837560
+-5.251 +1.501 0.0 0.0 42.61606580 -77.4837560
+-5.251 +3.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 +6.000 0.0 0.0 42.61606580 -77.4837560
+-5.251 +15.000 0.0 0.0 42.61606580 -77.4837560
+
+-6.749 -15.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -6.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -3.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 -1.501 0.0 0.0 54.77353420 -99.5882440
+-6.749 -1.500 0.0 0.0 56.15324314 -99.5882440
+-6.749 -1.499 0.0 0.0 57.53295207 -99.5882440
+-6.749 +1.499 0.0 0.0 57.53295207 -99.5882440
+-6.749 +1.500 0.0 0.0 56.15324314 -99.5882440
+-6.749 +1.501 0.0 0.0 54.77353420 -99.5882440
+-6.749 +3.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 +6.000 0.0 0.0 54.77353420 -99.5882440
+-6.749 +15.000 0.0 0.0 54.77353420 -99.5882440
+
+-6.750 -15.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -6.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -3.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 -1.501 0.0 0.0 54.78165000 -99.6030000
+-6.750 -1.500 0.0 0.0 55.47159927 -99.6030000
+-6.750 -1.499 0.0 0.0 56.16154855 -99.6030000
+-6.750 +1.499 0.0 0.0 56.16154855 -99.6030000
+-6.750 +1.500 0.0 0.0 55.47159927 -99.6030000
+-6.750 +1.501 0.0 0.0 54.78165000 -99.6030000
+-6.750 +3.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 +6.000 0.0 0.0 54.78165000 -99.6030000
+-6.750 +15.000 0.0 0.0 54.78165000 -99.6030000
+
+-6.751 -15.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -6.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -3.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.501 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.500 0.0 0.0 54.78976580 -99.6177560
+-6.751 -1.499 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.499 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.500 0.0 0.0 54.78976580 -99.6177560
+-6.751 +1.501 0.0 0.0 54.78976580 -99.6177560
+-6.751 +3.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 +6.000 0.0 0.0 54.78976580 -99.6177560
+-6.751 +15.000 0.0 0.0 54.78976580 -99.6177560
+
+-7.500 -15.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -6.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -3.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.501 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.500 0.0 0.0 60.86850000 -110.6700000
+-7.500 -1.499 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.499 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.500 0.0 0.0 60.86850000 -110.6700000
+-7.500 +1.501 0.0 0.0 60.86850000 -110.6700000
+-7.500 +3.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 +6.000 0.0 0.0 60.86850000 -110.6700000
+-7.500 +15.000 0.0 0.0 60.86850000 -110.6700000
More information about the CIG-COMMITS
mailing list