-
Notifications
You must be signed in to change notification settings - Fork 96
/
step05_onefault.cfg
88 lines (77 loc) · 2.91 KB
/
step05_onefault.cfg
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
[pylithapp.metadata]
# y
# ^
# |
# --> x
#
# ----------
# | |
# Ux= 0 | | Ux= 0
# | |
# | |
# ----------
# Uy=0
#
# Dirichlet boundary conditions
#
# boundary_yneg: Uy(x,-100*km) = 0
# boundary_xneg: Ux(-126*km,y) = 0
# boundary_xpos: Uy(+74*km,y) = 0
base = [pylithapp.cfg]
description = Reverse prescribed slip with zero displacement Dirichlet boundary condition.
authors = [Charles Williams]
keywords = [reverse slip, prescribed slip]
arguments = [step05_onefault.cfg, solver_onefault.cfg]
version = 1.0.0
pylith_version = [>=3.0, <4.0]
features = [
Static simulation,
pylith.materials.Elasticity,
pylith.materials.IsotropicLinearElasticity,
pylith.faults.FaultCohesiveKin,
pylith.faults.KinSrcStep,
spatialdata.spatialdb.UniformDB
]
# ----------------------------------------------------------------------
# Monitoring and parameter viewing.
# ----------------------------------------------------------------------
[pylithapp]
# You can view this file using the pylith_parameters application or
# online at https://geodynamics.github.io/pylith_parameters/.
dump_parameters.filename = output/step05_onefault-parameters.json
problem.progress_monitor.filename = output/step05_onefault-progress.txt
# Set the name of the problem that will be used to construct the
# output filenames. The default directory for output is 'output'.
problem.defaults.name = step05_onefault
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
# Solving the elasticity equation with a fault means the solution must include
# both displacement and Lagrange multiplier subfields.
[pylithapp.problem]
solution = pylith.problems.SolnDispLagrange
defaults.quadrature_order = 1
[pylithapp.problem.solution.subfields]
displacement.basis_order = 1
lagrange_fault.basis_order = 1
# ----------------------------------------------------------------------
# fault
# ----------------------------------------------------------------------
[pylithapp.problem]
interfaces = [fault]
[pylithapp.problem.interfaces.fault]
# The `label` and `label_value` correspond to the name and tag of the physical group
# for the fault in the Gmsh Python script. The `edge` and `edge_value` correspond to
# the name and tag of the physical group for the fault end in the Gmsh Python script.
label = fault
label_value = 20
edge = fault_end
edge_value = 21
observers.observer.data_fields = [slip]
[pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
# We have uniform slip, so we use a `UniformDB` spatial database.
db_auxiliary_field = spatialdata.spatialdb.UniformDB
db_auxiliary_field.description = Fault rupture auxiliary field spatial database
db_auxiliary_field.values = [initiation_time, final_slip_left_lateral, final_slip_opening]
db_auxiliary_field.data = [0.0*s, -2.0*m, 0.0*m]
# End of file