[cig-commits] r19573 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Thu Feb 2 09:44:27 PST 2012
Author: brad
Date: 2012-02-02 09:44:26 -0800 (Thu, 02 Feb 2012)
New Revision: 19573
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tractions.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tri3_200m.cfg
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/output/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tractions.spatialdb
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
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/quad4_050m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.cfg
Log:
Rework of tpv210-2d setup.
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/friction.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/friction.spatialdb 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/friction.spatialdb 2012-02-02 17:44:26 UTC (rev 19573)
@@ -1,16 +0,0 @@
-#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 = 2
- data-dim = 1
- space-dim = 2
- cs-data = cartesian {
- to-meters = 1.0e+3 // x, y are in km
- space-dim = 2
- }
-}
--7.4999 0.0 0.760 0.448 0.5 0.2
--7.5001 0.0 10000.0 0.448 0.5 1000.0
-
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou 2012-02-02 17:44:26 UTC (rev 19573)
@@ -1,7 +1,7 @@
# ----------------------------------------------------------------------
# Set vertex spacing with increasing spacing away from fault
# ----------------------------------------------------------------------
-#{bias_factor=1.015}
+#{bias_factor=1.02}
# Function for computing discretization size at end of curve with bias
#{fun_dx='dxStart*bias_factor**ceil( ln(1-curveL/dxStart*(1-bias_factor))/ln(bias_factor))'}
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-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/matprops.spatialdb 2012-02-02 17:44:26 UTC (rev 19573)
@@ -1,25 +1,14 @@
#SPATIAL.ascii 1
SimpleDB {
- // number of physical properties
num-values = 3
-
- // Names and units of physical properties
value-names = density vs vp
value-units = kg/m**3 m/s m/s
-
- // Number of locations where physical properties are given
- // 1 location -> data is uniform -> data dimension is 0
num-locs = 1
data-dim = 0
-
- // Problem is in 2-D
space-dim = 2
-
- // Coordinate system (2-D with coordinates in meters)
cs-data = cartesian {
to-meters = 1.0
space-dim = 2
}
}
-// Columns are x, y, density, vs, vp
-0.0 0.0 2700.0 3300.0 5716.0
+0.0 0.0 2700.0 3300.0 5716.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-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/pylithapp.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -7,19 +7,20 @@
# ----------------------------------------------------------------------
[pylithapp.journal.info]
timedependent = 1
-explicitlumped = 1
petsc = 1
+meshimporter = 1
meshiocubit = 1
+explicitlumped = 1
+explicitlumpedtri3 = 1
elasticityexplicit = 1
-meshimporter = 1
-#quadrature2d = 1
-#fiatlagrange = 1
+elasticityexplicittri3 = 1
faultcohesivedyn = 1
[pylithapp.journal.debug]
#pylithapp = 1
#problem = 1
-#explicit = 1
+#explicitlumped = 1
+#explicitlumpedtri3 = 1
# ----------------------------------------------------------------------
# mesh_generator
@@ -28,6 +29,11 @@
reader = pylith.meshio.MeshIOCubit
reorder_mesh = True
+distributor.partitioner = parmetis
+distributor.write_partition = True
+distributor.data_writer = pylith.meshio.DataWriterHDF5Mesh
+distributor.data_writer.filename = output/partition.h5
+
[pylithapp.mesh_generator.reader]
coordsys.space_dim = 2
@@ -37,44 +43,40 @@
[pylithapp.timedependent]
dimension = 2
-# Set interfaces to an array with 1 faults: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
formulation = pylith.problems.ExplicitLumped
-
-# Nondimensionalize problem using wave propagation parameters.
+formulation.norm_viscosity = 0.1
normalizer = spatialdata.units.NondimElasticDynamic
+normalizer.shear_wave_speed = 3333*m/s
-bc = [x_neg,x_pos,y_neg]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_neg = pylith.bc.AbsorbingDampers
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.y_neg = pylith.bc.AbsorbingDampers
-
[pylithapp.timedependent.formulation.time_step]
-total_time = 15.01*s
+#total_time = 15.01*s
+total_time = 0.0*s
dt = 0.01*s
+
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent]
-
materials = [elastic]
materials.elastic = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.elastic]
label = Elastic material
id = 1
+db_properties.label = Elastic properties
db_properties.iohandler.filename = matprops.spatialdb
+
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+bc = [x_neg,x_pos,y_neg]
+bc.x_neg = pylith.bc.AbsorbingDampers
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.y_neg = pylith.bc.AbsorbingDampers
+
[pylithapp.timedependent.bc.x_pos]
label = face_xpos
db.label = Absorbing BC +x
@@ -93,84 +95,53 @@
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+interfaces.fault = pylith.faults.FaultCohesiveDyn
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-fault = pylith.faults.FaultCohesiveDyn
-
[pylithapp.timedependent.interfaces.fault]
id = 10
label = fault
db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.label = Initial fault tractions
db_initial_tractions.query_type = linear
-db_initial_tractions.iohandler.filename = tractions.spatialdb
-db_initial_tractions.label = "Initial fault tractions"
friction = pylith.friction.SlipWeakening
+friction.label = Slip weakening
friction.db_properties = spatialdata.spatialdb.SimpleDB
-friction.db_properties.iohandler.filename = friction.spatialdb
friction.db_properties.label = Slip weakening
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
vertex_data_fields=[displacement,velocity]
output_freq = time_step
time_step = 0.04999*s
-writer.time_format = %06.3f
+writer = pylith.meshio.DataWriterHDF5Mesh
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
vertex_info_fields = [strike_dir,normal_dir,initial_traction]
vertex_data_fields = [slip,slip_rate,traction]
output_freq = time_step
time_step = 0.04999*s
-writer.time_format = %06.3f
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
-# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
+cell_info_fields = []
+cell_data_fields = []
cell_filter = pylith.meshio.CellFilterAvgMesh
output_freq = time_step
time_step = 15.0*s
-writer.time_format = %06.3f
+writer = pylith.meshio.DataWriterHDF5Mesh
# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
-pc_type = asm
-
-# Change the preconditioner settings.
-sub_pc_factor_shift_type = nonzero
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-15
-ksp_max_it = 100
-ksp_gmres_restart = 50
-snes_max_it = 200
-
-ksp_monitor = true
-ksp_view = true
-ksp_converged_reason = true
-
-snes_monitor = true
-snes_view = true
-snes_converged_reason = true
-
log_summary = true
-# Friction sensitivity solve
-friction_pc_type = asm
-friction_sub_pc_factor_shift_type = nonzero
-friction_ksp_max_it = 25
-friction_ksp_gmres_restart = 30
-
-friction_ksp_monitor = true
-friction_ksp_view = true
-friction_ksp_converged_reason = true
-
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou 2012-02-02 17:44:26 UTC (rev 19573)
@@ -14,6 +14,9 @@
surface all scheme pave
playback 'gradient.jou'
+# A to D
+curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -38,3 +41,4 @@
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.cfg 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -5,7 +5,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator.reader]
-filename = quad4_100m.exo
+filename = quad4_100m_refine.exo
# ----------------------------------------------------------------------
# problem
@@ -18,12 +18,12 @@
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-writer.filename = output/quad4_100m.vtk
+writer.filename = output/quad4_100m_refine.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/quad4_100m-fault.vtk
+writer.filename = output/quad4_100m_refine-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/quad4_100m-elastic.vtk
+writer.filename = output/quad4_100m_refine-elastic.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-02 17:44:26 UTC (rev 19573)
@@ -41,3 +41,4 @@
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.cfg 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -5,7 +5,7 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator.reader]
-filename = quad4_200m.exo
+filename = quad4_200m_refine.exo
# ----------------------------------------------------------------------
# problem
@@ -18,12 +18,12 @@
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-writer.filename = output/quad4_200m.vtk
+writer.filename = output/quad4_200m_refine.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/quad4_200m-fault.vtk
+writer.filename = output/quad4_200m_refine-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/quad4_200m-elastic.vtk
+writer.filename = output/quad4_200m_refine-elastic.vtk
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-02 17:44:26 UTC (rev 19573)
@@ -14,6 +14,9 @@
surface all scheme pave
playback 'gradient.jou'
+# A to D
+curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -37,3 +40,4 @@
# End of file
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -0,0 +1,13 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+db_initial_tractions.iohandler.filename = tpv10_tractions.spatialdb
+
+friction.db_properties.iohandler.filename = tpv10_friction.spatialdb
+friction.db_properties.query_type = nearest
+
+
+# End of file
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_friction.spatialdb (from rev 17829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/friction.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_friction.spatialdb 2012-02-02 17:44:26 UTC (rev 19573)
@@ -0,0 +1,16 @@
+#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 = 2
+ data-dim = 1
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+-7.499 0.0 0.760 0.448 0.5 0.2
+-7.500 0.0 10000.0 0.448 0.5 0.2
+
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tractions.spatialdb (from rev 17829, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tractions.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tractions.spatialdb 2012-02-02 17:44:26 UTC (rev 19573)
@@ -0,0 +1,19 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = traction-shear traction-normal
+ value-units = MPa 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
+ }
+}
+ 0.0 0.0 0.0 0.0
+-5.2499 0.0 42.60876158 -77.4675244
+-5.2500 0.0 59.51801330 -77.4690000
+-6.7500 0.0 76.46601710 -99.6030000
+-6.7501 0.0 54.78246158 -99.6044756
+-7.5000 0.0 60.86850000 -110.6700000
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tri3_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tri3_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tpv10_tri3_200m.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -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/tpv10_tri3_200m.h5
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/tpv10_tri3_200m-fault.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tpv10_tri3_200m-elastic.h5
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tractions.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tractions.spatialdb 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tractions.spatialdb 2012-02-02 17:44:26 UTC (rev 19573)
@@ -1,19 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 2
- value-names = traction-shear traction-normal
- value-units = MPa 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
- }
-}
-0.0 0.0 0.0 0.0
--5.249500000000000 -9.092400714332822 42.603892100000 -77.461622000000
--5.250000000000000 -9.093266739736606 59.518013300000 -77.469000000000
--6.750000000000000 -11.691342951089920 76.466017100000 -99.603000000000
--6.750500000000000 -11.692208976493704 54.785707900000 -99.610378000000
--7.500000000000000 -12.990381056766578 60.868500000000 -110.670000000000
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-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -2,27 +2,39 @@
[pylithapp]
# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+formulation = pylith.problems.ExplicitLumpedTri3
+
+
+# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent.materials.elastic]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
+
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.x_pos]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
+quadrature.cell.quad_order = 1
[pylithapp.timedependent.bc.x_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
+quadrature.cell.quad_order = 1
[pylithapp.timedependent.bc.y_neg]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
+quadrature.cell.quad_order = 1
+
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
@@ -30,3 +42,5 @@
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
+
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.cfg 2012-02-02 17:00:33 UTC (rev 19572)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/tri3_100m.cfg 2012-02-02 17:44:26 UTC (rev 19573)
@@ -5,19 +5,19 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator.reader]
-filename = tri3_100m.exo
+filename = tri3_100m_refine.exo
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
-writer.filename = output/tri3_100m.vtk
+writer.filename = output/tri3_100m_refine.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = output/tri3_100m-fault.vtk
+writer.filename = output/tri3_100m_refine-fault.vtk
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = output/tri3_100m-elastic.vtk
+writer.filename = output/tri3_100m_refine-elastic.vtk
More information about the CIG-COMMITS
mailing list