AIPS Reduction of AR603: VLA Obs. of 10jul06
============================================
mpr 17jul06-

***4cm finished***

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
  10jul06 is MJD 53926.0
  B configuration
  Integration time: 3.3sec for all

  Known problems:
    ***Dynamics: AR603-A-B-P,C-B-P x2 
    * Wind/API rms phase/clouds:
      10:58  1.6 m/s    2.8d    20% mixed clouds
      12:28  2.3 m/s   12.3d    sky clear
      13:59  0.6 m/s    2.3d    sky clear
    * ant. 13,14,16,18,24 stolen for EVLA
    * ant.  6 stowed due to subreflector problems
    * ant. 13 18 moved recently -- lousy positions
    * Ptg. failures: no ref. ptg. in this run

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

----
 4cm
----
*Analyzed 16aug06*

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

  The Sources
  -----------
  0137+331  flux cal.  3C48 
    elev yy
     4cm  uvra   0,  40   B: inner 6 ant./arm
     4cm  11:49:20-11:50:50 --> 0.49311 days  1.6mins

  0319+415  pol'n cal. aka 3C84
    elev yy
     4cm  uvra  0,  0 ; primary ; 20.70 Jy
     4cm  11:53:03-14:27:53 --> 0.55295 days  4.5mins

  02450+62281 aka 0244+624 phase cal. for 02405+61138 aka LS61303  1.34 degs.
    elev yy
     4cm  uvra  0,  0 ; primary ; 1.88 Jy
     4cm  13:02:10-14:14:50 --> 0.56840 days  8.9mins
  N0239+60 phase cal. for 02405+61138 aka LS61303 ???? degs.
    elev yy
     4cm  uvra  0,  0 ; ??? ; ???? Jy
     4cm  13:15:50-14:17:00 --> 0.57390 days  2.8mins
  02405+61138 aka LS61303  1.34 degs.
    elev yy
     4cm  12:59:50-14:15:20 --> 0.56830 days 20.3mins

  20077+40298 aka 2007+404  phase cal. for 20324+40575 aka CygX-3 4.7 degs.
    elev yy
     4cm  uvra  0,  0 ; primary ; 2.85 Jy
     4cm  11:03:50-11:44:10 --> 0.47371 days  4.8mins
  20324+40575 aka CygX-3  4.7 degs.
    elev yy
     4cm  11:01:20-11:45:40 --> 0.47524 days 16.2mins

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
    : ant. 12: high ph.diff. in latter half of one scan
    --> phase diff < 24.5 degrees

3. Elevation dependences: already done in FILLM!

4. SETJY 0137+331
  sour '0137+331','' ; opty 'calc'; bif 0 ; eif 0 ; freqid 1
  zerosp 0 ; calcode '' ; aparm 0
  SETJY1: BIF =  1 EIF =  2 /Range of IFs
  SETJY1: '0137+331        ' IF =  1 FLUX = 3.1637 (Jy calcd)
  SETJY1: '0137+331        ' IF =  2 FLUX = 3.1449 (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 '3c48' ; band 'X'

5b. CALIB phases: CL/1, FG/1 -> SN/1
  USE DOCAL=1, per Juan
  a. calsour '0137+331','' --> SN/1
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante  0
    flagver 1
    solint  0.1667 ; 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 '3c48_X' ; 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
    -->   672 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0059 +-  0.0011

  b. calsour '0319+415','' $ --> SN/2
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  $ [avoids closure errors]
    ante 0
    flagver 1
    solint  0.1667 ; 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
    -->  1756 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0035 +-  0.0008

  c. calsour '02450+62281','20077+40298','' --> SN/3
    docal 1 ; gainuse 1
    uvra   0,   0
    wtuv 0.0  [avoids closure errors]
    ante 0
    flagver 1
    solint  0.1667 ; 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
    -->  7656 good sol'ns
      No closure errors
      Avg. closures on 0244+624 up to 1.6d
      CALIB1: Average closure rms =  0.0172 +-  0.0084

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 '0137+331','0319+415','02450+62281','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 '02450+62281',''
    calcode ' '
    sour 'N0239+60','02405+61138',''
    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 '0137+331','' --> 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

    in2na '3c48_X' ; 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
    -->    84 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0009 +-  0.0001

  b. calsour '0319+415','' $ --> 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 now OK
    -->   252 good sol'ns
      No closure errors
      CALIB1: Average closure rms =  0.0105 +-  0.0033

  c. calsour '02450+62281','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
    -->  1428 good sol'ns
      No closure errors
      Avg. closures on 0244+624 up to 1.2%/0.5d
      CALIB1: Average closure rms =  0.0014 +-  0.0004

8. GETJY SN/2
  calsour '0137+331',''
  sour '0319+415','02450+62281','20077+40298',''
  freqid 1
  snver 2
  --> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
      GETJY1: 20077+40298     :  0   A     1    2.27146 +/-   0.00567
      GETJY1:                              2    2.26860 +/-   0.00565
      GETJY1: 0319+415        :  0   B     1   15.86484 +/-   0.03410
      GETJY1:                              2   15.86946 +/-   0.03223
      GETJY1: 02450+62281     :  0   B     1    0.69100 +/-   0.00180
      GETJY1:                              2    0.69296 +/-   0.00182
      GETJY1: Appears to have ended successfully
      GETJY1: scamper      31DEC06 TST: Cpu=       0.0  Real=       0

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
  --> ant. yyy seem to be mis-pointing

9. CLCAL A&P: CL/2, SN/2 -> CL/3
  a. Gains for all calibrators 
    snver 2 ; gainver 2 ; gainuse 3
    calsour '0137+331','0319+415','02450+62281','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 '02450+62281',''
    calcode ' '
    sour 'N0239+60','02405+61138',''
    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. < 207 mJy
    --> phase < 12.5 degs.

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

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. ANBPL/weights
  docal 2 ; gainuse 3 ; flagver 1 ; nplot 9 ; bparm 0 17 0 ; doarray 1
  bif 1 ; eif 2 ; avgif -1 ; do3col 1 ; stokes '' ; ante 0; opcode 'alsi'
    --> weights 450-1400 except
      ant. 19: 350- 575
      ant. 28: 350- 850 -- also starts low, vs. others start high --
        cooling?

12. TVFLG RLLR CL/3, FG/1 -> FG/1
  docal 2 ; gainuse 3 ; flagver 1 ; sour '' ; docat -1
  dparm 0 ; dparm(6) 3.3333333 ; dparm(3) 1; STOKES 'RLLR'
  a. calcode '*'
    --> amp. diff. < 109 mJy

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

12b. TVFLG V CL/3, FG/1 -> FG/1
  docal 2 ; gainuse 3 ; flagver 1 ; sour '' ; docat -1
  dparm 0 ; dparm(6) 3.3333333 ; dparm(3) 1; STOKES 'V'
  a. calcode '*'
    --> amp. diff. < 835 mJy

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

13. TASAV -> XTASAV.1

14. Pol'n calibration using unpolarized source, per Greg Taylor
  a. PCAL
    N.B. pmodel=0 gives `Solution indeterminate'

    calsour '0319+415',''
    docal 2 ; gainuse 3 ; flagver 1
    pmodel 15.867,0,0,0
    solint 1/6 ; soltype 'appr'
    prtlev 1
    refant 10 $ (E8)
    -->
      IF 1:
        Largest RCP amp.: ant.  7    0.03399 +/- 0.00016
          No solutions in RCP for ant.:
             6 10 13 14 16 18 24 26 29
        Largest LCP amp.: ant. 22    0.02949 +/- 0.00017
          No solutions in LCP for ant.:
             6 13 14 16 18 24 26 29
      IF  2:
        Largest RCP amp.: ant.  7    0.04174 +/- 0.00018
          No solutions in RCP for ant.:
             6 10 13 14 16 18 24 26 29
        Largest LCP amp.: ant. 22    0.02647 +/- 0.00018
          No solutions in LCP for ant.:
             6 13 14 16 18 24 26 29

  b. RLDIF
    sour '0137+331',''
    docal 2 ; gainuse 3 ; flagver 1
    dopol 1
    uvra   0,  40
    wtuv 0.0 $ [avoids closure errors]
    ante  00
    -->
      IF 1:
        Ampscalar average of matrix  =    93.06(  0.123) sigma =  0.945
        Vector average of upper data =    96.73          sigma = 12.204
        Vector average of lower data =    90.11          sigma = 11.938
      IF 2:
        Ampscalar average of matrix  =   100.37(  0.236) sigma =  1.811
        Vector average of upper data =   104.07          sigma = 11.422
        Vector average of lower data =    97.29          sigma = 11.835
      outputs:
        AIPS 1: RLDIF:  Task to return Right - Left phase difference
        AIPS 1: Adverbs     Values                 Comments
        AIPS 1: ----------------------------------------------------------------
        AIPS 1: CLCORPRM     93.06      100.3729   R-L phase difference returned
        AIPS 1:            *rest 0

  c. CLCOR -> CL/4
    for i=1:2 ; clcorp(i)=-138-clcorp(i);end
    sour ''; ante 0 ; timer 0 ; bif 1 ; eif 2 ; freqid 1
    gainver 3 ; gainuse 0
    opcode 'polr'
    stokes 'l'

  d. RLDIF as above, with new GAINUSE --> -138.00 degs., as expected!


15. TASAV -> XTASAV.11

16. SPLIT -> * 926.XSPLIT
  docal 2 ; gainuse 4 ; dopol 1 ; flagver 1 ; outcla 'Xsplit' ; outseq 0
  sour '' ; calcode '-cal' ; DOUVCOMP -1

17. IMAGR
  a. N0239+60 926.XSPLIT.1 -> XBM001, XCL001  ref'd to 0244+644
    (1) Outseq 1: Robust -0.5, boxed CLEAN 
      timer 0 ; bif 0 ; eif 0 ; stokes '' ; ante 0 ; base 0
      flagver -1 
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.04 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,5,512,513
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 100,100 ; trc  900, 350
      -->  10080 vis.
        776 x 641 mas at  15.1 deg
        WTNOISE 1.205
        CLEANed 129.144 mJy 
        -4.928 to 128.3060 mJy/beam (10,10-1010,1010: -4.928, 128.31 mJy/bm;
          pk at  551, 513)
        rms  924.2 microJy/beam, from histogram

        Fine detection, obvious phase errors
        TVMAX pixr -5.4e-3,5.4e-3
          AIPS 1: Map pixel position:   511.45  513.30
          AIPS 1: Skypos: RA 02 39 37.95524  DEC 60 05 41.9231
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  128.683 Milli JY/BEAM
        JMFIT: blc 430,432 ; trc 587,596
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 430,432 ; trc 587,596
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 430,432 ; trc 587,596

        --> 128.3060 +/- 0.924 mJy; int. mJy; TVSTAT int. 131.47 mJy
            [pk @ 511,513]                                  [2006 pix.]

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

    (1b) Outse 1: Robust 0.5, Stokes QUV
      -->  10080 vis.
        807 x 674 mas at  11.6 deg
        WTNOISE 1.044
       Q: NADA         rms  119.8 microJy/beam
       U: -0.713 mJy   rms  122.7 microJy/beam
         clbox -1,6,537,479
       V: -0.073 mJy   rms  122.7 microJy/beam
         clbox -1,2,514,504, -1,2,512,522 (one pos., one neg.)

    (1c) CALIB phase-only, SOLINT 1/6 -- > XPCAL .1
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver 2
      clr2na ; ncomp 0; nmap 0 ; 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.13,0,0
      uvra 0 ; solint 1/6 ; soltype 'l1' ; solmode 'p' ; outcla 'xpcal' 
        solcon 0
      -->  1344 good sol'ns
        CALIB1: Average closure rms =  0.0983 +-  0.0086

    (11) XPCAL.1 -> Outseq 11: Robust -0.5, boxed CLEAN 
      timer 0 ; bif 0 ; eif 0 ; stokes '' ; ante 0 ; base 0
      flagver -1 
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.04 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 10 ; clbox 0
        clbox(1, 1)~  -1      6    513    513
        clbox(1, 2)~ 529    467    543    489
        clbox(1, 3)~ 541    441    561    475
        clbox(1, 4)~ 549    435    571    459
        clbox(1, 5)~ 569    387    593    423
        clbox(1, 6)~ 519    483    539    509
        clbox(1, 7)~ 545    389    573    429
        clbox(1, 8)~ 587    323    615    357
        clbox(1, 9)~ 525    117    689    293
        clbox(1,10)~ 467    535    507    589
      do3d 1 ; overlap 2
      niter  2000
      "noise" box: blc 100,100 ; trc  350, 900
      -->  10080 (lost 0) vis.
        776 x 641 mas at  15.1 deg
        WTNOISE 1.205
        CLEANed 147.793 mJy 
        -0.421 to 132.0744 mJy/beam (10,10-1010,1010: -0.421, 132.07 mJy/bm;
          pk at  512, 513)
        rms  100.5 microJy/beam, from histogram

        Much nicer
        TVMAX pixr -0.35e-3,0.35e-3
          AIPS 1: Map pixel position:   512.00  513.00
          AIPS 1: Skypos: RA 02 39 37.95228  DEC 60 05 41.9110
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  132.068 Milli JY/BEAM

          AIPS 1: Map pixel position:   550.01  450.48
          AIPS 1: Skypos: RA 02 39 37.74898  DEC 60 05 39.4102
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =    1.319 Milli JY/BEAM
          
          AIPS 1: Map pixel position:   576.12  397.32
          AIPS 1: Skypos: RA 02 39 37.60935  DEC 60 05 37.2838
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =    1.322 Milli JY/BEAM
          
          AIPS 1: Map pixel position:   598.69  338.63
          AIPS 1: Skypos: RA 02 39 37.48865  DEC 60 05 34.9362
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  643.743 Micro JY/BEAM
          
          AIPS 1: Map pixel position:   487.54  552.13
          AIPS 1: Skypos: RA 02 39 38.08314  DEC 60 05 43.4763
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =  565.200 Micro JY/BEAM
        JMFIT: blc 430,432 ; trc 587,596
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 430,432 ; trc 587,596
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 430,432 ; trc 587,596

        --> 132.0744 +/- 0.1005 mJy; int. mJy; TVSTAT int. 132.58 mJy
            [pk @ 512,513]                                  [3098 pix.]

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

  b. LSI61303 926.XSPLIT.1 -> XBM001, XCL001  1.3 degs. 
    (0) Outseq XXX: The Big Picture -- skipped based on 6cm
      
    (1) Outseq 1: Robust -0.5, boxed CLEAN
      FLAGVER -1
      timer 0 ; bif 0 ; eif 0 ; stokes '' ; ante 0 ; base 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.04 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,5,512,488
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 100,100; trc  300, 900
      -->  73060 vis.
        754 x 567 mas at  15.7 deg
        WTNOISE 1.353
        CLEANed  66.160 mJy 
        -0.897 to 66.0387 mJy/beam (10,10-1010,1010: -0.896, 66.039 mJy/beam;
          peak at 512,488)
        rms 166.1 microJy/beam, from histogram

        Fine detection; some phase errors
        TVMAX: pixr -0.83e-3,0.83e-3
          AIPS 1: Map pixel position:   511.88  488.20
          AIPS 1: Skypos: RA 02 40 31.66457  DEC 61 13 45.5999
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =   66.063 Milli JY/BEAM
        JMFIT: blc 444,423 ; trc 580,555
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 444,423 ; trc 580,555
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 444,423 ; trc 580,555

        --> 66.0387 +/- 0.1661 mJy/bm; int. mJy; TVSTAT int. 65.830 mJy
            [pk @ 512,488]                                  [2019 pix.]

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

    (1b) Outse 1: Robust 0.5, Stokes QUV
      -->  73040 vis.
        790 x 628 mas at  14.5 deg
        WTNOISE 1.075
       Q:  nada       rms  47.5 microJy/beam
       U:  nada       rms  49.8 microJy/beam
       V:  nada       rms  54.7 microJy/beam

    (1c) CALIB phase-only, SOLINT 1/6 -- > XPCAL .1
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver 2
      clr2na ; ncomp 0; nmap 0 ; outna ''; outseq 0
      refant 10 $ (E8)
      aparm 0 ; aparm(6) 2 ; minamper 0 ; minphser 0 ; cparm 0 0 0 0 0
      snver 0
      smodel 0.066,0,0
      uvra 0 ; solint 1/6 ; soltype 'l1' ; solmode 'p' ; outcla 'xpcal' 
        solcon 0
      -->  9912 good sol'ns
        CALIB1: Average closure rms =  0.2032 +-  0.0246

    (11) XPCAL.1 -> Outseq 11: Robust -0.5, boxed CLEAN 
      FLAGVER -1
      timer 0 ; bif 0 ; eif 0 ; stokes '' ; ante 0 ; base 0
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255
      cells 0.04 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 1 ; clbox -1,5,512,513
      do3d 1 ; overlap 2
      niter   200
      "noise" box: blc 100,100; trc  300, 900
      -->  73060 (lost 0) vis.
        754 x 567 mas at  15.7 deg
        WTNOISE 1.353
        CLEANed  66.937 mJy 
        -0.391 to 66.9645 mJy/beam (10,10-1010,1010: -0.391, 66.964 mJy/beam;
          peak at 512,513)
        rms  67.7 microJy/beam, from histogram

        Much better; a hint of variability
        TVMAX: pixr -0.39e-3,0.39e-3
          AIPS 1: Map pixel position:   511.99  512.99
          AIPS 1: Skypos: RA 02 40 31.66394  DEC 61 13 46.5918
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =   66.960 Milli JY/BEAM
        JMFIT: blc 433,430 ; trc 591,600
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 433,430 ; trc 591,600
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 433,430 ; trc 591,600

        --> 66.9645 +/- 0.0677 mJy/bm; int. mJy; TVSTAT int. 65.830 mJy
            [pk @ 512,513]                                  [2019 pix.]

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

     (11b) UVPLT 
       bparm 11  9 0 ; bparm(8) 100
         --> timer 0 13 00 0, 0 13 20 0 is 68.3 mJy
             timer 0 13 58 0, 0 14 15 0 is 65.5 mJy
       bparm 11 10 0 ; bparm(8) 100
         --> steady at +/- 0.3 mJy

  c. CYGX-3   926.XSPLIT.1 -> XBM001, XCL001  4.7 degs. 
    (0) Outseq XXX: The Big Picture -- skipped based on earlier stuff
      
    (1) Outseq 1: Robust -0.5, careful boxed CLEAN
      FLAGVER -1
      timer 0 ; bif 0 ; eif 0 ; ante 0 ; base 0; stokes ''
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255 $ REVISED BELOW
      cells 0.02 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 3 ; clbox 0
        clbox(1,1)~  -1      5    523    512
        clbox(1,2)~ 508    482    539    544
        clbox(1,3)~ 507    466    533    489
      do3d 1 ; overlap 2
      niter   600
      $ Try for more careful cleaning:
      gain 0.02 ; minpa 513; factor -0.3 ; cmethod 'dft'
      "noise" box: blc 100,100 ; trc  400,900
      -->  55540 vis.
        757 x 594 mas at -69.9 deg
        WTNOISE 1.329
        CLEANed  11.415 mJy
        -0.333 to   7.5822 mJy/beam (10,10-1010,1010: -0.258, 7.5822 mJy/beam;
          peak at 523,511)
        rms  69.8 microJy/beam, from histogram

        TVMAX pixr -0.28e-3,0.28e-3
          AIPS 1: Map pixel position:   521.34  496.54
          AIPS 1: Skypos: RA 20 32 25.77441  DEC 40 57 27.6879
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =    6.023 Milli JY/BEAM
        JMFIT: blc 367,348 ; trc 683,682
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 367,348 ; trc 683,682
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 367,348 ; trc 683,682

        -->  7.5822 +/- 0.0698 mJy/bm; int. mJy; TVSTAT int. 11.652 mJy
            [pk @ 523,511]                                  [10552 pix.]

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

    (1b) CCRES ICL001.1 --> IRS001.1
      bmaj 0.05 ; bmin 0.05; bpa 0
      TVMAX
        AIPS 1: Map pixel position:   523.00  517.00
        AIPS 1: Skypos: RA 20 32 25.77148  DEC 40 57 28.0969
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    1.052 Milli JY/BEAM
        ...N1
        
        AIPS 1: Map pixel position:   523.89  534.15
        AIPS 1: Skypos: RA 20 32 25.76991  DEC 40 57 28.4400
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    1.865 Milli JY/BEAM
        ...N2
        
        AIPS 1: Map pixel position:   522.96  511.72
        AIPS 1: Skypos: RA 20 32 25.77155  DEC 40 57 27.9915
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    3.482 Milli JY/BEAM
        ...peak -- center?
        
        AIPS 1: Map pixel position:   522.91  510.85
        AIPS 1: Skypos: RA 20 32 25.77164  DEC 40 57 27.9740
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    2.789 Milli JY/BEAM
        ...S1

        AIPS 1: Map pixel position:   522.65  493.01
        AIPS 1: Skypos: RA 20 32 25.77209  DEC 40 57 27.6172
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    1.472 Milli JY/BEAM
        ...S2
        
    (1c) Outse 1: Robust 0.5, Stokes QUV
      -->  55540 vis.
        832 x 653 mas at -77.2 deg
        WTNOISE 1.086
       Q:  -0.721 mJy  rms  47.7 microJy/beam
         niter 200; clbox  -1,5,526,507
       U:  nada        rms  47.6 microJy/beam
       V:  nada        rms  46.1 microJy/beam

    (1d) CALIB phase-only, SOLINT 2 -- > XPCAL .1
      calsour '' ; wtuv 0.01 ; docal -1 ; flagver -1
      clr2na 
        in2na inna ; in2cla 'xcl001'; in2se 1; in2di indi
        ncomp 1e6,0; nmap 1 ; inver 0
      outna ''; outseq 0
      refant 10 $ (E8)
      aparm 0 ; aparm(6) 2 ; minamper 0 ; minphser 0 ; cparm 0 0 0 0 0
      snver 0
      smodel 0
      uvra 0 ; solint 2 ; soltype 'l1' ; solmode 'p' ; outcla 'xpcal' 
        solcon 0
      aparm(3) 1 $ avg. CPs
      aparm(5) 0 $ DON'T avg. CPs
      aparm(7) 4 $ SNR >= 4
      cparm(2) 1 $ normalize gains
      -->   982 good sol'ns
             26 failed
        CALIB1: Average closure rms =  0.4746 +-  0.1413

    (11) XPCAL.1 -> Outseq 11: Robust -0.5, boxed CLEAN 
      FLAGVER -1
      timer 0 ; bif 0 ; eif 0 ; ante 0 ; base 0; stokes ''
      uvwtfn '' ; uvsize 0 ; robust -0.5 ; uvbox 0
      uvra 0, 0 ; uvtap 0
      minpa 255 $ REVISED BELOW
      cells 0.02 ; imsize 1024 ; fldsize 0
      nfield 1
        rash  0 ; decsh 0
      boxfile ''
      oboxfile 'MPR1:tmp.oboxes
      nbox 3 ; clbox 0
        clbox(1,1)~  -1      5    523    512
        clbox(1,2)~ 505    484    543    545
        clbox(1,3)~ 509    463    535    491
      do3d 1 ; overlap 2
      niter  1000
      $ Try for more careful cleaning:
      gain 0.02 ; minpa 513; factor -0.3 ; cmethod 'dft'
      "noise" box: blc 100,100 ; trc  400,900
      -->  55540 (lost 0) vis.
        757 x 594 mas at -69.9 deg
        WTNOISE 1.329
        CLEANed  11.927 mJy
        -0.274 to   7.7661 mJy/beam (10,10-1010,1010: -0.206, 7.7661 mJy/beam;
          peak at 523,511)
        rms  63.8 microJy/beam, from histogram

        Much better

        TVMAX pixr -0.20e-3,0.20e-3
          AIPS 1: Map pixel position:   522.79  511.24
          AIPS 1: Skypos: RA 20 32 25.77185  DEC 40 57 27.9818
          AIPS 1: Skypos: IPOL  8460.100 MHZ
          AIPS 1: PEAK =    7.768 Milli JY/BEAM
        JMFIT: blc 367,348 ; trc 683,682
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid -1,-1,-1,  1, 1, 1
          niter 2000
        JMFIT: blc 367,348 ; trc 683,682
          ngauss 2 ; ctype 0 3 0
          gmax 0 ; gpos 0 ; gwid 0
          dowid  1, 1, 1,  1, 1, 1
          niter 2000
        IMEAN: blc 367,348 ; trc 683,682

        -->  7.7661 +/- 0.0638 mJy/bm; int. mJy; TVSTAT int. 11.933 mJy
            [pk @ 523,511]                                  [11452 pix.]

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

    (11b) CCRES ICL001.1 --> IRS001.11
      bmaj 0.05 ; bmin 0.05; bpa 0
      TVMAX
        AIPS 1: Map pixel position:   524.04  531.46
        AIPS 1: Skypos: RA 20 32 25.76964  DEC 40 57 28.3863
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    1.177 Milli JY/BEAM
        ...N2
        
        AIPS 1: Map pixel position:   522.83  511.76
        AIPS 1: Skypos: RA 20 32 25.77177  DEC 40 57 27.9921
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =    3.438 Milli JY/BEAM
        ...peak -- center???
        
        AIPS 1: Map pixel position:   522.78  498.05
        AIPS 1: Skypos: RA 20 32 25.77187  DEC 40 57 27.7179
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =  817.398 Micro JY/BEAM
        ...S2
        
        AIPS 1: Map pixel position:   523.36  476.96
        AIPS 1: Skypos: RA 20 32 25.77085  DEC 40 57 27.2962
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =  395.690 Micro JY/BEAM
        ...S3
        
        AIPS 1: Map pixel position:   522.29  509.41
        AIPS 1: Skypos: RA 20 32 25.77273  DEC 40 57 27.9451
        AIPS 1: Skypos: IPOL  8460.100 MHZ
        AIPS 1: PEAK =  234.099 Micro JY/BEAM
        ...S1

------------------
================================================================================
================================================================================