samplerate = 200 * u.MHz
fedge = 200. * u.MHz
fedge_at_top = True
fref = 150. * u.MHz # ref. freq. for dispersion measure
# convert time to UTC; dates given in EDT
time0 = Time(date_dict[psr], scale='utc') + 4*u.hr
verbose = 'very'
do_waterfall = True
do_foldspec = True
dedisperse = 'incoherent'
with multifile(seq_file, raw_files) as fh1:
if not isinstance(phasepol, Polynomial):
from astropy.utils.data import get_pkg_data_filename
from pulsar.predictor import Polyco
polyco_file = get_pkg_data_filename(phasepol)
polyco = Polyco(polyco_file)
phasepol = polyco.phasepol(time0, rphase='fraction', t0=time0,
time_unit=u.second, convert=True)
nskip = int(round(
((Time('2013-07-25T22:15:00', scale='utc') - time0) /
(recsize * 2 / samplerate)).to(u.dimensionless_unscaled)))
if verbose:
print("Using start time {0} and phase polynomial {1}"
.format(time0, phasepol))
samplerate = 200 * u.MHz
fedge = 200. * u.MHz
fedge_at_top = True
fref = 150. * u.MHz # ref. freq. for dispersion measure
# convert time to UTC; dates given in EDT
time0 = Time(date_dict[psr], scale='utc') + 4 * u.hr
verbose = 'very'
do_waterfall = True
do_foldspec = True
dedisperse = 'incoherent'
with multifile(seq_file, raw_files, comm=comm) as fh1:
if not isinstance(phasepol, Polynomial):
from astropy.utils.data import get_pkg_data_filename
from pulsar.predictor import Polyco
polyco_file = get_pkg_data_filename(phasepol)
polyco = Polyco(polyco_file)
phasepol = polyco.phasepol(time0,
rphase='fraction',
t0=time0,
time_unit=u.second,
convert=True)
nskip = int(
round(
((Time('2013-07-25T22:15:00', scale='utc') - time0) /
samplerate = 200 * u.MHz
fedge = 200. * u.MHz
fedge_at_top = True
fref = 150. * u.MHz # ref. freq. for dispersion measure
# convert time to UTC; dates given in EDT
time0 = Time(date_dict[psr], scale='utc') + 4 * u.hr
verbose = 'very'
do_waterfall = True
do_foldspec = True
dedisperse = 'incoherent'
with multifile(seq_file, raw_files) as fh1:
if not isinstance(phasepol, Polynomial):
from astropy.utils.data import get_pkg_data_filename
from pulsar.predictor import Polyco
polyco_file = get_pkg_data_filename(phasepol)
polyco = Polyco(polyco_file)
phasepol = polyco.phasepol(time0,
rphase='fraction',
t0=time0,
time_unit=u.second,
convert=True)
nskip = int(
round(
((Time('2013-07-25T22:15:00', scale='utc') - time0) /
samplerate = 200 * u.MHz
fedge = 200. * u.MHz
fedge_at_top = True
fref = 150. * u.MHz # ref. freq. for dispersion measure
# convert time to UTC; dates given in EDT
time0 = Time(date_dict[psr], scale='utc') + 4*u.hr
verbose = 'very'
do_waterfall = True
do_foldspec = True
dedisperse = 'incoherent'
with multifile(seq_file, raw_files, comm=comm) as fh1:
if not isinstance(phasepol, Polynomial):
from astropy.utils.data import get_pkg_data_filename
from pulsar.predictor import Polyco
polyco_file = get_pkg_data_filename(phasepol)
polyco = Polyco(polyco_file)
phasepol = polyco.phasepol(time0, rphase='fraction', t0=time0,
time_unit=u.second, convert=True)
nskip = int(round(
((Time('2013-07-25T22:15:00', scale='utc') - time0) /
(recsize * 2 / samplerate)).to(u.dimensionless_unscaled)))
if verbose:
print("Using start time {0} and phase polynomial {1}"
.format(time0, phasepol))