i=i+1
ws,0
xyouts,real_part(ib[i]),imaginary(ib[i]),'o',/device,alignment=0.5,chars=1.5
ws,2
plot,w*.1,v[*,real_part(ib[i]),imaginary(ib[i])],xtit='!7k!5 (nm)',chars=1.5,ytit='V/I',lin=2,thick=3,yr=[-.4,.4]
for ii=-3,3,2 do begin for jj=-3,3 do oplot,w*.1,v[*,real_part(ib[i])+ii,imaginary(ib[i])+jj],lin=1
oplot,w*.1,vhin[*,real_part(ib[i]),imaginary(ib[i])],thick=3
print,i,ib[i]