# -*- Python -*-

# If this file is given, the settings in this file will override all
# others except for parameters specified directly on the command line.

[pylithapp]

# ----------------------------------------------------------------------
# problem
# Specify the problem settings.
# This is a time-dependent problem so we use that facility.
# ----------------------------------------------------------------------
[pylithapp.timedependent]

# We want an implicit formulation.
formulation = pylith.problems.Implicit

bc = [East,West,North,South,z_neg]

# Set gravity field (default is None)
#gravity_field = spatialdata.spatialdb.GravityField

[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]
#output = [domain]

# Set subdomain component to OutputSolnSubset (subset of domain).
output.subdomain = pylith.meshio.OutputSolnSubset

# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
# -x face
[pylithapp.timedependent.bc.West]
bc_dof = [0,1]
label = West
db_initial.label = Dirichlet BC on West

# -y face
[pylithapp.timedependent.bc.South]
bc_dof = [0,1]
label = South
db_initial.label = Dirichlet BC on South

# +x face
[pylithapp.timedependent.bc.East]
bc_dof = [0,1]
label = East
db_initial.label = Dirichlet BC on East

# +y face
[pylithapp.timedependent.bc.North]
bc_dof = [0,1]
label = North
db_initial.label = Dirichlet BC on North

# -z face
[pylithapp.timedependent.bc.z_neg]
bc_dof = [0,1,2]
label = z_neg
db_initial.label = Dirichlet BC on -z

# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
# Give basename for VTK output of solution over domain.
[pylithapp.problem.formulation.output.domain]
output_freq = time_step
time_step = 30.0*day
# We are using HDF5 output so we must change the default writer.
#writer = pylith.meshio.DataWriterHDF5Mesh
writer.filename = Test1.vtk
writer.time_format = %04.1f ; Time stamp will be xxx.x with the "." removed.
writer.time_constant = 1*day ; Normalize time stamp by years

# Give basename for VTK domain output of solution over ground surface.
[pylithapp.problem.formulation.output.subdomain]
label = z_pos ; Name of nodeset for subdomain
output_freq = time_step
time_step = 30.0*day
#writer = pylith.meshio.DataWriterHDF5SubMesh
writer.filename = Test1-groundsurf.vtk
writer.time_format = %04.1f
writer.time_constant = 1*day

# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.EL_BIFASE.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
output_freq = time_step
time_step = 30.0*day
cell_info_fields = []
cell_data_fields = [] 
#writer = pylith.meshio.DataWriterHDF5Mesh
#writer.filename = ../../../data/bjha/Test1-statevars.vtk
#writer.time_format = %04.1f ; Time stamp will be xxx.x with the "." removed.
#writer.time_constant = 10*day ; Normalize time stamp by years

[pylithapp.timedependent.materials.EL_MONO.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
output_freq = time_step
time_step = 30.0*day
cell_info_fields = [] 
cell_data_fields = [] 
#writer = pylith.meshio.DataWriterHDF5Mesh
#writer.filename = ../../../data/bjha/Test1-statevars.vtk
#writer.time_format = %04.1f ; Time stamp will be xxx.x with the "." removed.
#writer.time_constant = 10*day ; Normalize time stamp by years