f=midigui() ; Select a photometry file spawn,getenv('vltiCbin')+'/oir1dCompressData "'+f+'" '+getenv('prismhmask')+' out' d=oirgetdata('out') indext=where(d.tartyp2 eq 'T ') indexs=where(d.tartyp2 eq 'S ') plot,total(d(indext).data1-d(indexs).data1,2)