Flip sign of _y in celestial set_aux.  It was a remnant of old code that is confusing now.
Add Rachel's elliptic function