The increasing use of irregularly shaped, off-centre fields in advanced treatment techniques, particularly intensity modulated radiation therapy, has strained the limits of conventional, broad-beam dose calculation algorithms. More recent models, such as kernel-based pencil beams and Monte Carlo methods, are accurate but suffer from the time needed for calculations and from the lack of clearly established methods for determining the parameters needed to match calculations with the particular dosimetric characteristics of an individual machine. This paper presents the implementation of a model that uses an extended source model to calculate the variation of fluence at the patient surface for any arbitrarily shaped field. It uses a macropencil beam model to calculate phantom scatter. Both head scatter and phantom scatter models use exponential functions fit to a series of measurements to determine the model's parameters. The means by which the model can be implemented in a clinical setting using standard dosimetric equipment is presented. Results for two separate machines and three energies are presented. Comparisons with measurements for a set of regular and irregular fields demonstrate the accuracy of the model for conventional, conformal and intensity modulated treatments. For rectangular and irregular fields at depths up to 20 cm, the accuracy was better than < or =1.5%, compared with errors of up to 7.5% with a standard algorithm. For a 20-step intensity modulated field, the accuracy was 3.4% compared with 18% with the conventional algorithm. The advantages of this model for IMRT are discussed.