The old behaviour: averaging over a maximum number of instantaneous
eddies, becomes an option.

We do not use bisect because `max_delta` should be small, around 6,
and most often there should be only one time step between
instantaneous eddies, so the first guess at distance `max_delta`
should be close to the result. Moreover, we would need to define a
date function if we used bisect.