timestep_DufortFrankel

Module

Description

$Id$

DuFort-Frankel time stepping routine. As is well known, this proceedure has disadvantages, but it is well suited to deal with the time step constrain near the center in spherical polar coordinates.

Quick access

Variables:

a21, a31, a32, advance_substeps, after_substep_sum_time, b1, b2, b3, b4, beta_hat, c21, c31, c32, c41, c42, c43, calc_pencils_sub_cycle, dt_alpha_ts, dt_beta_hat, dt_decrease, dt_increase, dt_major, e1, e2, e3, e4, errcon, errdf, errmax, errmaxs, f1, farr, farraymin, gam, grow, itter, maxtry, ode, ode_timestep_first, ode_timestep_second, pde_energy_only, rk3, rkck, rkck3, rkl_coeff, safety, shrnk, split_update, stiff, substeps, swap, update_after_substep

Routines:

initialize_timestep(), time_step()

Needed modules

Variables

Subroutines and functions

subroutine  timestep/time_step(f, df, p)

Runge Kutta advance, accurate to order itorder. At the moment, itorder can be 1, 2, or 3.

2-apr-01/axel: coded

14-sep-01/axel: moved itorder to cdata

Parameters:
Use :

mpicomm, cdata, equ, particles_main, borderprofiles (border_quenching()), interstellar (calc_snr_damp_int()), shear (advance_shear()), special (special_after_timestep())

Call to:

pde(), border_quenching(), special_after_timestep()