Step 6: Error 8#

Error Message#

Listing 290 Output when running Step 6.#
 1$ pylith step06_twofaults.cfg
 2
 3-- many lines not shown --
 4
 5 >> /software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/problems/TimeDependent.py:132:run
 6 -- timedependent(info)
 7 -- Solving problem.
 80 TS dt 0.2 time -0.2
 9    0 SNES Function norm 1.933709845425e-02
10      Linear solve converged due to CONVERGED_ATOL iterations 23
11    1 SNES Function norm 9.419694087448e-13
12    Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 1
131 TS dt 0.2 time 0.
14    0 SNES Function norm 9.419697323909e-13
15    Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 0
162 TS dt 0.2 time 0.2
17    0 SNES Function norm 8.112736904304e-03
18      Linear solve converged due to CONVERGED_ATOL iterations 22
19    1 SNES Function norm 1.563260631101e-12
20    Nonlinear solve converged due to CONVERGED_FNORM_ABS iterations 1
213 TS dt 0.2 time 0.4
22 >> /software/unix/py3.12-venv/pylith-debug/lib/python3.12/site-packages/pylith/problems/Problem.py:199:finalize
23 -- timedependent(info)
24 -- 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#

Listing 291 Correct error in step06_twofaults.cfg.#
[pylithapp.problem.interfaces.fault]
...
db_auxiliary_field.query_type = linear