# Define how many time steps are used. The cadence right now has a probelmatically large effect on the outcome. This needs
#  to be smoothed out.
period = 6 * 3600
orbits = 50
tmax = period * orbits
cadence = 30
times = np.arange(0, tmax, cadence)