Time average separate scans

If you have observed the source with multiple source/reference cycles you will want to scan-average the quotient spectra together.

 ASAP>av = q.average_time()

If for some you want to average multiple sets of scantables together you can:

 ASAP>av = average_time(q1, q2, q3)

The default is to use integration time weighting. The alternative is to use none, variance, Tsys weighting, Tsys & integration time or median averaging.

 ASAP>av = average_time(q, weight='tintsys')

To use variance based weighting, you need to supply a mask saying which channel range you want it to calculate the variance from.

 ASAP>msk = scans.create_mask([200,400],[600,800])
 ASAP>av = average_time(scans, mask=msk, weight='var')

If you have not observed your data with Doppler tracking (or run freq_align explicitly) you should align the data in frequency before averaging.

 ASAP>av = scans.average_time(align=True)

Note that, if needed, you should run gain_el and opacity before you average the data in time (§ & ).