**** SCRIPT for AW602 ******
$  Fill in the data
$  REMEMBER to <setenv AIPSPP `pwd`> in directory of data files

TASK 'FILLM'
   default
   doconcat 1
   vlaobs = 'AW602'
   band = 'Q'
$$   next means include data within 1000 khz in one IF
   cparm(7) 1000.0
   cparm(8) 0.25
   outn 'AW602'; outs 1; outd 1
   infile 'AIPSPP:AW602_B030502.xp
   ncount 3
   go;wait
   infile 'AIPSPP:AW602_C030508.xp
   ncount 2
   go;wait

$  Copy the two FQ files into separate aips data bases

TASK 'UVCOP'
   default
   inn 'AW602'; incl 'LINE'; ins 1; indi 1
   outn 'AW602'; outcl 'freq'; outs 1; outd 1
   freqid = 1; qual = 1
   go; wait
   outs 2
   freqid = 2; qual = 2
   go; wait

TASK 'INDXR'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   cparm 5 5 0
   go; wait
   ins 2
   go; wait

$  Set the flux density of 3C147

TASK 'SETJY'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   optype 'calc'
   go; wait
   ins 2
   go;wait

$  Clip points above 2 Jy

TASK 'CLIPM'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   aparm 2.0, 2.0
   go;wait
   ins 2
   go;wait

$  Various flags found during previous reductions

TASK 'UVFLG'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   ante 0 0; timer = 0 21 27 10 0 21 48 30; wait;go
   ante 10 0; timer = 6 0 0 0 7 0 0 0; wait;go
   ante 22 0; timer = 0 20 06 00 0 20 35 00; wait;go
   ante 0 0; timer 0 0 0 0 10; bchan 1; echan 3;wait;go
   ante 0 0; timer 0 0 0 0 10; bchan 58; echan 63;wait;go
   ins 2
   ante 0 0; timer = 0 21 27 10 0 21 48 30; wait;go
   ante 10 0; timer = 6 0 0 0 7 0 0 0; wait;go
   ante 22 0; timer = 0 20 06 00 0 20 35 00; wait;go
   ante 0 0; timer 0 0 0 0 10; bchan 1; echan 3;wait;go
   ante 0 0; timer 0 0 0 0 10; bchan 58; echan 63;wait;go

$  Do the initial G calibration on three calibrators

TASK 'CALIB'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   cals '0319+415' '0336+323' '0542+498'  
   docal 1
   flagver 0
   refant 15
   solint 1
   snver 1
   go;wait
   ins 2
   go;wait

$  Apply this calibration

TASK 'CLCAL'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   flagver 1
   refant 15
   gainver 1;gainuse 2;snver 1;
   go;wait
   ins 2
   go;wait

$  Determine the bandpass using 3C84

TASK 'BPASS'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   refant 15
   cals(1) '0319+415'
   docal 1; gainuse 2  
   flagver 1
   solint -1
   bpass(5) = -1
   go;wait
   ins 2
   go;wait

$  Now, do G calibration using Bandpass

TASK 'CALIB'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   cals '0319+415' '0336+323' '0542+498'  
   docal 1; gainuse 2
   doband 1; bpver 1
   flagver 0
   refant 15
   solint 1
   snver 2
   go;wait
   ins 2
   go;wait

$  Get flux density of 3C84 and 0336

TASK 'GETJY'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   source '0319+415' '0336+323'
   cals(1) '0542+498'
   snver 2
   go;wait
   ins 2
   go;wait

$  Apply this final calibration to data

TASK 'CLCAL'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   flagver 1
   refant 15
   gainver 2;gainuse 3;snver 2;
   go; wait

$  Image and clean 0336 in two data bases

TASK 'IMAGR'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   source(1) '0336+323'
   qual 1
   docal 1; gainuse 0
   flagver 1
   doband 1
   bpver 1
   nchav 63
   outn source(1)
   outs 0; outd 1
   in2n inn; in2cl 'uvwork'; in2d 2; in2s 0
   cellsize 0.5
   imsize 512
   uvtaper 120
   uvwtfn 'na'
   niter 200
   dotv -1
   go;wait
   ins 2; qual 2
   go; wait

$  Image all target fields for both frequencies
$  dirty images only, combine three channels for image

TASK 'IMAGR'
   default
   inn 'AW602'; incl 'freq'; ins 1; indi 1
   docal 1; gainuse 0
   flagver 1
   doband 1
   bpver 1
   bchan 4; echan 57
   nchav 3
   chinc 3
   outd 1
   in2n inn; in2cl 'uvwork'; in2d 2; in2s 0
   cellsize 0.5
   imsize 512
   uvtaper 120
   uvwtfn 'na'
   niter 00
   dotv -1
   source(1) '03288+31150'
   outn source(1)
   ins 2; outs 2; qual 2
   go;wait
   source(1) '03288+31149'
   outn source(1)
   ins 2; outs 2; qual 2
   go;wait
   source(1) '03289+31148'
   outn source(1)
   ins 2; outs 2; qual 2
   go;wait
   source(1) '03289+31147'
   outn source(1)
   ins 2; outs 2; qual 2
   go;wait
   source(1) '03289+31146'
   outn source(1)
   ins 2; outs 2; qual 2
   go;wait
   source(1) '03289+31146'
   outn source(1)
   ins 1; outs 1; qual 1
   go;wait
   source(1) '03290+31145'
   outn source(1)
   ins 1; outs 1; qual 1
   go;wait
   source(1) '03290+31144'
   outn source(1)
   ins 1; outs 1; qual 1
   go;wait
   source(1) '03290+31143'
   outn source(1)
   ins 1; outs 1; qual 1
   go;wait
   source(1) '03291+31142'
   outn source(1)
   ins 1; outs 1; qual 1
   go;wait