Two-dimensional visualizations at t = 1000 for simulations with α_P = 0, 1, and 2, from left to right. The top panels display the full computational domain, while the bottom panels show a zoomed-in view of the inner region. Dashed lines, with increasing opacity, represent the slow magnetosonic (SM), Alfvén, and fast magnetosonic (FM) critical surfaces, respectively. White lines indicate magnetic-field lines anchored at R = 2, 10, 20, 30, and up to 90. The background color map shows the density in arbitrary units.