Non radial oblateness force and more precise planet data

This is a `NOTE(phl)`, here it looks like the code is actively not trying :-p