[cig-commits] r12883 - short/2.5D/benchmarks/savageprescott/parameters
willic3 at geodynamics.org
willic3 at geodynamics.org
Mon Sep 15 13:23:39 PDT 2008
Author: willic3
Date: 2008-09-15 13:23:39 -0700 (Mon, 15 Sep 2008)
New Revision: 12883
Added:
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh
Removed:
short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750.cfg
Modified:
short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
Log:
Added shell scripts to run the examples, fixed the uniform resolution case
to reflect the revised resolution, and added PETSc option to monitor
singular values.
Modified: short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg 2008-09-15 20:07:52 UTC (rev 12882)
+++ short/2.5D/benchmarks/savageprescott/parameters/pylithapp.cfg 2008-09-15 20:23:39 UTC (rev 12883)
@@ -78,5 +78,6 @@
ksp_view = true
log_summary = true
ksp_max_it = 100
-ksp_gmres_restart = 50
+ksp_gmres_restart = 70
+ksp_monitor_singular_value = true
# start_in_debugger = true
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh 2008-09-15 20:23:39 UTC (rev 12883)
@@ -0,0 +1 @@
+pylith spbm_hex8_graded_1000m.cfg timedep.cfg
Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_1000m.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh 2008-09-15 20:23:39 UTC (rev 12883)
@@ -0,0 +1 @@
+pylith spbm_hex8_graded_3000m.cfg timedep.cfg
Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_graded_3000m.sh
___________________________________________________________________
Name: svn:executable
+ *
Deleted: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750.cfg 2008-09-15 20:07:52 UTC (rev 12882)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750.cfg 2008-09-15 20:23:39 UTC (rev 12883)
@@ -1,127 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the boundary and
-# interface conditions. The general parameters are specificed in the
-# pylithapp.cfg file which PyLith reads by default. Information controlling
-# time stepping is contained in the file timedep.cfg.
-#
-# To run the simulation:
-# pylith spbm_hex8_unif_8000m.cfg timedep.cfg
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-# Change the mesh importer to the CUBIT importer.
-importer = pylith.meshio.MeshIOCubit
-
-[pylithapp.mesh_generator.importer]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_unif_8000m.exo
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
-bc = [x_pos,x_neg,z_neg]
-
-# Set interfaces to an array of 1 fault: 'fault'.
-interfaces = [fault]
-
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-
-# We change spatial database for the rate of change from FixedDOFDB
-# (which has zero velocity) to UniformDB (which allows nonzero
-# velocities).
-
-# Set the parameters for the desired boundary conditions.
-#
-# The label corresponds to the nodeset ID in CUBIT.
-
-[pylithapp.timedependent.bc.x_pos]
-fixed_dof = [0, 1]
-label = 11
-db.label = Dirichlet BC on +x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on +x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of -1 cm/year.
-rate_db.data = [0.0,-3.1688087814e-10]
-
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = 12
-db.label = Dirichlet BC on -x
-rate_db = spatialdata.spatialdb.UniformDB
-rate_db.label = Dirichlet rate BC on -x
-rate_db.values = [dof-0,dof-1]
-# This rate (in m/s) gives a value of 1 cm/year.
-rate_db.data = [0.0,3.1688087814e-10]
-
-[pylithapp.timedependent.bc.z_neg]
-fixed_dof = [2]
-label = 16
-db.label = Dirichlet BC on -z
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces]
-
-# Set the parameters for the fault interface conditions.
-
-[pylithapp.timedependent.interfaces.fault]
-label = 10
-quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-mat_db.iohandler.filename = mat_elastic.spatialdb
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_unif_8000m.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK domain output of solution over ground surface.
-[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
-writer.filename = ../results/spbm_hex8_unif_8000m-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_unif_8000m-fault.vtk
-time_format = %04.0f
-time_constant = 1.0*year
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_unif_8000m-statevars-elastic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.maxwell.output]
-cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = ../results/spbm_hex8_unif_8000m-statevars-maxwell.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.cfg 2008-09-15 20:23:39 UTC (rev 12883)
@@ -0,0 +1,127 @@
+# -*- Python -*-
+[pylithapp]
+
+# This is not a self-contained simulation configuration file. This
+# file only specifies parameters specific to the boundary and
+# interface conditions. The general parameters are specificed in the
+# pylithapp.cfg file which PyLith reads by default. Information controlling
+# time stepping is contained in the file timedep.cfg.
+#
+# To run the simulation:
+# pylith spbm_hex8_unif_6750m.cfg timedep.cfg
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+# Change the mesh importer to the CUBIT importer.
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+# Set filename of mesh to import.
+filename = ../meshes/spbm_hex8_unif_6750m.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
+
+# Set interfaces to an array of 1 fault: 'fault'.
+interfaces = [fault]
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+
+# We change spatial database for the rate of change from FixedDOFDB
+# (which has zero velocity) to UniformDB (which allows nonzero
+# velocities).
+
+# Set the parameters for the desired boundary conditions.
+#
+# The label corresponds to the nodeset ID in CUBIT.
+
+[pylithapp.timedependent.bc.x_pos]
+fixed_dof = [0, 1]
+label = 11
+db.label = Dirichlet BC on +x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on +x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of -1 cm/year.
+rate_db.data = [0.0,-3.1688087814e-10]
+
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = 12
+db.label = Dirichlet BC on -x
+rate_db = spatialdata.spatialdb.UniformDB
+rate_db.label = Dirichlet rate BC on -x
+rate_db.values = [dof-0,dof-1]
+# This rate (in m/s) gives a value of 1 cm/year.
+rate_db.data = [0.0,3.1688087814e-10]
+
+[pylithapp.timedependent.bc.z_neg]
+fixed_dof = [2]
+label = 16
+db.label = Dirichlet BC on -z
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+
+# Set the parameters for the fault interface conditions.
+
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+mat_db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/spbm_hex8_unif_6750m.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK domain output of solution over ground surface.
+[pylithapp.problem.formulation.output.subdomain]
+label = 17 ; nodeset for subdomain
+writer.filename = ../results/spbm_hex8_unif_6750m-groundsurf.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Give basename for vtk fault rupture output.
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/spbm_hex8_unif_6750m-fault.vtk
+time_format = %04.0f
+time_constant = 1.0*year
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif_6750m-statevars-elastic.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.maxwell.output]
+cell_filter = pylith.meshio.CellFilterAvg
+writer.filename = ../results/spbm_hex8_unif_6750m-statevars-maxwell.vtk
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
Added: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh
===================================================================
--- short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh (rev 0)
+++ short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh 2008-09-15 20:23:39 UTC (rev 12883)
@@ -0,0 +1 @@
+pylith spbm_hex8_unif_6750m.cfg timedep.cfg
Property changes on: short/2.5D/benchmarks/savageprescott/parameters/spbm_hex8_unif_6750m.sh
___________________________________________________________________
Name: svn:executable
+ *
More information about the cig-commits
mailing list