A turbulence closure of one-and-a-half statistical order, where forecast equations are retained for mean wind and temperature, for wind and temperature variance, and for molecular dissipation of wind and temperature. Other unknowns such as fluxes or variables of higher statistical order are approximated or parameterized as a function of the mean quantities, variances E, and dissipation rates ε.