pro fivekeof

bmin = -100
bmax = 100
mag = bytarr(1024,1024)
in = '/synoptic/eof/fd_M_96m_01d/fd_M_96m_01d.004021/fd_M_96m_01d.4021.0006.fits'
dat = rfits(in, head=hd, /scale)
mag = bytscl(dat, min=bmin, max=bmax)

save_dev = !d.name
set_plot, 'z'
tv2, 1024, 1024, /init
erase
tv2, mag
out=tvrd()
outfil = '/home/jeneen/5kmag.gif'
write_gif, outfil, out

device,/close

set_plot, save_dev

end