Top: velocity profiles as functions of the latitudinal coordinate, θ, normalised to the local sound speed. Bottom: mean magnetic field profiles normalised to the vertical magnetic field in the disk mid-plane, ⟨B_z⟩(r, θ = π/2) as functions of the latitudinal coordinate. The grey zones correspond to the turbulent regions in Fig. 3. All profiles are radially averaged between r₁ = 8 and r₂ = 14 after being normalised.