Realistic three-dimensional magneto-thermal simulations of magnetars with strong, large-scale toroidal magnetic fields accurately describe the observed light curves of 10 out of 19 magnetars in quiescence and allow their rotational orientation to be further constrained.
