#!/usr/bin/env ParselTongue

from AIPS import AIPS
from Wizardry.AIPSData import AIPSUVData
#from AIPSData import AIPSUVData

 
from pylab import plot, show
from math import atan2
 
AIPS.userno = 702
 
data = AIPSUVData('VC384B', 'UVDATA', 1, 8)
phase = []
time = []
for visibility in data:
    if visibility.baseline == [1, 2]:
        re = visibility.visibility[0, 0, 0, 0]
        im = visibility.visibility[0, 0, 0, 1]
        phase.append(atan2(re, im))
        time.append(visibility.time)
 
plot(time, phase, '.')
show()