AIPS Reduction of AM858: VLA Obs. of 17mar06
============================================
mpr 21mar06-22mar06

***finished first round***

Still to do:

THE BASICS:
  On-line disk FILLM w/ doweight=1 ; douvcomp=-1 ; cparm(8)=1/6; bparm=0
    ==> set up for calibration of weights (docal=2), 10s CL table entries,
      and DEFAULT GAIN/OPACITY CORRECTIONS
  17mar06 is MJD 53811.0
  A configuration
  Integration time: 3.3sec for all

  Known problems:
    * Wind/API rms phase/clouds:
      13:31  0.3 m/s    1.6d    70% stratiform clouds
    * ant. 14,16 positions recently updated
    * ant. 13,18 moved recently -- lousy positions
    ***EVLA ant. 13,14,16,18 included***
    * EVLA ant. 18 working only on IFs AC (IF 1) at 4cm
    * 13:35-14:02  correlator errors in quad 4, array 2 -- had to re-seat
      multiplier card
    * Ptg. failures:
      13:51 1924+156  all ok!
      14:15 1331+305  all ok!
      14:22 2007+404  mostly bad (all but 16!)

                      Location Of VLA Antennas
 
                               N72 ( 1)
                               N56 (27)
                               N48 (20)
                               N40 ( 2)
                               N32 ( 7)
                               N24 (26)
                               N16 (14)*
                               N8  ( 9)
                               N6  (18)*
                           ( 6) W8   E8  (10)
                         ( 8) W16      E16 (17)
                       (23) W24          E24 (16)*
                     (11) W32              E32 (19)
                   (  )                      E40 (15)
                 ( 4) W48                      E48 (12)
               ( 5) W56                          E56 ( 3)
             (21) W64                              E64 (25)
           (22) W72                                  E72 (28)
                          EVLA:MP  (13)
                          VLA: OUT (24)
                          VPT: OUT (29)
* => EVLA ANTENNA
MP=> Master Pad

  EVLA evaluation:
    4cm: no data on ant. 18 except 1L
      ant. 18 IF 1L phase is very odd: with 3.3s records, we get
        one value, no data, phase jumps by 100d, then repeats
        amp. is basically steady, and consistent with others
      all great, except ant. 13 last scan (0 14 10 - 0 14 15)
      raw amps. all fine eccept ant. 13 vanishes in last scan (factor 50)
    1.3cm:  no ant. 13,18
      otherwise great phases
      raw amps.: ant. 14-VLA fine
                 ant. 16-VLA and EVLA low by 1.4-1.6

===============================================================================

----
 4cm
----
*Analyzed 22mar06*

0. IF 1: 8435.1 GHz  50 MHz
   IF 2: 8485.1 GHz  50 MHz

  The Sources
  -----------
  1331+305  flux cal.  3C286
    elev yy
     4cm  uvra  50, 300   A: inner 3 ant./arm; ABC: reduce flux by 1%
     4cm  14:11:40-14:13:20 --> 0.59201 days  1.7mins

  19247+15407 aka 1924+156 phase cal. for 19152+10567 aka G1915+105  5.3 degs
    elev yy
     4cm  uvra  0,  0 ; secondary ; 0.44 Jy
     4cm  13:36:03-13:49:10 --> 0.57135 days  5.7mins
  19152+10567 aka G1915+105  5.3 degs.
    elev yy
     4cm  13:34:20-13:48:30 --> 0.57049 days  9.2mins

0. EVLA tests
  a. VPLOT phase vs. time
    docal -1 ; calcode '*' ; nplot 9; bparm 0,2 ; dotv 1
      flagver -1; do3col 1
      ante 13, 14, 16, 18,  1,  2,  3, 0
      base ante
    stokes 'rr'; bif 0; eif 0
    stokes 'll'; bif 0; eif 0
    --> no data on ant. 18 except 1L
      ant. 18 IF 1L phase is very odd: with 3.3s records, we get
        one value, no data, phase jumps by 100d, then repeats
        amp. is basically steady, and consistent with others
      all great, except ant. 13 last scan (0 14 10 - 0 14 15)
      raw amps. all fine eccept ant. 13 vanishes in last scan (factor 50)

1. Skipped PCAL transfer for now (no recent solutions yet)

2. TVFLG -> FG/1
  source '' ; docal -1 ; flagver 1 ; docat -1 ; calcode '*' ; dparm(6) 3.33333
  dparm(3) 1
  : EVLA ant. 13: random phases on 1331+305
  : EVLA ant. 18: random phases on IF 1L -- seems to bop up and down
    --> phase diff < 19.2 degrees on 1331+305 (ant. 21, 22 -- end of W arm)
                   < 12.3 degrees on 1924+156

3. SETJY 1331+305
  sour '1331+305','' ; opty 'calc'; bif 0 ; eif 0 ; freqid 1
  zerosp 0 ; calcode '' ; aparm 0
  SETJY1: BIF =  1 EIF =  2 /Range of IFs
  SETJY1: '1331+305        ' IF =  1 FLUX = 5.2158 (Jy calcd)
  SETJY1: '1331+305        ' IF =  2 FLUX = 5.1947 (Jy calcd)
  SETJY1: / Using (1999.2) VLA or Reynolds (1934-638) coefficients
  SETJY1: Appears to have ended successfully
  SETJY1: scamper      31DEC06 TST: Cpu=       0.0  Real=       0

4a. CALRD
  object '3c286' ; band 'X'

4b. CALIB CL/1, FG/1 -> SN/1-2
  USE DOCAL=1, per Juan
  a. calsour '1331+305','' --> SN/1
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante  0
    flagver 1
    solint 0 ; solmode ''
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5

    in2na '3c286_X' ; in2cla 'model' ; in2se 1 ; inver 1 ; ncomp 1e5
    nmap 1
    cmethod 'dft'

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are OK
    -->   100 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0008 +-  0.0001

  b. calsour '19247+15407','' --> SN/2
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante 0
    flagver 1
    solint 0 ; solmode ''
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are OK
    -->   612 good sol'ns
      No closure errors
      Avg. closures up to 1.1%/0.5d
      CALIB1: Average closure rms =  0.0012 +-  0.0002

5. GETJY SN/1-2
  calsour '1331+305',''
  sour '19247+15407',''
  freqid 1
  snver 0
  ante -13,-14,-16,-18,0
  --> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
      GETJY1: 19247+15407     :  0   A     1    0.78879 +/-   0.00128
      GETJY1:                              2    0.78819 +/-   0.00124
      GETJY1: Appears to have ended successfully
      GETJY1: scamper      31DEC06 TST: Cpu=       0.0  Real=       0

5b. EVLA checks
  a. LISTR/GAIN
    --> amp. fine
    --> phase wrap on 13 (80d in 14mins)

6. CLCAL CL/1, SN/1-2 -> CL/2
  a. Gains for all calibrators 
    snver 0 ; gainver 1 ; gainuse 2
    calsour '1331+305','19247+15407',''
    calcode ' '
    sour calsour
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'self'
    dobtween -1
    cutoff 10
    --> CL/2

  b. snver 0 ; gainver 1 ; gainuse 2
    calsour '19247+15407',''
    calcode ' '
    sour '19152+10567',''
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'simp'
    cutoff 10
    --> CL/2

7. TVFLG CL/2, FG/1 -> FG/1
  docal 2 ; gainuse 2 ; flagver 1 ; sour '' ; docat -1
  dparm 0 ; dparm(6) 3.3333333 ; dparm(3) 1
  a. calcode '*'
    --> amp. diff. < 132 mJy (EVLA ant. 14 LCP higher than the rest)
    --> phase < 12 degs.

  b. calcode '-cal'
    --> amp. diff. <  95 mJy

8. SNPLT
  inext 'ty' ; bif 0 ; eif 0 ; stokes '' ; xinc 1
  nplots 9 ; opcode 'alsi'
  a. TSYS nominal sensitivities
    opty 'tsys' ; xaxis 0
    SKIPPED FOR NOW
    --> spiffy: 0.yy-0.yy mostly
      ant. yy: 

  b. TANT
    opty 'tant' ; xaxis 2
    SKIPPED FOR NOW
    --> yy-yyy K, for elev. yy to yy, except
      ant. yy: yyy-yyy K

8c. VPLOT/weights
  docal 2 ; gainuse 2 ; flagver 1 ; nplot 9 ; bparm 0 16 0 ; doarray 1
  bif 1 ; eif 2 ; avgif -1 ; do3col 1 ; stokes 'RR' ; ante 17 0
    --> weights 650-1100, except
      ant. 12: 580-770
      EVLA ant. 13: RCP 950-980 (missing 3c286, which is lower weight for
        all others)
                    LCP 830 IF 1L; 940 IF 2L
                  
      EVLA ant. 14: RCP 640-700 IF 1R, 560-590 IF 2R
                    LCP 300-430, as expected
      EVLA ant. 16: 550-900
      ant. 19: 450-650

9. Pol'n calibration -> postponed until PCAL works and/or we can transfer

10. TASAV -> XTASAV.1

11. SPLIT -> * 811.XSPLIT
  docal 2 ; gainuse 2 ; DOPOL -1 ; flagver 1 ; outcla 'Xsplit' ; outseq 0
  sour '' ; calcode '-cal' ; DOUVCOMP -1

12. IMAGR
  a. GRS1915  811.XSPLIT -> XBM001, XCL001  5.3degs.
    (-1) UVFLG EVLA ants. -> FG/1
      flagver 1 ; reason 'evla ant'; opcode 'flag' ; dohist 1
      ante  13, 14, 16, 18,0

    (0) Outseq XXX: The Big Picture -- skipped, based on B config.
      
    (1) Outseq 1: Robust -0.5, boxed CLEAN 
      FLAGVER 1
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.02 ; imsize 2048 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 6 ; clbox 0
        clbox(1,1)~  -1      3   1026   1022
        clbox(1,2)~1010    995   1018   1005
        clbox(1,3)~1013   1003   1024   1015
        clbox(1,4)~1019   1013   1027   1023
        clbox(1,5)~1023   1024   1031   1030
        clbox(1,6)~1025   1029   1031   1035
      do3d 1 ; overlap 2
      niter  1000
      "noise" box: blc 600,590 ; trc 1450,940
      -->  25491 vis.
        221 x 193 mas at   6.8 deg
        WTNOISE 1.311
        CLEANed  13.121 mJy
        -0.456 to 11.882 mJy/beam (20,20-2020,2020: -0.418,11.882 mJy/beam;
          peak at 1026,1022)
        rms  80.3 microJy/beam, from histogram

        Fine detection
        TVMAX pixr -0.36e-3,0.36e-3
          AIPS 1: Map pixel position:  1025.60 1022.06
          AIPS 1: Skypos: RA 19 15 11.54752  DEC 10 56 44.7082
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =   11.933 Milli JY/BEAM
          
          AIPS 1: Map pixel position:  1017.98 1005.52
          AIPS 1: Skypos: RA 19 15 11.55788  DEC 10 56 44.3773
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  469.418 Micro JY/BEAM
          
          AIPS 1: Map pixel position:  1022.39 1010.57
          AIPS 1: Skypos: RA 19 15 11.55188  DEC 10 56 44.4783
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  925.177 Micro JY/BEAM
        JMFIT: blc  978, 967 ; trc 1074,1077
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gwid 0 ; gpos 0 
          dowid -1,-1,-1, 1, 1, 1
          niter 2000
        JMFIT: blc  978, 967 ; trc 1074,1077
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gwid 0 ; gpos 0 
          dowid  1, 1, 1, 1, 1, 1
          niter 2000
        IMEAN: blc  978, 967 ; trc 1074,1077

        --> 11.882 +/- 0.0803 mJy/bm; int. mJy; TVSTAT int. 13.016 mJy
            [pk @ 1026,1022]                                    [1114 pix.]

          Corr. for JMFIT zero: +/- 0.0803 mJy/bm; int. mJy;
            TVSTAT int. mJy

    (2) Outseq 2: Robust 0.5, boxed CLEAN 
      FLAGVER 1
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust 0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.02 ; imsize 2048 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 4 ; clbox 0
        clbox(1,1)~  -1      3   1026   1022
        clbox(1,2)~1012    999   1019   1007
        clbox(1,3)~1016   1008   1022   1016
        clbox(1,4)~  -1      4   1030   1034  $ this is very iffy
      do3d 1 ; overlap 2
      niter  1000
      "noise" box: blc 600,590 ; trc 1450,940
      -->  25491 vis.
        247 x 210 mas at  -6.4 deg
        WTNOISE 1.081
        CLEANed  13.245 mJy
        -0.411 to 11.899 mJy/beam (20,20-2020,2020: -0.351,11.899 mJy/beam;
          peak at 1026,1022)
        rms  67.3 microJy/beam, from histogram

        Lovely detection
        TVMAX pixr -0.30e-3,0.30e-3
          AIPS 1: Map pixel position:  1025.60 1022.05
          AIPS 1: Skypos: RA 19 15 11.54753  DEC 10 56 44.7080
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =   11.944 Milli JY/BEAM

          AIPS 1: Map pixel position:  1022.45 1011.58
          AIPS 1: Skypos: RA 19 15 11.55181  DEC 10 56 44.4986
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =    1.978 Milli JY/BEAM
          
          AIPS 1: Map pixel position:  1016.26 1003.50
          AIPS 1: Skypos: RA 19 15 11.56021  DEC 10 56 44.3371
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  556.273 Micro JY/BEAM
        JMFIT: blc  978, 967 ; trc 1074,1077
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gwid 0 ; gpos 0 
          dowid -1,-1,-1, 1, 1, 1
          niter 2000
        JMFIT: blc  978, 967 ; trc 1074,1077
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gwid 0 ; gpos 0 
          dowid  1, 1, 1, 1, 1, 1
          niter 2000
        IMEAN: blc  978, 967 ; trc 1074,1077

        --> 11.899 +/- 0.0673 mJy/bm; int. mJy; TVSTAT int. 13.045 mJy
            [pk @ 1026,1022]                                    [1301 pix.]

          Corr. for JMFIT zero: +/- -1.0673 mJy/bm; int. mJy;
            TVSTAT int. mJy

================================================================================

-----
1.3cm
-----
*Analyzed 22mar06*

0. IF 1: 22485.1 GHz
   IF 2: 22435.1 GHz 

  The Sources
  -----------
  1331+305  flux cal.  3C286
    elev yy
   1.3cm  uvra   0, 185   inner 2 ant./arm
   1.3cm  14:16:20-14:17:20 --> 0.59502 days  1.1mins

  19247+15407 aka 1924+156 phase cal. for 19152+10567 aka G1915+105  5.3 degs.
    elev yy
   1.3cm  uvra  0,  0 ; secondary ; 0.40-0.65 Jy
   1.3cm  13:54:03-14:07:30 --> 0.58387 days  6.0mins
  19152+10567 aka G1915+105  5.3 degs.
    elev yy
   1.3cm  13:52:20-14:09:00 --> 0.58386 days 10.8mins

  20077+40298 aka 2007+404 phase cal. for 20324+40575 aka CygX-3  4.7 degs.
    elev yy
   1.3cm  uvra  0,  0 ; secondary ; 2.3-1.1  Jy
   1.3cm  14:25:23-14:31:13 --> 0.60301 days  2.9mins
  20324+40575 aka CygX-3  4.7 degs.
    elev yy
   1.3cm  14:23:40-14:30:20 --> 0.60213 days  4.7mins

0. EVLA tests
  a. VPLOT phase vs. time
    docal -1 ; calcode '*' ; nplot 9; bparm 0,2 ; dotv 1
      flagver -1; do3col 1
      ante 13, 14, 16, 18,  1,  2,  3, 0
      base ante
    stokes 'rr'; bif 0; eif 0
    stokes 'll'; bif 0; eif 0
    
    -->  no ant. 13,18
      otherwise great phases
      raw amps.: ant. 14-VLA fine
                 ant. 16-VLA and EVLA low by 1.4-1.6

1. Skipped PCAL transfer for now (no recent solutions yet)

2. TVFLG -> FG/1
  source '' ; docal -1 ; flagver 1 ; docat -1 ; calcode '*' ; dparm(6) 3.33333
  dparm(3) 1
    --> phase diff < 23.9 degrees on 1924+156
                   < 31.2 degrees on 1331+305
                   < 24.0 degrees on 2007+404

3. Elevation dependences: already done in FILLM!

4. SETJY 1331+305
  sour '1331+305','' ; opty 'calc'; bif 0 ; eif 0 ; freqid 1
  zerosp 0 ; calcode '' ; aparm 0
  SETJY1: BIF =  1 EIF =  2 /Range of IFs
  SETJY1: '1331+305        ' IF =  1 FLUX = 2.5170 (Jy calcd)
  SETJY1: '1331+305        ' IF =  2 FLUX = 2.5214 (Jy calcd)
  SETJY1: / Using (1999.2) VLA or Reynolds (1934-638) coefficients
  SETJY1: Appears to have ended successfully
  SETJY1: scamper      31DEC06 TST: Cpu=       0.0  Real=       0

5a. CALRD
  object '3c286' ; band 'K'

5b. CALIB phases: CL/1, FG/1 -> SN/1
  USE DOCAL=1, per Juan
  a. calsour '1331+305','' --> SN/1
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante  0
    flagver 1
    solint 1/3 ; solmode 'p'
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5

    in2na '3c286_K' ; in2cla 'model' ; in2se 1 ; inver 1 ; ncomp 1e5
    nmap 1
    cmethod 'dft'
    snver 1

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are OK
    -->   287 good sol'ns
            1 failed
      No closure errors
      CALIB1: Average closure rms =  0.0150 +-  0.0033

  b. calsour '19247+15407','20077+40298','' --> SN/2
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante 0
    flagver 1
    solint 1/3 ; solmode 'p'
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5
    snver 1

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are OK
    -->  1800 good sol'ns
      No closure errors
      Avg. closures up to 1.4d on 1924+156
      CALIB1: Average closure rms =  0.0210 +-  0.0071

5c. SNPLT to check phases: SN/1
  inext 'sn' ; inver 1 ; pixr -180,180 ; nplot 9 ; xinc 1
  opty 'phas' ; opcode 'alsi' ; do3col 1 ; dotv 1
  --> looks ok! 

6. CLCAL phases: CL/1, SN/1 -> CL/2
  a. Gains for all calibrators 
    snver 1 ; gainver 1 ; gainuse 2
    calsour '1331+305','19247+15407','20077+40298',''
    calcode ' '
    sour calsour
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'self'
    dobtween -1
    cutoff 10
    --> CL/2

  b. snver 1 ; gainver 1 ; gainuse 2
    calsour '19247+15407',''
    calcode ' '
    sour '19152+10567',''
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'simp'
    cutoff 10
    --> CL/2

  c. snver 1 ; gainver 1 ; gainuse 2
    calsour '20077+40298',''
    calcode ' '
    sour '20324+40575',''
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'simp'
    cutoff 10
    --> CL/2

7. CALIB amplitudes: CL/2, FG/1 -> SN/2
  USE DOCAL=1, per Juan
  a. calsour '1331+305','' --> SN/1
    docal 1 ; gainuse 2
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante  0
    flagver 1
    solint 0 ; solmode ''
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5

    in2na '3c286_K' ; in2cla 'model' ; in2se 1 ; inver 1 ; ncomp 1e5
    nmap 1
    cmethod 'dft'
    snver 2

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are now OK
    -->   100 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0011 +-  0.0000

  b. calsour '19247+15407','20077+40298','' --> SN/2
    docal 1 ; gainuse 2
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante 0
    flagver 1
    solint 0 ; solmode ''
    soltype 'l1'
    cparm 0 ; cparm(5) 0
    aparm 0 ; aparm(6) 2 ; minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
    refant 10 (E8)
    doflag -2.5
    snver 2

    NOTE that I use VECTOR averaging (cparm(5)=0), since phases are OK
    -->   900 good sol'ns
      No closure errors
      Avg. closures up to 1.6%/0.9d on 1924+156
      CALIB1: Average closure rms =  0.0017 +-  0.0001

8. GETJY SN/2
  calsour '1331+305',''
  sour '19247+15407','20077+40298',''
  freqid 1
  snver 2
  ante -13,-14,-16,-18,0
  --> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
      GETJY1: 19247+15407     :  0   A     1    0.75241 +/-   0.00472
      GETJY1:                              2    0.75228 +/-   0.00487
      GETJY1: 20077+40298     :  0   A     1    1.87302 +/-   0.01704
      GETJY1:                              2    1.87386 +/-   0.01707
      GETJY1: Appears to have ended successfully
      GETJY1: scamper      31DEC06 TST: Cpu=       0.0  Real=       0

8b. EVLA checks
  a. LISTR/GAIN
    EVLA ant. 14 gains higher on 2007+404 -- elevation effect?
      also high gains in LCP
    EVLA ant. 16 much higher gains than the rest, and even more so on 2007+404

    No sign of mis-pointing.  Hmm.

8b. SNPLT to check amplitudes: SN/2
  inext 'sn' ; inver 2 ; pixr 0 ; nplot 9 ; xinc 1
  opty 'amp' ; opcode 'alsi' ; do3col 1 ; dotv 1
  --> OK except as above

9. CLCAL A&P: CL/2, SN/2 -> CL/3
  a. Gains for all calibrators 
    snver 2 ; gainver 2 ; gainuse 3
    calsour '1331+305','19247+15407','20077+40298',''
    calcode ' '
    sour calsour
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'self'
    dobtween -1
    cutoff 10
    --> CL/2

  b. snver 2 ; gainver 2 ; gainuse 3
    calsour '19247+15407',''
    calcode ' '
    sour '19152+10567',''
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'simp'
    cutoff 10
    --> CL/2

  c. snver 2 ; gainver 2 ; gainuse 3
    calsour '20077+40298',''
    calcode ' '
    sour '20324+40575',''
    soucode ' '
    freqid 1
    refant 10 (E8)
    interpol 'simp'
    cutoff 10
    --> CL/2

10. TVFLG CL/3, FG/1 -> FG/1
  docal 2 ; gainuse 3 ; flagver 1 ; sour '' ; docat -1
  dparm 0 ; dparm(6) 3.3333333 ; dparm(3) 1
  a. calcode '*'
    --> amp. diff. < 302 mJy
    --> phase < 24.1 degs.
      ...EVLA ant. 16 noisier than the rest

  b. calcode '-cal'
    : IF 1L huge amp. on most baselines (up to 762 Jy) 14:03:30-13:03:40
      inclusive -- just after correlator error report
    --> amp. diff. < 311 mJy
    --> ph. diff. < 40.7d on CygX-3

11. SNPLT
  inext 'ty' ; bif 0 ; eif 0 ; stokes '' ; xinc 1
  nplots 9 ; opcode 'alsi'
  a. TSYS nominal sensitivities
    opty 'tsys' ; xaxis 0
    SKIPPED FOR NOW
    --> spiffy: 0.yy-0.yy mostly
      ant. yy: 

  b. TANT
    opty 'tant' ; xaxis 2
    SKIPPED FOR NOW
    --> yy-yyy K, for elev. yy to yy, except
      ant. yy: yyy-yyy K

11c. VPLOT/weights
  docal 2 ; gainuse 3 ; flagver 1 ; nplot 9 ; bparm 0 16 0 ; doarray 1
  bif 1 ; eif 2 ; avgif -1 ; do3col 1 ; stokes 'LL' ; ante 20 0
    --> weights  70-300 except
      EVLA ant. 14:  60-140 --> lower than VLA -- all IFs comparable
      EVLA ant. 16:  50-140 --> lower than VLA -- all IFs comparable

12. Pol'n calibration -> postponed until PCAL works and/or we can transfer

13. TASAV -> KTASAV.1

14. SPLIT -> * 811.KSPLIT
  docal 2 ; gainuse 3 ; DOPOL -1 ; flagver 1 ; outcla 'Ksplit' ; outseq 0
  sour '' ; calcode '-cal' ; DOUVCOMP -1

15. IMAGR
  a. GRS1915  811.KSPLIT -> KBM001, KCL001  5.3degs.
    (-1) UVFLG EVLA ants. -> FG/1
      flagver 1 ; reason 'evla ant'; opcode 'flag' ; dohist 1
      ante  13, 14, 16, 18,0

    (0) Outseq XXX: The Big Picture -- skipped, based on B config.

    (1) Outseq 1: NA, boxed CLEAN 
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn 'na' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,5,2052,2044
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  34508 vis.
       110.2 x  86.8 mas at  -3.9 deg
        WTNOISE 1.000
        CLEANed   9.004 mJy
        -0.627 to  8.649 mJy/beam (40,40-4050,4050: -0.593, 8.649 mJy/beam;
          peak at 2051,2043)
        rms 119.4 microJy/beam, from histogram

        OK detection
        TVMAX pixr -0.60e-3,0.60e-3
          AIPS 1: Map pixel position:  2051.41 2043.25
          AIPS 1: Skypos: RA 19 15 11.5473821  DEC 10 56 44.709513
          AIPS 1: Skypos: IPOL  22460.100 MHZ
          AIPS 1: PEAK =    8.704 Milli JY/BEAM
        JMFIT: blc 2016,2006 ; trc 2087,2080
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2016,2006 ; trc 2087,2080
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2016,2006 ; trc 2087,2080

        --> 8.6487 +/- 0.1194 mJy/bm; int. mJy; TVSTAT int. 9.2946 mJy
            [pk @ 2051,2043]                                   [622 pix.]

          Corr. for JMFIT zero: +/- 0.1194 mJy/bm; int. mJy;
            TVSTAT int. mJy

    (2) Outseq 2: Robust 0.5
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust 0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,5,2051,2043
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  34508 vis.
        92.1 x  75.3 mas at  -1.3 deg
        WTNOISE 1.077
        CLEANed   9.001 mJy
        -0.715 to  8.554 mJy/beam (40,40-4050,4050: -0.663, 8.554 mJy/beam;
          peak at 2051,2043)
        rms 124.2 microJy/beam, from histogram

        OK detection
        TVMAX pixr -0.58e-3,0.58e-3
        JMFIT: blc 2016,2006 ; trc 2087,2080
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2016,2006 ; trc 2087,2080
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2016,2006 ; trc 2087,2080

        --> 8.5541 +/- 0.1242 mJy/bm; int. mJy; TVSTAT int. 9.1122 mJy
            [pk @ 2051,2043]                                   [464 pix.]

          Corr. for JMFIT zero: +/- 0.1242 mJy/bm; int. mJy;
            TVSTAT int. mJy

  b. CYGX-3   811.KSPLIT -> KBM001, KCL001  4.7degs.
    (-1) UVFLG EVLA ants. -> FG/1
      flagver 1 ; reason 'evla ant'; opcode 'flag' ; dohist 1
      ante  13, 14, 16, 18,0

    (0) Outseq XXX: The Big Picture -- skipped, based on B config.

    (1) Outseq 1: Robust -0.5, boxed CLEAN 
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,4,2070,2041
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  17304 vis.
        84.0 x  67.0 mas at   24.6 deg
        WTNOISE 1.328
        CLEANed 671.572 mJy
        -37.901 to 605.266 mJy/beam (40,40-4050,4050: -37.901,605.27 mJy/beam;
          peak at 2070,2040)
        rms 4560.1 microJy/beam, from histogram

        Fine detection, needs self-cal
        TVMAX pixr -38e-3,38e-3
          AIPS 1: Map pixel position:  2069.75 2040.39
          AIPS 1: Skypos: RA 20 32 25.7716990  DEC 40 57 27.930865
          AIPS 1: Skypos: IPOL  22460.100 MHZ
          AIPS 1: PEAK =  608.989 Milli JY/BEAM
        JMFIT: blc 2035,2003 ; trc 2105,2078
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2035,2003 ; trc 2105,2078
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2035,2003 ; trc 2105,2078

        --> 605.266 +/- 4.560 mJy/bm; int. mJy; TVSTAT int. 672.94 mJy
            [pk @ 2070,2040]                                   [425 pix.]

          Corr. for JMFIT zero: +/- 4.560 mJy/bm; int. mJy;
            TVSTAT int. mJy

    (1b) CALIB phase-only, dot, SOLINT 1/3 -- > KPCAL .1
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver 1
      clr2na ; nmap 0 ; ncomp 0 ; flux -1e6 ; cmethod 'dft'
      outna ''; outseq 0
      refant 10 $ E8
      aparm 0 ; aparm(6) 2 ; minamper 0 ; minphser 100 ; cparm 0 0 0 100 0
      snver 0
      smodel 0.7,0
      uvra 0 ; solint 1/3 ; soltype 'l1' ; solmode 'p' ; outcla 'Kpcal' 
        solcon 0
      aparm(3) 0 ; aparm(5) 0 $ DO NOT average IFs/pol'ns
      -->  1196 good sol'ns
        CALIB1: Average closure rms =  0.0282 +-  0.0042

    (11) KPCAL.1 -> Outseq 11: Robust -0.5, uvra 0,  boxed CLEAN 
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,4,2048,2049
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  17304 vis.
        84.0 x  67.0 mas at   24.6 deg
        WTNOISE 1.328
        CLEANed 726.982 mJy
         -7.162 to 650.445 mJy/beam (40,40-4050,4050: -7.172, 650.44 mJy/beam;
          peak at 2048,2049)
        rms  776.5 microJy/beam, from histogram

        Much better.  Still hints of N/S extension -- must track this down!!
        TVMAX pixr -6.9e-3,6.9e-3
          AIPS 1: Map pixel position:  2048.00 2049.01
          AIPS 1: Skypos: RA 20 32 25.7908996  DEC 40 57 28.017112
          AIPS 1: Skypos: IPOL  22460.100 MHZ
          AIPS 1: PEAK =  649.851 Milli JY/BEAM
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2012,2011 ; trc 2083,2088

        --> 650.445 +/- 0.7765 mJy/bm; int. mJy; TVSTAT int. 725.69 mJy
            [pk @ 2048,2049]                                   [642 pix.]

          Corr. for JMFIT zero: +/- 0.7765 mJy/bm; int. mJy;
            TVSTAT int. mJy

    (11b) KPCAL.1: CALIB A&P, dot, SOLINT 1 -- > KAPCAL .1
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver 1
      clr2na ; nmap 0 ; ncomp 0 ; flux -1e6 ; cmethod 'dft'
      outna ''; outseq 0
      refant 10 $ E8
      aparm 0 ; aparm(6) 2 ; minamper 0 ; minphser 100 ; cparm 0 0 0 100 0
      snver 0
      smodel 0.725,0
      uvra 0 ; solint 1 ; soltype 'l1' ; solmode 'a&p' ; outcla 'Kapcal' 
        solcon 0
      aparm(3) 0 ; aparm(5) 0 $ DO NOT average IFs/pol'ns
      cparm(2) 1 $ normalize gains to unity to avoid flux wander
      -->   552 good sol'ns
        CALIB1: All selected antennas are used in the gain normalization
        CALIB1: RPOL, IF= 1 The average gain over these antennas is  1.048E+00
        CALIB1: RPOL, IF= 2 The average gain over these antennas is  1.048E+00
        CALIB1: LPOL, IF= 1 The average gain over these antennas is  1.047E+00
        CALIB1: LPOL, IF= 2 The average gain over these antennas is  1.048E+00
        CALIB1: The global average gain which is actually
        CALIB1: applied to all data, is   1.048E+00

        CALIB1: Average closure rms =  0.0694 +-  0.0061

    (21) KAPCAL.1 -> Outseq 21: Robust -0.5, uvra 0,  boxed CLEAN 
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,4,2048,2049
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  17304 vis.
        86.2 x  67.1 mas at   21.6 deg
        WTNOISE 1.366
        CLEANed 684.318 mJy
         -7.135 to 654.250 mJy/beam (40,40-4050,4050: -7.135, 654.25 mJy/beam;
          peak at 2048,2049)
        rms 1097.0 microJy/beam, from histogram

        Worse than phase-only.  Hmmm....
        TVMAX pixr -8.7e-3,8.7e-3
          AIPS 1: Map pixel position:  2048.00 2049.01
          AIPS 1: Skypos: RA 20 32 25.7908991  DEC 40 57 28.017083
          AIPS 1: Skypos: IPOL  22460.100 MHZ
          AIPS 1: PEAK =  653.612 Milli JY/BEAM
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2012,2011 ; trc 2083,2088

        --> 654.250 +/- 1.097  mJy/bm; int. mJy; TVSTAT int. 693.31 mJy
            [pk @ 2048,2049]                                   [706 pix.]

          Corr. for JMFIT zero: +/- 1.097  mJy/bm; int. mJy;
            TVSTAT int. mJy

    (11c) KPCAL.1: CALIB A&P, CLEAN model, SOLINT 1 -- > KAPCAL .2
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver 1
      clr2na ; nmap 0 ; ncomp 0 ; flux -1e6 ; cmethod 'dft'
      in2na inna; in2cla 'icl001'; in2se 11 ; ncomp 1e6,0; nmap 1
      outna ''; outseq 0
      refant 10 $ E8
      aparm 0 ; aparm(6) 2 ; minamper 0 ; minphser 100 ; cparm 0 0 0 100 0
      snver 0
      smodel 0
      uvra 0 ; solint 1 ; soltype 'l1' ; solmode 'a&p' ; outcla 'Kapcal' 
        solcon 0
      aparm(3) 0 ; aparm(5) 0 $ DO NOT average IFs/pol'ns
      cparm(2) 1 $ normalize gains to unity to avoid flux wander
      -->   552 good sol'ns
        CALIB1: All selected antennas are used in the gain normalization
        CALIB1: RPOL, IF= 1 The average gain over these antennas is  1.002E+00
        CALIB1: RPOL, IF= 2 The average gain over these antennas is  1.001E+00
        CALIB1: LPOL, IF= 1 The average gain over these antennas is  1.001E+00
        CALIB1: LPOL, IF= 2 The average gain over these antennas is  1.001E+00
        CALIB1: The global average gain which is actually
        CALIB1: applied to all data, is   1.001E+00

        CALIB1: Average closure rms =  0.0344 +-  0.0121

    (31) KAPCAL.2 -> Outseq 31: Robust -0.5, uvra 0,  boxed CLEAN 
      FLAGVER 1 $ no EVLA antennas
      timer 0 ; bif 0 ; eif 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.01 ; imsize 4096 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,4,2048,2049
        added clbox(1,2)~ 2045,2043,2051,2056 during CLEAN !!
      do3d 1 ; overlap 2
      niter   500
      "noise" box: blc 1230,2290 ; trc 2700,3000
      -->  17304 vis.
        84.0 x  66.9 mas at   24.7 deg
        WTNOISE 1.329
        CLEANed 727.520 mJy
         -1.799 to 647.365 mJy/beam (40,40-4050,4050: -1.799, 647.37 mJy/beam;
          peak at 2048,2049)
        rms  303.0 microJy/beam, from histogram

        MUCH better
        TVMAX pixr -1.8e-3,1.8e-3
          AIPS 1: Map pixel position:  2048.00 2049.01
          AIPS 1: Skypos: RA 20 32 25.7908991  DEC 40 57 28.017125
          AIPS 1: Skypos: IPOL  22460.100 MHZ
          AIPS 1: PEAK =  646.779 Milli JY/BEAM
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
          ...no obvious residuals
        JMFIT: blc 2012,2011 ; trc 2083,2088
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 2012,2011 ; trc 2083,2088

        --> 647.365 +/- 0.3030 mJy/bm; int. mJy; TVSTAT int. 727.97 mJy
            [pk @ 2048,2049]                                   [729 pix.]

          Corr. for JMFIT zero: +/- 0.3030 mJy/bm; int. mJy;
            TVSTAT int. mJy

================================================================================
================================================================================