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