All spectra with mode ID modeid will be summed by SUM_BDA with the following command: IDL> data_out = sum_bda(index, data, modeid, nsum)
