A Constraint-Preserving Finite Difference Method for the Damped Wave Map Equation to the Sphere