Heat-assisted magnetic recording (HAMR) technology represents the most promising candidate to replace the current perpendicular recording paradigm to achieve higher storage densities. To better understand HAMR dynamics in granular media we need to describe accurately the magnetization dynamics up to temperatures close to the Curie point. To this end we propose a multiscale approach based on the Landau-Lifshitz-Bloch (LLB) equation of motion parametrized using atomistic calculations. The LLB formalism describes the magnetization dynamics at finite temperature and allows us to efficiently simulate large system sizes and long time scales. Atomistic simulations provide the required temperature dependent input quantities for the LLB equation, such as the equilibrium magnetization and the anisotropy and can be used to capture the detailed magnetization dynamics. The multiscale approach makes it possible to overcome the computational limitations of atomistic models in dealing with large systems, such as a recording track, while incorporating the basic physics of the HAMR process. We investigate the magnetization dynamics of a single FePt grain as a function of the properties of the temperature profile and applied field and test the LLB results against atomistic calculations. Our results prove the appropriateness and potential of the approach proposed here where the granular model is able to reproduce the atomistic simulations and capture the main properties of a HAMR medium.