Left: 2D visualization at t = 1000 for the simulation performed over the interval [0, π]. 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, up to 130. The background colormap shows the density in arbitrary units. Upper right: Latitudinal profiles of the magnetic field components B_R, B_θ, B_ϕ (normalized to B_z at the disk midplane), and components of the velocity v_R, v_θ, v_ϕ (normalized to the Keplerian velocity v_K at the disk midplane), measured at R = 2.0 and t = 1000. The central gray-shaded region corresponds to the disk, the white region to the disk wind, and the outer gray-shaded regions to the spine. Lower right: Profiles of the four MHD invariants κ, ω_*, λ, e along the magnetic field line anchored at R = 2.0, measured at t = 1000. The shaded area indicates the resistive MHD region, while the vertical dashed lines mark the positions of the SM, A, and FM critical surfaces.