As my understanding progressed, the code became a mess. Thought I’d need Cartesian xyz, then spherical, then orbits. Eventually, I use great circles, with the particles’ velocities updated using a Coriolis acceleration approximation.
Click on the planet (which rotates very fast) and watch the particles move under (hopefully a fairly accurate version of) the Coriolis effect.
Thanks!