I have finally uploaded a small Python program for plotting certain aspects of the Weierstrass Elliptic P function. It can be found in my github repository. These plots visualize the behaviour of the Weierstrass elliptic P function on the standard square torus. This is something I have studied in the context of my doctoral thesis (see p. 39, ff.).

The first video visualizes the global behaviour of the Weierstrass P function. At each point in time, what you see is the image of a circle of the standards square torus of the form Cy = (x, y)|0 < =x < =1 under the Weierstrass P function on the Riemann sphere.

The second video visualizes the behaviour of the Weierstrass P function along the circle C0.

And finally, the third video visualizes the behaviour of the Weierstrass P function along the circle C0.5.