Developed for prediction of turbulent flux in second order closure equation, turbulent flux can be calculated by relation between turbulent kinetic energy (TKE) and dissipation equation. Four cases and 20 experiments show the length scale depends on the value of input data. Stability function increases significantly. Also, a covariance temperature has influence on stability function for heat (Sh) and  The turbulent kinetic energy has the strongest influence with heat flux, momentum flux in x- and y-directions, and dissipation.

planetary boundary layer, turbulence, WRF.