In subroutine transit_factor there is a typo with respect to the units of the variable au. It is given in meters but the comment reads km.