Finishing Relaxation

It is possible to start a dynamical run after relaxation finishes without starting the code. The parameter ``TRELOFF'' gives the time at which relaxation finishes, and if this is less than the time at which the run itself finishes, ``TF'', the code will switch to a dynamical mode, turning on the effects of AV and radiation reaction if the user indicated they should be used.

For a single star configuration the process is extremely simple. All velocities are set to zero and the run is relaunched. For the corotating binary case, All velocities are set such that

$$\vec{v}=\vec{\Omega}\times\vec{r},$$ (14)

where $\vec{\Omega}$ is taken to point in the z-direction. Note that this approach eliminates all spurious velocities which had yet to damp out.