This paper is a study on an accurate and computationally efficient method for reconstructing pointwise power distributions from coarse-mesh nodal calculations. The modern nodal codes can calculate global reactor power shapes and criticality very efficiently and accurately. But inherent in the nodal procedures, there is inevitable loss of information on local heterogeneous quantities. In this study, an improved form function method which reflects the exponential transition of the thermal flux near the assembly surface is developed for the reconstruction of the heterogeneous fluxes. Use of the new form function method in several pressurized water reactor (PWR) benchmark problems reduces the maximum errors in the reconstructed thermal flux to those in the reconstructed fast flux. Even for assemblies adjacent to the steel baffle in realistic PWR cores, use of this method also results in improved pointwise power reconstruction.