import spharpy
import numpy as np
coords = spharpy.samplings.equal_area(n_max=0, n_points=500)
data = np.sin(coords.azimuth) * np.cos(coords.elevation)
spharpy.plot.balloon_wireframe(coords, data, cmap_encoding='phase')