(A) Steady angular velocity component (see Eq. (8)) as extracted from a hydrodynamical EULAG simulation, plotted as isocontours of rotational frequency. (B) and (C) Total angular velocity at epoch of Grand Maximum and Minimum, respectively, as produced by our non-kinematic α²Ω mean-field dynamo model (specifically, simulation 4 in Table 1). The interface between the convection zone and underlying radiative zone is indicated by the white dashed line.