Stationary differential rotation δΩ(r, θ)/Ω₀ (first and third panels) and norm of the total meridional velocity field (second and fourth panels) with its associated streamlines in black, in the viscous anelastic regime. The parameters of the simulation shown in the two left panels are E = 10⁻⁵, P_r(N₀/Ω₀)² = 10⁴, Re_c = 10⁻², and ρ_i/ρ₀ = 20.9 (run 5.2 in Table 7), and in the two right panelsE = 10⁻⁴, P_r(N₀/Ω₀)² = 10⁴, Re_c = 10, and ρ_i/ρ₀ = 20.9 (run 4.3 in Table 7). In the first and third panels, the coloured contours are normalised to the top value Ω₀.