xyz_2d＝array([sin(theta_2d)*sin(phi_2d),