So you want a contour plot of the defect density. That's pretty essential information - how could anyone know you didn't mean a contour plot of the electron density, or electrostatic potential, or something else...
Yes, you can do this kind of study using ATKClassical, since ATK contains the essential ingredients needed: the Tersoff potential by de Brito Mota, the capability to do MD (also the melt and quench preparation step), and to compute the local coordination of atoms, i.e. produce the data that can be turned into the contour plot.
But it's not a single script to do all this, it's a big scientific study.