[cig-commits] r18955 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters
willic3 at geodynamics.org
willic3 at geodynamics.org
Tue Sep 20 21:41:52 PDT 2011
Author: willic3
Date: 2011-09-20 21:41:52 -0700 (Tue, 20 Sep 2011)
New Revision: 18955
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/asm.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
Log:
More cleaning up of S&P example.
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/asm.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/asm.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/asm.cfg 2011-09-21 04:41:52 UTC (rev 18955)
@@ -0,0 +1,10 @@
+# -*- Python -*-
+[pylithapp]
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+# Set the solver options.
+
+[pylithapp.petsc]
+pc_type = asm
+sub_pc_factor_shift_type = nonzero
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg 2011-09-21 04:13:05 UTC (rev 18954)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/hex8_lock20_tau0-1.cfg 2011-09-21 04:41:52 UTC (rev 18955)
@@ -2,53 +2,40 @@
[pylithapp]
# This is not a self-contained simulation configuration file. This
-# file only specifies parameters specific to the mesh used.
+# file only specifies material information and output file names for a
+# particular mesh and set of parameters.
# The general parameters are specificed in the pylithapp.cfg file which
-# PyLith reads by default. Information controlling time-dependent fault
+# PyLith reads by default. Mesh information is contained in hex8.cfg,
+# fault locking information is contained in fault_lockxx.cfg, and solution
+# parameters are contained in asm.cfg and fieldsplit.cfg.
# behavior is contained in the file timedep.cfg.
#
# To run the simulation:
-# pylith spbm_hex8_graded_20km.cfg timedep.cfg
+# pylith hex8.cfg fault_lock20.cfg fieldsplit.cfg hex8_lock20_tau0-1.cfg
# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator.reader]
-# Set filename of mesh to import.
-filename = ../meshes/spbm_hex8_graded_20km.exo
-
-# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials]
-elastic.quadrature.cell = pylith.feassemble.FIATLagrange
-elastic.quadrature.cell.dimension = 3
-maxwell.quadrature.cell = pylith.feassemble.FIATLagrange
-maxwell.quadrature.cell.dimension = 3
+[pylithapp.timedependent.materials.maxwell]
+db_properties.iohandler.filename = mat_maxwell_tau0-1.spatialdb
+db_properties.label = Viscoelastic properties yielding tau0 = 1
# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.fault]
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Give basename for VTK files.
[pylithapp.problem.formulation.output.domain.writer]
-filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km.vtk
+filename = ../results/hex8_lock20_tau0-1/hex8_lock20_tau0-1.h5
[pylithapp.problem.formulation.output.subdomain]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-groundsurf.vtk
+writer.filename = ../results/hex8_lock20_tau0-1/hex8_lock20_tau0-1-groundsurf.h5
[pylithapp.timedependent.interfaces.fault.output.writer]
-filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-fault.vtk
+filename = ../results/hex8_lock20_tau0-1/hex8_lock20_tau0-1-fault.h5
[pylithapp.timedependent.materials.elastic.output]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-elastic.vtk
+writer.filename = ../results/hex8_lock20_tau0-1/hex8_lock20_tau0-1-elastic.h5
[pylithapp.timedependent.materials.maxwell.output]
-writer.filename = ../results/spbm_hex8_graded_20km/spbm_hex8_graded_20km-maxwell.vtk
+writer.filename = ../results/hex8_lock20_tau0-1/hex8_lock20_tau0-1-maxwell.h5
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg 2011-09-21 04:13:05 UTC (rev 18954)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/parameters/pylithapp.cfg 2011-09-21 04:41:52 UTC (rev 18955)
@@ -38,7 +38,6 @@
[pylithapp.mesh_generator]
# Change the mesh reader to the CUBIT reader.
reader = pylith.meshio.MeshIOCubit
-reader.use_nodeset_names = False
reorder_mesh = true
[pylithapp.mesh_generator.reader]
@@ -117,11 +116,14 @@
# Apply a rate of -1 cm/year.
db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
+# Fix z-displacements on bottom boundary.
[pylithapp.timedependent.bc.z_neg]
bc_dof = [2]
label = face_zneg
db_initial.label = Dirichlet BC on -z
+# Fix x-displacements on y-boundaries to prevent movement perpendicular to
+# the fault.
[pylithapp.timedependent.bc.y_pos]
bc_dof = [0]
label = face_ypos_nofault
@@ -143,30 +145,45 @@
# output
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.domain.writer]
-time_format = %04.0f
-time_constant = 1.0*year
+[pylithapp.problem.formulation.output.domain]
+output_freq = time_step
+time_step = 50*year
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
# Give basename for VTK domain output of solution over ground surface.
[pylithapp.problem.formulation.output.subdomain]
label = face_zpos ; nodeset for subdomain
vertex_data_fields = [displacement,velocity]
+output_freq = time_step
+time_step = 10*year
+writer = pylith.meshio.DataWriterHDF5SubMesh
writer.time_format = %04.0f
writer.time_constant = 1.0*year
# Give basename for vtk fault rupture output.
-[pylithapp.timedependent.interfaces.fault.output.writer]
-time_format = %04.0f
-time_constant = 1.0*year
+[pylithapp.timedependent.interfaces.fault.output]
+output_freq = time_step
+time_step = 10*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
+output_freq = time_step
+time_step = 50*year
cell_filter = pylith.meshio.CellFilterAvgMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
writer.time_format = %04.0f
writer.time_constant = 1.0*year
[pylithapp.timedependent.materials.maxwell.output]
+output_freq = time_step
+time_step = 50*year
cell_filter = pylith.meshio.CellFilterAvgMesh
+writer = pylith.meshio.DataWriterHDF5Mesh
writer.time_format = %04.0f
writer.time_constant = 1.0*year
@@ -176,8 +193,6 @@
# Set the solver options.
[pylithapp.petsc]
-pc_type = asm
-sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 1000
More information about the CIG-COMMITS
mailing list