[cig-commits] r19603 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Wed Feb 8 17:15:51 PST 2012
Author: brad
Date: 2012-02-08 17:15:50 -0800 (Wed, 08 Feb 2012)
New Revision: 19603
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_state.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_stress.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_050m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_xface_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_yface_tractions.spatialdb
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/matprops.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg
Log:
Started working on tpv13.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/matprops.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/matprops.spatialdb 2012-02-08 06:56:57 UTC (rev 19602)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/matprops.spatialdb 2012-02-09 01:15:50 UTC (rev 19603)
@@ -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 = 2
@@ -11,4 +11,4 @@
space-dim = 2
}
}
-0.0 0.0 2700.0 3300.0 5716.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-2d/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg 2012-02-08 06:56:57 UTC (rev 19602)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -15,6 +15,7 @@
elasticityexplicit = 1
elasticityexplicittri3 = 1
faultcohesivedyn = 1
+fiatsimplex = 1
[pylithapp.journal.debug]
#pylithapp = 1
@@ -102,9 +103,6 @@
id = 10
label = fault
-db_initial_tractions = spatialdata.spatialdb.SimpleDB
-db_initial_tractions.label = Initial fault tractions
-
friction = pylith.friction.SlipWeakening
friction.label = Slip weakening
@@ -121,7 +119,7 @@
writer = pylith.meshio.DataWriterHDF5Mesh
[pylithapp.timedependent.interfaces.fault.output]
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+vertex_info_fields = [strike_dir,normal_dir]
vertex_data_fields = [slip,slip_rate,traction]
output_freq = time_step
time_step = 0.04999*s
@@ -139,7 +137,7 @@
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
-log_summary = true
+#log_summary = true
# End of file
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-08 06:56:57 UTC (rev 19602)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -5,6 +5,8 @@
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
db_initial_tractions.iohandler.filename = tpv10_tractions.spatialdb
db_initial_tractions.query_type = linear
@@ -12,4 +14,11 @@
friction.db_properties.query_type = nearest
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+
+
# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg 2012-02-08 06:56:57 UTC (rev 19602)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -12,6 +12,8 @@
# 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
@@ -19,4 +21,11 @@
friction.db_properties.query_type = nearest
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+
+
# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,72 @@
+[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, -1.0, 0.0]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 8.00001*s
+
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+#[pylithapp.timedependent.materials]
+#elastic = pylith.materials.DruckerPragerPlaneStrain
+
+[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
+
+#db_initial_state = spatialdata.spatialdb.SimpleDB
+#db_initial_state.label = Initial state variables
+#db_initial_state.iohandler.filename = tpv13_initial_state.spatialdb
+#db_initial_state.query_type = linear
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_neg,x_pos,y_neg,x_neg_grav,x_pos_grav,y_neg_grav]
+bc.x_neg_grav = pylith.bc.Neumann
+bc.x_pos_grav = pylith.bc.Neumann
+bc.y_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_xface_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_xface_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_yface_tractions.spatialdb
+db_initial.query_type = linear
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+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/tpv13_initial_state.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_state.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_state.spatialdb 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,18 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 5
+ value-names = stress-zz-initial plastic-strain-xx plastic-strain-yy plastic-strain-xy plastic-strain-zz
+ value-units = MPa None None None None
+ num-locs = 4
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 2
+ }
+}
+0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 -11.95115 -134.35684 0.0 0.0 0.0 0.0
+0.0 -11.95116 -199.10633 0.0 0.0 0.0 0.0
+0.0 -80.00000 -1332.8 0.0 0.0 0.0 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_stress.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_stress.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_stress.spatialdb 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,16 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = stress-xx stress-yy stress-xy
+ value-units = MPa MPa MPa
+ num-locs = 2
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 2
+ }
+}
+0.0 0.0 0.0 0.0 0.0
+0.0 -80.00000 -465.94688 -1332.8 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,23 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg_grav]
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_050m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_050m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_050m.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -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/tpv13_tri3_050m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_tri3_050m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_tri3_050m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_100m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_100m.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -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/tpv13_tri3_100m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_tri3_100m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_tri3_100m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_tri3_200m.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -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.01*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv13_tri3_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_tri3_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_tri3_200m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_xface_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_xface_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_xface_tractions.spatialdb 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,16 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear traction-normal change-start-time
+ value-units = MPa MPa year
+ num-locs = 2
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 2
+ }
+}
+0.0 0.0 0.0 0.0 0.0
+0.0 -80.00000 0.0 -465.94688 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_yface_tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_yface_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_yface_tractions.spatialdb 2012-02-09 01:15:50 UTC (rev 19603)
@@ -0,0 +1,16 @@
+// -*- C++ -*-
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear traction-normal change-start-time
+ value-units = MPa MPa year
+ num-locs = 2
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3
+ space-dim = 2
+ }
+}
+0.0 0.0 0.0 0.0 0.0
+0.0 -80.00000 0.0 -1332.8 0.0
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg 2012-02-08 06:56:57 UTC (rev 19602)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg 2012-02-09 01:15:50 UTC (rev 19603)
@@ -7,7 +7,6 @@
[pylithapp.timedependent]
formulation = pylith.problems.ExplicitLumpedTri3
-
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
@@ -22,17 +21,17 @@
[pylithapp.timedependent.bc.x_pos]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.x_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
# ----------------------------------------------------------------------
More information about the CIG-COMMITS
mailing list