added y_harmonics, a method to calculate spherical harmonics for given
angles theta and phi, and added the ability to calculate coefficients of
decomposition of the isosurface shape onto a sphericla harmonic basis.
updated run_cif_to_surface to calculate the shape coefficients