A multi-layer zone fire growth model is developed to predict the vertical distributions of the temperature in a single room. The fire room volume is divided into a number of horizontal layers, in which the temperature and other physical properties are assumed to be uniform. The principal equations for each laminated horizontal layer are derived from the conservation equations of mass and energy. The implemented fire sub-models are introduced, including the combustion, fluid flow and heat transfer models. Combined with these sub-models, the zone equations for the gas temperature of each layer are solved by Runge-Kutta method for each time step. The results of the sample calculations compare well with the results of experiments conducted by Steckler et al.