#                 '/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)
Ejemplo n.º 2
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:
Ejemplo n.º 3
0
    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:
Ejemplo n.º 4
0
    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)
Ejemplo n.º 5
0
        [
            "/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)