[pylithapp]

# This is not a self-contained simulation configuration file. This
# file only specifies the general parameters common to the simulations
# in this directory.

# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
# Turn on some journals to show progress.
[pylithapp.journal.info]
timedependent = 1
implicit = 1
petsc = 1
solverlinear = 1
meshiocubit = 1
implicitelasticity = 1
faultcohesivekin = 1
fiatlagrange = 1
pylithapp = 1
materials = 1

# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
#debug = 1   ; uncomment to get very verbose mesh information

# Change the default mesh reader to the CUBIT reader.
reader = pylith.meshio.MeshIOCubit

# Optimize ordering of mesh cells and vertices using reverse
# Cuthill-KcKee algorithm.
#reorder_mesh = True

[pylithapp.mesh_generator.reader]
# Set filename of mesh to import.
filename = mesh/box_tet4_1000m.exo

# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent]

# Set materials to an array of 2 materials:
#   'upper_crust' and 'lower_crust'
materials = [upper_crust,lower_crust]
# Uncomment the following for step1.cfg
materials.lower_crust = pylith.materials.MaxwellIsotropic3D

[pylithapp.timedependent.materials.upper_crust]
label = Upper crust material
id = 1
db_properties.label = Properties for upper crust
#db_properties.iohandler.filename = spatialdb/mat_elastic.spatialdb
db_properties.iohandler.filename = spatialdb/mat_elastic_vector.spatialdb
# For hex cells
#quadrature.cell = pylith.feassemble.FIATLagrange
# For tet cells
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.dimension = 3

[pylithapp.timedependent.materials.lower_crust]
label = Lower crust material
id = 2
db_properties.label = Properties for lower crust
#db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
# Uncomment the following for step1.cfg
db_properties.iohandler.filename = spatialdb/mat_maxwell_vector.spatialdb
#db_properties.iohandler.filename = spatialdb/mat_elastic_vector.spatialdb
# For hex cells
#quadrature.cell = pylith.feassemble.FIATLagrange
# For tet cells
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.dimension = 3

# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
# Set the solver options.
[pylithapp.petsc]
#malloc = True
#malloc_debug = True
#malloc_dump = True

# Preconditioner settings.
pc_type = asm
sub_pc_factor_shift_type = nonzero

# Convergence parameters.
#ksp_rtol = 1.0e-10
#ksp_atol = 1.0e-20
#ksp_max_it = 300
ksp_rtol = 1.0e-8
ksp_atol = 1.0e-12
ksp_max_it = 300
ksp_gmres_restart = 50

# Linear solver monitoring options.
ksp_monitor = true
ksp_view = true
ksp_converged_reason = true

# Nonlinear solver monitoring options.
#snes_rtol = 1.0e-8
#snes_atol = 1.0e-12
#snes_max_it = 100
#snes_monitor = true
#linesearch_monitor = true
#snes_view = true
#snes_converged_reason = true

# PETSc summary -- useful for performance information.
log_summary = true

# Uncomment to launch gdb when starting PyLith.
# start_in_debugger = true

# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Names of output files are set in stepX.cfg. We consolidate all of the
# output settings that are common to all of the simulations here.

[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

# Domain
[pylithapp.problem.formulation.output.domain]
#output_freq = time_step
#time_step = 5.0*year
writer = pylith.meshio.DataWriterHDF5Mesh

# Ground surface
[pylithapp.problem.formulation.output.subdomain]
#output_freq = time_step
#time_step = 5.0*year
label = face_zpos ; Name of CUBIT nodeset for ground surface.
writer = pylith.meshio.DataWriterHDF5SubMesh

# Materials
[pylithapp.timedependent.materials.upper_crust.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
#output_freq = time_step
#time_step = 9.99999*year
writer = pylith.meshio.DataWriterHDF5Mesh

[pylithapp.timedependent.materials.lower_crust.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
#output_freq = time_step
#time_step = 9.99999*year
writer = pylith.meshio.DataWriterHDF5Mesh