Cross-sectional cuts of density (top left), temperature (top right), current density squared (bottom left), and vorticity squared (bottom right) at the end of the ideal MHD simulation, $\overline{t}=820$. Temperature is normalized with respect to the initial value at r = 0. Logarithmic scale is used for each variable except for the density to optimize visualization. The cross-sectional cuts are done at the tube center, z = 0, except for current density squared which is done at z = R. The complete temporal evolution is available as an online movie.