reads,'image3',a,nx,ny
tvwindow,0,xs=nx*5,ys=ny*5
tvscl,congrid(alog(a>(1.e-4*max(a))),5*nx,5*ny)
grey=(findgen(5*nx)/(5*nx-1))#replicate(1,20)
tvscl,alog(grey>(1.e-4*max(grey)))
bar=replicate(1,5/0.043)#replicate(1,4)*255
tv,bar,20,40
xyouts,15,50,'1 arcsec',/dev,size=2.8
psprint,'plate1'
stop
end