With the importance of permanent disposal of high-level radioactive waste (HLW) generated in Korea, the deep geological disposal system based on the KBS-3 type is being developed. Since the deep geological repository must provide the long-term isolation of HLW from the surface environment and normal habitats for humans, plants, and animals, it is essential to assess the longterm performance of the disposal facility considering thermal-hydraulic-mechanical-chemical (TH- M-C) evolution. Decay heat dissipated from HLW contained in the canister causes an increase in temperature in the adjacent area. The requirement for the maximum temperature is established in consideration of the possibility of bentonite degradation. Therefore, when designing the repository, the temperature in the region of interest should be identified in detail through the thermal evolution assessment to ensure that the design requirement is satisfied. In the thermal evolution analysis, it is needed to evaluate the temperature distribution over the entire area of the disposal panel to consider the heat generated from both a single canister and adjacent canisters. Computational fluid dynamics (CFD) codes are widely used for detailed temperature analysis but are limited to simulating a wide range. Accordingly, in this study, we developed an analytical solution-based program for efficiently calculating the temperature distribution throughout the deposition panel, which is based on threedimensional heat conduction equations. The code developed can assess the temperature distribution of engineered and natural barrier systems. Principal parameters to be inputted are as follows: (a) geometry of the panel (e.g. width, length, height, spacing between canisters), (b) geometry of the canister (e.g. diameter, height), (c) thermal properties of bentonite and host-rock, (d) initial conditions (e.g. residual heat, temperature), and (e) time information (e.g. canister emplacement rate, time-interval, period). Through the calculation for the conceptual problem of a deposition panel capable of accommodating 900 (i.e. 30×30) canisters, it was confirmed that the program can adequately predict when and where the maximum temperature will occur. It is expected that the overall temperature distribution within the panel can be obtained by the evaluation of the entire region using this program reflecting the detailed design of the repository to be developed in the future. In addition, the thermal evolution analysis considering the influence of other canisters can be performed by applying the results as boundary conditions in the CFD analysis.