The bathymetry time step is used now for the earthquake phase time steps
Changed divergence correction to be more diffusive - more stable but allows unphysical saturation in certain situations