On simulation of local fluxes in molecular junctions