EVLA Data reduction, (AR708_5), on J123707 CO following notes 
from Jeff Wagg

on aida125
===============================
UVLOD
-First download the data from archive in FITS format (uvfits)

-Read the data in (takes a while!)
default 'FITLD'; douvcomp -1; 
datain '/disk1/hodge/RAWData/AR708_5Ka.uvfits 
outn 'AR708J1237';

-> read in without error
-> user number 708, for the project code
--------------------------------
# UVSRT (to put data in time order, which is necessary for INDXR)
# (Also takes a while)
# 
# task 'UVSRT'; sort 'T'; outn ''; outcl ''; outdi 1; rotate 0; defer 0;
# go;
# 
-------------------------------
INDXR (to create NX and CL tables)

task 'indxr'; infi ''; prtlev 0; 
cparm 1 10 0.5; bparm 0

*I used diff values than Jeff b/c I'm guessing they've fixed the online
opacities now. But what he used for bparm is below
# cparm 1 10 0.5; bparm 20 1 0 0 0 0 0 0 0, -1

-> I got the error:
INDXR1: Warning:   0.165 percent of times not found in WX table
INDXR1: This may be a normal amount

--------------------------------------
PRTAN 
get central antenna:

task 'PRTAN'; docrt 124; invers 1;

refant: 2 
------------------------
LISTR
get summary of observations

task 'LISTR'; optype 'SCAN'; sources ''; calcode ''; docrt 124

Source = J123707+6214

CALIBRATORS:

primary calib: J1331+3030 = 3C286

secondary cal: J1302+5748

1302+578   J2000  A 13h02m52.465282s  57d48'37.609420"  Aug01  CJ2    
1300+580   B1950  A 13h00m47.131400s  58d04'43.188000"
-----------------------------------------------------
BAND        A B C D    FLUX(Jy)    UVMIN(kL)  UVMAX(kL)
=====================================================
 20cm    L  S S X X       0.31            15         visplot
  6cm    C  S S S S       0.52                       visplot
3.7cm    X  S S S S       0.89
  2cm    U  W W W W       0.4
0.7cm    Q  W W W W       0.4 

--------------------------------------
VLANT
-to update estimates of antenna positions

task 'vlant'; default

-> made CL 2
-----------------------------
POSSM (on the flux calibrator)
-Plot phase vs. freq channel to check for phase slopes/wraps due to bad 
delays in the EVLA data 

task 'possm'; default; doband -1; docal 1; gainuse 0;
bchan 0; echan 0
sources 'J1331+3030' ''
calcod ''; antenna 2 0
aparm 1 1 0 0, -200 200
codetype 'PHAS'
stokes 'RR'
bif 1; eif 1
nplots 9; dotv 1

-> antennas 11,13,14,16,17,18,19,26 have no fringes

----------------------------
UVFLG
-flag antennas with no fringes

task 'uvflg'; default; antenn 11 13 14 16 17 18 19 26
basel 0; sour ''; timer 0; stokes ''; bif 0; eif 0
opcod 'flag'; aparm 0; bchan 0; echan 0; reason ''

-> made FG 1
-----------------------------
PRTUV
print some uv data to get integration times on calibrator/source:
(which we will need for FRING)

task 'PRTUV'; sources ''; cparm 0; dparm 0; cparm(9)=204; docrt 124

   calib: 1s
  source: 1s
---------------------------
FRING
-to fix the delays using the bandpass calibrator

task 'fring'; default
calsour 'J1331+3030' '';
docal 1; gainuse 1; flagver 0;
refant 2
dparm 1 1000 100 1 0 0 0 0 1 0
timer 0 10 51 0 0 10 52 0
solint 1.2

-> Made SN table 1

-------------------------------
CLCAL
-To take the correct delays from the SN Table 1 and apply them to 
the CL table made by INDXR, producing a new CL table)

task 'clcal'; default; calcod ''
source ''; calsour 'J1331+3030' ''; samptyp 'box'; bparm 24 24
interpol ''; snver 1; invers 1; gainver 0; gainuse 0;
refant 2

-> Made CL table 3
-----------------------------
POSSM
-Check that the phases are now flat

task 'possm'; default; doband -1; 
sources 'J1331+3030' ''
calcod ''; antenna 2 0
aparm 1 1 0 0 0 -200 200
codetype 'PHAS'
stokes 'RR'
bif 1; eif 1
nplots 9; dotv 1
docal 1; gainuse 3

-> flat

---------------------------------------
UVFLG
-add flags for shadowing

task 'uvflg'; default;
aparm(5)=25; opcod 'FLAG'

-> wrote 62 flags
------------------------------------
QUACK
-flag first few secs of data

task 'quack'; opcod 'BEG'; aparm 0 0.08333; sourc ''

-> *It said it flagged 11.33 min of data in 68 scans, but when I look at 
listr, opcod 'list', flagver=1, nothing is different, and there are
no new entries in the FG table (using PRTAB)...will have to flag by hand
below

--------------------------------------
TVLFG
-look for data that needs to be flagged

default 'tvflg'; docal 1; doband -1; bchan 1; echan 64; 
nchav 4; bif 1; eif 2; calcod ''; timer 0; stokes ''
sources 'J1331+3030' 'J1302+5748';
flagver 1; outfgver 1; docat -1; dohist -1
gainuse 0; 
dparm(6)=1

*When loading, gave error:
TVFLG1: WARNING: SUBARRAY  1 SOME ANTENNAS LACK POSITION INFO

*and also:
TVFLG1: GRIDTB WARNING: using ECHAN   67 not input ECHAN   64

-> flagged:
anten 1 ~10:36, IF1
anten 3, ~10:20, IF 1,2
anten 5, ~9:38, IF 1,2
anten 6, ~6:47, IF1
anten 15, ~8:15, IF1,2
anten 23, ~6:05, IF1
anten 24, ~6:05, IF1,2
anten 25, ~7:12, IF1,2
anten 27, ~8:52, IF1
I also flagged the ends of two scans near the end which looked bright

UVFLG
-Also flagged the following with uvflg, which I noticed in a previous
iteration:
anten 22 0; basel 0; timer 0; bif 1; eif 0; stokes ''; sour ''; aparm 0; reason ''
anten 6 0; basel 23 0; bif 2; eif 2; timer 0; stokes 'rl'; bchan 62; echan 64; 

--------------------------------------
SETJY
-to set flux density of primary calibrator

task 'SETJY'; sources 'J1331+3030' ''; optype 'CALC'; freqid=1; calcode '';
aparm 0,0; bif 0; eif 0

SETJY1: 'J1331+3030      ' IF =  1 FLUX = 1.9181 (Jy calcd)
SETJY1: 'J1331+3030      ' IF =  2 FLUX = 1.9125 (Jy calcd)

-> Why are these different than Jeff's values???? 
He got:
SETJY1: 'J1331+3030 ' IF = 1 FLUX = 1.8732 (Jy calcd)
SETJY1: 'J1331+3030 ' IF = 2 FLUX = 1.8679 (Jy calcd)

---------------------------
CALIB
*change refant

-Delete SN1 (from delay cal since it was already applied)
task 'extdest'; inext 'sn'; invers 1

Task 'CALIB'; Calsour 'J1331+3030' ''; Qual -1; Calcode '';
Selband -1; Selfreq -1; Freqid 1; Timerang 0; ichansel 10 54 1 0;
Dofit 0; Antuse 0; Subarray 0; Uvrange 0 0; Wtuv 0; 
Docalib 1; Gainuse 0; flagver 0; Doband -1; Bpver -1; Smooth 0;
In2name ''; In2class ''; In2seq 0; In2disk 0; Invers 0; Ncomp 0; 
Flux 0; Nmaps 0; Cmethod ''; Cmodel ''; Smodel 0; Outname ''; 
Outclass ''; Outseq 0; Outdisk 0; Solint 0;  doflag 0
Aparm 0; Soltype ''; Solmode 'A&P'; Solcon 0; Minamper 0;
Minphser 0; cparm 0; Snver 0; Antwt 0; Gainerr 0; baddisk 0
uvrange 0 ; anten 0; refant 2

-> made SN table 1 
CALIB1: RPOL, IF= 1 The average gain over these antennas is  9.359E+00
CALIB1: RPOL, IF= 2 The average gain over these antennas is  1.231E+01
CALIB1: LPOL, IF= 1 The average gain over these antennas is  9.774E+00
CALIB1: LPOL, IF= 2 The average gain over these antennas is  1.243E+01
CALIB1: Found           68 good solutions
CALIB1: Average closure rms =  0.00009 +-  0.00002

tget calib; calsour 'J1302+5748','';

-> made SN table 2 
CALIB1: RPOL, IF= 1 The average gain over these antennas is  9.545E+00
CALIB1: RPOL, IF= 2 The average gain over these antennas is  1.262E+01
CALIB1: LPOL, IF= 1 The average gain over these antennas is  9.930E+00
CALIB1: LPOL, IF= 2 The average gain over these antennas is  1.259E+01
CALIB1: Found         2472 good solutions
CALIB1: Failed on        4 solutions
CALIB1: Average closure rms =  0.00033 +-  0.00007

-------------------------------
GETJY
get flux of secondary calib.

-Then get new fluxes for phase calibrator:
task 'GETJY'; Sources 'J1302+5748',''; Soucode '' ; 
Calsour 'J1331+3030',''; Qual -1; Calcode ''; Bif 0; Eif 0; Timerang 0;
Antennas 0; Subarray 0 ; Selband -1; Selfreq -1; Freqid 1; Snver 0

GETJY1: J1302+5748      :  0         1    0.35819 +/-   0.00504  0.00672
GETJY1:                              2    0.35822 +/-   0.00510  0.00710

-> WHOA, twice what Jeff got?? He got:
GETJY1: J1302+5748 : 0 1 0.18532 +/- 0.00264 0.00343
GETJY1:                2 0.18479 +/- 0.00261 0.00352

(Since my setjy results were also diff than his, I tried an experiment 
and set my values to his, but my getjy results are still the same.)

---------------------------------------
CLCAL   
-interpolate gain table

Task 'CLCAL'; Sources ''; Soucode ''; 
Calsour 'J1302+5748',''; Qual -1; Calcode ''; Timerang 0; Subarray 0;
Antennas 0; Selband -1; Selfreq -1; Freqid 1; Opcode ''; 
Interpol '2PT'; Cutoff 0; Samptype ''; Bparm 0;
Snver 0; inver 0; gainver 0; gainuse 0; refant 2

-> made CL table 4
---------------------------------------
LISTR
-Look at CL table to verify calibration of 2ndry calibrator(s)

task 'listr'; sources 'J1302+5748' ''; optyp 'matx'; docalib 1; 
inext 'cl'; gainuse 0; dparm  5, 1, 0; stokes ''

-> Phases are beautiful, but amplitudes are horrible.....