Title: A 3.3 Million Year Simulation of the Geodynamo using MAG Summary: These images show output from a MAG calculation simulating the long-term behavior of the geodynamo driven by compositional buoyancy derived from solidification of the Earth's inner core. Model Parameters: Rayleigh number=8x10^5; Ekman number=6.5x10^-3; Prandtl number=1; magnetic Prandtl number=20; heat source=-1 (chemical convection); truncation l(max)=32. Fig magSSC1: Time series of rms fluid velocity (magnetic Reynolds number units), rms internal magnetic field (Elsasser number units), rms CMB dipole (Elsasser number units) and dipole tilt. Time axis is dipole decay times, approximately 20 kyr in the core. Shaded region marks Reverse-to-Normal polarity change. Fig magSSC2: Detailed time series during the polarity change shown in Fig magSSC1. Times a and b correspond to the images shown in Figs magSSC3 and magSSC4, respectively. Fig magSSC3: Dynamo structure prior to the polarity reversal (time a in Fig magSSC2). Left images are on the core-mantle boundary; right images are in the outer core volume. Continental outlines are for reference. Fig magSSC4: Dynamo structure after the polarity reversal (time b in Fig magSSC2). Left images are on the core-mantle boundary; right images are in the outer core volume. Continental outlines are for reference. Fig magSSC5: Geomagnetic dipole paths and dipole statistics from the simulation.