Abstract. Parameterization of turbulent fluxes under stably stratified conditions has always been a challenge. Current surface fluxes calculation schemes either need iterations or suffer low accuracy. In this paper, a non-iterative scheme is proposed to approach the classic iterative computation results using multiple regressions. It can be applied to the full range of roughness status 10 ≤ z/z₀ ≤ 10⁵ and −0.5 ≤ log (z₀/z_0h) ≤ 30 under stable conditions 0 < Ri_B ≤ 2.5. The maximum (average) relative errors for the turbulent transfer coefficients for momentum and sensible heat are 12% (1%) and 9% (1%), respectively.