[cig-commits] r12411 - short/3D/PyLith/trunk/examples/3d/hex8

[brad at geodynamics.org]

Author: brad
Date: 2008-07-14 10:33:35 -0700 (Mon, 14 Jul 2008)
New Revision: 12411

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
Log:
Switched to using adaptive time stepping in 3d/hex8 savageprescott example.

Modified: short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg	2008-07-14 17:26:18 UTC (rev 12410)
+++ short/3D/PyLith/trunk/examples/3d/hex8/savageprescott.cfg	2008-07-14 17:33:35 UTC (rev 12411)
@@ -13,6 +13,8 @@
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent]
+# Change initial guess for solution from vector of zeros to previous
+# solution.
 formulation.solver.initial_guess_zero = False
 
 # Set bc to an array of 3 boundary conditions: 'x_pos','x_neg', and 'z_neg'.
@@ -21,11 +23,16 @@
 # Set interfaces to an array of 1 fault: 'fault'.
 interfaces = [fault]
 
+# Change time stepping algorithm from uniform time step, to adaptive
+# time stepping.
+implicit.time_step = pylith.problems.TimeStepAdapt
+
 [pylithapp.timedependent.implicit.time_step]
 # Define the total time for the simulation and the default time step size.
 # Overwrite values in pylithapp.cfg
 total_time = 700.0*year ; total time of simulation
-dt = 10.0*year ; time step
+max_dt = 50.0*year ; time step
+stability_factor = 1.0 ; use time step equal to stable value from materials
 
 [pylithapp.timedependent.implicit]
 # Set the output to an array of 2 output managers.