Comparison of the saturation regimes in the simulation with the local HSE BC of Zingale et al. (2002) (opaque lines) and in the 2dDf1e–2 run with the BCs described in Sect. 3.3 (transparent lines). Top: Volume-averaged time evolution of the available potential energy density ⟨E_A⟩_𝒱, and of the different contributions to the kinetic ⟨E_K⟩_𝒱 volume energy density. Bottom: Volume-averaged time evolution of the APE injection rate ⟨ε_i⟩_𝒱, APE thermal dissipation rate −⟨ε_κ⟩_𝒱, and the (opposite amount of) reversible buoyancy work −⟨𝒜⟩_𝒱 during the same run. Dashed lines are indicative of negative quantities whose sign has been switched for the purpose of visualisation.