i=i+1
ws,0
xyouts,real_part(ir[i]),imaginary(ir[i]),'+',/device,alignment=0.5,chars=1.5
ws,2
plot,w*.1,v[*,real_part(ir[i]),imaginary(ir[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(ir[i])+ii,imaginary(ir[i])+jj],lin=1
oplot,w*.1,vhin[*,real_part(ir[i]),imaginary(ir[i])],thick=3
print,i,ir[i]