We investigate β-ensembles with real-analytic potential V in the intermediate temperature regime, 1 >>β>>1/N. We obtain global and mesoscopic CLTs after recentering by a refined equilibrium measure. Our result reveals a transition between random matrix/Poisson regimes in the variance of the fluctuations depending on the mesoscopic scale and the parameter α=βN. Our main technique consists in studying loop equations in this intermediate temperature regime and the proofs involve the construction of an appropriate re-centering which refines the usual equilibrium measure and two optimal local law (one for each regime).