# Plot amp v time for calibrator:
uvplt vis=0228+673 line=wide,1,2 device=/xw axis=time,amp

# Plot phase v time for calibrator:
uvplt vis=0228+673 line=wide,1,2 device=/xw axis=time,phase

	--> Ant. 3 and Ant. 9 not so hot.

# Phase calibrate calibrator.
selfcal vis=0228+673 line=wide,1,2 interval=20 refant=1

# Look at it.  Make a map for kicks.
uvplt vis=0228+673 line=wide,1,2 device=/xw axis=time,phase
uvplt vis=0228+673 line=wide,1,2 device=/xw axis=time,phase average=5
invert vis=0228+673 map=0228+673.mp beam=0228+673.bm line=wide,1,2 cell=2 imsize=128
cgdisp in=0228+673.mp device=/xw
cgdisp in=0228+673.bm device=/xw

# Apply gains from calibrator to source.
gpcopy vis=0228+673 out=IC10

# Make dirty map of source.  Look at it.
invert vis=IC10 options=systemp line=velocity,40,-350,1,1 map=junk.mp beam=junk.bm cell=2 imsize=128
cgdisp in=junk.mp device=/xw options=3value csize=1,1,1,1
cgdisp in=junk.mp device=/xw options=3value csize=1,1,1,1 'region=images(18,21)'

# Clean and restore
clean map=junk.mp beam=junk.bm out=junk_c.mp
cgdisp in=junk_c.mp device=/xw options=3value csize=1,1,1,1
restor model=junk_c.mp beam=junk.bm map=junk.mp mode=clean out=junk_cr.mp
cgdisp in=junk_cr.mp device=/xw

# Make a velocity map; spectra...
velplot in=junk_cr.mp device=/xw

# Now try combining both tracks.  First, bring over night 1's data and call
# it IC10_1.  Well that was easy.  So, it's already phase calibrated; jump
# right to invert, using both data sets.  Maybe this worked?

invert vis=IC10,IC10_1 options=systemp line=velocity,40,-350,1,1 map=supah.mp beam=supah.bm cell=2 imsize=128
cgdisp in=supah.mp device=/xw options=3value csize=1,1,1,1
clean map=supah.mp beam=supah.bm out=supah_c.mp
restor model=supah_c.mp beam=supah.bm map=supah.mp mode=clean out=supah_cr.mp
cgdisp in=supah_cr.mp device=/xw options=3value csize=1,1,1,1 'region=images(18,18)'