do_waterfall = True
    do_foldspec = True
    coherent_dedispersion = False

    # with open(file1, 'rb') as fh1:
    # with multifile('/mnt/aro/hdd2_node7/algonquin/sequence.2013-07-24T15:06:16.3.dat',
    #                ['/mnt/aro/hdd2_node7/algonquin/raw_voltage.2013-07-24T15:06:16.0.dat',
    #                 '/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:
    do_waterfall = True
    do_foldspec = True
    coherent_dedispersion = False

    # with open(file1, 'rb') as fh1:
    # with multifile('/mnt/aro/hdd2_node7/algonquin/sequence.2013-07-24T15:06:16.3.dat',
    #                ['/mnt/aro/hdd2_node7/algonquin/raw_voltage.2013-07-24T15:06:16.0.dat',
    #                 '/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_node7/sequence.2013-07-24T16:38:38.3.dat',
                   ['/mnt/data-pen1/pen/njones/VLBI_july212013/hdd2_node7/raw_voltage.2013-07-24T16:38:38.0.dat',
                    '/mnt/data-pen1/pen/njones/VLBI_july212013/hdd1_node7/raw_voltage.2013-07-24T16:38:38.1.dat',
                    '/mnt/data-pen1/pen/njones/VLBI_july212013/hdd3_node7/raw_voltage.2013-07-24T16:38:38.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: