f=midigui()
tag='hd107446'
oir1dcompressdata,f,getenv('prismhmask'),tag
d=getoirdata(tag+'.compressed.fits')
plot,d(500:900).data1(100)

plot,d(500:900).data1(100)-d(500:900).data2(100)