! Started logfile: VR422201_X.log on Wed Apr  2 00:30:39 2008
observe VR422201_X.edt
! Reading UV FITS file: VR422201_X.edt
! AN table 1: 88 integrations on 136 of 136 possible baselines.
! AN table 2: 111 integrations on 136 of 136 possible baselines.
! AN table 3: 131 integrations on 136 of 136 possible baselines.
! Apparent sampling: 0.287455 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: VR422201
! 
! 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 1214 lines of history.
! 
! Reading 51604 visibilities.
![@gonzo6-0.1 VR422201_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 = 6.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 1019 telescope corrections were flagged in sub-array 1.
!  A total of 1245 telescope corrections were flagged in sub-array 2.
!  A total of 1501 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 2.
!  A total of 1029 telescope corrections were flagged in sub-array 1.
!  A total of 1224 telescope corrections were flagged in sub-array 2.
!  A total of 1498 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 3.
!  A total of 1013 telescope corrections were flagged in sub-array 1.
!  A total of 1140 telescope corrections were flagged in sub-array 2.
!  A total of 1489 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 4.
!  A total of 1071 telescope corrections were flagged in sub-array 1.
!  A total of 1176 telescope corrections were flagged in sub-array 2.
!  A total of 1550 telescope corrections were flagged in sub-array 3.
! 
! Fit before self-cal, rms=0.768295Jy  sigma=0.013908
! Fit after  self-cal, rms=0.766847Jy  sigma=0.013897
! 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.5313 mas, bmaj=0.8017 mas, bpa=-5.77 degrees
! Estimated noise=575.372 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.963014 Jy
! Total flux subtracted in 50 components = 0.963014 Jy
! Clean residual min=-0.044871 max=0.272057 Jy/beam
! Clean residual mean=0.000320 rms=0.013539 Jy/beam
! Combined flux in latest and established models = 0.963014 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 3 model components to the UV plane model.
! The established model now contains 3 components and 0.963014 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.787797Jy  sigma=0.014261
! Fit after  self-cal, rms=0.787790Jy  sigma=0.014260
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.276612 Jy
! Total flux subtracted in 50 components = 0.276612 Jy
! Clean residual min=-0.037176 max=0.221117 Jy/beam
! Clean residual mean=0.000268 rms=0.011093 Jy/beam
! Combined flux in latest and established models = 1.23963 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 7 model components to the UV plane model.
! The established model now contains 9 components and 1.23963 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.653538Jy  sigma=0.011828
! Fit after  self-cal, rms=0.653135Jy  sigma=0.011824
! Inverting map 
! Added new window around map position (-0.5, -2.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.211344 Jy
! Total flux subtracted in 50 components = 0.211344 Jy
! Clean residual min=-0.031284 max=0.116736 Jy/beam
! Clean residual mean=0.000225 rms=0.009351 Jy/beam
! Combined flux in latest and established models = 1.45097 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 22 components and 1.45097 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.554852Jy  sigma=0.010091
! Fit after  self-cal, rms=0.546850Jy  sigma=0.010023
! Inverting map 
! Added new window around map position (-0.5, -1.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.153952 Jy
! Total flux subtracted in 50 components = 0.153952 Jy
! Clean residual min=-0.026051 max=0.086589 Jy/beam
! Clean residual mean=0.000197 rms=0.007945 Jy/beam
! Combined flux in latest and established models = 1.60492 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 30 components and 1.60492 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.481723Jy  sigma=0.008820
! Fit after  self-cal, rms=0.478273Jy  sigma=0.008795
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.121727 Jy
! Total flux subtracted in 50 components = 0.121727 Jy
! Clean residual min=-0.022785 max=0.073449 Jy/beam
! Clean residual mean=0.000173 rms=0.007004 Jy/beam
! Combined flux in latest and established models = 1.72665 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 49 components and 1.72665 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.429853Jy  sigma=0.007879
! Fit after  self-cal, rms=0.427723Jy  sigma=0.007864
! Inverting map 
! Added new window around map position (0.1, -2.4).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.10113 Jy
! Total flux subtracted in 50 components = 0.10113 Jy
! Clean residual min=-0.020275 max=0.060792 Jy/beam
! Clean residual mean=0.000154 rms=0.006290 Jy/beam
! Combined flux in latest and established models = 1.82778 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 19 model components to the UV plane model.
! The established model now contains 57 components and 1.82778 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.390136Jy  sigma=0.007142
! Fit after  self-cal, rms=0.388555Jy  sigma=0.007132
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0856272 Jy
! Total flux subtracted in 50 components = 0.0856272 Jy
! Clean residual min=-0.018648 max=0.053339 Jy/beam
! Clean residual mean=0.000137 rms=0.005725 Jy/beam
! Combined flux in latest and established models = 1.9134 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 22 model components to the UV plane model.
! The established model now contains 63 components and 1.9134 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.358392Jy  sigma=0.006537
! Fit after  self-cal, rms=0.357226Jy  sigma=0.006530
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0728924 Jy
! Total flux subtracted in 50 components = 0.0728924 Jy
! Clean residual min=-0.017423 max=0.047143 Jy/beam
! Clean residual mean=0.000123 rms=0.005277 Jy/beam
! Combined flux in latest and established models = 1.9863 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 22 model components to the UV plane model.
! The established model now contains 71 components and 1.9863 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.333067Jy  sigma=0.006040
! Fit after  self-cal, rms=0.332280Jy  sigma=0.006035
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0627954 Jy
! Total flux subtracted in 50 components = 0.0627954 Jy
! Clean residual min=-0.016434 max=0.042589 Jy/beam
! Clean residual mean=0.000110 rms=0.004923 Jy/beam
! Combined flux in latest and established models = 2.04909 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 86 components and 2.04909 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.312828Jy  sigma=0.005631
! Fit after  self-cal, rms=0.312160Jy  sigma=0.005627
! Inverting map 
! Added new window around map position (-0.5, -3.4).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0557753 Jy
! Total flux subtracted in 50 components = 0.0557753 Jy
! Clean residual min=-0.016064 max=0.038508 Jy/beam
! Clean residual mean=0.000099 rms=0.004628 Jy/beam
! Combined flux in latest and established models = 2.10487 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 97 components and 2.10487 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.295979Jy  sigma=0.005284
! Fit after  self-cal, rms=0.295377Jy  sigma=0.005281
! Inverting map 
! Added new window around map position (0.5, 1.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0496135 Jy
! Total flux subtracted in 50 components = 0.0496135 Jy
! Clean residual min=-0.015323 max=0.032118 Jy/beam
! Clean residual mean=0.000089 rms=0.004376 Jy/beam
! Combined flux in latest and established models = 2.15448 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 108 components and 2.15448 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.281675Jy  sigma=0.004984
! Fit after  self-cal, rms=0.281369Jy  sigma=0.004982
! Inverting map 
! Added new window around map position (0.1, -3.4).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0441582 Jy
! Total flux subtracted in 50 components = 0.0441582 Jy
! Clean residual min=-0.014940 max=0.030229 Jy/beam
! Clean residual mean=0.000080 rms=0.004174 Jy/beam
! Combined flux in latest and established models = 2.19864 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 118 components and 2.19864 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.269940Jy  sigma=0.004727
! Fit after  self-cal, rms=0.269605Jy  sigma=0.004726
! Inverting map 
! Added new window around map position (0.7, -3.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0399974 Jy
! Total flux subtracted in 50 components = 0.0399974 Jy
! Clean residual min=-0.014365 max=0.026604 Jy/beam
! Clean residual mean=0.000072 rms=0.003996 Jy/beam
! Combined flux in latest and established models = 2.23864 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 128 components and 2.23864 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.259886Jy  sigma=0.004503
! Fit after  self-cal, rms=0.259566Jy  sigma=0.004502
! Inverting map 
! Added new window around map position (-1.1, -3.2).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0362129 Jy
! Total flux subtracted in 50 components = 0.0362129 Jy
! Clean residual min=-0.013997 max=0.023553 Jy/beam
! Clean residual mean=0.000065 rms=0.003843 Jy/beam
! Combined flux in latest and established models = 2.27485 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 139 components and 2.27485 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.251148Jy  sigma=0.004306
! Fit after  self-cal, rms=0.250827Jy  sigma=0.004305
! Inverting map 
! Added new window around map position (-1.1, -2.2).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0329094 Jy
! Total flux subtracted in 50 components = 0.0329094 Jy
! Clean residual min=-0.013566 max=0.022060 Jy/beam
! Clean residual mean=0.000058 rms=0.003709 Jy/beam
! Combined flux in latest and established models = 2.30776 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 35 model components to the UV plane model.
! The established model now contains 151 components and 2.30776 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.243534Jy  sigma=0.004131
! Fit after  self-cal, rms=0.243223Jy  sigma=0.004130
! Inverting map 
! Added new window around map position (0.2, -4.5).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0301121 Jy
! Total flux subtracted in 50 components = 0.0301121 Jy
! Clean residual min=-0.013085 max=0.021645 Jy/beam
! Clean residual mean=0.000052 rms=0.003591 Jy/beam
! Combined flux in latest and established models = 2.33787 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 30 model components to the UV plane model.
! The established model now contains 159 components and 2.33787 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.236991Jy  sigma=0.003979
! Fit after  self-cal, rms=0.236726Jy  sigma=0.003978
! Inverting map 
! Added new window around map position (-1.1, -4.8).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0279251 Jy
! Total flux subtracted in 50 components = 0.0279251 Jy
! Clean residual min=-0.012959 max=0.020560 Jy/beam
! Clean residual mean=0.000046 rms=0.003486 Jy/beam
! Combined flux in latest and established models = 2.3658 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 33 model components to the UV plane model.
! The established model now contains 170 components and 2.3658 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.231346Jy  sigma=0.003847
! Fit after  self-cal, rms=0.231110Jy  sigma=0.003846
! 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.6715 mas, bmaj=1.091 mas, bpa=0.5483 degrees
! Estimated noise=268.412 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.0834851 Jy
! Component: 100  -  total flux cleaned = 0.145449 Jy
! Total flux subtracted in 100 components = 0.145449 Jy
! Clean residual min=-0.015309 max=0.042632 Jy/beam
! Clean residual mean=0.002325 rms=0.009043 Jy/beam
! Combined flux in latest and established models = 2.51125 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 73 model components to the UV plane model.
! The established model now contains 226 components and 2.51125 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.208671Jy  sigma=0.003205
! Fit after  self-cal, rms=0.207631Jy  sigma=0.003193
! 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 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        1.02     GC        1.03     HN        1.01     
!   KK        1.00*    KP        0.99     LA        0.98     MC        0.92     
!   MK        1.00*    NL        1.13     NY        0.99     ON        0.99     
!   OV        1.00*    PT        0.95     SC        0.98     WF        0.93     
!   WZ        1.09     
! 
!  A total of 21 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        0.98     GC        1.00*    HN        0.95     
!   KK        1.00*    KP        0.91     LA        0.90     MC        1.00*    
!   MK        1.00     NL        1.16     NY        0.92     ON        1.20     
!   OV        0.98     PT        0.88     SC        1.04     WF        0.92     
!   WZ        2.63     
! 
!  A total of 12 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00     FD        1.02     GC        1.07     HN        0.98     
!   KK        1.06     KP        0.93     LA        0.95     MC        1.07     
!   MK        0.98     NL        1.05     NY        0.95     ON        1.00     
!   OV        1.01     PT        0.92     SC        0.99     WF        0.92     
!   WZ        1.14     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.01     FD        1.02     GC        1.12     HN        1.00     
!   KK        1.00*    KP        0.99     LA        1.00     MC        0.93     
!   MK        1.00*    NL        1.13     NY        1.25     ON        1.00     
!   OV        1.00*    PT        0.96     SC        1.01     WF        0.92     
!   WZ        1.01     
! 
!  A total of 30 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        0.97     GC        1.00*    HN        0.95     
!   KK        1.00*    KP        0.92     LA        0.91     MC        1.00*    
!   MK        0.99     NL        1.12     NY        0.94     ON        1.18     
!   OV        0.96     PT        0.91     SC        1.04     WF        0.94     
!   WZ        2.10     
! 
!  A total of 15 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00     FD        1.01     GC        1.13     HN        0.99     
!   KK        1.05     KP        0.97     LA        0.95     MC        1.06     
!   MK        0.97     NL        1.06     NY        0.96     ON        1.00     
!   OV        1.00     PT        0.94     SC        0.99     WF        0.94     
!   WZ        1.05     
! 
! 
! Correcting IF 3.
!  A total of 9 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.03     FD        1.03     GC        1.30     HN        1.01     
!   KK        1.00*    KP        1.01     LA        0.98     MC        0.93     
!   MK        1.00*    NL        1.14     NY        0.98     ON        0.99     
!   OV        1.00*    PT        0.96     SC        0.96     WF        1.01     
!   WZ        1.04     
! 
!  A total of 24 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.04     FD        1.00     GC        1.00*    HN        0.96     
!   KK        1.01     KP        0.94     LA        0.90     MC        1.00*    
!   MK        0.98     NL        1.11     NY        0.95     ON        1.18     
!   OV        0.98     PT        0.92     SC        1.05     WF        1.02     
!   WZ        1.95     
! 
!  A total of 18 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.01     FD        1.00     GC        1.11     HN        0.98     
!   KK        1.02     KP        0.96     LA        0.95     MC        1.07     
!   MK        0.97     NL        1.02     NY        0.96     ON        1.02     
!   OV        1.00     PT        0.94     SC        0.97     WF        0.98     
!   WZ        1.05     
! 
! 
! Correcting IF 4.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.01     FD        1.03     GC        1.00     HN        1.03     
!   KK        1.00*    KP        1.00     LA        1.00     MC        0.93     
!   MK        1.00*    NL        1.10     NY        1.24     ON        0.97     
!   OV        1.00*    PT        0.98     SC        0.98     WF        1.06     
!   WZ        1.00*    
! 
!  A total of 27 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        0.99     GC        1.00*    HN        0.95     
!   KK        1.03     KP        0.93     LA        0.92     MC        1.00*    
!   MK        1.00     NL        1.15     NY        0.98     ON        1.20     
!   OV        0.98     PT        0.95     SC        1.12     WF        1.10     
!   WZ        1.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.98     FD        1.01     GC        1.10     HN        0.97     
!   KK        1.03     KP        0.96     LA        0.94     MC        1.07     
!   MK        0.98     NL        1.02     NY        0.95     ON        1.04     
!   OV        1.00     PT        0.95     SC        1.02     WF        1.08     
!   WZ        1.00*    
! 
! 
! Fit before self-cal, rms=0.207631Jy  sigma=0.003193
! Fit after  self-cal, rms=0.152621Jy  sigma=0.001924
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.6803 mas, bmaj=1.12 mas, bpa=1.076 degrees
! Estimated noise=264.3 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.0263903 Jy
! Component: 100  -  total flux cleaned = 0.0356599 Jy
! Total flux subtracted in 100 components = 0.0356599 Jy
! Clean residual min=-0.009619 max=0.016948 Jy/beam
! Clean residual mean=0.001406 rms=0.002421 Jy/beam
! Combined flux in latest and established models = 2.54691 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 23 model components to the UV plane model.
! The established model now contains 249 components and 2.54691 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.148680Jy  sigma=0.001813
! Fit after  self-cal, rms=0.148042Jy  sigma=0.001806
! Inverting map 
! Added new window around map position (-0.2, -6.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0155016 Jy
! Component: 100  -  total flux cleaned = 0.0220044 Jy
! Total flux subtracted in 100 components = 0.0220044 Jy
! Clean residual min=-0.007557 max=0.013206 Jy/beam
! Clean residual mean=0.001228 rms=0.002061 Jy/beam
! Combined flux in latest and established models = 2.56891 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 39 model components to the UV plane model.
! The established model now contains 282 components and 2.56891 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.146082Jy  sigma=0.001754
! Fit after  self-cal, rms=0.145771Jy  sigma=0.001751
! Inverting map 
! Added new window around map position (1.5, -5.9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0124946 Jy
! Component: 100  -  total flux cleaned = 0.016664 Jy
! Total flux subtracted in 100 components = 0.016664 Jy
! Clean residual min=-0.006324 max=0.010461 Jy/beam
! Clean residual mean=0.001098 rms=0.001853 Jy/beam
! Combined flux in latest and established models = 2.58557 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 301 components and 2.58557 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.144370Jy  sigma=0.001718
! Fit after  self-cal, rms=0.144164Jy  sigma=0.001716
! Inverting map 
! Added new window around map position (4, 5.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00667542 Jy
! Component: 100  -  total flux cleaned = 0.0112561 Jy
! Total flux subtracted in 100 components = 0.0112561 Jy
! Clean residual min=-0.005613 max=0.009025 Jy/beam
! Clean residual mean=0.001006 rms=0.001747 Jy/beam
! Combined flux in latest and established models = 2.59683 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 314 components and 2.59683 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.143242Jy  sigma=0.001697
! Fit after  self-cal, rms=0.143140Jy  sigma=0.001697
! Inverting map 
! Added new window around map position (6.2, -13.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00507103 Jy
! Component: 100  -  total flux cleaned = 0.00806608 Jy
! Total flux subtracted in 100 components = 0.00806608 Jy
! Clean residual min=-0.005541 max=0.008829 Jy/beam
! Clean residual mean=0.000948 rms=0.001669 Jy/beam
! Combined flux in latest and established models = 2.6049 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 40 model components to the UV plane model.
! The established model now contains 331 components and 2.6049 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.142436Jy  sigma=0.001683
! Fit after  self-cal, rms=0.142362Jy  sigma=0.001683
! Inverting map 
! Added new window around map position (-7.2, 1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00601672 Jy
! Component: 100  -  total flux cleaned = 0.00924248 Jy
! Total flux subtracted in 100 components = 0.00924248 Jy
! Clean residual min=-0.005326 max=0.008204 Jy/beam
! Clean residual mean=0.000876 rms=0.001604 Jy/beam
! Combined flux in latest and established models = 2.61414 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 35 model components to the UV plane model.
! The established model now contains 343 components and 2.61414 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.141769Jy  sigma=0.001671
! Fit after  self-cal, rms=0.141711Jy  sigma=0.001671
! Inverting map 
! Added new window around map position (2, -9.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0054656 Jy
! Component: 100  -  total flux cleaned = 0.00771396 Jy
! Total flux subtracted in 100 components = 0.00771396 Jy
! Clean residual min=-0.005871 max=0.007690 Jy/beam
! Clean residual mean=0.000815 rms=0.001561 Jy/beam
! Combined flux in latest and established models = 2.62185 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 28 model components to the UV plane model.
! The established model now contains 354 components and 2.62185 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.141170Jy  sigma=0.001661
! Fit after  self-cal, rms=0.141124Jy  sigma=0.001661
! 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)
! Added new window around map position (-0.2, -18.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00731058 Jy
! Component: 100  -  total flux cleaned = 0.0106031 Jy
! Total flux subtracted in 100 components = 0.0106031 Jy
! Clean residual min=-0.005630 max=0.007088 Jy/beam
! Clean residual mean=0.000735 rms=0.001495 Jy/beam
! Combined flux in latest and established models = 2.63246 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 37 model components to the UV plane model.
! The established model now contains 373 components and 2.63246 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.140610Jy  sigma=0.001652
! Fit after  self-cal, rms=0.123825Jy  sigma=0.001439
! Inverting map and beam 
! Estimated beam: bmin=0.6794 mas, bmaj=1.12 mas, bpa=1.193 degrees
! Estimated noise=263.561 mJy/beam.
! Added new window around map position (1, -4.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00929896 Jy
! Component: 100  -  total flux cleaned = 0.0139948 Jy
! Total flux subtracted in 100 components = 0.0139948 Jy
! Clean residual min=-0.004283 max=0.007025 Jy/beam
! Clean residual mean=0.000722 rms=0.001157 Jy/beam
! Combined flux in latest and established models = 2.64645 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 400 components and 2.64645 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.123141Jy  sigma=0.001422
! Fit after  self-cal, rms=0.122582Jy  sigma=0.001416
! Inverting map and beam 
! Estimated beam: bmin=0.6793 mas, bmaj=1.12 mas, bpa=1.198 degrees
! Estimated noise=263.585 mJy/beam.
! Added new window around map position (3, -5.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.008212 Jy
! Component: 100  -  total flux cleaned = 0.0122838 Jy
! Total flux subtracted in 100 components = 0.0122838 Jy
! Clean residual min=-0.003919 max=0.005337 Jy/beam
! Clean residual mean=0.000667 rms=0.001068 Jy/beam
! Combined flux in latest and established models = 2.65873 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 39 model components to the UV plane model.
! The established model now contains 425 components and 2.65873 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.122025Jy  sigma=0.001402
! Fit after  self-cal, rms=0.121755Jy  sigma=0.001399
! Inverting map and beam 
! Estimated beam: bmin=0.6792 mas, bmaj=1.12 mas, bpa=1.206 degrees
! Estimated noise=263.569 mJy/beam.
! Added new window around map position (-2.3, -4.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00613718 Jy
! Component: 100  -  total flux cleaned = 0.0100539 Jy
! Total flux subtracted in 100 components = 0.0100539 Jy
! Clean residual min=-0.003613 max=0.005025 Jy/beam
! Clean residual mean=0.000627 rms=0.001008 Jy/beam
! Combined flux in latest and established models = 2.66879 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 447 components and 2.66879 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.121331Jy  sigma=0.001389
! Fit after  self-cal, rms=0.121170Jy  sigma=0.001387
! Inverting map and beam 
! Estimated beam: bmin=0.6791 mas, bmaj=1.12 mas, bpa=1.218 degrees
! Estimated noise=263.563 mJy/beam.
! Added new window around map position (6.7, -19.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00527349 Jy
! Component: 100  -  total flux cleaned = 0.0088419 Jy
! Total flux subtracted in 100 components = 0.0088419 Jy
! Clean residual min=-0.003323 max=0.004728 Jy/beam
! Clean residual mean=0.000590 rms=0.000952 Jy/beam
! Combined flux in latest and established models = 2.67763 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 45 model components to the UV plane model.
! The established model now contains 463 components and 2.67763 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.120820Jy  sigma=0.001379
! Fit after  self-cal, rms=0.120695Jy  sigma=0.001378
! Inverting map and beam 
! Estimated beam: bmin=0.679 mas, bmaj=1.12 mas, bpa=1.23 degrees
! Estimated noise=263.547 mJy/beam.
! Added new window around map position (-8.3, -13.9).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00459035 Jy
! Component: 100  -  total flux cleaned = 0.0081136 Jy
! Total flux subtracted in 100 components = 0.0081136 Jy
! Clean residual min=-0.003064 max=0.004510 Jy/beam
! Clean residual mean=0.000548 rms=0.000910 Jy/beam
! Combined flux in latest and established models = 2.68574 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 51 model components to the UV plane model.
! The established model now contains 481 components and 2.68574 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.120380Jy  sigma=0.001371
! Fit after  self-cal, rms=0.120284Jy  sigma=0.001370
! Inverting map and beam 
! Estimated beam: bmin=0.6788 mas, bmaj=1.12 mas, bpa=1.241 degrees
! Estimated noise=263.525 mJy/beam.
! Added new window around map position (-4.9, -3.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00402254 Jy
! Component: 100  -  total flux cleaned = 0.00728743 Jy
! Total flux subtracted in 100 components = 0.00728743 Jy
! Clean residual min=-0.002834 max=0.004375 Jy/beam
! Clean residual mean=0.000514 rms=0.000867 Jy/beam
! Combined flux in latest and established models = 2.69303 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 45 model components to the UV plane model.
! The established model now contains 497 components and 2.69303 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.120024Jy  sigma=0.001364
! Fit after  self-cal, rms=0.119937Jy  sigma=0.001364
! Inverting map and beam 
! Estimated beam: bmin=0.6788 mas, bmaj=1.12 mas, bpa=1.251 degrees
! Estimated noise=263.495 mJy/beam.
! Added new window around map position (0.3, 2.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.003414 Jy
! Component: 100  -  total flux cleaned = 0.0066575 Jy
! Total flux subtracted in 100 components = 0.0066575 Jy
! Clean residual min=-0.002689 max=0.003942 Jy/beam
! Clean residual mean=0.000483 rms=0.000836 Jy/beam
! Combined flux in latest and established models = 2.69969 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 44 model components to the UV plane model.
! The established model now contains 510 components and 2.69969 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.119687Jy  sigma=0.001359
! Fit after  self-cal, rms=0.119615Jy  sigma=0.001358
! Inverting map and beam 
! Estimated beam: bmin=0.6787 mas, bmaj=1.12 mas, bpa=1.266 degrees
! Estimated noise=263.458 mJy/beam.
! Added new window around map position (8.6, -15.1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00337635 Jy
! Component: 100  -  total flux cleaned = 0.00657089 Jy
! Total flux subtracted in 100 components = 0.00657088 Jy
! Clean residual min=-0.002594 max=0.003554 Jy/beam
! Clean residual mean=0.000452 rms=0.000801 Jy/beam
! Combined flux in latest and established models = 2.70626 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 46 model components to the UV plane model.
! The established model now contains 526 components and 2.70626 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.119397Jy  sigma=0.001354
! Fit after  self-cal, rms=0.119331Jy  sigma=0.001353
! Inverting map and beam 
! Estimated beam: bmin=0.6786 mas, bmaj=1.12 mas, bpa=1.278 degrees
! Estimated noise=263.43 mJy/beam.
! Added new window around map position (0.8, -1.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00306386 Jy
! Component: 100  -  total flux cleaned = 0.00571569 Jy
! Total flux subtracted in 100 components = 0.00571569 Jy
! Clean residual min=-0.002473 max=0.003467 Jy/beam
! Clean residual mean=0.000419 rms=0.000773 Jy/beam
! Combined flux in latest and established models = 2.71197 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 50 model components to the UV plane model.
! The established model now contains 537 components and 2.71197 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.119154Jy  sigma=0.001349
! Fit after  self-cal, rms=0.119098Jy  sigma=0.001349
! Inverting map and beam 
! Estimated beam: bmin=0.6785 mas, bmaj=1.121 mas, bpa=1.29 degrees
! Estimated noise=263.411 mJy/beam.
! Added new window around map position (5, -10.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00322413 Jy
! Component: 100  -  total flux cleaned = 0.00564248 Jy
! Total flux subtracted in 100 components = 0.00564248 Jy
! Clean residual min=-0.002351 max=0.003379 Jy/beam
! Clean residual mean=0.000387 rms=0.000739 Jy/beam
! Combined flux in latest and established models = 2.71762 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 39 model components to the UV plane model.
! The established model now contains 551 components and 2.71762 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118932Jy  sigma=0.001345
! Fit after  self-cal, rms=0.118875Jy  sigma=0.001345
! Inverting map and beam 
! Estimated beam: bmin=0.6785 mas, bmaj=1.121 mas, bpa=1.301 degrees
! Estimated noise=263.393 mJy/beam.
! Added new window around map position (11.5, -13.3).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00298103 Jy
! Component: 100  -  total flux cleaned = 0.00541276 Jy
! Total flux subtracted in 100 components = 0.00541276 Jy
! Clean residual min=-0.002327 max=0.003114 Jy/beam
! Clean residual mean=0.000363 rms=0.000713 Jy/beam
! Combined flux in latest and established models = 2.72303 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 44 model components to the UV plane model.
! The established model now contains 560 components and 2.72303 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118726Jy  sigma=0.001342
! Fit after  self-cal, rms=0.118692Jy  sigma=0.001342
! Inverting map and beam 
! Estimated beam: bmin=0.6785 mas, bmaj=1.121 mas, bpa=1.31 degrees
! Estimated noise=263.378 mJy/beam.
! Added new window around map position (13.6, 15.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00279009 Jy
! Component: 100  -  total flux cleaned = 0.00521241 Jy
! Total flux subtracted in 100 components = 0.00521241 Jy
! Clean residual min=-0.002371 max=0.002864 Jy/beam
! Clean residual mean=0.000336 rms=0.000688 Jy/beam
! Combined flux in latest and established models = 2.72824 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 49 model components to the UV plane model.
! The established model now contains 572 components and 2.72824 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118560Jy  sigma=0.001339
! Fit after  self-cal, rms=0.118534Jy  sigma=0.001339
! Inverting map and beam 
! Estimated beam: bmin=0.6784 mas, bmaj=1.121 mas, bpa=1.318 degrees
! Estimated noise=263.371 mJy/beam.
! Added new window around map position (4.3, -20.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00246294 Jy
! Component: 100  -  total flux cleaned = 0.00485433 Jy
! Total flux subtracted in 100 components = 0.00485433 Jy
! Clean residual min=-0.002342 max=0.002763 Jy/beam
! Clean residual mean=0.000319 rms=0.000662 Jy/beam
! Combined flux in latest and established models = 2.7331 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 49 model components to the UV plane model.
! The established model now contains 586 components and 2.7331 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118416Jy  sigma=0.001336
! Fit after  self-cal, rms=0.118385Jy  sigma=0.001336
! Inverting map and beam 
! Estimated beam: bmin=0.6784 mas, bmaj=1.121 mas, bpa=1.326 degrees
! Estimated noise=263.361 mJy/beam.
! Added new window around map position (3.3, -11.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00225787 Jy
! Component: 100  -  total flux cleaned = 0.00441431 Jy
! Total flux subtracted in 100 components = 0.00441431 Jy
! Clean residual min=-0.002283 max=0.002684 Jy/beam
! Clean residual mean=0.000301 rms=0.000642 Jy/beam
! Combined flux in latest and established models = 2.73751 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 47 model components to the UV plane model.
! The established model now contains 597 components and 2.73751 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118279Jy  sigma=0.001334
! Fit after  self-cal, rms=0.118257Jy  sigma=0.001334
! Inverting map and beam 
! Estimated beam: bmin=0.6784 mas, bmaj=1.121 mas, bpa=1.334 degrees
! Estimated noise=263.356 mJy/beam.
! Added new window around map position (-0.3, -11.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00230923 Jy
! Component: 100  -  total flux cleaned = 0.004293 Jy
! Total flux subtracted in 100 components = 0.004293 Jy
! Clean residual min=-0.002282 max=0.002505 Jy/beam
! Clean residual mean=0.000276 rms=0.000624 Jy/beam
! Combined flux in latest and established models = 2.7418 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 617 components and 2.7418 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118155Jy  sigma=0.001331
! Fit after  self-cal, rms=0.118138Jy  sigma=0.001331
! Inverting map and beam 
! Estimated beam: bmin=0.6784 mas, bmaj=1.121 mas, bpa=1.34 degrees
! Estimated noise=263.354 mJy/beam.
! Added new window around map position (8.5, -12.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00226359 Jy
! Component: 100  -  total flux cleaned = 0.00414937 Jy
! Total flux subtracted in 100 components = 0.00414937 Jy
! Clean residual min=-0.002202 max=0.002449 Jy/beam
! Clean residual mean=0.000259 rms=0.000605 Jy/beam
! Combined flux in latest and established models = 2.74595 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 631 components and 2.74595 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118045Jy  sigma=0.001330
! Fit after  self-cal, rms=0.118034Jy  sigma=0.001330
! Inverting map and beam 
! Estimated beam: bmin=0.6784 mas, bmaj=1.121 mas, bpa=1.346 degrees
! Estimated noise=263.352 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.0018178 Jy
! Component: 100  -  total flux cleaned = 0.00358352 Jy
! Total flux subtracted in 100 components = 0.00358352 Jy
! Clean residual min=-0.002192 max=0.002367 Jy/beam
! Clean residual mean=0.000243 rms=0.000591 Jy/beam
! Combined flux in latest and established models = 2.74954 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 60 model components to the UV plane model.
! The established model now contains 644 components and 2.74954 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.117957Jy  sigma=0.001328
! Fit after  self-cal, rms=0.117627Jy  sigma=0.001323
! Inverting map 
print "***********************************"
! *********************************** 
print "*** Finished Amplitude Self-Cal ***"
! *** Finished Amplitude Self-Cal *** 
print "***********************************"
! *********************************** 
save VR422201_X
! Writing UV FITS file: VR422201_X.uvf
! Writing 644 model components to file: VR422201_X.mod
! wwins: Wrote 35 windows to VR422201_X.win
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6784 x 1.121 at 1.346 degrees (North through East)
! Clean map  min=-0.0020615  max=1.2649 Jy/beam
! Writing clean map to FITS file: VR422201_X.fits
! Writing difmap environment to: VR422201_X.par
print "****************************"
! **************************** 
print "*** Now, make some plots ***"
! *** Now, make some plots *** 
print "****************************"
! **************************** 
![@gonzo_plot VR422201_X]
device VR422201_Xuv.ps/vps,1,2
! Attempting to open device: 'VR422201_Xuv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device VR422201_Xim.ps/vps
! Attempting to open device: 'VR422201_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.6784 mas, bmaj=1.121 mas, bpa=1.346 degrees
! Estimated noise=263.352 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6784 x 1.121 at 1.346 degrees (North through East)
! Clean map  min=-0.0020615  max=1.2649 Jy/beam
![Exited script file: gonzo_plot]
![Exited script file: gonzo6-0.1]
quit
! Quitting program
! Log file VR422201_X.log closed on Wed Apr  2 00:31:29 2008