[Kwant] how to plot the Charge density with multiple orbitals per site

[Xiangru Kong]

Hi all,

I want to plot the the Charge density using kwant.plot method, I want to know 
how to plot the  Charge density in a system with multiple orbitals per site.

Here is some pesuo codes:

.....
syst = kwant.Builder()
.....
model= kwant.wraparound.wraparound(syst).finalized()
ham=model.hamiltonian_submatrix

(eigen,evecs) = la.eigh(ham)

def site_size(i):
    return 2 * wf[i] / wf.max()

kwant.plot(model, site_size=site_size, site_color=(0, 0, 1, 
0.3),fig_size=(10,10));
*********
In my codes, I just found len(model.sites)=6, that is the number of sites in 
the unit cell, however, due to there are multiple orbitals per site, the 
len(evecs[:,i])=22 orbitals. So I want to know the real infortion when using 
kwant.plot? Like what I just posted, it only plot evecs[:6,i]**2 functions?

Thanks
Xiangru