Many NWP and climate models use prognostic turbulence kinetic energy (TKE) equation. TKE dissipation ratio is parameterized via suitable ratio of TKE and so-called integral length scale for turbulence, which is weighted by a certain numerical coefficient. The scientific community does not have unique agreement on formulation of integral turbulence scales. Therefore, in literature we find various forms of those integral scales. However, usually they are derived from autocorrelation functions and Fourier spectra of wind speed components. These are the so-called 1/e and zero crossing scales and the scale at which a normalized Fourier spectra achieve their maximum value. The aim of this study is to evaluate integral turbulence scales and the related coefficient, for the first time for bora flows. For the analysis we used dana collected with „WindMaster“ ultrasonic anemometers at heights 2, 5 and 10 m, close to new Maslenica bridge, with a frequency sampling of 20 Hz. In the period from 09 October 2015 to 09 October 2016 recorded were 48 bora events with cumulative duration of 1206 hours. Such a large dataset gives us an opportunity to evaluate integral turbulence scales and the related numerical coefficient to test the parametrization of the TKE dissipation rate for the atmospheric surface boundary layer in bora flows. So-called 1/e scale proved to be the best estimate for the integral turbulence scale. At the end, based on the analyzed dataset, it has been observed a significant disagreement of the numerical coefficient for bora flows compared with values used in the literature, proposed by e.g. Cuxart et al. (2006). This leaves us space for further studies in future.