We present a model for the scaling of mixing in weakly rotating stratified flows characterized by their Rossby, Froude and Reynolds numbers Ro, F r, Re. It is based on quasiequipartition between kinetic and potential modes, sub-dominant vertical velocity w, and lessening of the energy transfer to small scales as measured by a dissipation efficiency β = V / D , with V the kinetic energy dissipation and D = u 3 rms /L int its dimensional expression, w, u rms the vertical and rms velocities, and L int the integral scale. We determine the domains of validity of such laws for a large numerical study of the unforced Boussinesq equations mostly on grids of 1024 3 points, with Ro/F r 2.5, and with 1600 Re ≈ 5.4 × 10 4 ; the Prandtl number is one, initial conditions are either isotropic and at large scale for the velocity, and zero for the temperature θ, or in geostrophic balance. Three regimes in Froude number, as for stratified flows, are observed: dominant waves, eddy-wave interactions and strong turbulence. A wave-turbulence balance for the transfer time τ tr = N τ 2 N L , with τ N L = L int /u rms the turn-over time and N the Brunt-Väisälä frequency, leads to β growing linearly with F r in the intermediate regime, with a saturation at β ≈ 0.3 or more, depending on initial conditions for larger Froude numbers. The Ellison scale is also found to scale linearly with F r. The flux Richardson number, with B f = N wθ the buoyancy flux, transitions for roughly the same parameter values as for β. These regimes for the present study are delimited bythe mixing efficiency, putting together the three relationships of the model allows for the prediction of the scaling Γ f ∼ F r −2 ∼ R −1 B in the low and intermediate regimes for high Re, whereas for higher Froude numbers, Γ f ∼ R −1/2 B, a scaling already found in observations: as turbulence strengthens, β ∼ 1, w ≈ u rms , and smaller buoyancy fluxes altogether correspond to a decoupling of velocity and temperature fluctuations, the latter becoming passive.