The Discrete Ordinates Interpolation Method (DOIM) is tested in three-dimensional enclosures. The radiative transfer equation (RTE) is solved for a linear source term and the DOIM is formulated for a gray medium. Several interpolation methods can be applied to the DOIM scheme. Among them, the interpolation method applicable to an unstructured grid system is discussed. In a regular hexahedron enclosure, radiative wall heat fluxes are calculated and compared with exact solutions. The enclosure has an absorbing, emitting and nonscattering medium and a constant temperature distribution. These results are obtained with varying optical depths (xD = 0.1, 1.0, 10.0). Also, the same calculations are performed in an irregular hexahedron enclosure. The DOIM is applied to an unstructured grid system as well as a structured grid system for the same regular hexahedron enclosure. They are compared with the exact solutions and the computational efficiencies are discussed. When compared with the analytic solutions, results of the DOIM are in good agreement for three-dimensional enclosures. Furthermore, the DOIM can be easily applied to the unstructured grid system, which proves the reliability and versatility of the DOIM