[cig-commits] r15173 - in short/3D/PyLith/trunk/examples/bar_shearwave: . hex8 quad4 tet4 tri3
brad at geodynamics.org
brad at geodynamics.org
Wed Jun 10 11:44:04 PDT 2009
Author: brad
Date: 2009-06-10 11:44:03 -0700 (Wed, 10 Jun 2009)
New Revision: 15173
Added:
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
Removed:
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/shearwave.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/shearwave.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/shearwave.cfg
Modified:
short/3D/PyLith/trunk/examples/bar_shearwave/README
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/README
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/README
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/Makefile.am
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/README
Log:
Renamed shearwave.cfg to pylithapp.cfg. Cleanup.
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/README 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/README 2009-06-10 18:44:03 UTC (rev 15173)
@@ -8,7 +8,7 @@
The problem can be run by entering
- pylith shearwave.cfg
+ pylith
in any of the subdirectories.
@@ -23,7 +23,7 @@
Each of the subdirectories contains the files listed below (along with
other files for the mesh).
- - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
- matprops.spatialdb: Material database of physical properties
- shearwave_slip.spatialdb: Spatial database with values for final slip
- shearwave_sliprate.spatialdb: Spatial database with values for the
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/Makefile.am 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/Makefile.am 2009-06-10 18:44:03 UTC (rev 15173)
@@ -19,7 +19,7 @@
geometry.jou \
matprops.spatialdb \
mesh_hex8_200m.jou \
- shearwave.cfg \
+ pylithapp.cfg \
shearwave_slip.spatialdb \
shearwave_risetime.spatialdb \
shearwave_sliptime.spatialdb
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/README 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/README 2009-06-10 18:44:03 UTC (rev 15173)
@@ -4,7 +4,7 @@
- geometry.jou: CUBIT journal file defining the geometry of the mesh
- mesh_hex8_200m.jou: CUBIT journal file for meshing the problem
- - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
- matprops.spatialdb: Material database of physical properties
- bar_hex8_200m.exo: Exodus file with mesh generated by CUBIT
- shearwave_slip.spatialdb: Spatial database with values for final slip
@@ -15,4 +15,4 @@
To run the problem, enter
- pylith shearwave.cfg
+ pylith
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg (from rev 15170, short/3D/PyLith/trunk/examples/bar_shearwave/hex8/shearwave.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -0,0 +1,215 @@
+# -*- Python -*-
+
+# These settings define an explicit time-stepping problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#explicit = 1
+#petsc = 1
+#solverlinear = 1
+meshiocubit = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = bar_hex8_200m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+
+# Change to an explicit time stepping formulation
+formulation = pylith.problems.Explicit
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+# Set interfaces to an array with 1 fault: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 1 material 'elastic'.
+materials = [elastic]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 12.0*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+
+# Label for material
+label = Elastic material
+
+# Nodeset id from CUBIT
+id = 1
+
+# Spatial database with physical properties for elastic material
+db_properties.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 3-D Lagrange cell with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing boundary condition on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D Lagrange cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing boundary condition on -x face of bar
+
+# Nodeset id from CUBIT
+label = 12
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D Lagrange cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+# Dirichlet boundary condition on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+fixed_dof = [0, 2]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10'). Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements. The default ID for a fault is 100. If you have a
+# fault in your mesh you must:
+# 1. If you create your own fault ID, make sure it does not conflict
+# with any of you material ID's.
+# 2. If you use the default fault ID, make sure that none of your
+# material ID's are equal to 100.
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# Specify the basis functions and quadrature:
+# 2-D Lagrange cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# Switch to Brune slip time function
+eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = shearwave_slip.spatialdb
+
+# Database specifying rise time.
+rise_time.iohandler.filename = shearwave_risetime.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+# Change the preconditioner settings (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+#ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# We set three output manager parameters for each type of output.
+# skip - Number of time steps to skip between output
+# writer.filename - Root name for VTK file
+# writer.time_format - Format of time stamp (decimal point is removed)
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+skip = 3
+writer.filename = output/shearwave.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 3
+writer.filename = output/shearwave-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+skip = 33
+writer.filename = output/shearwave-statevars.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvgMesh
Property changes on: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
___________________________________________________________________
Name: svn:mergeinfo
+
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/shearwave.cfg 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/shearwave.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -1,222 +0,0 @@
-# -*- Python -*-
-
-# The settings in this file (shearwave.cfg) will be read if it is
-# specified on the command line:
-# 'pylith shearwave.cfg'
-#
-# If this file is given, the settings in this file will override all
-# others except for parameters specified directly on the command line.
-#
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-#solverlinear = 1
-meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-filename = bar_hex8_200m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 3
-
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
-bc = [x_pos,x_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.elastic]
-
-# Label for material
-label = Elastic material
-
-# Nodeset id from CUBIT
-id = 1
-
-# Spatial database with physical properties for elastic material
-properties_db.iohandler.filename = matprops.spatialdb
-
-# Set the basis functions and quadrature:
-# 3-D Lagrange cell with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 3
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 2-D Lagrange cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 2-D Lagrange cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos]
-# Dirichlet boundary condition on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-fixed_dof = [0, 2]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# Specify the basis functions and quadrature:
-# 2-D Lagrange cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
-
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
-
-# Database specifying rise time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero =
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-10
-ksp_max_it = 50
-ksp_gmres_restart = 10
-
-#ksp_monitor = true
-#ksp_view = true
-log_summary = true
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-skip = 3
-writer.filename = output/shearwave.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.filename = output/shearwave-fault.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-skip = 33
-writer.filename = output/shearwave-statevars.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/Makefile.am 2009-06-10 18:44:03 UTC (rev 15173)
@@ -19,7 +19,7 @@
geometry.jou \
matprops.spatialdb \
mesh_quad4_200m.jou \
- shearwave.cfg \
+ pylithapp.cfg \
shearwave_slip.spatialdb \
shearwave_risetime.spatialdb \
shearwave_sliptime.spatialdb
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/README 2009-06-10 18:44:03 UTC (rev 15173)
@@ -4,7 +4,7 @@
- geometry.jou: CUBIT journal file defining the geometry of the mesh
- mesh_quad4_200m.jou: CUBIT journal file for meshing the problem
- - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
- matprops.spatialdb: Material database of physical properties
- bar_quad4_200m.exo: Exodus file with mesh generated by CUBIT
- shearwave_slip.spatialdb: Spatial database with values for final slip
@@ -15,4 +15,4 @@
To run the problem, enter
- pylith shearwave.cfg
+ pylith
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg (from rev 15170, short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -0,0 +1,219 @@
+# -*- Python -*-
+
+# These settings define an explicit time-stepping problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#explicit = 1
+#petsc = 1
+#solverlinear = 1
+meshiocubit = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = bar_quad4_200m.exo
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+# Change to an explicit time stepping formulation
+formulation = pylith.problems.Explicit
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+# Set interfaces to an array with 1 fault: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 1 material 'elastic'.
+materials = [elastic]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 12.0*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+
+# Label for material
+label = Elastic material
+
+# Nodeset id from CUBIT
+id = 1
+
+# Spatial database with physical properties for elastic material
+db_properties.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 2-D Lagrange cell with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing boundary condition on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing boundary condition on -x face of bar
+
+# Nodeset id from CUBIT
+label = 12
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+# Dirichlet boundary condition on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+fixed_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10'). Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements. The default ID for a fault is 100. If you have a
+# fault in your mesh you must:
+# 1. If you create your own fault ID, make sure it does not conflict
+# with any of you material ID's.
+# 2. If you use the default fault ID, make sure that none of your
+# material ID's are equal to 100.
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+# Switch to Brune slip time function
+eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = shearwave_slip.spatialdb
+
+# Database specifying rise_time.
+rise_time.iohandler.filename = shearwave_risetime.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+# Change the preconditioner settings (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+#ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# We set three output manager parameters for each type of output.
+# skip - Number of time steps to skip between output
+# writer.filename - Root name for VTK file
+# writer.time_format - Format of time stamp (decimal point is removed)
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+skip = 3
+writer.filename = output/shearwave.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 3
+writer.filename = output/shearwave-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+skip = 33
+writer.filename = output/shearwave-statevars.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvgMesh
Property changes on: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
___________________________________________________________________
Name: svn:mergeinfo
+
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave.cfg 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -1,226 +0,0 @@
-# -*- Python -*-
-
-# The settings in this file (shearwave.cfg) will be read if it is
-# specified on the command line:
-# 'pylith shearwave.cfg'
-#
-# If this file is given, the settings in this file will override all
-# others except for parameters specified directly on the command line.
-#
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-#solverlinear = 1
-meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-filename = bar_quad4_200m.exo
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 2
-
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
-bc = [x_pos,x_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-elastic = pylith.materials.ElasticPlaneStrain
-
-[pylithapp.timedependent.materials.elastic]
-
-# Label for material
-label = Elastic material
-
-# Nodeset id from CUBIT
-id = 1
-
-# Spatial database with physical properties for elastic material
-properties_db.iohandler.filename = matprops.spatialdb
-
-# Set the basis functions and quadrature:
-# 2-D Lagrange cell with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos]
-# Dirichlet boundary condition on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-fixed_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# Specify the basis functions and quadrature:
-# 1-D Lagrange cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
-
-# Database specifying rise_time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero =
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-10
-ksp_max_it = 50
-ksp_gmres_restart = 10
-
-#ksp_monitor = true
-#ksp_view = true
-log_summary = true
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-skip = 3
-writer.filename = output/shearwave.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.filename = output/shearwave-fault.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-skip = 33
-writer.filename = output/shearwave-statevars.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/Makefile.am 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/Makefile.am 2009-06-10 18:44:03 UTC (rev 15173)
@@ -19,7 +19,7 @@
bar_tet4_200m.pset \
matprops.spatialdb \
mesh_tet4_200m.lagrit \
- shearwave.cfg \
+ pylithapp.cfg \
shearwave_slip.spatialdb \
shearwave_risetime.spatialdb \
shearwave_sliptime.spatialdb
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/README 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/README 2009-06-10 18:44:03 UTC (rev 15173)
@@ -3,7 +3,7 @@
Files associated with the example problem:
- mesh_tet4_200m.lagrit: LaGriT script file for meshing the problem
- - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
- matprops.spatialdb: Material database of physical properties
- bar_tet4_200m.gmv: GMV file with mesh generated by LaGriT
- bar_tet4_200m.pset: Pset file with point sets for mesh generated by LaGriT
@@ -15,4 +15,4 @@
To run the problem, enter
- pylith shearwave.cfg
+ pylith
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg (from rev 15170, short/3D/PyLith/trunk/examples/bar_shearwave/tet4/shearwave.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -0,0 +1,241 @@
+# -*- Python -*-
+
+# These settings define an explicit time-stepping problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#explicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiolagrit = 1
+#homogeneous = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOLagrit
+
+[pylithapp.mesh_generator.reader]
+filename_gmv = bar_tet4_200m.gmv
+filename_pset = bar_tet4_200m.pset
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+
+# Change to an explicit time stepping formulation
+formulation = pylith.problems.Explicit
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+# Set interfaces to an array with 1 fault: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 2 materials 'pos' and 'neg'.
+materials = [pos,neg]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 12.0*s
+dt = 0.0333333333333333333*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.pos]
+
+# Label for material
+label = elastic material +x
+
+# Material id (itetclr) from LaGriT
+id = 1
+
+# Spatial database with physical properties for elastic material
+db_properties.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 3-D simplex cell with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.materials.neg]
+
+# Label for material
+label = elastic material -x
+
+# Material id (itetclr) from LaGriT
+id = 2
+
+# Spatial database with physical properties for elastic material
+db_properties.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 3-D simplex cell with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing boundary condition on +x face of bar
+
+# Pset name from LaGriT
+label = boundary_xp
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing boundary condition on +x face of bar
+
+# Pset name from LaGriT
+label = boundary_xm
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+# Dirichlet boundary condition on all vertices except fault vertices
+
+fixed_dof = [0, 2]
+
+# Pset name from LaGriT
+label = pts_notfault
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10'). Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements. The default ID for a fault is 100. If you have a
+# fault in your mesh you must:
+# 1. If you create your own fault ID, make sure it does not conflict
+# with any of you material ID's.
+# 2. If you use the default fault ID, make sure that none of your
+# material ID's are equal to 100.
+
+[pylithapp.timedependent.interfaces.fault]
+# Material id
+id = 100
+
+# Pset name from LaGriT
+label = fault
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+# Switch to Brune slip time function
+eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = shearwave_slip.spatialdb
+
+# Database specifying rise time.
+rise_time.iohandler.filename = shearwave_risetime.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+# Change the preconditioner settings (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+#ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# We set three output manager parameters for each type of output.
+# skip - Number of time steps to skip between output
+# writer.filename - Root name for VTK file
+# writer.time_format - Format of time stamp (decimal point is removed)
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+skip = 3
+writer.filename = output/shearwave.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 3
+writer.filename = output/shearwave-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.pos.output]
+skip = 33
+writer.filename = output/shearwave-statevars-pos.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvgMesh
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.neg.output]
+skip = 33
+writer.filename = output/shearwave-statevars-neg.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvgMesh
Property changes on: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
___________________________________________________________________
Name: svn:mergeinfo
+
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/shearwave.cfg 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/shearwave.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -1,248 +0,0 @@
-# -*- Python -*-
-
-# The settings in this file (shearwave.cfg) will be read if it is
-# specified on the command line:
-# 'pylith shearwave.cfg'
-#
-# If this file is given, the settings in this file will override all
-# others except for parameters specified directly on the command line.
-#
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-#solverlinear = 1
-#meshiolagrit = 1
-#homogeneous = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOLagrit
-
-[pylithapp.mesh_generator.reader]
-filename_gmv = bar_tet4_200m.gmv
-filename_pset = bar_tet4_200m.pset
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 3
-
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
-bc = [x_pos,x_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 2 materials 'pos' and 'neg'.
-materials = [pos,neg]
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.0333333333333333333*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.pos]
-
-# Label for material
-label = elastic material +x
-
-# Material id (itetclr) from LaGriT
-id = 1
-
-# Spatial database with physical properties for elastic material
-properties_db.iohandler.filename = matprops.spatialdb
-
-# Set the basis functions and quadrature:
-# 3-D simplex cell with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.materials.neg]
-
-# Label for material
-label = elastic material -x
-
-# Material id (itetclr) from LaGriT
-id = 2
-
-# Spatial database with physical properties for elastic material
-properties_db.iohandler.filename = matprops.spatialdb
-
-# Set the basis functions and quadrature:
-# 3-D simplex cell with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Pset name from LaGriT
-label = boundary_xp
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 2-D simplex cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on +x face of bar
-
-# Pset name from LaGriT
-label = boundary_xm
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 2-D simplex cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos]
-# Dirichlet boundary condition on all vertices except fault vertices
-
-fixed_dof = [0, 2]
-
-# Pset name from LaGriT
-label = pts_notfault
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-[pylithapp.timedependent.interfaces.fault]
-# Material id
-id = 100
-
-# Pset name from LaGriT
-label = fault
-
-# Specify the basis functions and quadrature:
-# 2-D simplex cell in 3-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
-
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
-
-# Database specifying rise time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero =
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-10
-ksp_max_it = 50
-ksp_gmres_restart = 10
-
-#ksp_monitor = true
-#ksp_view = true
-log_summary = true
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-skip = 3
-writer.filename = output/shearwave.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.filename = output/shearwave-fault.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.pos.output]
-skip = 33
-writer.filename = output/shearwave-statevars-pos.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.neg.output]
-skip = 33
-writer.filename = output/shearwave-statevars-neg.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/Makefile.am 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/Makefile.am 2009-06-10 18:44:03 UTC (rev 15173)
@@ -19,7 +19,7 @@
geometry.jou \
matprops.spatialdb \
mesh_tri3_200m.jou \
- shearwave.cfg \
+ pylithapp.cfg \
shearwave_slip.spatialdb \
shearwave_risetime.spatialdb \
shearwave_sliptime.spatialdb
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/README
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/README 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/README 2009-06-10 18:44:03 UTC (rev 15173)
@@ -4,7 +4,7 @@
- geometry.jou: CUBIT journal file defining the geometry of the mesh
- mesh_tri3_200m.jou: CUBIT journal file for meshing the problem
- - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+ - pylithapp.cfg: Parameter settings for the shear wave dislocation problem
- matprops.spatialdb: Material database of physical properties
- bar_tri3_200m.exo: Exodus file with mesh generated by CUBIT
- shearwave_slip.spatialdb: Spatial database with values for final slip
@@ -15,4 +15,4 @@
To run the problem, enter
- pylith shearwave.cfg
+ pylith
Copied: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg (from rev 15170, short/3D/PyLith/trunk/examples/bar_shearwave/tri3/shearwave.cfg)
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -0,0 +1,219 @@
+# -*- Python -*-
+
+# These settings define an explicit time-stepping problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#explicit = 1
+#petsc = 1
+#solverlinear = 1
+meshiocubit = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = bar_tri3_200m.exo
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+
+# Change to an explicit time stepping formulation
+formulation = pylith.problems.Explicit
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+# Set interfaces to an array with 1 fault: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 1 material 'elastic'.
+materials = [elastic]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 12.0*s
+dt = 0.0333333333333333*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials]
+elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+
+# Label for material
+label = Elastic material
+
+# Nodeset id from CUBIT
+id = 1
+
+# Spatial database with physical properties for elastic material
+db_properties.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 2-D simplex cell with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing boundary condition on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D simplex cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing boundary condition on -x face of bar
+
+# Nodeset id from CUBIT
+label = 12
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 1-D simplex cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+# Dirichlet boundary condition on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+fixed_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10'). Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements. The default ID for a fault is 100. If you have a
+# fault in your mesh you must:
+# 1. If you create your own fault ID, make sure it does not conflict
+# with any of you material ID's.
+# 2. If you use the default fault ID, make sure that none of your
+# material ID's are equal to 100.
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# Specify the basis functions and quadrature:
+# 1-D simplex cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+quadrature.cell.quad_order = 2
+
+# Switch to Brune slip time function
+eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = shearwave_slip.spatialdb
+
+# Database specifying rise time.
+rise_time.iohandler.filename = shearwave_risetime.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+# Change the preconditioner settings (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-10
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+#ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# We set three output manager parameters for each type of output.
+# skip - Number of time steps to skip between output
+# writer.filename - Root name for VTK file
+# writer.time_format - Format of time stamp (decimal point is removed)
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+skip = 3
+writer.filename = output/shearwave.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 3
+writer.filename = output/shearwave-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+skip = 33
+writer.filename = output/shearwave-statevars.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvgMesh
Property changes on: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
___________________________________________________________________
Name: svn:mergeinfo
+
Deleted: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/shearwave.cfg 2009-06-10 18:19:15 UTC (rev 15172)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/shearwave.cfg 2009-06-10 18:44:03 UTC (rev 15173)
@@ -1,226 +0,0 @@
-# -*- Python -*-
-
-# The settings in this file (shearwave.cfg) will be read if it is
-# specified on the command line:
-# 'pylith shearwave.cfg'
-#
-# If this file is given, the settings in this file will override all
-# others except for parameters specified directly on the command line.
-#
-# These settings define an explicit time-stepping problem using a fault with
-# kinematically-specified slip.
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-#explicit = 1
-#petsc = 1
-#solverlinear = 1
-meshiocubit = 1
-#explicitelasticity = 1
-#quadrature2d = 1
-#fiatlagrange = 1
-#faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-reader = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.reader]
-filename = bar_tri3_200m.exo
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-dimension = 2
-
-# Change to an explicit time stepping formulation
-formulation = pylith.problems.Explicit
-
-# Nondimensionalize problem using wave propagation parameters.
-normalizer = spatialdata.units.NondimElasticDynamic
-
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
-bc = [x_pos,x_neg,y_pos]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-# Set interfaces to an array with 1 fault: 'fault'.
-interfaces = [fault]
-
-# Set materials to an array with 1 material 'elastic'.
-materials = [elastic]
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.0333333333333333*s
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-elastic = pylith.materials.ElasticPlaneStrain
-
-[pylithapp.timedependent.materials.elastic]
-
-# Label for material
-label = Elastic material
-
-# Nodeset id from CUBIT
-id = 1
-
-# Spatial database with physical properties for elastic material
-properties_db.iohandler.filename = matprops.spatialdb
-
-# Set the basis functions and quadrature:
-# 2-D simplex cell with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing boundary condition on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing boundary condition on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos]
-# Dirichlet boundary condition on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-fixed_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 100
-
-# Nodeset id from CUBIT
-label = 10
-
-# Specify the basis functions and quadrature:
-# 1-D simplex cell in 2-D space with 2nd order quadrature
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-quadrature.cell.quad_order = 2
-
-# Switch to Brune slip time function
-eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
-
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-
-# Database specifying the final slip.
-slip.iohandler.filename = shearwave_slip.spatialdb
-
-# Database specifying rise time.
-rise_time.iohandler.filename = shearwave_risetime.spatialdb
-
-# Database specifying time at which slip begins at each point.
-slip_time.iohandler.filename = shearwave_sliptime.spatialdb
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-ksp_type = gmres
-pc_type = asm
-# Change the preconditioner settings (must turn off
-# shift_positive_definite and turn on shift_nonzero).
-sub_pc_factor_shift_positive_definite = 0
-sub_pc_factor_shift_nonzero =
-
-ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-10
-ksp_max_it = 50
-ksp_gmres_restart = 10
-
-#ksp_monitor = true
-#ksp_view = true
-log_summary = true
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# We set three output manager parameters for each type of output.
-# skip - Number of time steps to skip between output
-# writer.filename - Root name for VTK file
-# writer.time_format - Format of time stamp (decimal point is removed)
-
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-skip = 3
-writer.filename = output/shearwave.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-skip = 3
-writer.filename = output/shearwave-fault.vtk
-writer.time_format = %05.2f
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-skip = 33
-writer.filename = output/shearwave-statevars.vtk
-writer.time_format = %05.2f
-
-# Averate state variables over cell
-cell_filter = pylith.meshio.CellFilterAvgMesh
More information about the CIG-COMMITS
mailing list