Figure 2: Comparison of the orbital angular velocity versus time for different values of the Courant parameter . As can be seen, the simulations with = 0.25–0.5 result in stable runs that converge to the same value, implying that a smaller or equivalently a smaller is not necessary and would only use extra CPU time. For comparison, we plot a simulation with to show that the stability is lost for .