We develop a reduced model for hard-magnetic, thin, linear-elastic shells that can be actuated through an external magnetic field, with geometrically exact strain measures. Assuming a reduced kinematics based on the Kirchhoff-Love assumption, we derive a reduced twodimensional magneto-elastic energy that can be minimized through numerical analysis. In parallel, we simplify the reduced energy by expanding it up to the second order in the displacement field and provide a physical interpretation. Our theoretical analysis allows us to identify and interpret the two primary mechanisms dictating the magneto-elastic response: a combination of equivalent magnetic pressure and forces at the first order, and distributed magnetic torques at the second order. We contrast our reduced framework against a threedimensional nonlinear model by investigating three test cases involving the indentation and the pressure buckling of shells under magnetic loading. We find excellent agreement between the two approaches, thereby verifying our reduced model for shells undergoing nonlinear and non-axisymmetric deformations. We believe that our model for magneto-elastic shells will serve as a valuable tool for the rational design of magnetic structures, enriching the set of reduced magnetic models.