The Leipzig Wind Profile is generally known as a typical neutral planetary boundary layer flow. But it became clear from the present research that it was not completely neutral but weakly stable. We examined whether we could simulate the Leipzig Wind Profile by using a ($k-{varepsilon}$) turbulence model including the equation of potential temperature. By solving analytically the Second Moment Closure Model under the assumption of local equilibrium and under the condition of a stratified flow, we expressed the turbulent diffusion coefficients (both momentum and thermal) as functions of flux Richardson number. Our ($k-{varepsilon}$) turbulence model which included the equation of potential temperature and the turbulent diffusion coefficients varying with flux Richardson number reproduced the Leipzig Wind Profile.