It is better to use atan2 here, otherwise incorrect values will be produced. (The same for other methods). https://github.com/qsantos/spyce/blob/2026837180df3fecebb3219df4756cf93b4de103/spyce/coordinates.py#L41 Also all moons in solar system missing NorthPole values.