acf = spec.tocovdata(nr, nt, rate=1)