Comparison of the r and z velocity components in the initial set-up (thin solid line) with the quasi-stationary solutions in the numerical simulations in the HD (dot-dashed line) and the MHD (long-dashed line) cases, with Ω = 0.2Ω_ms. The meaning of the lines is the same as in Fig. 3. Left panels: values averaged over 10 stellar rotations, to smooth out the oscillations. Right panels: values at r = 15R_⋆. In the z-direction, an additional match is shown with squares in the cases with larger magnetic field. It has the same (analytically derived) functional dependence as shown with triangles for the smaller magnetic field, only the proportionality coefficient is different.