Numerical methods for simulating selective laser melting (SLM) have been widely carried out to understand the physical behaviors behind the process. Numerical simulation at the macroscale allows the relationship between input parameters (laser power, scanning speed, powder layer thickness, etc.) and output results (distortion, residual stress, etc.) to be investigated. However, the macroscale thermal models solved by the finite element method cannot predict the melt pool depth correctly as they ignore the effect of fluid flow in the melting pool, especially in the case of the presence of deep penetration. To remedy this limitation, an easy-implemented temperature-dependent heat source is proposed. This heat source can adjust its parameters during the simulation to compensate for these neglected thermal effects related to the fluid flow and keyhole, and the heat source’s parameters become fixed once the temperatures of the points of interest become stable. Contrary to the conventional heat source model, parameters of the proposed heat source do not require a calibration with experiments for each process parameter. The proposed model is validated by comparing its results with those of the anisotropic thermal conductivity method and experimental measurements.