# '/mnt/aro/hdd1_node7/algonquin/raw_voltage.2013-07-24T15:06:16.1.dat',
# '/mnt/aro/hdd3_node7/algonquin/raw_voltage.2013-07-24T15:06:16.2.dat']) as fh1:
#with multifile('/mnt/aro/hdd2_node7/algonquin/sequence.2013-07-27T16:55:17.3.dat',
# ['/mnt/aro/hdd2_node7/algonquin/raw_voltage.2013-07-27T16:55:17.0.dat',
# '/mnt/aro/hdd1_node7/algonquin/raw_voltage.2013-07-27T16:55:17.1.dat',
# '/mnt/aro/hdd3_node7/algonquin/raw_voltage.2013-07-27T16:55:17.2.dat']) as fh1:
with multifile('/mnt/data-pen1/pen/njones/VLBI_july212013/hdd2_node9/sequence.2013-07-25T18:14:20.3.dat',
['/mnt/data-pen1/pen/njones/VLBI_july212013/hdd2_node9/raw_voltage.2013-07-25T18:14:20.0.dat',
'/mnt/data-pen1/pen/njones/VLBI_july212013/hdd1_node9/raw_voltage.2013-07-25T18:14:20.1.dat',
'/mnt/data-pen1/pen/njones/VLBI_july212013/hdd3_node9/raw_voltage.2013-07-25T18:14:20.2.dat']) as fh1:
foldspec2, waterfall = fold(fh1, '4bit', samplerate,
fedge, fedge_at_top, nchan,
nt, ntint, nhead,
ngate, ntbin, ntw, dm, fref, phasepol,
coherent=coherent_dedispersion,
do_waterfall=do_waterfall,
do_foldspec=do_foldspec,
verbose=verbose, progress_interval=1)
if do_waterfall:
np.save("aro{}waterfall.npy".format(psr), waterfall)
if do_foldspec:
np.save("aro{}foldspec2.npy".format(psr), foldspec2)
f2 = foldspec2.copy()
f2[0] = 0.
foldspec1 = f2.sum(axis=2)
fluxes = foldspec1.sum(axis=0)
foldspec3 = f2.sum(axis=0)
dedisperse = "by-channel"
# with open(file1, 'rb') as fh1:
with multifile(seq_file, raw_files) as fh1:
foldspec2, waterfall = fold(
fh1,
"4bit",
samplerate,
fedge,
fedge_at_top,
nchan,
nt,
ntint,
nhead,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
dedisperse=dedisperse,
do_waterfall=do_waterfall,
do_foldspec=do_foldspec,
verbose=verbose,
progress_interval=1,
)
if do_waterfall:
np.save("aro{}waterfall.npy".format(psr), waterfall)
if do_foldspec:
with multifile(
'/scratch/aro/hdd2_node7/algonquin/sequence.0329_VLBI_june30.dat',
[
'/scratch/aro/hdd2_node7/algonquin/raw_voltage.0329_VLBI_june30.1.dat',
'/scratch/aro/hdd1_node7/algonquin/raw_voltage.0329_VLBI_june30.2.dat',
'/scratch/aro/hdd3_node7/algonquin/raw_voltage.0329_VLBI_june30.3.dat'
]) as fh1:
foldspec2, waterfall = fold(fh1,
np.int8,
samplerate,
fedge,
fedge_at_top,
nblock,
nt,
ntint,
nhead,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
do_waterfall=do_waterfall,
verbose=verbose,
progress_interval=1)
f2 = foldspec2.copy()
f2[0] = 0.
foldspec1 = f2.sum(axis=2)
fluxes = foldspec1.sum(axis=0)
foldspec3 = f2.sum(axis=0)
if igate is not None:
do_foldspec = True
dedisperse = 'by-channel'
#with open(file1, 'rb') as fh1:
with multifile(seq_file, raw_files) as fh1:
foldspec2, waterfall = fold(fh1,
'4bit',
samplerate,
fedge,
fedge_at_top,
nchan,
nt,
ntint,
nhead,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
dedisperse=dedisperse,
do_waterfall=do_waterfall,
do_foldspec=do_foldspec,
verbose=verbose,
progress_interval=1)
if do_waterfall:
np.save("aro{}waterfall.npy".format(psr), waterfall)
if do_foldspec:
np.save("aro{}foldspec2.npy".format(psr), foldspec2)
[
"/scratch/aro/hdd2_node7/algonquin/raw_voltage.0329_VLBI_june30.1.dat",
"/scratch/aro/hdd1_node7/algonquin/raw_voltage.0329_VLBI_june30.2.dat",
"/scratch/aro/hdd3_node7/algonquin/raw_voltage.0329_VLBI_june30.3.dat",
],
) as fh1:
foldspec2, waterfall = fold(
fh1,
np.int8,
samplerate,
fedge,
fedge_at_top,
nblock,
nt,
ntint,
nhead,
ngate,
ntbin,
ntw,
dm,
fref,
phasepol,
do_waterfall=do_waterfall,
verbose=verbose,
progress_interval=1,
)
f2 = foldspec2.copy()
f2[0] = 0.0
foldspec1 = f2.sum(axis=2)
fluxes = foldspec1.sum(axis=0)
foldspec3 = f2.sum(axis=0)