alpha0 = -k0 * np.sin(theta0) * np.cos(phi0)