Turbulence parameterization over complex terrain is a difficult task due to lack of observational data and thorough theoretical understanding. Based on limited observational results and theoretical analysis, we have developed a turbulence parameterization scheme for shear dominated atmospheric conditions. In this scheme, the mean wind from a diagnostic wind model, the three-dimensional wind field (3DWF) model, is used to compute the atmospheric shear deformation and the observed potential temperature profile is used to compute the Brunt-Vaisala frequency. A turbulent diffusivity parameter is then related to model grid length and local shear and buoyancy conditions. The scheme also considers the ground surface morphology, such as buildings and forest canopy, that enhances the turbulence production.