! Started logfile: 3C274_X.log on Thu Aug 14 14:31:21 2008
observe 3C274_X.edt
! Reading UV FITS file: 3C274_X.edt
! AN table 1: 75 integrations on 171 of 171 possible baselines.
! AN table 2: 90 integrations on 171 of 171 possible baselines.
! Apparent sampling: 0.119121 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 3C274
! 
! 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  8.40997e+09        8e+06          1        8e+06
!  02        2  8.47997e+09        8e+06          1        8e+06
!  03        3  8.79497e+09        8e+06          1        8e+06
!  04        4  8.89997e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1526 lines of history.
! 
! Reading 13444 visibilities.
![@gonzo5-0.1 3C274_X]
print "***********************************"
! *********************************** 
print "*** Using gonzo instruction set ***"
! *** Using gonzo instruction set *** 
print "***********************************"
! *********************************** 
float field_size
field_size = 1024
float field_cell
field_cell = 0.1
integer clean_niter
clean_niter = 50
float clean_gain
clean_gain = 0.03
float dynam
dynam = 5.0
float win_mult
win_mult = 1.8
float solint
solint = 3.6/6.0
logical doflt
doflt = true
logical doamp
float new_peak
float flux_cutoff
mapsize field_size,field_cell
! Map grid = 1024x1024 pixels with 0.100x0.100 milli-arcsec cellsize.
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
startmod
! Applying default point source starting model.
! Performing phase self-cal
! Adding 1 model components to the UV plane model.
! The established model now contains 1 components and 1 Jy
! 
! Correcting IF 1.
!  A total of 973 telescope corrections were flagged in sub-array 1.
!  A total of 1436 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 976 telescope corrections were flagged in sub-array 1.
!  A total of 1436 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 976 telescope corrections were flagged in sub-array 1.
!  A total of 1448 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 973 telescope corrections were flagged in sub-array 1.
!  A total of 1442 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.586320Jy  sigma=8.038358
! Fit after  self-cal, rms=0.581928Jy  sigma=8.006648
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
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.
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.5416 mas, bmaj=1.034 mas, bpa=1.345 degrees
! Estimated noise=0.807874 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (0, 0).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.51406 Jy
! Total flux subtracted in 50 components = 0.51406 Jy
! Clean residual min=-0.037727 max=0.203176 Jy/beam
! Clean residual mean=-0.000069 rms=0.012810 Jy/beam
! Combined flux in latest and established models = 0.51406 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 8 model components to the UV plane model.
! The established model now contains 8 components and 0.51406 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.505166Jy  sigma=7.722543
! Fit after  self-cal, rms=0.504210Jy  sigma=7.717612
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.220622 Jy
! Total flux subtracted in 50 components = 0.220622 Jy
! Clean residual min=-0.027578 max=0.159176 Jy/beam
! Clean residual mean=-0.000027 rms=0.009404 Jy/beam
! Combined flux in latest and established models = 0.734681 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 14 model components to the UV plane model.
! The established model now contains 17 components and 0.734681 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.399271Jy  sigma=5.879650
! Fit after  self-cal, rms=0.398004Jy  sigma=5.870026
! Inverting map 
! Added new window around map position (-0.8, 0.2).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.151424 Jy
! Total flux subtracted in 50 components = 0.151424 Jy
! Clean residual min=-0.024835 max=0.094909 Jy/beam
! Clean residual mean=0.000004 rms=0.007291 Jy/beam
! Combined flux in latest and established models = 0.886106 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 15 model components to the UV plane model.
! The established model now contains 30 components and 0.886106 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.333843Jy  sigma=4.778609
! Fit after  self-cal, rms=0.331197Jy  sigma=4.737997
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0992841 Jy
! Total flux subtracted in 50 components = 0.0992841 Jy
! Clean residual min=-0.022454 max=0.070575 Jy/beam
! Clean residual mean=0.000025 rms=0.006072 Jy/beam
! Combined flux in latest and established models = 0.98539 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 18 model components to the UV plane model.
! The established model now contains 44 components and 0.98539 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.296726Jy  sigma=4.131732
! Fit after  self-cal, rms=0.295781Jy  sigma=4.116323
! Inverting map 
! Added new window around map position (-1.4, 0.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0761476 Jy
! Total flux subtracted in 50 components = 0.0761476 Jy
! Clean residual min=-0.020685 max=0.056744 Jy/beam
! Clean residual mean=0.000042 rms=0.005279 Jy/beam
! Combined flux in latest and established models = 1.06154 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 25 model components to the UV plane model.
! The established model now contains 59 components and 1.06154 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.272344Jy  sigma=3.702069
! Fit after  self-cal, rms=0.271934Jy  sigma=3.693521
! Inverting map 
! Added new window around map position (0.8, -0.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0627183 Jy
! Total flux subtracted in 50 components = 0.0627183 Jy
! Clean residual min=-0.018813 max=0.038506 Jy/beam
! Clean residual mean=0.000055 rms=0.004681 Jy/beam
! Combined flux in latest and established models = 1.12426 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 21 model components to the UV plane model.
! The established model now contains 73 components and 1.12426 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.254541Jy  sigma=3.384286
! Fit after  self-cal, rms=0.254310Jy  sigma=3.382125
! Inverting map 
! Added new window around map position (-2.1, 0.2).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0510059 Jy
! Total flux subtracted in 50 components = 0.0510059 Jy
! Clean residual min=-0.017674 max=0.030663 Jy/beam
! Clean residual mean=0.000065 rms=0.004261 Jy/beam
! Combined flux in latest and established models = 1.17526 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 26 model components to the UV plane model.
! The established model now contains 91 components and 1.17526 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.241676Jy  sigma=3.161283
! Fit after  self-cal, rms=0.241498Jy  sigma=3.158952
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0427029 Jy
! Total flux subtracted in 50 components = 0.0427029 Jy
! Clean residual min=-0.016871 max=0.025967 Jy/beam
! Clean residual mean=0.000074 rms=0.003941 Jy/beam
! Combined flux in latest and established models = 1.21797 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 25 model components to the UV plane model.
! The established model now contains 101 components and 1.21797 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.232018Jy  sigma=2.995866
! Fit after  self-cal, rms=0.231912Jy  sigma=2.994359
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0361699 Jy
! Total flux subtracted in 50 components = 0.0361699 Jy
! Clean residual min=-0.016097 max=0.022135 Jy/beam
! Clean residual mean=0.000081 rms=0.003698 Jy/beam
! Combined flux in latest and established models = 1.25414 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 27 model components to the UV plane model.
! The established model now contains 110 components and 1.25414 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.224634Jy  sigma=2.872524
! Fit after  self-cal, rms=0.224572Jy  sigma=2.871356
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0309927 Jy
! Total flux subtracted in 50 components = 0.0309927 Jy
! Clean residual min=-0.015455 max=0.019117 Jy/beam
! Clean residual mean=0.000087 rms=0.003512 Jy/beam
! Combined flux in latest and established models = 1.28513 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 27 model components to the UV plane model.
! The established model now contains 119 components and 1.28513 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.218894Jy  sigma=2.779227
! Fit after  self-cal, rms=0.218858Jy  sigma=2.778341
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0268218 Jy
! Total flux subtracted in 50 components = 0.0268218 Jy
! Clean residual min=-0.014770 max=0.017064 Jy/beam
! Clean residual mean=0.000092 rms=0.003368 Jy/beam
! Combined flux in latest and established models = 1.31195 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 29 model components to the UV plane model.
! The established model now contains 128 components and 1.31195 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.214344Jy  sigma=2.707913
! Fit after  self-cal, rms=0.214326Jy  sigma=2.707209
! Inverting map 
! Added new window around map position (-2.7, -0.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.023493 Jy
! Total flux subtracted in 50 components = 0.023493 Jy
! Clean residual min=-0.014205 max=0.016338 Jy/beam
! Clean residual mean=0.000097 rms=0.003253 Jy/beam
! Combined flux in latest and established models = 1.33544 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 27 model components to the UV plane model.
! The established model now contains 133 components and 1.33544 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.210736Jy  sigma=2.652920
! Fit after  self-cal, rms=0.210729Jy  sigma=2.652387
! Inverting map 
! Added new window around map position (-3, 1.6).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0210184 Jy
! Total flux subtracted in 50 components = 0.0210184 Jy
! Clean residual min=-0.013660 max=0.015425 Jy/beam
! Clean residual mean=0.000100 rms=0.003157 Jy/beam
! Combined flux in latest and established models = 1.35646 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 31 model components to the UV plane model.
! The established model now contains 144 components and 1.35646 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.207775Jy  sigma=2.608456
! Fit after  self-cal, rms=0.207753Jy  sigma=2.607995
! Inverting map 
print "****************************************"
! **************************************** 
print "*** Finished Uniform Weighting CLEAN ***"
! *** Finished Uniform Weighting CLEAN *** 
print "****************************************"
! **************************************** 
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.
win_mult = win_mult * 1.6
clean_niter = clean_niter * 2
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.6094 mas, bmaj=1.188 mas, bpa=-0.802 degrees
! Estimated noise=0.548448 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0387275 Jy
! Component: 100  -  total flux cleaned = 0.0680226 Jy
! Total flux subtracted in 100 components = 0.0680226 Jy
! Clean residual min=-0.014268 max=0.022068 Jy/beam
! Clean residual mean=0.001560 rms=0.004277 Jy/beam
! Combined flux in latest and established models = 1.42448 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 58 model components to the UV plane model.
! The established model now contains 187 components and 1.42448 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.199754Jy  sigma=2.501726
! Fit after  self-cal, rms=0.199670Jy  sigma=2.496789
! Inverting map 
! Added new window around map position (-7.6, 1.4).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0265021 Jy
! Component: 100  -  total flux cleaned = 0.0477608 Jy
! Total flux subtracted in 100 components = 0.0477608 Jy
! Clean residual min=-0.014321 max=0.018701 Jy/beam
! Clean residual mean=0.001479 rms=0.003854 Jy/beam
! Combined flux in latest and established models = 1.47224 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 49 model components to the UV plane model.
! The established model now contains 208 components and 1.47224 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.195329Jy  sigma=2.446732
! Fit after  self-cal, rms=0.195307Jy  sigma=2.443305
! Inverting map 
! Added new window around map position (-5.9, 1.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.021062 Jy
! Component: 100  -  total flux cleaned = 0.0379213 Jy
! Total flux subtracted in 100 components = 0.0379213 Jy
! Clean residual min=-0.013815 max=0.017030 Jy/beam
! Clean residual mean=0.001410 rms=0.003579 Jy/beam
! Combined flux in latest and established models = 1.51017 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 47 model components to the UV plane model.
! The established model now contains 231 components and 1.51017 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.192411Jy  sigma=2.412887
! Fit after  self-cal, rms=0.192333Jy  sigma=2.410278
! Inverting map 
! Added new window around map position (-10.4, 3.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0178938 Jy
! Component: 100  -  total flux cleaned = 0.032337 Jy
! Total flux subtracted in 100 components = 0.032337 Jy
! Clean residual min=-0.013408 max=0.013709 Jy/beam
! Clean residual mean=0.001332 rms=0.003378 Jy/beam
! Combined flux in latest and established models = 1.5425 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 55 model components to the UV plane model.
! The established model now contains 258 components and 1.5425 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.190237Jy  sigma=2.380624
! Fit after  self-cal, rms=0.190226Jy  sigma=2.379287
! Inverting map 
print "****************************************"
! **************************************** 
print "*** Finished Natural Weighting CLEAN ***"
! *** Finished Natural Weighting CLEAN *** 
print "****************************************"
! **************************************** 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 9 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.06     FD        0.90     GC        1.01     HH        1.00*    
!   HN        1.03     KK        1.11     KP        1.11     LA        0.87     
!   MA        1.02     MK        1.09     NL        1.14     NY        0.97     
!   OV        0.94     PT        0.91     SC        1.21     TS        1.13     
!   WF        0.90     WZ        1.00*    GG        1.00*    
! 
!  A total of 48 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.02     FD        0.86     GC        1.11     HH        1.19     
!   HN        1.01     KK        1.13     KP        1.07     LA        0.89     
!   MA        0.88     MK        1.37     NL        1.02     NY        0.90     
!   OV        0.96     PT        0.90     SC        1.12     TS        1.00*    
!   WF        1.02     WZ        0.73     GG        1.00*    
! 
! 
! Correcting IF 2.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.08     FD        0.98     GC        1.02     HH        1.00*    
!   HN        0.92     KK        1.11     KP        0.97     LA        0.88     
!   MA        1.00     MK        0.98     NL        0.96     NY        0.99     
!   OV        1.00     PT        1.08     SC        1.02     TS        0.91     
!   WF        0.94     WZ        1.00*    GG        1.00*    
! 
!  A total of 48 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.06     FD        0.89     GC        1.15     HH        1.13     
!   HN        0.97     KK        1.19     KP        0.96     LA        0.83     
!   MA        0.85     MK        1.32     NL        0.95     NY        0.89     
!   OV        1.03     PT        1.07     SC        0.96     TS        1.00*    
!   WF        1.05     WZ        0.88     GG        1.00*    
! 
! 
! Correcting IF 3.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.04     FD        1.02     GC        1.05     HH        1.00*    
!   HN        0.80     KK        1.07     KP        0.87     LA        1.05     
!   MA        1.01     MK        1.08     NL        1.06     NY        1.06     
!   OV        0.89     PT        1.00     SC        0.91     TS        0.92     
!   WF        0.86     WZ        1.00*    GG        1.00*    
! 
!  A total of 42 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.98     FD        0.93     GC        4.49     HH        1.10     
!   HN        0.80     KK        1.09     KP        0.87     LA        0.99     
!   MA        0.89     MK        1.34     NL        1.01     NY        0.89     
!   OV        0.94     PT        0.97     SC        0.87     TS        1.00*    
!   WF        0.96     WZ        0.88     GG        1.00*    
! 
! 
! Correcting IF 4.
!  A total of 9 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.02     FD        1.07     GC        0.99     HH        1.00*    
!   HN        0.83     KK        1.06     KP        0.88     LA        1.21     
!   MA        0.86     MK        1.06     NL        1.08     NY        0.96     
!   OV        1.00     PT        1.04     SC        0.93     TS        0.91     
!   WF        0.89     WZ        1.00*    GG        1.00*    
! 
!  A total of 42 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00     FD        1.02     GC        1.04     HH        1.38     
!   HN        0.84     KK        1.04     KP        0.88     LA        1.14     
!   MA        0.72     MK        1.40     NL        1.01     NY        0.86     
!   OV        1.01     PT        1.01     SC        0.92     TS        1.00*    
!   WF        1.01     WZ        0.85     GG        1.00*    
! 
! 
! Fit before self-cal, rms=0.190226Jy  sigma=2.379287
! Fit after  self-cal, rms=0.149138Jy  sigma=1.615517
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.6069 mas, bmaj=1.231 mas, bpa=-1.832 degrees
! Estimated noise=0.542599 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (-13.3, 3.1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0141092 Jy
! Component: 100  -  total flux cleaned = 0.0226328 Jy
! Total flux subtracted in 100 components = 0.0226328 Jy
! Clean residual min=-0.007197 max=0.011312 Jy/beam
! Clean residual mean=0.000413 rms=0.002084 Jy/beam
! Combined flux in latest and established models = 1.56514 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 25 model components to the UV plane model.
! The established model now contains 283 components and 1.56514 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147656Jy  sigma=1.591685
! Fit after  self-cal, rms=0.147457Jy  sigma=1.589428
! Inverting map 
! Added new window around map position (-3.8, 1.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0106801 Jy
! Component: 100  -  total flux cleaned = 0.0168218 Jy
! Total flux subtracted in 100 components = 0.0168218 Jy
! Clean residual min=-0.007202 max=0.008277 Jy/beam
! Clean residual mean=0.000363 rms=0.001978 Jy/beam
! Combined flux in latest and established models = 1.58196 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 34 model components to the UV plane model.
! The established model now contains 306 components and 1.58196 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.146374Jy  sigma=1.572591
! Fit after  self-cal, rms=0.146246Jy  sigma=1.571270
! Inverting map 
! Added new window around map position (-12.5, 5.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00795382 Jy
! Component: 100  -  total flux cleaned = 0.0128981 Jy
! Total flux subtracted in 100 components = 0.0128981 Jy
! Clean residual min=-0.007049 max=0.008208 Jy/beam
! Clean residual mean=0.000315 rms=0.001893 Jy/beam
! Combined flux in latest and established models = 1.59486 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 38 model components to the UV plane model.
! The established model now contains 328 components and 1.59486 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.145502Jy  sigma=1.561775
! Fit after  self-cal, rms=0.145430Jy  sigma=1.561164
! Inverting map 
! Added new window around map position (-20.4, 9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00908638 Jy
! Component: 100  -  total flux cleaned = 0.0157984 Jy
! Total flux subtracted in 100 components = 0.0157984 Jy
! Clean residual min=-0.006756 max=0.006667 Jy/beam
! Clean residual mean=0.000254 rms=0.001785 Jy/beam
! Combined flux in latest and established models = 1.61065 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 44 model components to the UV plane model.
! The established model now contains 354 components and 1.61065 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.144780Jy  sigma=1.550390
! Fit after  self-cal, rms=0.144722Jy  sigma=1.549813
! Inverting map 
dynam = dynam - 1.0
doamp = true
solint = solint * (6.0/3.6)
flux_cutoff = imstat(rms) * dynam
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00561547 Jy
! Component: 100  -  total flux cleaned = 0.00982235 Jy
! Total flux subtracted in 100 components = 0.00982235 Jy
! Clean residual min=-0.006642 max=0.006544 Jy/beam
! Clean residual mean=0.000218 rms=0.001727 Jy/beam
! Combined flux in latest and established models = 1.62048 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 52 model components to the UV plane model.
! The established model now contains 373 components and 1.62048 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.144278Jy  sigma=1.544110
! Fit after  self-cal, rms=0.143627Jy  sigma=1.419164
! Inverting map and beam 
! Estimated beam: bmin=0.6075 mas, bmaj=1.235 mas, bpa=-1.937 degrees
! Estimated noise=0.542631 mJy/beam.
! Added new window around map position (-18, 4.1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00800817 Jy
! Component: 100  -  total flux cleaned = 0.013546 Jy
! Total flux subtracted in 100 components = 0.013546 Jy
! Clean residual min=-0.004192 max=0.005952 Jy/beam
! Clean residual mean=0.000393 rms=0.001417 Jy/beam
! Combined flux in latest and established models = 1.63402 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 21 model components to the UV plane model.
! The established model now contains 388 components and 1.63402 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.143030Jy  sigma=1.410715
! Fit after  self-cal, rms=0.144464Jy  sigma=1.404501
! Inverting map and beam 
! Estimated beam: bmin=0.607 mas, bmaj=1.234 mas, bpa=-1.927 degrees
! Estimated noise=0.542402 mJy/beam.
! Added new window around map position (-7.1, 3.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00624569 Jy
! Component: 100  -  total flux cleaned = 0.0108519 Jy
! Total flux subtracted in 100 components = 0.0108519 Jy
! Clean residual min=-0.004039 max=0.005699 Jy/beam
! Clean residual mean=0.000357 rms=0.001351 Jy/beam
! Combined flux in latest and established models = 1.64487 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 38 model components to the UV plane model.
! The established model now contains 409 components and 1.64487 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.144037Jy  sigma=1.398506
! Fit after  self-cal, rms=0.145609Jy  sigma=1.397225
! Inverting map and beam 
! Estimated beam: bmin=0.6069 mas, bmaj=1.234 mas, bpa=-1.928 degrees
! Estimated noise=0.542473 mJy/beam.
! Added new window around map position (-22.5, 9.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00646674 Jy
! Component: 100  -  total flux cleaned = 0.0113602 Jy
! Total flux subtracted in 100 components = 0.0113602 Jy
! Clean residual min=-0.004044 max=0.005240 Jy/beam
! Clean residual mean=0.000324 rms=0.001295 Jy/beam
! Combined flux in latest and established models = 1.65623 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 42 model components to the UV plane model.
! The established model now contains 428 components and 1.65623 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.145226Jy  sigma=1.391605
! Fit after  self-cal, rms=0.146349Jy  sigma=1.391139
! Inverting map and beam 
! Estimated beam: bmin=0.6071 mas, bmaj=1.233 mas, bpa=-1.923 degrees
! Estimated noise=0.542647 mJy/beam.
! Added new window around map position (-11.6, 3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00600698 Jy
! Component: 100  -  total flux cleaned = 0.0103489 Jy
! Total flux subtracted in 100 components = 0.0103489 Jy
! Clean residual min=-0.003976 max=0.004912 Jy/beam
! Clean residual mean=0.000303 rms=0.001245 Jy/beam
! Combined flux in latest and established models = 1.66658 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 40 model components to the UV plane model.
! The established model now contains 442 components and 1.66658 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.146010Jy  sigma=1.385962
! Fit after  self-cal, rms=0.146840Jy  sigma=1.385662
! Inverting map and beam 
! Estimated beam: bmin=0.6071 mas, bmaj=1.234 mas, bpa=-1.945 degrees
! Estimated noise=0.542648 mJy/beam.
! Added new window around map position (12.6, 19.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00483289 Jy
! Component: 100  -  total flux cleaned = 0.00886637 Jy
! Total flux subtracted in 100 components = 0.00886637 Jy
! Clean residual min=-0.003833 max=0.004711 Jy/beam
! Clean residual mean=0.000288 rms=0.001201 Jy/beam
! Combined flux in latest and established models = 1.67545 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 54 model components to the UV plane model.
! The established model now contains 462 components and 1.67545 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.146549Jy  sigma=1.381713
! Fit after  self-cal, rms=0.147178Jy  sigma=1.381467
! Inverting map and beam 
! Estimated beam: bmin=0.6071 mas, bmaj=1.235 mas, bpa=-1.948 degrees
! Estimated noise=0.542505 mJy/beam.
! Added new window around map position (-16.3, 4).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00460753 Jy
! Component: 100  -  total flux cleaned = 0.00807831 Jy
! Total flux subtracted in 100 components = 0.00807831 Jy
! Clean residual min=-0.003667 max=0.004480 Jy/beam
! Clean residual mean=0.000271 rms=0.001164 Jy/beam
! Combined flux in latest and established models = 1.68353 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 48 model components to the UV plane model.
! The established model now contains 477 components and 1.68353 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.146903Jy  sigma=1.377776
! Fit after  self-cal, rms=0.147294Jy  sigma=1.377516
! Inverting map and beam 
! Estimated beam: bmin=0.6069 mas, bmaj=1.235 mas, bpa=-1.953 degrees
! Estimated noise=0.542273 mJy/beam.
! Added new window around map position (-20.3, 4.9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00436164 Jy
! Component: 100  -  total flux cleaned = 0.00769777 Jy
! Total flux subtracted in 100 components = 0.00769777 Jy
! Clean residual min=-0.003630 max=0.004217 Jy/beam
! Clean residual mean=0.000258 rms=0.001128 Jy/beam
! Combined flux in latest and established models = 1.69123 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 52 model components to the UV plane model.
! The established model now contains 495 components and 1.69123 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147021Jy  sigma=1.374163
! Fit after  self-cal, rms=0.147329Jy  sigma=1.373872
! Inverting map and beam 
! Estimated beam: bmin=0.6067 mas, bmaj=1.236 mas, bpa=-1.957 degrees
! Estimated noise=0.542027 mJy/beam.
! Added new window around map position (-9, 3.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00351249 Jy
! Component: 100  -  total flux cleaned = 0.00654382 Jy
! Total flux subtracted in 100 components = 0.00654382 Jy
! Clean residual min=-0.003535 max=0.004076 Jy/beam
! Clean residual mean=0.000246 rms=0.001097 Jy/beam
! Combined flux in latest and established models = 1.69777 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 48 model components to the UV plane model.
! The established model now contains 509 components and 1.69777 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147094Jy  sigma=1.371045
! Fit after  self-cal, rms=0.147397Jy  sigma=1.370761
! Inverting map and beam 
! Estimated beam: bmin=0.6064 mas, bmaj=1.237 mas, bpa=-1.97 degrees
! Estimated noise=0.541786 mJy/beam.
! Added new window around map position (-18.3, 7.9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00378083 Jy
! Component: 100  -  total flux cleaned = 0.00670658 Jy
! Total flux subtracted in 100 components = 0.00670658 Jy
! Clean residual min=-0.003496 max=0.003927 Jy/beam
! Clean residual mean=0.000233 rms=0.001062 Jy/beam
! Combined flux in latest and established models = 1.70448 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 56 model components to the UV plane model.
! The established model now contains 532 components and 1.70448 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147185Jy  sigma=1.367919
! Fit after  self-cal, rms=0.147354Jy  sigma=1.367639
! Inverting map and beam 
! Estimated beam: bmin=0.6061 mas, bmaj=1.238 mas, bpa=-1.979 degrees
! Estimated noise=0.541527 mJy/beam.
! Added new window around map position (-18.4, 12).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0032019 Jy
! Component: 100  -  total flux cleaned = 0.00616074 Jy
! Total flux subtracted in 100 components = 0.00616074 Jy
! Clean residual min=-0.003456 max=0.003790 Jy/beam
! Clean residual mean=0.000220 rms=0.001031 Jy/beam
! Combined flux in latest and established models = 1.71064 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 53 model components to the UV plane model.
! The established model now contains 545 components and 1.71064 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147159Jy  sigma=1.365137
! Fit after  self-cal, rms=0.147360Jy  sigma=1.364969
! Inverting map and beam 
! Estimated beam: bmin=0.6059 mas, bmaj=1.239 mas, bpa=-1.996 degrees
! Estimated noise=0.541259 mJy/beam.
! Added new window around map position (-0.3, -13.9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00287449 Jy
! Component: 100  -  total flux cleaned = 0.005563 Jy
! Total flux subtracted in 100 components = 0.005563 Jy
! Clean residual min=-0.003490 max=0.003558 Jy/beam
! Clean residual mean=0.000210 rms=0.001003 Jy/beam
! Combined flux in latest and established models = 1.7162 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 52 model components to the UV plane model.
! The established model now contains 556 components and 1.7162 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147180Jy  sigma=1.362781
! Fit after  self-cal, rms=0.147363Jy  sigma=1.362659
! Inverting map and beam 
! Estimated beam: bmin=0.6057 mas, bmaj=1.24 mas, bpa=-2.01 degrees
! Estimated noise=0.540969 mJy/beam.
! Added new window around map position (2.7, -12).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00275067 Jy
! Component: 100  -  total flux cleaned = 0.00559798 Jy
! Total flux subtracted in 100 components = 0.00559798 Jy
! Clean residual min=-0.003492 max=0.003480 Jy/beam
! Clean residual mean=0.000195 rms=0.000977 Jy/beam
! Combined flux in latest and established models = 1.7218 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 57 model components to the UV plane model.
! The established model now contains 570 components and 1.7218 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147199Jy  sigma=1.360688
! Fit after  self-cal, rms=0.147361Jy  sigma=1.360551
! Inverting map and beam 
! Estimated beam: bmin=0.6055 mas, bmaj=1.241 mas, bpa=-2.021 degrees
! Estimated noise=0.540682 mJy/beam.
doamp = false
solint = solint / (6.0/3.6)
flux_cutoff = imstat(rms) * dynam
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00181061 Jy
! Component: 100  -  total flux cleaned = 0.00451898 Jy
! Total flux subtracted in 100 components = 0.00451898 Jy
! Clean residual min=-0.003479 max=0.003438 Jy/beam
! Clean residual mean=0.000185 rms=0.000957 Jy/beam
! Combined flux in latest and established models = 1.72632 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 69 model components to the UV plane model.
! The established model now contains 584 components and 1.72632 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.147216Jy  sigma=1.359041
! Fit after  self-cal, rms=0.146828Jy  sigma=1.355175
! Inverting map 
print "***********************************"
! *********************************** 
print "*** Finished Amplitude Self-Cal ***"
! *** Finished Amplitude Self-Cal *** 
print "***********************************"
! *********************************** 
save 3C274_X
! Writing UV FITS file: 3C274_X.uvf
! Writing 584 model components to file: 3C274_X.mod
! wwins: Wrote 26 windows to 3C274_X.win
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6055 x 1.241 at -2.021 degrees (North through East)
! Clean map  min=-0.0033239  max=0.59567 Jy/beam
! Writing clean map to FITS file: 3C274_X.fits
! Writing difmap environment to: 3C274_X.par
print "****************************"
! **************************** 
print "*** Now, make some plots ***"
! *** Now, make some plots *** 
print "****************************"
! **************************** 
![@gonzo_plot 3C274_X]
device 3C274_Xuv.ps/vps,1,2
! Attempting to open device: '3C274_Xuv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device 3C274_Xim.ps/vps
! Attempting to open device: '3C274_Xim.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=0.6055 mas, bmaj=1.241 mas, bpa=-2.021 degrees
! Estimated noise=0.540682 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6055 x 1.241 at -2.021 degrees (North through East)
! Clean map  min=-0.0033239  max=0.59567 Jy/beam
![Exited script file: gonzo_plot]
![Exited script file: gonzo5-0.1]
quit
! Quitting program
! Log file 3C274_X.log closed on Thu Aug 14 14:31:54 2008