The NOVAS C geo_posvel() function is imprecise for geodetic observer locations on Earth, since it does not account for polar offsets x_p, y_p. As such the geocentric GCRS observer position and velocity vectors are accurate at the arc-seconds level only.