[cig-commits] r19596 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Tue Feb 7 12:46:54 PST 2012
Author: brad
Date: 2012-02-07 12:46:53 -0800 (Tue, 07 Feb 2012)
New Revision: 19596
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_050m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_050m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_200m.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_faultinfo.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg
Log:
Added cfg files for tpv12.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_faultinfo.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_faultinfo.py 2012-02-07 18:48:16 UTC (rev 19595)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/plot_faultinfo.py 2012-02-07 20:46:53 UTC (rev 19596)
@@ -8,7 +8,7 @@
# ======================================================================
#
-sim = "tpv11"
+sim = "tpv12"
cell = "tri3"
dx = 200
@@ -29,20 +29,26 @@
dipDist = vertices[:,1]/sin(60.0*pi/180.0)
# Expected tractions
-mask = numpy.bitwise_and(dipDist <= -10.5e+3, dipDist >= -13.5e+3)
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":
+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
-pylab.plot(shearE, dipDist, 'ko',
+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[:,0], dipDist, 'rx',
- normalE, dipDist, 'ko',
+ normalE, dipDist, 'k*',
traction[:,1], dipDist, 'rx')
pylab.show()
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg 2012-02-07 18:48:16 UTC (rev 19595)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -104,7 +104,6 @@
db_initial_tractions = spatialdata.spatialdb.SimpleDB
db_initial_tractions.label = Initial fault tractions
-db_initial_tractions.query_type = linear
friction = pylith.friction.SlipWeakening
friction.label = Slip weakening
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_offfaultdata.py 2012-02-07 18:48:16 UTC (rev 19595)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_offfaultdata.py 2012-02-07 20:46:53 UTC (rev 19596)
@@ -8,7 +8,7 @@
# ----------------------------------------------------------------------
#
-sim = "tpv10"
+sim = "tpv12"
cell = "quad4"
dx = 50
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_onfaultdata.py 2012-02-07 18:48:16 UTC (rev 19595)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_onfaultdata.py 2012-02-07 20:46:53 UTC (rev 19596)
@@ -8,9 +8,9 @@
# ----------------------------------------------------------------------
#
-sim = "tpv11"
-cell = "tri3"
-dx = 200
+sim = "tpv12"
+cell = "quad4"
+dx = 50
inputRoot = "output/%s_%s_%03dm-fault" % (sim,cell,dx)
outdir = "scecfiles/%s_%s_%03dm/" % (sim,cell,dx)
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg 2012-02-07 18:48:16 UTC (rev 19595)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -6,6 +6,7 @@
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
db_initial_tractions.iohandler.filename = tpv10_tractions.spatialdb
+db_initial_tractions.query_type = linear
friction.db_properties.iohandler.filename = tpv10_friction.spatialdb
friction.db_properties.query_type = nearest
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,22 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 8.00001*s
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions.iohandler.filename = tpv12_tractions.spatialdb
+db_initial_tractions.query_type = linear
+
+friction.db_properties.iohandler.filename = tpv12_friction.spatialdb
+friction.db_properties.query_type = nearest
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_friction.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_friction.spatialdb 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,19 @@
+#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 = 6
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+-5.2499 0.0 0.70 0.10 0.5 0.2
+-5.2500 0.0 0.54 0.10 0.5 0.2
+-6.7500 0.0 0.54 0.10 0.5 0.2
+-6.7501 0.0 0.70 0.10 0.5 0.2
+-7.4999 0.0 0.70 0.10 0.5 0.2
+-7.5000 0.0 10000.0 0.10 0.5 0.2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_050m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_050m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_050m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_050m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation.norm_viscosity = 0.4
+normalizer.wave_period = 0.15*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.0025*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_quad4_050m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_quad4_050m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_quad4_050m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_100m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_100m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_100m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation.norm_viscosity = 0.4
+normalizer.wave_period = 0.3*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.005*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_quad4_100m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_quad4_100m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_quad4_100m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,30 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation.norm_viscosity = 0.4
+normalizer.wave_period = 0.6*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_quad4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_quad4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_quad4_200m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tractions.spatialdb 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,17 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = traction-shear traction-normal
+ value-units = MPa MPa
+ num-locs = 4
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+ 0.0 0.0 0.0 0.0
+ 0.0 -11.95115 56.07458 -101.98215
+ 0.0 -11.95116 0.0 -199.10633
+ 0.0 -20.00000 0.0 -333.2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_050m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_050m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_050m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_050m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.15*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.005*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_tri3_050m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_tri3_050m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_tri3_050m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_100m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_100m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_100m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.3*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_tri3_100m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_tri3_100m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_tri3_100m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_tri3_200m.cfg 2012-02-07 20:46:53 UTC (rev 19596)
@@ -0,0 +1,29 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tri3_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.6*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.02*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_tri3_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_tri3_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_tri3_200m-elastic.h5
More information about the CIG-COMMITS
mailing list