Add the qualifier 'N_disk_crossings' to the function 'orbit_calculator'

With the new qualifier 'N_disk_crossings', it is possible to stopp
the computation of orbits at an arbitrary disk crossing event rather
than only at the first one.