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))