Created: May 09, 2017 - JHZ
Updated: July 10, 2017 - JHZ

##########################################################################
Note:
Testing data:  from JVLA observations of NRAO530 (calibrator) & Sgr A (target) -

Processes:
Time-dependent & antenna-based delay correction derived from NRAO530 and apply to Sgr A -
##########################################################################
    Split day-1 data into 14 sets (scan-based)

    in Miriad:
    set target=nrao530_ad1
    uvaver vis=$target.uv options=nocal,nopass select='time('5:50,6:00')' \
           out=$target.uv_1
    echo '1 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('6:00,6:20')' \
           out=$target.uv_2
    echo '2 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('6:20,6:40')' \
           out=$target.uv_3
    echo '3 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('6:40,6:50')' \
           out=$target.uv_4
    echo '4 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('7:00,7:10')' \
           out=$target.uv_5
    echo '5 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('7:10,7:30')' \
           out=$target.uv_6
    echo '6 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('7:30,7:50')' \
           out=$target.uv_7
    echo '7 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('7:50,8:00')' \
           out=$target.uv_8
    echo '8 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('8:00,8:20')' \
           out=$target.uv_9
    echo '9 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('8:20,8:40')' \
           out=$target.uv_10
    echo '10 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('8:40,9:00')' \
           out=$target.uv_11
    echo '11 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('9:00,9:10')' \
           out=$target.uv_12
    echo '12 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('9:10,9:30')' \
           out=$target.uv_13
    echo '13 next'
    $<
    uvaver vis=$target.uv options=nocal,nopass select='time('9:30,9:50')' \
           out=$target.uv_14
    echo '14 end'

    Compute 14 sets of BP solutions in time sequence for each of 16 spectral windows

    set maxstep=14
    set i=0
    while ($i < $maxstep)
        @ i++
    smamfcal vis=$target.uv_$i select='source('$sour')' edge=5,5 \
         refant=1 interval=5
    #     options=msmooth smooth=3,1,0.9 bcut=10
    $<
    GPPLT:
    smagpplt vis=$target.uv_$i device=/xs \
            yaxis=amp \
            options=bandpass  \
            nxy=1,10 yrange=0.9,1.1
    $<
    smagpplt vis=$target.uv_$i device=/xs \
            yaxis=phase \
            options=bandpass  \
            nxy=1,10 yrange=-5,5
    echo 'target.uv_'$i
    end
    exit

    Interpolate the BP solutions and apply to the Sgr A data

    \rm -r $target.uv.bp
    smatbpass bpfile=$target.uv nfiles=14 vis=sgra_ad1.uv out=sgra_ad1.uv.bp \
            bptime=5
    echo 'done with time-dependent bpass/delay corrections!'
    exit