Due to a mistake when migrating from NOVAS C 3.1, the polar wobble corrections set via cel_pole() were not applied to the true obliquity and equation of equinoxes in e_tilt(). This only affects the old, now disfavored, way of incorporating Earth orientation parameters. The more modern way is via wobble(), which is used by cel2ter() / ter2cel() as well as by cirs_to_itrs() / itrs_to_cirs() and tod_to_itrs() and itr_to_tod().

To give stronger preference to the favored (more modern) method, the cel_pole() function and the global PSI_COR / EPS_COR variables will be deprecated going forward.