[cig-commits] r19850 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210
brad at geodynamics.org
brad at geodynamics.org
Thu Mar 22 11:46:34 PDT 2012
Author: brad
Date: 2012-03-22 11:46:34 -0700 (Thu, 22 Mar 2012)
New Revision: 19850
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tet4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_bottom_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_initial_stress.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_lateral_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4_200m.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/matprops.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py
Log:
Added files for tpv12 and tpv13.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/matprops.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/matprops.spatialdb 2012-03-22 15:19:16 UTC (rev 19849)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/matprops.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -1,8 +1,8 @@
#SPATIAL.ascii 1
SimpleDB {
- num-values = 3
- value-names = density vs vp
- value-units = kg/m**3 m/s m/s
+ num-values = 6
+ value-names = density vs vp cohesion friction-angle dilatation-angle
+ value-units = kg/m**3 m/s m/s MPa degree degree
num-locs = 1
data-dim = 0
space-dim = 3
@@ -11,4 +11,4 @@
space-dim = 3
}
}
-0.0 0.0 0.0 2700.0 3300.0 5716.0
+0.0 0.0 0.0 2700.0 3300.0 5716.0 5.0 40.36453657309736 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py 2012-03-22 15:19:16 UTC (rev 19849)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/plot_faultinfo.py 2012-03-22 18:46:34 UTC (rev 19850)
@@ -8,7 +8,7 @@
# ======================================================================
#
-sim = "tpv11"
+sim = "tpv12"
cell = "tet4"
dx = 200
@@ -31,20 +31,20 @@
# 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
+ 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
+ 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
+ shearE = mask*0.549847*normalE + ~mask*0.0
pylab.plot(shearE, dipDist, 'k*',
traction[:,1], dipDist, 'rx',
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py 2012-03-22 15:19:16 UTC (rev 19849)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_offfaultdata.py 2012-03-22 18:46:34 UTC (rev 19850)
@@ -5,53 +5,36 @@
# Brad T. Aagaard
# U.S. Geological Survey
#
-# <LicenseText>
-#
# ----------------------------------------------------------------------
#
+sim = "tpv10"
cell = "tet4"
-dx = 200
-dt = 0.05
+dx = 100
-# outputRoot = "output/%s_%3dm_%s" % (cell,dx,"refine")
-# outdir = "scecfiles/%s_%3dm_%s/" % (cell,dx,"refine")
+inputRoot = "output/%s_%s_%03dm" % (sim,cell,dx)
+outdir = "scecfiles/%s_%s_%03dm/" % (sim,cell,dx)
-outputRoot = "output/%s_%3dm" % (cell,dx)
-outdir = "scecfiles/%s_%3dm/" % (cell,dx)
-
# ----------------------------------------------------------------------
-
+import tables
import numpy
import time
-from pylith.utils.VTKDataReader import VTKDataReader
-
# ----------------------------------------------------------------------
-timestamps = numpy.arange(50,15001,50)
-targets = numpy.array([[-3000.0, 0.0, -0.0],
- [-2000.0, 0.0, -0.0],
- [-1000.0, 0.0, -0.0],
- [+1000.0, 0.0, -0.0],
- [+2000.0, 0.0, -0.0],
- [+3000.0, 0.0, -0.0],
- [-1000.0, 0.0,-300.0],
- [ -500.0, 0.0,-300.0],
- [ +500.0, 0.0,-300.0],
- [+1000.0, 0.0,-300.0],
- [-3000.0, +12000.0, -0.0],
- [+3000.0, +12000.0, -0.0]])
+targets = numpy.array([[-3000.0, 0.0, -0.0],
+ [-2000.0, 0.0, -0.0],
+ [-1000.0, 0.0, -0.0],
+ [+1000.0, 0.0, -0.0],
+ [+2000.0, 0.0, -0.0],
+ [+3000.0, 0.0, -0.0],
+ [-3000.0, +12000.0, -0.0],
+ [+3000.0, +12000.0, -0.0]])
-reader = VTKDataReader()
tolerance = 1.0e-6
-# Get vertices and find indices of target locations
-#filename = "%s-off_fault_t%05d.vtk" % (outputRoot,timestamps[0])
-filename = "%s_t%05d.vtk" % (outputRoot,timestamps[0])
-data = reader.read(filename)
-
-vertices = numpy.array(data['vertices'])
+h5 = tables.openFile("%s.h5" % inputRoot, 'r')
+vertices = h5.root.geometry.vertices[:]
ntargets = targets.shape[0]
indices = []
for target in targets:
@@ -59,33 +42,32 @@
(vertices[:,1]-target[1])**2 +
(vertices[:,2]-target[2])**2 )**0.5
min = numpy.min(dist)
- indices.append(numpy.where(dist <= min+tolerance)[0])
+ indices.append(numpy.where(dist <= min+tolerance)[0][0])
-print "Indices", indices
-print "Coordinates of selected points:",vertices[indices,:]
+print "Indices: ", indices
+print "Coordinates of selected points:\n",vertices[indices,:]
-# Extract values
-nsteps = timestamps.shape[0]
-disp = numpy.zeros((nsteps,ntargets,3)) # 3-D array (time, targets, components)
-vel = numpy.zeros((nsteps,ntargets,3))
-itime = 0
-for timestamp in timestamps:
- #filename = "%s-off_fault_t%05d.vtk" % (outputRoot,timestamp)
- filename = "%s_t%05d.vtk" % (outputRoot,timestamp)
- data = reader.read(filename)
- fields = data['vertex_fields']
- disp[itime,0:ntargets,:] = fields['displacement'][indices,:].squeeze()
- vel[itime,0:ntargets,:] = fields['velocity'][indices,:].squeeze()
- itime += 1
+# Get datasets
+disp = h5.root.vertex_fields.displacement[:]
+vel = h5.root.vertex_fields.velocity[:]
+timeStamps = h5.root.time[:].ravel()
+nsteps = timeStamps.shape[0]
+dt = timeStamps[1] - timeStamps[0]
+h5.close()
+
+# Extract locations
+disp = disp[:,indices,:]
+vel = vel[:,indices,:]
+
# Write data
headerA = \
- "# problem = TPV210\n" + \
- "# author = Surendra Somala, Brad Aagaard\n" + \
+ "# problem = %s\n" % sim.upper() + \
+ "# author = Brad Aagaard\n" + \
"# date = %s\n" % (time.asctime()) + \
"# code = PyLith\n" + \
- "# code_version = 1.5.0a\n" + \
+ "# code_version = 1.7.0a (scecdynrup branch)\n" + \
"# element_size = %s\n" % dx
headerB = \
"# Time series in 7 columns of E14.6:\n" + \
@@ -106,11 +88,9 @@
locName = "%+04dst%03ddp%03d"
lengthScale = 1000.0
-timeScale = 1000.0
-dip = -targets[:,2] / lengthScale
-strike = targets[:,1] / lengthScale
body = -targets[:,0] / lengthScale
-time = timestamps / timeScale
+strike = targets[:,1] / lengthScale
+dip = -targets[:,2] / lengthScale
for iloc in xrange(ntargets):
pt = locName % (round(10*body[iloc]),
@@ -125,7 +105,7 @@
strike[iloc],
dip[iloc]))
fout.write(headerB)
- data = numpy.transpose((time,
+ data = numpy.transpose((timeStamps,
disp[:,iloc,1],
vel[:,iloc,1],
-disp[:,iloc,2],
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py 2012-03-22 15:19:16 UTC (rev 19849)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_onfaultdata.py 2012-03-22 18:46:34 UTC (rev 19850)
@@ -5,14 +5,12 @@
# Brad T. Aagaard
# U.S. Geological Survey
#
-# <LicenseText>
-#
# ----------------------------------------------------------------------
#
-sim = "tpv11"
+sim = "tpv10"
cell = "tet4"
-dx = 200
+dx = 100
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/tabulate_ruptime.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py 2012-03-22 15:19:16 UTC (rev 19849)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tabulate_ruptime.py 2012-03-22 18:46:34 UTC (rev 19850)
@@ -8,9 +8,9 @@
# ----------------------------------------------------------------------
#
-sim = "tpv11"
+sim = "tpv10"
cell = "tet4"
-dx = 200
+dx = 100
inputRoot = "output/%s_%s_%03dm-fault" % (sim,cell,dx)
outdir = "scecfiles/%s_%s_%03dm/" % (sim,cell,dx)
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12.cfg 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,30 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 8.00001*s
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
+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
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,dip_dir,normal_dir,initial_traction,static_coefficient]
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_friction.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_friction.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,100 @@
+#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 = 80
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
+-5.249 -15.000 0.0 10000.0 0.10 0.5 0.2
+-5.249 -14.999 0.0 0.70 0.10 0.5 0.2
+-5.249 -1.501 0.0 0.70 0.10 0.5 0.2
+-5.249 -1.500 0.0 0.70 0.10 0.5 0.2
+-5.249 -1.499 0.0 0.70 0.10 0.5 0.2
+-5.249 +1.499 0.0 0.70 0.10 0.5 0.2
+-5.249 +1.500 0.0 0.70 0.10 0.5 0.2
+-5.249 +1.501 0.0 0.70 0.10 0.5 0.2
+-5.249 +14.999 0.0 0.70 0.10 0.5 0.2
+-5.249 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-5.250 -15.000 0.0 10000.0 0.10 0.5 0.2
+-5.250 -14.999 0.0 0.70 0.10 0.5 0.2
+-5.250 -1.501 0.0 0.70 0.10 0.5 0.2
+-5.250 -1.500 0.0 0.62 0.10 0.5 0.2
+-5.250 -1.499 0.0 0.66 0.10 0.5 0.2
+-5.250 +1.499 0.0 0.66 0.10 0.5 0.2
+-5.250 +1.500 0.0 0.62 0.10 0.5 0.2
+-5.250 +1.501 0.0 0.70 0.10 0.5 0.2
+-5.250 +14.999 0.0 0.70 0.10 0.5 0.2
+-5.250 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-5.251 -15.000 0.0 10000.0 0.10 0.5 0.2
+-5.251 -14.999 0.0 0.70 0.10 0.5 0.2
+-5.251 -1.501 0.0 0.70 0.10 0.5 0.2
+-5.251 -1.500 0.0 0.62 0.10 0.5 0.2
+-5.251 -1.499 0.0 0.54 0.10 0.5 0.2
+-5.251 +1.499 0.0 0.54 0.10 0.5 0.2
+-5.251 +1.500 0.0 0.62 0.10 0.5 0.2
+-5.251 +1.501 0.0 0.70 0.10 0.5 0.2
+-5.251 +14.999 0.0 0.70 0.10 0.5 0.2
+-5.251 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-6.749 -15.000 0.0 10000.0 0.10 0.5 0.2
+-6.749 -14.999 0.0 0.70 0.10 0.5 0.2
+-6.749 -1.501 0.0 0.70 0.10 0.5 0.2
+-6.749 -1.500 0.0 0.62 0.10 0.5 0.2
+-6.749 -1.499 0.0 0.54 0.10 0.5 0.2
+-6.749 +1.499 0.0 0.54 0.10 0.5 0.2
+-6.749 +1.500 0.0 0.62 0.10 0.5 0.2
+-6.749 +1.501 0.0 0.70 0.10 0.5 0.2
+-6.749 +14.999 0.0 0.70 0.10 0.5 0.2
+-6.749 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-6.750 -15.000 0.0 10000.0 0.10 0.5 0.2
+-6.750 -14.999 0.0 0.70 0.10 0.5 0.2
+-6.750 -1.501 0.0 0.70 0.10 0.5 0.2
+-6.750 -1.500 0.0 0.66 0.10 0.5 0.2
+-6.750 -1.499 0.0 0.62 0.10 0.5 0.2
+-6.750 +1.499 0.0 0.62 0.10 0.5 0.2
+-6.750 +1.500 0.0 0.66 0.10 0.5 0.2
+-6.750 +1.501 0.0 0.70 0.10 0.5 0.2
+-6.750 +14.999 0.0 0.70 0.10 0.5 0.2
+-6.750 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-6.751 -15.000 0.0 10000.0 0.10 0.5 0.2
+-6.751 -14.999 0.0 0.70 0.10 0.5 0.2
+-6.751 -1.501 0.0 0.70 0.10 0.5 0.2
+-6.751 -1.500 0.0 0.70 0.10 0.5 0.2
+-6.751 -1.499 0.0 0.70 0.10 0.5 0.2
+-6.751 +1.499 0.0 0.70 0.10 0.5 0.2
+-6.751 +1.500 0.0 0.70 0.10 0.5 0.2
+-6.751 +1.501 0.0 0.70 0.10 0.5 0.2
+-6.751 +14.999 0.0 0.70 0.10 0.5 0.2
+-6.751 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-7.499 -15.000 0.0 10000.0 0.10 0.5 0.2
+-7.499 -14.999 0.0 0.70 0.10 0.5 0.2
+-7.499 -1.501 0.0 0.70 0.10 0.5 0.2
+-7.499 -1.500 0.0 0.70 0.10 0.5 0.2
+-7.499 -1.499 0.0 0.70 0.10 0.5 0.2
+-7.499 +1.499 0.0 0.70 0.10 0.5 0.2
+-7.499 +1.500 0.0 0.70 0.10 0.5 0.2
+-7.499 +1.501 0.0 0.70 0.10 0.5 0.2
+-7.499 +14.999 0.0 0.70 0.10 0.5 0.2
+-7.499 +15.000 0.0 10000.0 0.10 0.5 0.2
+
+-7.500 -15.000 0.0 10000.0 0.10 0.5 0.2
+-7.500 -14.999 0.0 10000.0 0.10 0.5 0.2
+-7.500 -1.501 0.0 10000.0 0.10 0.5 0.2
+-7.500 -1.500 0.0 10000.0 0.10 0.5 0.2
+-7.500 -1.499 0.0 10000.0 0.10 0.5 0.2
+-7.500 +1.499 0.0 10000.0 0.10 0.5 0.2
+-7.500 +1.500 0.0 10000.0 0.10 0.5 0.2
+-7.500 +1.501 0.0 10000.0 0.10 0.5 0.2
+-7.500 +14.999 0.0 10000.0 0.10 0.5 0.2
+-7.500 +15.000 0.0 10000.0 0.10 0.5 0.2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tet4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tet4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tet4_200m.cfg 2012-03-22 18:46:34 UTC (rev 19850)
@@ -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.0125*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv12_tet4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv12_tet4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv12_tet4_200m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv12_tractions.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,18 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear-leftlateral traction-shear-updip traction-normal
+ value-units = MPa MPa MPa
+ num-locs = 5
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
+ 0.0 0.0 0.0 0.0 0.0 0.0
+ 0.0 0.0 -11.950 0.0 -56.06918 -101.97234
+ 0.0 0.0 -11.951 0.0 -28.03694 -150.54227
+ 0.0 0.0 -11.952 0.0 0.0 -199.12032
+ 0.0 0.0 -20.000 0.0 0.0 -333.2
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13.cfg 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,95 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+gravity_field = spatialdata.spatialdb.GravityField
+# Adjust gravity to compensate for fluid pressure g = (rho-rhoH20)/rho*9.8
+gravity_field.acceleration = 6.17037037037037*m/s**2
+gravity_field.gravity_dir = [0.0, 0.0, -1.0]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.00001*s
+
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.DruckerPrager3D
+elastic.allow_tensile_yield = True
+
+[pylithapp.timedependent.materials.elastic]
+db_initial_stress = spatialdata.spatialdb.SimpleDB
+db_initial_stress.label = Initial stress
+db_initial_stress.iohandler.filename = tpv13_initial_stress.spatialdb
+db_initial_stress.query_type = linear
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_neg,x_pos,y_neg,y_pos,z_neg,x_neg_grav,x_pos_grav,y_neg_grav,y_pos_grav,z_neg_grav]
+bc.x_neg_grav = pylith.bc.Neumann
+bc.x_pos_grav = pylith.bc.Neumann
+bc.y_neg_grav = pylith.bc.Neumann
+bc.y_pos_grav = pylith.bc.Neumann
+bc.z_neg_grav = pylith.bc.Neumann
+
+[pylithapp.timedependent.bc.x_pos_grav]
+label = face_xpos
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Boundary tractions +x
+db_initial.iohandler.filename = tpv13_lateral_tractions.spatialdb
+db_initial.query_type = linear
+
+[pylithapp.timedependent.bc.x_neg_grav]
+label = face_xneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Boundary tractions -x
+db_initial.iohandler.filename = tpv13_lateral_tractions.spatialdb
+db_initial.query_type = linear
+
+[pylithapp.timedependent.bc.y_pos_grav]
+label = face_ypos
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Boundary tractions +y
+db_initial.iohandler.filename = tpv13_lateral_tractions.spatialdb
+db_initial.query_type = linear
+
+[pylithapp.timedependent.bc.y_neg_grav]
+label = face_yneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Boundary tractions -y
+db_initial.iohandler.filename = tpv13_lateral_tractions.spatialdb
+db_initial.query_type = linear
+
+[pylithapp.timedependent.bc.z_neg_grav]
+label = face_zneg
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Boundary tractions -z
+db_initial.iohandler.filename = tpv13_bottom_tractions.spatialdb
+db_initial.query_type = linear
+up_dir = [0, -1, 0]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+friction.db_properties.iohandler.filename = tpv12_friction.spatialdb
+friction.db_properties.query_type = nearest
+
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,dip_dir,normal_dir]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+#[pylithapp.timedependent.materials.elastic.output]
+#cell_info_fields = [density,mu,lambda,alpha_yield,beta,alpha_flow]
+#cell_data_fields = [total_strain,stress,plastic_strain]
+#time_step = 0.05*s
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_bottom_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_bottom_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_bottom_tractions.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,16 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = traction-shear-horiz traction-shear-vert traction-normal change-start-time
+ value-units = MPa MPa MPa year
+ num-locs = 2
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 0.0 -80.00000 0.0 0.0 -1332.8 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_initial_stress.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_initial_stress.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,19 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 6
+ value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
+ value-units = MPa MPa MPa MPa MPa MPa
+ num-locs = 5
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 0.0 -11.950 -69.60082 -69.60082 -199.08700 0.0 0.0 0.0
+0.0 0.0 -11.951 -134.35515 -134.35515 -199.10366 0.0 0.0 0.0
+0.0 0.0 -11.952 -199.12032 -199.12032 -199.12032 0.0 0.0 0.0
+0.0 0.0 -80.000 -1332.8 -1332.8 -1332.8 0.0 0.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_lateral_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_lateral_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_lateral_tractions.spatialdb 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,19 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = traction-shear-horiz traction-shear-vert traction-normal change-start-time
+ value-units = MPa MPa MPa year
+ num-locs = 5
+ data-dim = 1
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 3
+ }
+}
+0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 0.0 -11.950 0.0 0.0 -69.60082 0.0
+0.0 0.0 -11.951 0.0 0.0 -134.35515 0.0
+0.0 0.0 -11.952 0.0 0.0 -199.12032 0.0
+0.0 0.0 -80.00000 0.0 0.0 -1332.8 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4.cfg 2012-03-22 18:46:34 UTC (rev 19850)
@@ -0,0 +1,32 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/tpv13_tet4_200m.cfg 2012-03-22 18:46:34 UTC (rev 19850)
@@ -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.0125*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv13_tet4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_tet4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_tet4_200m-elastic.h5
More information about the CIG-COMMITS
mailing list