I wanted to treat it as a surface of revolution, and then simply use
Mod[x[t],2Pi] to have each geodesic loop back onto this single portion of the surface. However, as you can see, there are some strange artifacts on graph when I use
Mod. If I treat the surface like normal, then there are no artifacts, so I know it's not part of the calculation.
How should I go about getting these to disappear? Ideally, I'd like to plot geodesics for large times, and watch the path. However every time the path leaves one end, it produces one of these artifacts, which clutters up the picture considerably.