Step 6: Error 8#
Error Message#
1$ pylith step06_twofaults.cfg
2
3 -- Solving problem.
40 TS dt 0.01 time 0.
5 0 SNES Function norm 1.975141284264e-02
6 Linear solve converged due to CONVERGED_ATOL iterations 408
7 1 SNES Function norm 7.771687291628e-13
8 Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 1
91 TS dt 0.01 time 0.01
10 >> /software/baagaard/py38-venv/pylith-debug/lib/python3.8/site-packages/pylith/problems/Problem.py:201:finalize
11 -- timedependent(info)
12 -- Finalizing problem.
Troubleshooting Strategy#
The simulation ran without errors.
In visualizing the output we notice the slip distribution contains a sharp transition from 0 m to 2.0 m; we intended to prescribe slip that is uniform above y=-20 km and tapers linearly to 0 at y=-30 km.
We load the JSON parameter file into the PyLith Parameter Viewer and find that we are using the default query_type
of nearest
for the earthquake rupture parameters.
For our intended piecewise linear variation in slip, we need to use linear
for the query_type
.
Resolution#
[pylithapp.problem.interfaces.fault]
...
db_auxiliary_field.query_type = linear