nx=128 & ny=128
grey=(findgen(nx)/(nx-1))#replicate(1,6)
tvscl,alog(grey>(1.e-4*max(grey)))
bar1=replicate(1,2./0.043)#replicate(1,2)*255 
bar2=replicate(1,1./0.043)#replicate(1,2)*255
bar8=replicate(1,1.6/0.043)#replicate(1,2)*255
;
tvwindow,0,xs=4*nx,ys=6*ny
;
reads,'image1sim',sim1,nx1,ny1
reads,'image1'   ,dat1,nx1,ny1
dat1=dat1>0.
scale2=255./max(sim1) & scale1=255./max(dat1) & scale=min([scale1,scale2])
;
tv,(congrid(dat1*scale,nx,ny)<255),0*nx,5*ny
tv,grey*255,0*nx,5*ny & tv,bar2,0*nx+5,5*ny+10
xyouts,0*nx+8,5*ny+15,'2.',/dev,size=1.
xyouts,0*nx+8,5*ny+15,'2"0',/dev,size=1.
;
tv,(congrid(dat1*scale*4,nx,ny)<255),1*nx,5*ny
tv,(grey*255*4)<255,1*nx,5*ny & tv,bar2,1*nx+5,5*ny+10
xyouts,1*nx+8,5*ny+15,'2.',/dev,size=1.
xyouts,1*nx+8,5*ny+15,'2"0',/dev,size=1.
;
tv,(congrid(sim1*scale,nx,ny)<255),2*nx,5*ny
tv,grey*255,2*nx,5*ny & tv,bar2,2*nx+5,5*ny+10
xyouts,2*nx+8,5*ny+15,'2.',/dev,size=1.
xyouts,2*nx+8,5*ny+15,'2"0',/dev,size=1.
;
tv,(congrid(sim1*scale*4,nx,ny)<255),3*nx,5*ny
tv,(grey*255*4)<255,3*nx,5*ny & tv,bar2,3*nx+5,5*ny+10
xyouts,3*nx+8,5*ny+15,'2.',/dev,size=1.
xyouts,3*nx+8,5*ny+15,'2"0',/dev,size=1.
;
xyouts,2*nx-50,6*ny-15,'PC5 547nm Focus -312',/dev,size=1.
;
; IMAGE 2
;
reads,'image2sim',sim2,nx,ny
reads,'image2'   ,dat2,nx,ny
dat2=dat2>0.
scale2=255./max(sim2) & scale1=255./max(dat2) & scale=min([scale1,scale2])
;
tv,(dat2*scale)<255,0*nx,4*ny
tv,grey*255,0*nx,4*ny & tv,bar1,0*nx+5,4*ny+10
xyouts,0*nx+20,4*ny+15,'2.',/dev,size=1.
xyouts,0*nx+20,4*ny+15,'2"0',/dev,size=1.
;
tv,(dat2*scale*4)<255,1*nx,4*ny
tv,(grey*255*4)<255,1*nx,4*ny & tv,bar1,1*nx+5,4*ny+10
xyouts,1*nx+20,4*ny+15,'2.',/dev,size=1.
xyouts,1*nx+20,4*ny+15,'2"0',/dev,size=1.
;
tv,(sim2*scale)<255,2*nx,4*ny
tv,grey*255,2*nx,4*ny & tv,bar1,2*nx+5,4*ny+10
xyouts,2*nx+20,4*ny+15,'2.',/dev,size=1.
xyouts,2*nx+20,4*ny+15,'2"0',/dev,size=1.
;
tv,(sim2*scale*4)<255,3*nx,4*ny
tv,(grey*255*4)<255,3*nx,4*ny & tv,bar1,3*nx+5,4*ny+10
xyouts,3*nx+20,4*ny+15,'2.',/dev,size=1.
xyouts,3*nx+20,4*ny+15,'2"0',/dev,size=1.
;
xyouts,2*nx-50,5*ny-15,'PC5 547nm Focus +342',/dev,size=1.
;
; IMAGE 3
;
reads,'image3sim',sim3,nx,ny
reads,'image3'   ,dat3,nx,ny
dat3=dat3>0.
mn=nx/2-nx/16 &mx=nx/2+nx/16
scale2=255./max(sim3) & scale1=255./max(dat3) & scale=min([scale1,scale2])
;
tv,(congrid(dat3(mn:mx,mn:mx),nx,ny)*scale1)<255,0*nx,3*ny
tv,grey*255,0*nx,3*ny & tv,bar8,0*nx+5,3*ny+10
xyouts,0*nx+15,3*ny+15,'0"2',/dev,size=1.
xyouts,0*nx+15,3*ny+15,'0.',/dev,size=1.
;
tv,(dat3*scale*50                          )<255 ,1*nx,3*ny
tv,(grey*255*50)<255,1*nx,3*ny & tv,bar1,1*nx+5,3*ny+10
xyouts,1*nx+20,3*ny+15,'2"0',/dev,size=1.
xyouts,1*nx+20,3*ny+15,'2.',/dev,size=1.
;
tv,(congrid(sim3(mn:mx,mn:mx),nx,ny)*scale2)<255 ,2*nx,3*ny
tv,grey*255,2*nx,3*ny & tv,bar8,2*nx+5,3*ny+10
xyouts,2*nx+15,3*ny+15,'0"2',/dev,size=1.
xyouts,2*nx+15,3*ny+15,'0.',/dev,size=1.
;
tv,(sim3*scale*50                          )<255 ,3*nx,3*ny
tv,(grey*255*50)<255,3*nx,3*ny & tv,bar1,3*nx+5,3*ny+10
xyouts,3*nx+20,3*ny+15,'2"0',/dev,size=1.
xyouts,3*nx+20,3*ny+15,'2.',/dev,size=1.
;
xyouts,2*nx-50,4*ny-15,'PC6 486nm Focus -10',/dev,size=1.
;
; IMAGE 4
;
reads,'image4sim',sim4,nx,ny
reads,'image4'   ,dat4,nx,ny
dat4=dat4>0.
mn=nx/2-nx/16 &mx=nx/2+nx/16
scale2=255./max(sim4) & scale1=255./max(dat4) & scale=min([scale1,scale2])
;
tv,(congrid(dat4(mn-1:mx-1,mn-2:mx-2),nx,ny)*scale1)<255 ,0*nx,2*ny
tv,grey*255,0*nx,2*ny & tv,bar8,0*nx+5,2*ny+10
xyouts,0*nx+15,2*ny+15,'0"2',/dev,size=1.
xyouts,0*nx+15,2*ny+15,'0.',/dev,size=1.
;
tv,(dat4*scale*50                                  )<255 ,1*nx,2*ny
tv,(grey*255*50)<255,1*nx,2*ny & tv,bar1,1*nx+5,2*ny+10
xyouts,1*nx+20,2*ny+15,'2"0',/dev,size=1.
xyouts,1*nx+20,2*ny+15,'2.',/dev,size=1.
;
tv,(congrid(sim4(mn:mx,mn:mx),nx,ny)*scale2        )<255 ,2*nx,2*ny
tv,grey*255,2*nx,2*ny & tv,bar8,2*nx+5,2*ny+10
xyouts,2*nx+15,2*ny+15,'0"2',/dev,size=1.
xyouts,2*nx+15,2*ny+15,'0.',/dev,size=1.
;
tv,(sim4*scale*50                                  )<255 ,3*nx,2*ny
tv,(grey*255*50)<255,3*nx,2*ny & tv,bar1,3*nx+5,2*ny+10
xyouts,3*nx+20,2*ny+15,'2"0',/dev,size=1.
xyouts,3*nx+20,2*ny+15,'2.',/dev,size=1.
;
xyouts,2*nx-50,3*ny-15,'PC6 230nm Focus -10',/dev,size=1.
;
; IMAGE 5
;
reads,'image5sim',sim5,nx5,ny5
reads,'image5n'   ,dat5,nx5,ny5
dat5s=smooth(dat5,2)>0.
mn=nx5/2-nx5/16 &mx=nx5/2+nx5/16
scale2=255./max(sim5) & scale1=255./max(dat5) & scale=min([scale1,scale2])
;
tv,(congrid(dat5(mn:mx,mn-1:mx-1),nx,ny)*scale1)<255 ,0*nx,1*ny
tv,grey*255,0*nx,1*ny & tv,bar8,0*nx+5,1*ny+10
xyouts,0*nx+15,1*ny+15,'0"2',/dev,size=1.
xyouts,0*nx+15,1*ny+15,'0.',/dev,size=1.
;
tv,(congrid(dat5s*scale*50,nx,ny)               )<255 ,1*nx,1*ny
tv,(grey*255*50)<255,1*nx,1*ny & tv,bar1,1*nx+5,1*ny+10
xyouts,1*nx+20,1*ny+15,'2"0',/dev,size=1.
xyouts,1*nx+20,1*ny+15,'2.',/dev,size=1.
;
tv,(congrid(sim5(mn:mx,mn:mx),nx,ny)*scale2    )<255 ,2*nx,1*ny
tv,grey*255,2*nx,1*ny & tv,bar8,2*nx+5,1*ny+10
xyouts,2*nx+15,1*ny+15,'0"2',/dev,size=1.
xyouts,2*nx+15,1*ny+15,'0.',/dev,size=1.
;
tv,(congrid(sim5*scale*50,nx,ny)               )<255 ,3*nx,1*ny
tv,(grey*255*50)<255,3*nx,1*ny & tv,bar1,3*nx+5,1*ny+10
xyouts,3*nx+20,1*ny+15,'2"0',/dev,size=1.
xyouts,3*nx+20,1*ny+15,'2.',/dev,size=1.
;
xyouts,2*nx-50,2*ny-15,'F96 486nm Focus -10',/dev,size=1.
;
; IMAGE 6
;
reads,'image6sim',sim6,nx6,ny6
reads,'image6'   ,dat6,nx6,ny6
dat6s=smooth(dat6,2)>0.
mn=nx6/2-nx6/16 &mx=nx6/2+nx6/16
scale2=255./max(sim6) & scale1=255./max(dat6) & scale=min([scale1,scale2])
;
tv,(congrid(dat6(mn-1:mx-1,mn-4:mx-4),nx,ny)*scale1)<255 ,0*nx,0*ny
tv,grey*255,0*nx,0*ny & tv,bar8,0*nx+5,0*ny+10
xyouts,0*nx+15,0*ny+15,'0"2',/dev,size=1.
xyouts,0*nx+15,0*ny+15,'0.',/dev,size=1.
;
tv,(congrid(dat6s*scale*50,nx,ny)           )<255 ,1*nx,0*ny
tv,(grey*255*50)<255,1*nx,0*ny & tv,bar1,1*nx+5,0*ny+10
xyouts,1*nx+20,0*ny+15,'2"0',/dev,size=1.
xyouts,1*nx+20,0*ny+15,'2.',/dev,size=1.
;
tv,(congrid(sim6(mn:mx,mn:mx),nx,ny)*scale2)<255 ,2*nx,0*ny
tv,grey*255,2*nx,0*ny & tv,bar8,2*nx+5,0*ny+10
xyouts,2*nx+15,0*ny+15,'0"2',/dev,size=1.
xyouts,2*nx+15,0*ny+15,'0.',/dev,size=1.
;
tv,(congrid(sim6*scale*50,nx,ny)           )<255 ,3*nx,0*ny
tv,(grey*255*50)<255,3*nx,0*ny & tv,bar1,3*nx+5,0*ny+10
xyouts,3*nx+20,0*ny+15,'2"0',/dev,size=1.
xyouts,3*nx+20,0*ny+15,'2.',/dev,size=1.
;
xyouts,2*nx-50,1*ny-15,'F96 120nm Focus -10',/dev,size=1.
;
; Vertical bars
;
tv,replicate(128,1,128*6),1*128-1,0
for j=0,5 do tv,replicate(200,1,105),2*128-1,j*128
tv,replicate(128,1,128*6),3*128-1,0
psprint,'plate2'
stop
end