image = intarr(64,64,29)
;
;filnam = '/home/soi/data/nov664/lev0/
filnam = '/home/soi/data/nov667/lev0/
for i=0,28 do begin
    infil = filnam + string(i, format='(i4.4)') + '.fits'
    img = rfits(infil)
    image(0,0,i) = rebin(img, 64, 64) 
end
;
for i=0,28 do tvscl,rebin(image(*,*,i), 128, 128), i
stop
mtmtune,image,m1,m2,trans

end