Hello,

I am trying to run a simple post-seismic relaxation model that I created with power-law viscoelastic lower crust/mantle. Below are the first few lines in my .cfg file (step02_postseis_powerlaw.cfg) where I define the time step:

[pylithapp]

dump_parameters.filename = output/step02-postseis-powerlaw.json

problem.progress_monitor.filename = output/step02-postseis-powerlaw-progress.txt

# ----------------------------------------------------------------------

# problem

# ----------------------------------------------------------------------

[pylithapp.problem]

# For this problem we must switch to a nonlinear solver.

implicit.solver = pylith.problems.SolverNonlinear

[pylithapp.problem.formulation]

time_step = pylith.problems.TimeStepAdapt ; Change the time step algorithm

[pylithapp.problem.formulation.time_step]

total_time = 20.0*year
max_dt = 1.0*year

adapt_skip = 10 ; Default value

stability_factor = 2.0 ; Default value

However, when executing the scripts (pylith step02_mat_powerlaw.cfg step02_postseis_powerlaw.cfg solver_fieldsplit.cfg), it appears that the time step is way too small compared to the 20-year total time (in the order of second):

/Users/Alpine/PyLith/pylith-2.2.2-darwin-10.13.6/lib/python2.7/site-packages/pylith/problems/TimeDependent.py:195:run

â€“ timedependent(info)

â€“ Preparing to advance solution from time t=0s to t=2.75653s.

I wonder why pylith is assigning such a small time step and how can I increase it (to maybe 0.5 year)?

I have tried using a uniform time step but it ended up in errors:

RuntimeError: Current nondimensionalized time step of 5.0000e-03 exceeds the nondimensionalized stable time step of 1.7470e-09.

application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0

Thanks in advance for any help,

-Guo Cheng