import glob
from asap import *

files = glob.glob("skydips/*.rpf")
taus = skydip(files)
scan = scanatable("mydata.rpf")
q = scan.auto_quotient()
b = q.auto_polybaseline(order=2, insitu=False)
av = b.average_time(weight="tsys")
avp = av.average_pol()
corrected = avp.opacity(taus, insitu=False)