Results of the mean flows from the hydrodynamic calculation. (a) Differential rotation ⟨Ω⟩=Ω₀ + ⟨v_ϕ⟩/(r sin θ). The black dotted curve shows the location of the base of the convection zone at r = 0.71 R_⊙, below which the stratification is subadiabatic. (b) Meridional circulation ⟨v_θ⟩. The black solid and dashed curves show contours of the streamfunction Ψ defined by ρ₀v_m = ∇ × (Ψe_ϕ) where v_m = ⟨v_r⟩e_r + ⟨v_θ⟩e_θ. The meridional circulation is counter-clockwise (clockwise) in the northern (southern) hemisphere, i.e., the flow is poleward (equatorward) at the surface (base).