A vertical two dimensional model for thermo stratified flows is proposed, which solves the velocity and pressure by a fractional step method, and Large Eddy Simulation model is used to calculate the Reynolds stress and turbulent diffusion effect.