Fix a bug in the function 'orbit_calculator'

There was a bug in the function that prevented it from passing the
model-dependent Sun-Galactic-center distance to the sub-routine
'cel2gal'. Consequently, the default value (8.4 kpc) was always used.