The functions read_flags
and index_good_frames
(written by Roy van Boekel) do the following:
flags=read_flags('tag')
(don't add flag.fits
!)
gd = oirgetdata('tag.groupdelay.fits')
goodframes = index_good_frames(gd,flags)
MIA+EWS> help,goodframes GOODFRAMES INT = Array[16000] MIA+EWS> print,goodframes[0:10] 0 0 1 1 1 1 1 1 0 0 1 MIA+EWS> print,total(goodframes) 9951.00