In radiotherapy, accurately calculated dose distributions of irregularly shaped photon beams are needed. In this paper, an algorithm is presented which enables the calculation of dose distributions due to irregular fields using pencil beam kernels derived from simple basic beam data usually measured on treatment units, i.e. central axis depth-dose curves and profiles. The only extra data that are needed, and are not currently measured, is the phantom scatter factor curve at the reference depth. The algorithm has been developed as an extension to a previously developed algorithm for rectangular fields which is based on the Milan-Bentley storage model. In the case of an irregular field, the depth dose and the boundary function are computed by convolution of a field intensity function with pencil beam kernels. The depth dose is computed by using a 'scatter' kernel, which is derived from the stored depth-dose curves and from the phantom scatter factor curve. The boundary function is computed by using a 'boundary' kernel, which is derived from the boundary profile of a number of large square fields. Because of the simplicity of the data used and the underlying concepts, which for instance do not separate the head scatter from the primary beam, this algorithm presents some shortcomings. On the other hand, this simplicity is also of great advantage and the inaccuracy is acceptable for most clinical situations.