Thermally driven slope flows are ubiquitous in nature and play a major role in regulating local microclimates in valleys, glaciers and ice-sheets. They control in large part surface momentum, heat and moisture fluxes, and their effects must therefore be accounted for in weather prediction and climate models. Due to the interplay between shear and buoyancy in generating and destroying turbulence, the thermal and hydrodynamic boundary layers that characterizes slope flows are very shallow, with turbulent motions of reduced size, when compared to those populating neutral and convective boundary layers, posing serious difficulties in terms of computational resolution and modeling requirements. Not surprisingly, therefore, considerable effort has been devoted in designing parameterizations, for use in larger scale models. In this direction goes the current work, where we use direct numerical simulations (DNS) to determine the vertical profiles of momentum and buoyancy eddy diffusivities, for pure slope flows over smooth surfaces. Mean flow profiles from one-dimensional models are then compared against DNS results at various sloping angles and results will be discussed.