How lammps calculates single-atom heat flux

Dear all,
I want to obtain a two-dimensional local heat flux, which can be achieved by a single-atom heat flux distribution. How can I calculate the single-atom heat flux in lammps?

Thanks!

What can be done with LAMMPS is described in the documentation. If you want to compute something that is not documented, then you will either have to write new C++ code that implements that functionality, or find and use a different MD code that has this feature, or try to compute this property in post-processing, if possible.