I am trying to model a single spark discharge for electric discharge machining process in abaqus using DFLUX subroutine. The heat flux for my model is Gaussian distributed for an axisymmetric model. The problem is I am not getting any temperature profile after running the subroutine. Can anyone please help if there is any error in my subroutine or any problem in my model. I am attaching the DFLUX subroutine. Thank you.