! Started logfile: difmap.log on Thu Nov  8 11:25:04 2007
observe CTA26_S.edt
! Reading UV FITS file: CTA26_S.edt
! AN table 1: 17 integrations on 105 of 105 possible baselines.
! AN table 2: 88 integrations on 105 of 105 possible baselines.
! AN table 3: 11 integrations on 105 of 105 possible baselines.
! AN table 4: 75 integrations on 105 of 105 possible baselines.
! Apparent sampling: 0.13438 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: CTA26
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  2.23697e+09        8e+06          1        8e+06
!  02        2  2.26697e+09        8e+06          1        8e+06
!  03        3  2.35697e+09        8e+06          1        8e+06
!  04        4  2.37697e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 970 lines of history.
! 
! Reading 10780 visibilities.
select
! Selecting polarization: RR,  channels: 1..4
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
uvweight 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
mapsize 1024,0.4
! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize.
peakwin 0.6
! Inverting map and beam 
! Estimated beam: bmin=1.84 mas, bmaj=3.271 mas, bpa=13.09 degrees
! Estimated noise=1393.03 mJy/beam.
! Added new window around map position (0, 0).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 1.60542 Jy
! Component: 100  -  total flux cleaned = 1.85037 Jy
! Total flux subtracted in 100 components = 1.85037 Jy
! Clean residual min=-0.145454 max=0.499473 Jy/beam
! Clean residual mean=0.000011 rms=0.036952 Jy/beam
! Combined flux in latest and established models = 1.85037 Jy
selfcal
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 9 components and 1.85037 Jy
! 
! Correcting IF 1.
!  A total of 120 telescope corrections were flagged in sub-array 1.
!  A total of 998 telescope corrections were flagged in sub-array 2.
!  A total of 110 telescope corrections were flagged in sub-array 3.
!  A total of 859 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 2.
!  A total of 120 telescope corrections were flagged in sub-array 1.
!  A total of 998 telescope corrections were flagged in sub-array 2.
!  A total of 110 telescope corrections were flagged in sub-array 3.
!  A total of 805 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
!  A total of 120 telescope corrections were flagged in sub-array 1.
!  A total of 998 telescope corrections were flagged in sub-array 2.
!  A total of 110 telescope corrections were flagged in sub-array 3.
!  A total of 856 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 4.
!  A total of 120 telescope corrections were flagged in sub-array 1.
!  A total of 998 telescope corrections were flagged in sub-array 2.
!  A total of 110 telescope corrections were flagged in sub-array 3.
!  A total of 856 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=0.763046Jy  sigma=0.006632
! Fit after  self-cal, rms=0.762162Jy  sigma=0.006657
clean
! Inverting map and beam 
! Estimated beam: bmin=1.846 mas, bmaj=3.439 mas, bpa=14.32 degrees
! Estimated noise=1399.81 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0224314 Jy
! Component: 100  -  total flux cleaned = 0.0325432 Jy
! Total flux subtracted in 100 components = 0.0325432 Jy
! Clean residual min=-0.139457 max=0.492375 Jy/beam
! Clean residual mean=-0.000003 rms=0.036521 Jy/beam
! Combined flux in latest and established models = 1.88291 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 1.88291 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.742614Jy  sigma=0.006490
! Fit after  self-cal, rms=0.742741Jy  sigma=0.006490
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00974552 Jy
! Component: 100  -  total flux cleaned = 0.0191284 Jy
! Total flux subtracted in 100 components = 0.0191284 Jy
! Clean residual min=-0.135748 max=0.483214 Jy/beam
! Clean residual mean=-0.000012 rms=0.035813 Jy/beam
! Combined flux in latest and established models = 1.90204 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 1.90204 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.727610Jy  sigma=0.006367
! Fit after  self-cal, rms=0.727807Jy  sigma=0.006367
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00912639 Jy
! Component: 100  -  total flux cleaned = 0.0181819 Jy
! Total flux subtracted in 100 components = 0.0181819 Jy
! Clean residual min=-0.132547 max=0.474960 Jy/beam
! Clean residual mean=-0.000021 rms=0.035163 Jy/beam
! Combined flux in latest and established models = 1.92022 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 1.92022 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.715104Jy  sigma=0.006254
! Fit after  self-cal, rms=0.715287Jy  sigma=0.006254
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00850936 Jy
! Component: 100  -  total flux cleaned = 0.0168987 Jy
! Total flux subtracted in 100 components = 0.0168987 Jy
! Clean residual min=-0.129676 max=0.467216 Jy/beam
! Clean residual mean=-0.000029 rms=0.034565 Jy/beam
! Combined flux in latest and established models = 1.93712 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 1.93712 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.704685Jy  sigma=0.006150
! Fit after  self-cal, rms=0.704870Jy  sigma=0.006150
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00816087 Jy
! Component: 100  -  total flux cleaned = 0.0157473 Jy
! Total flux subtracted in 100 components = 0.0157473 Jy
! Clean residual min=-0.127700 max=0.459847 Jy/beam
! Clean residual mean=-0.000036 rms=0.034015 Jy/beam
! Combined flux in latest and established models = 1.95287 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 1.95287 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.696104Jy  sigma=0.006056
! Fit after  self-cal, rms=0.696270Jy  sigma=0.006056
mapl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 0.6
! Added new window around map position (-1.6, -0.8).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
uvweight 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean
! Inverting map and beam 
! Estimated beam: bmin=2.45 mas, bmaj=4.481 mas, bpa=13.44 degrees
! Estimated noise=789.051 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.224712 Jy
! Component: 100  -  total flux cleaned = 0.227651 Jy
! Total flux subtracted in 100 components = 0.227651 Jy
! Clean residual min=-0.116716 max=0.334363 Jy/beam
! Clean residual mean=-0.000374 rms=0.034836 Jy/beam
! Combined flux in latest and established models = 2.18052 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 9 components and 2.18052 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.678795Jy  sigma=0.005653
! Fit after  self-cal, rms=0.678960Jy  sigma=0.005652
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00589117 Jy
! Component: 100  -  total flux cleaned = 0.00865485 Jy
! Total flux subtracted in 100 components = 0.00865485 Jy
! Clean residual min=-0.115157 max=0.327958 Jy/beam
! Clean residual mean=-0.000381 rms=0.034399 Jy/beam
! Combined flux in latest and established models = 2.18918 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 9 components and 2.18918 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.673657Jy  sigma=0.005599
! Fit after  self-cal, rms=0.673777Jy  sigma=0.005598
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0030361 Jy
! Component: 100  -  total flux cleaned = 0.00584557 Jy
! Total flux subtracted in 100 components = 0.00584557 Jy
! Clean residual min=-0.113274 max=0.323287 Jy/beam
! Clean residual mean=-0.000385 rms=0.034038 Jy/beam
! Combined flux in latest and established models = 2.19502 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 9 components and 2.19502 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.669051Jy  sigma=0.005547
! Fit after  self-cal, rms=0.669150Jy  sigma=0.005547
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00292186 Jy
! Component: 100  -  total flux cleaned = 0.00543604 Jy
! Total flux subtracted in 100 components = 0.00543604 Jy
! Clean residual min=-0.111511 max=0.319079 Jy/beam
! Clean residual mean=-0.000389 rms=0.033701 Jy/beam
! Combined flux in latest and established models = 2.20046 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 9 components and 2.20046 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.665090Jy  sigma=0.005499
! Fit after  self-cal, rms=0.665159Jy  sigma=0.005498
loglevs 4
! The new contour levels are:
!  -4 4 8 16 32 64
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.45 x 4.481 at 13.44 degrees (North through East)
! Clean map  min=-0.10258  max=1.9067 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00280317 Jy
! Component: 100  -  total flux cleaned = 0.00551446 Jy
! Total flux subtracted in 100 components = 0.00551446 Jy
! Clean residual min=-0.109860 max=0.314534 Jy/beam
! Clean residual mean=-0.000394 rms=0.033373 Jy/beam
! Combined flux in latest and established models = 2.20597 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 9 components and 2.20597 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.661681Jy  sigma=0.005452
! Fit after  self-cal, rms=0.661712Jy  sigma=0.005452
peakwin 0.6
! Inverting map 
! Added new window around map position (-2, -1.2).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
mapl
! 
! Move the cursor into the plot window and press 'H' for help
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 15 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.04     HH        2.41     KB        1.00*    KP        1.00     
!   LA        0.99     MK        1.00*    NL        0.86     NY        1.00*    
!   OV        1.00*    PT        0.98     SC        1.01     TC        1.16     
!   WF        1.08     WZ        2.38     BR        1.00*    
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.97     HH        1.00*    KB        1.00*    KP        0.64     
!   LA        0.98     MK        1.49     NL        1.03     NY        1.00*    
!   OV        0.84     PT        0.96     SC        1.18     TC        0.79     
!   WF        1.00     WZ        1.00*    BR        0.89     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.02     HH        1.00*    KB        2.74     KP        0.90     
!   LA        0.93     MK        0.85     NL        1.02     NY        1.00*    
!   OV        0.97     PT        0.96     SC        1.00*    TC        1.00*    
!   WF        1.00*    WZ        1.00*    BR        0.96     
! 
!  A total of 171 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   FD        0.88     HH        2.18     KB        1.00*    KP        1.64     
!   LA        0.90     MK        1.00*    NL        0.95     NY        1.00*    
!   OV        1.00*    PT        1.17     SC        1.00*    TC        1.00*    
!   WF        1.20     WZ        3.81     BR        1.00*    
! 
! 
! Correcting IF 2.
!  A total of 15 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.02     HH        2.51     KB        1.00*    KP        0.83     
!   LA        0.79     MK        1.00*    NL        0.93     NY        1.00*    
!   OV        1.00*    PT        1.02     SC        0.89     TC        1.04     
!   WF        1.03     WZ        2.61     BR        1.00*    
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.97     HH        1.00*    KB        1.00*    KP        0.76     
!   LA        0.96     MK        1.49     NL        1.04     NY        1.00*    
!   OV        0.83     PT        0.94     SC        1.08     TC        0.77     
!   WF        1.03     WZ        1.00*    BR        0.90     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.01     HH        1.00*    KB        2.81     KP        0.66     
!   LA        0.97     MK        0.87     NL        1.05     NY        1.00*    
!   OV        0.98     PT        0.95     SC        1.00*    TC        1.00*    
!   WF        1.00*    WZ        1.00*    BR        0.96     
! 
!  A total of 60 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   FD        0.91     HH        1.65     KB        1.00*    KP        1.42     
!   LA        0.92     MK        1.00*    NL        0.97     NY        1.00*    
!   OV        1.00*    PT        1.16     SC        1.00*    TC        0.99     
!   WF        1.15     WZ        2.64     BR        1.00*    
! 
! 
! Correcting IF 3.
!  A total of 15 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.01     HH        2.73     KB        1.00*    KP        0.88     
!   LA        0.88     MK        1.00*    NL        0.94     NY        1.00*    
!   OV        1.00*    PT        1.02     SC        0.89     TC        1.09     
!   WF        1.06     WZ        3.66     BR        1.00*    
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.99     HH        1.00*    KB        1.00*    KP        0.70     
!   LA        0.97     MK        1.61     NL        0.99     NY        1.00*    
!   OV        0.85     PT        0.96     SC        1.11     TC        0.74     
!   WF        1.10     WZ        1.00*    BR        0.92     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.02     HH        1.00*    KB        3.40     KP        0.79     
!   LA        0.95     MK        0.94     NL        1.05     NY        1.00*    
!   OV        0.95     PT        0.97     SC        1.00*    TC        1.00*    
!   WF        1.00*    WZ        1.00*    BR        0.95     
! 
!  A total of 174 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   FD        0.90     HH        2.46     KB        1.00*    KP        1.51     
!   LA        0.93     MK        1.00*    NL        0.96     NY        1.00*    
!   OV        1.00*    PT        1.16     SC        1.00*    TC        1.00*    
!   WF        1.38     WZ        5.26     BR        1.00*    
! 
! 
! Correcting IF 4.
!  A total of 15 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.03     HH        2.62     KB        1.00*    KP        0.88     
!   LA        0.84     MK        1.00*    NL        0.96     NY        1.00*    
!   OV        1.00*    PT        1.03     SC        0.91     TC        1.01     
!   WF        1.08     WZ        3.57     BR        1.00*    
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.99     HH        1.00*    KB        1.00*    KP        0.73     
!   LA        0.97     MK        1.66     NL        0.99     NY        1.00*    
!   OV        0.85     PT        0.97     SC        1.09     TC        0.72     
!   WF        1.07     WZ        1.00*    BR        0.91     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.02     HH        1.00*    KB        4.53     KP        0.78     
!   LA        0.95     MK        0.96     NL        1.06     NY        1.00*    
!   OV        0.95     PT        0.96     SC        1.00*    TC        1.00*    
!   WF        1.00*    WZ        1.00*    BR        0.93     
! 
!  A total of 174 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   FD        0.91     HH        2.21     KB        1.00*    KP        1.46     
!   LA        0.92     MK        1.00*    NL        0.95     NY        1.00*    
!   OV        1.00*    PT        1.16     SC        1.00*    TC        1.00*    
!   WF        1.46     WZ        5.38     BR        1.00*    
! 
! 
! Fit before self-cal, rms=0.661712Jy  sigma=0.005452
! Fit after  self-cal, rms=1.373746Jy  sigma=0.002180
mapl
! Inverting map and beam 
! Estimated beam: bmin=4.033 mas, bmaj=5.93 mas, bpa=16.61 degrees
! Estimated noise=795.409 mJy/beam.
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -0.000898617 Jy
! Component: 100  -  total flux cleaned = -0.000914968 Jy
! Total flux subtracted in 100 components = -0.000914968 Jy
! Clean residual min=-0.017890 max=0.032124 Jy/beam
! Clean residual mean=0.000023 rms=0.005636 Jy/beam
! Combined flux in latest and established models = 2.20506 Jy
selfcal
! Performing phase self-cal
! Adding 2 model components to the UV plane model.
! The established model now contains 11 components and 2.20506 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.372418Jy  sigma=0.002178
! Fit after  self-cal, rms=1.378804Jy  sigma=0.002175
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.033 x 5.93 at 16.61 degrees (North through East)
! Clean map  min=-0.016447  max=2.0673 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
loglevs 2
! The new contour levels are:
!  -2 2 4 8 16 32 64
mapl cln
! 
! Move the cursor into the plot window and press 'H' for help
loglevs 1
! The new contour levels are:
!  -1 1 2 4 8 16 32 64
mapl cln
! 
! Move the cursor into the plot window and press 'H' for help
loglevs 0.7
! The new contour levels are:
!  -0.7 0.7 1.4 2.8 5.6 11.2 22.4 44.8 89.6
mapl cln
! 
! Move the cursor into the plot window and press 'H' for help
radp
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 365 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=4.038 mas, bmaj=5.931 mas, bpa=16.56 degrees
! Estimated noise=797.045 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.000626254 Jy
! Component: 100  -  total flux cleaned = 0.000686353 Jy
! Total flux subtracted in 100 components = 0.000686353 Jy
! Clean residual min=-0.015953 max=0.031829 Jy/beam
! Clean residual mean=0.000017 rms=0.005479 Jy/beam
! Combined flux in latest and established models = 2.20574 Jy
selfcal
! Performing phase self-cal
! Adding 3 model components to the UV plane model.
! The established model now contains 11 components and 2.20574 Jy
! 
! Correcting IF 1.
!  A total of 4 telescope corrections were flagged in sub-array 2.
!  A total of 1 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 2.
!  A total of 21 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
!  A total of 1 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 4.
!  A total of 1 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=0.486154Jy  sigma=0.001836
! Fit after  self-cal, rms=0.473588Jy  sigma=0.001831
mapl cln
! Inverting map and beam 
! Estimated beam: bmin=4.038 mas, bmaj=5.934 mas, bpa=16.58 degrees
! Estimated noise=797.087 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.038 x 5.934 at 16.58 degrees (North through East)
! Clean map  min=-0.015535  max=2.0682 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 0.6
! Inverting map 
! Added new window around map position (4.8, 4.8).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.02047 Jy
! Component: 100  -  total flux cleaned = 0.0188228 Jy
! Total flux subtracted in 100 components = 0.0188228 Jy
! Clean residual min=-0.016064 max=0.033085 Jy/beam
! Clean residual mean=0.000009 rms=0.005114 Jy/beam
! Combined flux in latest and established models = 2.22457 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 15 components and 2.22457 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.469679Jy  sigma=0.001796
! Fit after  self-cal, rms=0.469475Jy  sigma=0.001791
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0029547 Jy
! Component: 100  -  total flux cleaned = 0.00155543 Jy
! Total flux subtracted in 100 components = 0.00155543 Jy
! Clean residual min=-0.017289 max=0.031512 Jy/beam
! Clean residual mean=0.000008 rms=0.004958 Jy/beam
! Combined flux in latest and established models = 2.22612 Jy
selfcal
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 18 components and 2.22612 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.467566Jy  sigma=0.001777
! Fit after  self-cal, rms=0.467538Jy  sigma=0.001776
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0002851 Jy
! Component: 100  -  total flux cleaned = -0.000382368 Jy
! Total flux subtracted in 100 components = -0.000382368 Jy
! Clean residual min=-0.017289 max=0.030935 Jy/beam
! Clean residual mean=0.000008 rms=0.004902 Jy/beam
! Combined flux in latest and established models = 2.22574 Jy
selfcal
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 19 components and 2.22574 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.466354Jy  sigma=0.001766
! Fit after  self-cal, rms=0.466334Jy  sigma=0.001766
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.038 x 5.934 at 16.58 degrees (North through East)
! Clean map  min=-0.0161  max=2.0528 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
mapl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
selfcal true,true,120
! Performing amp+phase self-cal over 120 minute time intervals
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 2.
!  A total of 3 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 2.
!  A total of 63 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=0.466334Jy  sigma=0.001766
! Fit after  self-cal, rms=0.519589Jy  sigma=0.001679
mapl
! Inverting map and beam 
! Estimated beam: bmin=4.149 mas, bmaj=6.124 mas, bpa=16.55 degrees
! Estimated noise=799.273 mJy/beam.
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00986529 Jy
! Component: 100  -  total flux cleaned = 0.00875435 Jy
! Total flux subtracted in 100 components = 0.00875435 Jy
! Clean residual min=-0.014941 max=0.032925 Jy/beam
! Clean residual mean=0.000011 rms=0.004647 Jy/beam
! Combined flux in latest and established models = 2.23449 Jy
selfcal
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 20 components and 2.23449 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.519404Jy  sigma=0.001665
! Fit after  self-cal, rms=0.519702Jy  sigma=0.001664
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00228519 Jy
! Component: 100  -  total flux cleaned = 0.00135051 Jy
! Total flux subtracted in 100 components = 0.00135051 Jy
! Clean residual min=-0.015822 max=0.032036 Jy/beam
! Clean residual mean=0.000010 rms=0.004547 Jy/beam
! Combined flux in latest and established models = 2.23584 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 20 components and 2.23584 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.519882Jy  sigma=0.001657
! Fit after  self-cal, rms=0.519920Jy  sigma=0.001657
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.000130308 Jy
! Component: 100  -  total flux cleaned = -0.000293965 Jy
! Total flux subtracted in 100 components = -0.000293965 Jy
! Clean residual min=-0.015685 max=0.031587 Jy/beam
! Clean residual mean=0.000009 rms=0.004495 Jy/beam
! Combined flux in latest and established models = 2.23555 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 20 components and 2.23555 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.520147Jy  sigma=0.001652
! Fit after  self-cal, rms=0.520131Jy  sigma=0.001652
radp
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 2 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=4.149 mas, bmaj=6.124 mas, bpa=16.55 degrees
! Estimated noise=799.274 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -0.000388039 Jy
! Component: 100  -  total flux cleaned = -0.000784122 Jy
! Total flux subtracted in 100 components = -0.000784122 Jy
! Clean residual min=-0.015318 max=0.031318 Jy/beam
! Clean residual mean=0.000009 rms=0.004457 Jy/beam
! Combined flux in latest and established models = 2.23477 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 20 components and 2.23477 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.518714Jy  sigma=0.001647
! Fit after  self-cal, rms=0.518536Jy  sigma=0.001647
peakwin 0.6
! Inverting map 
! Added new window around map position (-3.2, 55.2).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -0.000370851 Jy
! Component: 100  -  total flux cleaned = -0.000742776 Jy
! Total flux subtracted in 100 components = -0.000742776 Jy
! Clean residual min=-0.014940 max=0.031140 Jy/beam
! Clean residual mean=0.000009 rms=0.004428 Jy/beam
! Combined flux in latest and established models = 2.23402 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 20 components and 2.23402 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.518774Jy  sigma=0.001644
! Fit after  self-cal, rms=0.518725Jy  sigma=0.001644
selfcal true,true,10
! Performing amp+phase self-cal over 10 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.518725Jy  sigma=0.001644
! Fit after  self-cal, rms=0.522928Jy  sigma=0.001604
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 19 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=4.106 mas, bmaj=6.155 mas, bpa=16.93 degrees
! Estimated noise=799.138 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00771704 Jy
! Component: 100  -  total flux cleaned = 0.00735061 Jy
! Total flux subtracted in 100 components = 0.00735061 Jy
! Clean residual min=-0.013152 max=0.032528 Jy/beam
! Clean residual mean=0.000001 rms=0.004348 Jy/beam
! Combined flux in latest and established models = 2.24137 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 21 components and 2.24137 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
!  A total of 6 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.506604Jy  sigma=0.001594
! Fit after  self-cal, rms=0.500776Jy  sigma=0.001593
clean
! Inverting map and beam 
! Estimated beam: bmin=4.106 mas, bmaj=6.155 mas, bpa=16.93 degrees
! Estimated noise=799.139 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0016289 Jy
! Component: 100  -  total flux cleaned = 0.00129229 Jy
! Total flux subtracted in 100 components = 0.00129229 Jy
! Clean residual min=-0.013193 max=0.031892 Jy/beam
! Clean residual mean=0.000000 rms=0.004292 Jy/beam
! Combined flux in latest and established models = 2.24267 Jy
selfcal
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 21 components and 2.24267 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.500843Jy  sigma=0.001590
! Fit after  self-cal, rms=0.500900Jy  sigma=0.001590
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 8.23979e-05 Jy
! Component: 100  -  total flux cleaned = -0.000195157 Jy
! Total flux subtracted in 100 components = -0.000195157 Jy
! Clean residual min=-0.013296 max=0.031588 Jy/beam
! Clean residual mean=0.000000 rms=0.004263 Jy/beam
! Combined flux in latest and established models = 2.24247 Jy
selfcal
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 21 components and 2.24247 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.501025Jy  sigma=0.001587
! Fit after  self-cal, rms=0.501039Jy  sigma=0.001587
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.106 x 6.155 at 16.93 degrees (North through East)
! Clean map  min=-0.018072  max=2.0413 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 0.6
! Inverting map 
! Added new window around map position (-3.2, 55.6).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
radp
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 254 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=4.103 mas, bmaj=6.147 mas, bpa=16.96 degrees
! Estimated noise=801.592 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -0.00231892 Jy
! Component: 100  -  total flux cleaned = -0.00260994 Jy
! Total flux subtracted in 100 components = -0.00260994 Jy
! Clean residual min=-0.012782 max=0.029369 Jy/beam
! Clean residual mean=0.000025 rms=0.004012 Jy/beam
! Combined flux in latest and established models = 2.23986 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 23 components and 2.23986 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
!  A total of 2 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=0.438967Jy  sigma=0.001536
! Fit after  self-cal, rms=0.428655Jy  sigma=0.001531
clean
! Inverting map and beam 
! Estimated beam: bmin=4.104 mas, bmaj=6.148 mas, bpa=16.93 degrees
! Estimated noise=801.613 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.000436125 Jy
! Component: 100  -  total flux cleaned = 0.000174268 Jy
! Total flux subtracted in 100 components = 0.000174268 Jy
! Clean residual min=-0.012892 max=0.029155 Jy/beam
! Clean residual mean=0.000025 rms=0.003994 Jy/beam
! Combined flux in latest and established models = 2.24004 Jy
selfcal
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 23 components and 2.24004 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.428590Jy  sigma=0.001529
! Fit after  self-cal, rms=0.428589Jy  sigma=0.001529
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.104 x 6.148 at 16.93 degrees (North through East)
! Clean map  min=-0.016611  max=2.0372 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
save CTA26_S
! Writing UV FITS file: CTA26_S.uvf
! Writing 23 model components to file: CTA26_S.mod
! wwins: Wrote 2 windows to CTA26_S.win
! Writing clean map to FITS file: CTA26_S.fits
! Writing difmap environment to: CTA26_S.par
![@gonzo_plot CTA26_S]
device CTA26_Suv.ps/vps,1,2
! Attempting to open device: 'CTA26_Suv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device CTA26_Sim.ps/vps
! Attempting to open device: 'CTA26_Sim.ps/vps'
uvweight 0,-1
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
loglevs -1,2048,2
! The new contour levels are:
!  -1 1 2 4 8 16 32 64 128 256 512 1024 2048
cmul = imstat(rms)*3
! Inverting map and beam 
! Estimated beam: bmin=2.917 mas, bmaj=4.345 mas, bpa=29.6 degrees
! Estimated noise=939.769 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.917 x 4.345 at 29.6 degrees (North through East)
! Clean map  min=-0.023505  max=1.9003 Jy/beam
![Exited script file: gonzo_plot]
quit
! Quitting program
! Log file difmap.log closed on Thu Nov  8 11:34:03 2007