verbose = True
do_waterfall = True
foldspecs = []
waterfalls = []
for P in range(20):
file1 = file_fmt.format(S=0, P=P)
file2 = file_fmt.format(S=1, P=P)
fbottom = fwidth*(563.+P*20)
f2, wf = fold(file1, file2, '>f4',
fbottom, fwidth, nchan,
nt, ntint, nskip,
ngate, ntbin, ntw, dm, fref, phasepol,
do_waterfall=do_waterfall,
verbose=verbose, progress_interval=1)
np.save("lofar{}foldspec2{}.npy".format(psr, P), f2)
foldspecs.append(f2)
waterfalls.append(wf)
foldspec2 = np.concatenate(foldspecs, axis=0)
waterfall = np.concatenate(waterfalls, axis=0)
np.save("lofar{}foldspec2.npy".format(psr), foldspec2)
np.save("lofar{}waterfall.npy".format(psr), waterfall)
f2 = foldspec2.copy()
foldspec1 = f2.sum(axis=2)
for P in range(20):
file1 = file_fmt.format(S=0, P=P)
file2 = file_fmt.format(S=1, P=P)
fbottom = fwidth * (563. + P * 20)
f2, wf = fold(file1,
file2,
'>f4',
fbottom,
fwidth,
nchan,
nt,
ntint,
nskip,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
coherent=coherent,
do_waterfall=do_waterfall,
verbose=verbose,
progress_interval=1)
np.save("lofar{}foldspec2{}.npy".format(psr, P), f2)
foldspecs.append(f2)
waterfalls.append(wf)
foldspec2 = np.concatenate(foldspecs, axis=0)
file1 = file_fmt.format(S=0, P=P)
file2 = file_fmt.format(S=1, P=P)
fbottom = fwidth * (563.0 + P * 20)
f2, ic, wf = fold(
file1,
file2,
">f4",
fbottom,
fwidth,
nchan,
nt,
ntint,
nskip,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
coherent=coherent,
do_waterfall=do_waterfall,
verbose=verbose,
progress_interval=1,
)
np.save("lofar{}foldspec2{}.npy".format(psr, P), f2)
foldspecs.append(f2)
icounts.append(ic)
waterfalls.append(wf)