Differential rotation between the inner equator and the outer sphere normalised to the outer value of the rotation rate Ω₀ in the Eddington–Sweet regime, plotted as a function of P_r(N₀/Ω₀)²R_o/E. The Rossby number R_o is varied to encompass the linear and non-linear regimes. The Ekman number is fixed to E = 10⁻⁶. The green curve corresponds to simulations with P_r(N₀/Ω₀)² = 10⁻¹ (runs 1.4, 3.1, 4.1, 5.1, and 6.1 in Table 4) and the blue curve to simulations with P_r(N₀/Ω₀)² = 10⁻² (runs 2.2, 3.2, 4.2, and 5.2 in the same table).