casa6

            if line.strip().endswith('/'):

                break

            continue

        return

​

    def get_spw_ranges(self):

        # Get ranges for each spectral window

        try:

            tb.open(f'{self.vis}/SPECTRAL_WINDOW')

            for spw in self.spws:

                chan_freq = tb.getcol('CHAN_FREQ')[:,spw]

                min_freq = min(chan_freq) / (1e6)

                max_freq = max(chan_freq) / (1e6)

                self.spw_ranges[spw] = (min_freq, max_freq)

            tb.close()

        except Exception as e:

            print("ERROR FAILED TO GET SPW FREQ RANGES")

​

    def check_which_spw(self, freq):

        in_spw = []

        for spw, ranges in self.spw_ranges.items():

            if ranges[0] >= freq[0] and ranges[1] <= freq[1]:

                in_spw.append(spw)

        # If there is no spw that this belongs to default to what we were setting this value to before

        if not in_spw:

            print("WARNING: NO CORRESPONDING SPW FOUND FOR FREQ {}".format(freq))

            return [-1]

        else:

            return in_spw

​

    def process_gc(self, infp, keys):

        antenna_name = find_antenna(keys[0], self.gain_keys)

        if not antenna_name:

            print('Antenna missing from GAIN group')

            skip_values(infp)

            return

        try:

            antenna = get_antenna_index(antenna_name, self.ant_names)

        except:

            print('Antenna %s not present in FITS-IDI file' % antenna_name)

            skip_values(infp)

            return

        keys = dict(keys[0])

​

        try:

​

            self.get_spw_ranges()

            spw_ids = self.check_which_spw(keys['FREQ'])

        except Exception as e:

            spw_ids = [-1]

            print('ANTAB HAS NO FREQ INFORMATION DEFAULTING TO FREQ ID -1')

​

        dpfu = {}

        try:

            dpfu['R'] = keys['DPFU'][0]