Hi Barış,
you can for instance make the lead_onsite different from onsite and then print syst[lat(i,j)] for different regions to check. 
Regards,
Ousmane

Le sam. 6 juil. 2019 à 22:43, Barış Özgüler <ozguler@wisc.edu> a écrit :
Hi Ousmane,

How can I make sure that there is no overlap between the lead and normal region? Lead_shape() returns -L <= x <= L and we want the lead to satisfy x <= 0. If we don't use color_site(site), then we would get the attached plot. The black region not only represents the normal region but also represents the lead. However, the lead may also cover the region of x > 0 as the function Lead_shape() suggests. Is there another way to check what part of the shape is scattering region, what part is lead?

Thanks,

-- 
A. Barış Özgüler

PhD Candidate in Physics

https://home.physics.wisc.edu/ozguler/

University of Wisconsin–Madison
Madison, Wisconsin 53706, USA