Convert to use astropy erfa, but not yet arrays
Cleaner numpy version