Exemplo n.º 1
0
    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))
Exemplo n.º 2
0
    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) /
Exemplo n.º 3
0
    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) /
Exemplo n.º 4
0
    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))