pro myzhao1013
device,retain=2,decompose=0
window,0,xs=1024,ys=1024
file=file_search('G:\1013\-0.5_+0.9\*.fit',count=n0)
for i=0,n0-1 do begin
m0=swap_endian(readfits(file(i),mh))
wavelength=sxpar(mh,'WAVELNTH')
print,wavelength
m1=congrid(m0,1024,1024)
tvscl,m1
stop
wait,1
endfor
end