With the localizing gradient enhancement, a damage model for quasi-brittle materials is able to achieve regularized softening responses, with localized damage profiles corresponding to the development of macroscopic cracks, to resolve the numerical spurious effects induced by the conventional gradient enhancement. The typical implementation strategy for a gradient enhanced model is to solve the system of governing equations simultaneously. Focusing on the finite element (FE) package Abaqus, a user element subroutine is required to define the finite elements with additional degrees of freedom for the nonlocal field. Moreover, with user elements, additional effort is required to visualize the numerical results. To an inexperienced engineer/researcher, these requirements can be challenging. In this paper, a simple implementation of the localizing gradient damage model is elaborated. By utilizing the in-built coupled thermo-mechanical elements in Abaqus, the user only needs to define the material constitutive laws, as well as the sensitivity terms with respect to the field variables. Post-processing of results can be done directly in Abaqus. The applicability and ease of implementation are demonstrated via several examples, including those that utilize the Abaqus features of element deletion, contact between surfaces, as well as the incorporation of cohesive elements. Sample files can be downloaded from https://github.com/leonghien/Localizing-Gradient-Damage-with-UMAT-UMATHT