The TKE (Turbulent Kinetic Energy) closure scheme with 1(1)/(2) order precision and the method of reducing vertical layers are used to decrease computing amount. Modelling results agree well with observations. Compared with results of 3 order closure schemes, there is no obvious reduction of precision. This scheme provides the possibility of inducting and improving the precision of PBL parameterization in three dimensional atmospheric models.