Reynolds stresses associated with the extracted modes in our simulations summed over m = 1 − 39. Panels a–c: ζ_z-symmetric columnar convective modes, the n = 0 equatorial Rossby modes, and the mixed Rossby modes, respectively. Panel d: the total Reynolds stresses (including other modes and small-scale convection) is shown. Upper and lower panel: respectively corresponds to ρ₀⟨v_rv_ϕ⟩ and ρ₀⟨v_θv_ϕ⟩.