A fast interpolation algorithm using the partition of unity method on the sphere