© 2019 Elsevier B.V. Plasma heat flux in the peripheral plasma reaches the first wall (FW) along a magnetic field line, and it sometimes causes several MW/m2 orders of magnitude high heat flux concentration at narrow region, such as the edge of the blanket module. Thus, to assess and to reduce the heat load is key issue in the DEMO design activity. In this research, a new heat load analysis code is introduced based on the e-folding model. In this code, the decay length λ is changed depending on wall connection length of magnetic field line (defined as the length of magnetic field line from the wall to the wall), and parallel heat flux q// is calculated in each flux tube. This code can calculate the FW heat load simulating actual blanket module shapes. The 0.23 MW/m2 peak heat load at inner midplane in the case of the ideal FW (without gaps between blanket modules) is increased to 22 MW/m2 at toroidal module edge in the case of box shaped module. To shadow the edge and reduce such peak heat load, toroidal, and poloidal roof shaping is applied. Required roof height is analyzed from this code calculation. After shaping, peak heat load is reduced to 1.2 MW/m2. This value is under the allowable value 1.5 MW/m2, and in this case, surface temperature is also less than allowable temperature.