[cig-commits] r19732 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Tue Mar 6 14:17:07 PST 2012
Author: brad
Date: 2012-03-06 14:17:07 -0800 (Tue, 06 Mar 2012)
New Revision: 19732
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_050m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_200m.cfg
Modified:
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/tpv12_quad4_200m.cfg
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_xface_tractions.spatialdb
Log:
Revert to sharp transition in stress field for tpv13. Added quad4 cases.
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-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_offfaultdata.py 2012-03-06 22:17:07 UTC (rev 19732)
@@ -8,7 +8,7 @@
# ----------------------------------------------------------------------
#
-sim = "tpv12"
+sim = "tpv13"
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-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tabulate_onfaultdata.py 2012-03-06 22:17:07 UTC (rev 19732)
@@ -9,8 +9,8 @@
#
sim = "tpv13"
-cell = "tri3"
-dx = 100
+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/tpv12_quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg 2012-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv12_quad4_200m.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg 2012-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -73,10 +73,10 @@
# ----------------------------------------------------------------------
# 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
+#[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
Modified: 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 2012-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_state.spatialdb 2012-03-06 22:17:07 UTC (rev 19732)
@@ -12,8 +12,8 @@
space-dim = 2
}
}
-0.0 0.0 0.0 0.0 0.0 0.0 0.0
-0.0 -11.950 -134.34391 0.0 0.0 0.0 0.0
-0.0 -11.951 -166.72940 0.0 0.0 0.0 0.0
-0.0 -11.952 -199.12032 0.0 0.0 0.0 0.0
-0.0 -80.000 -1332.8 0.0 0.0 0.0 0.0
+0.0 0.0 0.0 0.0 0.0 0.0 0.0
+0.0 -11.950 -134.34391 0.0 0.0 0.0 0.0
+0.0 -11.951 -166.72940 0.0 0.0 0.0 0.0
+0.0 -11.952 -199.12032 0.0 0.0 0.0 0.0
+0.0 -80.000 -1332.8 0.0 0.0 0.0 0.0
Modified: 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 2012-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_initial_stress.spatialdb 2012-03-06 22:17:07 UTC (rev 19732)
@@ -16,4 +16,4 @@
0.0 -11.950 -69.60082 -199.08700 0.0
0.0 -11.951 -134.35515 -199.10366 0.0
0.0 -11.952 -199.12032 -199.12032 0.0
-0.0 -80.00000 -1332.8 -1332.8 0.0
+0.0 -80.000 -1332.8 -1332.8 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -0,0 +1,23 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos_grav]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg_grav]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg_grav]
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_050m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_050m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_050m.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -0,0 +1,28 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_050m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.15*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.0025*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv13_quad4_050m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_quad4_050m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_quad4_050m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_100m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_100m.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -0,0 +1,28 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_100m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.normalizer]
+wave_period = 0.3*s
+
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.005*s
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/tpv13_quad4_100m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_quad4_100m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_quad4_100m-elastic.h5
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_quad4_200m.cfg 2012-03-06 22:17:07 UTC (rev 19732)
@@ -0,0 +1,28 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_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_quad4_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv13_quad4_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv13_quad4_200m-elastic.h5
Modified: 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 2012-03-06 19:55:02 UTC (rev 19731)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv13_xface_tractions.spatialdb 2012-03-06 22:17:07 UTC (rev 19732)
@@ -4,7 +4,7 @@
num-values = 3
value-names = traction-shear traction-normal change-start-time
value-units = MPa MPa year
- num-locs = 4
+ num-locs = 5
data-dim = 1
space-dim = 2
cs-data = cartesian {
@@ -13,6 +13,7 @@
}
}
0.0 0.0 0.0 0.0 0.0
-0.0 -11.95115 0.0 -69.60751 0.0
-0.0 -11.95116 0.0 -199.10633 0.0
+0.0 -11.950 0.0 -69.60082 0.0
+0.0 -11.951 0.0 -134.35515 0.0
+0.0 -11.952 0.0 -199.12032 0.0
0.0 -80.00000 0.0 -1332.8 0.0
More information about the CIG-COMMITS
mailing list