Address invalid value encountered in generating urchin mesh

Specifically this MR will address:

• Raise exception on floating point error for mesh generation.
• Use more robust spherical harmonic function from pyfmmlib when available, instead of the one in scipy.