The returned rate in novas_helio_dist() is actually the magnitude of the SSB based motion, not the radial rate of recession (as it should have been).