sp3d = Ga('e_r e_th e_ph', g=[1, r**2, r**2 * sin(th)**2],