scipy.special.sph_harm is deprecated

[hoyer-a]

sph_harm will be removed in scipy 1.17.
The new sh-function is sph_harm_y

[mberz]

Thanks for catching this. It seems the primary change is the order of n and m and the change of coordinate variables.
Additionally, the function supports returning the derivative, which may be interesting to use when calculating the gradient with respect to phi

I guess in this case it could be beneficial to implement a switch for different scipy versions and release a patch before moving to spharpy 1.0.0

[mberz]

There are now also functions which return the spherical harmonics for multiple degrees and orders, same for the Legendre functions. This could be quite useful to avoid loops when building the basis matrix.

[xefonon]

This is great, I don't understand from the documentation if the derivatives are with respect to both theta and phi.