spherical_harmonics = spharpy.SphericalHarmonics(
    n_max=n_max, coordinates=sampling,
    basis_type='real', normalization='SN3D', condon_shortley=False
)