* Rotated SH functions to center zonal harmonics around y
* Implemented projection using polynomials