Interesting, does the mesh have to be of spherical topology or does it work on other meshes? What is your method for constraining?

No, it works with any given mesh from processing or imported.. I constrained the agents on the mesh by a simple routine of "correcting" at each iteration of the algorithm the position of the agent according to its closest mesh quad. There must be a better, in terms of computational power management, way to do this though...