Differential rotation in the viscous regime scaled by Re_c as a function of radius at the latitude θ = π/8. The analytical solution Eq. (48) (black dashed lines) is compared with the numerical results. Low Reynolds number runs (1.1–4.1 in Table 3, with Re_c = 10⁻² − 10⁻¹) are plotted with violet dotted lines and are all very close to the low Reynolds number version of the analytical solution, Eq. (49). Simulations at Re_c = 1 (runs 4.2 and 5.1 in Table 3) and Re_c = 10 (run 4.3 Table 3) are also in agreement with the analytical solution.