! Started logfile: 0710+439_S.log on Wed Aug 20 11:30:14 2008
observe 0710+439_S.edt
! Reading UV FITS file: 0710+439_S.edt
! AN table 1: 188 integrations on 171 of 171 possible baselines.
! AN table 2: 72 integrations on 171 of 171 possible baselines.
! Apparent sampling: 0.167139 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 0710+439
! 
! 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.22497e+09        8e+06          1        8e+06
!  02        2  2.24497e+09        8e+06          1        8e+06
!  03        3  2.33497e+09        8e+06          1        8e+06
!  04        4  2.36497e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1615 lines of history.
! 
! Reading 29724 visibilities.
![@gonzo7-0.4 0710+439_S]
print "***********************************"
! *********************************** 
print "*** Using gonzo instruction set ***"
! *** Using gonzo instruction set *** 
print "***********************************"
! *********************************** 
float field_size
field_size = 1024
float field_cell
field_cell = 0.4
integer clean_niter
clean_niter = 50
float clean_gain
clean_gain = 0.03
float dynam
dynam = 7.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.400x0.400 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 2465 telescope corrections were flagged in sub-array 1.
!  A total of 764 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 2465 telescope corrections were flagged in sub-array 1.
!  A total of 758 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 2465 telescope corrections were flagged in sub-array 1.
!  A total of 761 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 2465 telescope corrections were flagged in sub-array 1.
!  A total of 758 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.635842Jy  sigma=7.333903
! Fit after  self-cal, rms=0.634901Jy  sigma=7.328083
! 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=2.115 mas, bmaj=3.155 mas, bpa=15.04 degrees
! Estimated noise=0.82008 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.456073 Jy
! Total flux subtracted in 50 components = 0.456073 Jy
! Clean residual min=-0.063016 max=0.281729 Jy/beam
! Clean residual mean=0.000656 rms=0.018249 Jy/beam
! Combined flux in latest and established models = 0.456073 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 5 model components to the UV plane model.
! The established model now contains 5 components and 0.456073 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.600135Jy  sigma=7.407013
! Fit after  self-cal, rms=0.600161Jy  sigma=7.406868
! Inverting map 
! Added new window around map position (0.4, -8.4).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.247631 Jy
! Total flux subtracted in 50 components = 0.247631 Jy
! Clean residual min=-0.048594 max=0.115868 Jy/beam
! Clean residual mean=0.000535 rms=0.013894 Jy/beam
! Combined flux in latest and established models = 0.703704 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 12 components and 0.703704 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.490956Jy  sigma=6.010281
! Fit after  self-cal, rms=0.479523Jy  sigma=5.872569
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.1571 Jy
! Total flux subtracted in 50 components = 0.1571 Jy
! Clean residual min=-0.040980 max=0.122040 Jy/beam
! Clean residual mean=0.000468 rms=0.011555 Jy/beam
! Combined flux in latest and established models = 0.860804 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 21 components and 0.860804 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.420369Jy  sigma=5.040968
! Fit after  self-cal, rms=0.418262Jy  sigma=5.016365
! Inverting map 
! Added new window around map position (-0.4, -23.6).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.128528 Jy
! Total flux subtracted in 50 components = 0.128528 Jy
! Clean residual min=-0.036032 max=0.071052 Jy/beam
! Clean residual mean=0.000384 rms=0.009839 Jy/beam
! Combined flux in latest and established models = 0.989333 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 33 components and 0.989333 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.377464Jy  sigma=4.438926
! Fit after  self-cal, rms=0.375152Jy  sigma=4.416626
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0991422 Jy
! Total flux subtracted in 50 components = 0.0991422 Jy
! Clean residual min=-0.032169 max=0.059481 Jy/beam
! Clean residual mean=0.000332 rms=0.008588 Jy/beam
! Combined flux in latest and established models = 1.08847 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 40 components and 1.08847 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.345637Jy  sigma=3.966092
! Fit after  self-cal, rms=0.344930Jy  sigma=3.959877
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0813656 Jy
! Total flux subtracted in 50 components = 0.0813656 Jy
! Clean residual min=-0.028883 max=0.051625 Jy/beam
! Clean residual mean=0.000289 rms=0.007636 Jy/beam
! Combined flux in latest and established models = 1.16984 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 16 model components to the UV plane model.
! The established model now contains 45 components and 1.16984 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.322523Jy  sigma=3.611119
! Fit after  self-cal, rms=0.322294Jy  sigma=3.608174
! 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=2.706 mas, bmaj=3.986 mas, bpa=12.69 degrees
! Estimated noise=0.399314 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.8, -4).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.12183 Jy
! Component: 100  -  total flux cleaned = 0.206812 Jy
! Total flux subtracted in 100 components = 0.206812 Jy
! Clean residual min=-0.021075 max=0.048832 Jy/beam
! Clean residual mean=0.002936 rms=0.008250 Jy/beam
! Combined flux in latest and established models = 1.37665 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 53 model components to the UV plane model.
! The established model now contains 83 components and 1.37665 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276579Jy  sigma=2.864755
! Fit after  self-cal, rms=0.276433Jy  sigma=2.855776
! 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 102 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.90     FD        1.00     GC        1.07     HH        1.00*    
!   HN        0.87     KK        0.71     KP        0.86     LA        1.01     
!   MA        0.81     MK        0.95     NL        0.90     NY        1.11     
!   OV        0.90     PT        0.95     SC        0.80     TS        1.00*    
!   WF        1.00     WZ        1.07     GG        1.00*    
! 
!  A total of 9 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.87     FD        1.04     GC        0.97     HH        1.00*    
!   HN        0.88     KK        1.00*    KP        0.84     LA        1.00     
!   MA        3.28     MK        1.00*    NL        0.87     NY        1.21     
!   OV        1.00*    PT        0.98     SC        0.89     TS        1.00*    
!   WF        0.92     WZ        1.04     GG        1.00*    
! 
! 
! Correcting IF 2.
!  A total of 102 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.90     FD        0.93     GC        0.99     HH        1.00*    
!   HN        0.95     KK        0.69     KP        0.98     LA        1.03     
!   MA        1.64     MK        1.01     NL        1.00     NY        1.06     
!   OV        0.85     PT        0.81     SC        0.94     TS        1.00*    
!   WF        1.01     WZ        1.00     GG        1.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.86     FD        0.96     GC        0.94     HH        1.00*    
!   HN        0.94     KK        1.00*    KP        0.96     LA        1.01     
!   MA        1.31     MK        1.00*    NL        1.01     NY        1.19     
!   OV        1.00*    PT        0.83     SC        1.10     TS        1.00*    
!   WF        0.91     WZ        0.94     GG        1.00*    
! 
! 
! Correcting IF 3.
!  A total of 102 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.96     FD        0.93     GC        1.00     HH        1.00*    
!   HN        1.20     KK        0.69     KP        1.11     LA        0.87     
!   MA        0.88     MK        0.98     NL        0.96     NY        1.08     
!   OV        0.96     PT        0.89     SC        1.07     TS        1.00*    
!   WF        1.09     WZ        1.01     GG        1.00*    
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.92     FD        0.94     GC        0.60     HH        1.00*    
!   HN        1.19     KK        1.00*    KP        1.05     LA        0.85     
!   MA        4.50     MK        1.00*    NL        0.92     NY        1.15     
!   OV        1.00*    PT        0.87     SC        1.15     TS        1.00*    
!   WF        0.97     WZ        0.97     GG        1.00*    
! 
! 
! Correcting IF 4.
!  A total of 102 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.95     FD        0.81     GC        1.01     HH        1.00*    
!   HN        1.17     KK        0.75     KP        1.10     LA        0.74     
!   MA        1.07     MK        0.96     NL        0.97     NY        1.09     
!   OV        0.88     PT        0.85     SC        1.05     TS        1.00*    
!   WF        1.11     WZ        1.07     GG        1.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.94     FD        0.86     GC        0.71     HH        1.00*    
!   HN        1.14     KK        1.00*    KP        1.03     LA        0.72     
!   MA        1.18     MK        1.00*    NL        0.95     NY        1.17     
!   OV        1.00*    PT        0.88     SC        1.20     TS        1.00*    
!   WF        0.98     WZ        1.02     GG        1.00*    
! 
! 
! Fit before self-cal, rms=0.276433Jy  sigma=2.855776
! Fit after  self-cal, rms=0.255524Jy  sigma=2.019463
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.807 mas, bmaj=4.054 mas, bpa=10.99 degrees
! Estimated noise=0.375864 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.0225767 Jy
! Component: 100  -  total flux cleaned = 0.0315856 Jy
! Total flux subtracted in 100 components = 0.0315856 Jy
! Clean residual min=-0.011145 max=0.019393 Jy/beam
! Clean residual mean=0.001143 rms=0.002815 Jy/beam
! Combined flux in latest and established models = 1.40824 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 16 model components to the UV plane model.
! The established model now contains 98 components and 1.40824 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.251873Jy  sigma=1.958978
! Fit after  self-cal, rms=0.251923Jy  sigma=1.951453
! Inverting map 
! Added new window around map position (-0.4, 6.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0144633 Jy
! Component: 100  -  total flux cleaned = 0.0227708 Jy
! Total flux subtracted in 100 components = 0.0227708 Jy
! Clean residual min=-0.008872 max=0.013756 Jy/beam
! Clean residual mean=0.000860 rms=0.002359 Jy/beam
! Combined flux in latest and established models = 1.43101 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 105 components and 1.43101 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.250307Jy  sigma=1.916561
! Fit after  self-cal, rms=0.250323Jy  sigma=1.915379
! 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 (-31.2, -9.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0111541 Jy
! Component: 100  -  total flux cleaned = 0.0164379 Jy
! Total flux subtracted in 100 components = 0.0164379 Jy
! Clean residual min=-0.008158 max=0.011034 Jy/beam
! Clean residual mean=0.000656 rms=0.002093 Jy/beam
! Combined flux in latest and established models = 1.44745 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 30 model components to the UV plane model.
! The established model now contains 127 components and 1.44745 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.249383Jy  sigma=1.898784
! Fit after  self-cal, rms=0.258104Jy  sigma=1.825990
! Inverting map and beam 
! Estimated beam: bmin=2.807 mas, bmaj=4.077 mas, bpa=11.43 degrees
! Estimated noise=0.374897 mJy/beam.
! Added new window around map position (-6, 5.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00922945 Jy
! Component: 100  -  total flux cleaned = 0.013906 Jy
! Total flux subtracted in 100 components = 0.013906 Jy
! Clean residual min=-0.007103 max=0.008666 Jy/beam
! Clean residual mean=0.000530 rms=0.001600 Jy/beam
! Combined flux in latest and established models = 1.46135 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 25 model components to the UV plane model.
! The established model now contains 147 components and 1.46135 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.257311Jy  sigma=1.813896
! Fit after  self-cal, rms=0.257520Jy  sigma=1.810357
! Inverting map and beam 
! Estimated beam: bmin=2.802 mas, bmaj=4.074 mas, bpa=11.52 degrees
! Estimated noise=0.373973 mJy/beam.
! Added new window around map position (6.4, -22.8).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00702794 Jy
! Component: 100  -  total flux cleaned = 0.011251 Jy
! Total flux subtracted in 100 components = 0.011251 Jy
! Clean residual min=-0.006456 max=0.007240 Jy/beam
! Clean residual mean=0.000431 rms=0.001379 Jy/beam
! Combined flux in latest and established models = 1.4726 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 22 model components to the UV plane model.
! The established model now contains 157 components and 1.4726 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.257028Jy  sigma=1.801531
! Fit after  self-cal, rms=0.261067Jy  sigma=1.801998
! Inverting map and beam 
! Estimated beam: bmin=2.799 mas, bmaj=4.071 mas, bpa=11.55 degrees
! Estimated noise=0.373389 mJy/beam.
! Added new window around map position (29.2, -9.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00519756 Jy
! Component: 100  -  total flux cleaned = 0.00903028 Jy
! Total flux subtracted in 100 components = 0.00903028 Jy
! Clean residual min=-0.006116 max=0.006619 Jy/beam
! Clean residual mean=0.000367 rms=0.001237 Jy/beam
! Combined flux in latest and established models = 1.48163 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 23 model components to the UV plane model.
! The established model now contains 165 components and 1.48163 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.260717Jy  sigma=1.796243
! Fit after  self-cal, rms=0.264747Jy  sigma=1.797093
! Inverting map and beam 
! Estimated beam: bmin=2.798 mas, bmaj=4.07 mas, bpa=11.56 degrees
! Estimated noise=0.372801 mJy/beam.
! Added new window around map position (-11.2, 16.4).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00472721 Jy
! Component: 100  -  total flux cleaned = 0.00757022 Jy
! Total flux subtracted in 100 components = 0.00757022 Jy
! Clean residual min=-0.005702 max=0.005909 Jy/beam
! Clean residual mean=0.000306 rms=0.001139 Jy/beam
! Combined flux in latest and established models = 1.4892 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 26 model components to the UV plane model.
! The established model now contains 177 components and 1.4892 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.264474Jy  sigma=1.792825
! Fit after  self-cal, rms=0.268268Jy  sigma=1.793512
! Inverting map and beam 
! Estimated beam: bmin=2.797 mas, bmaj=4.067 mas, bpa=11.57 degrees
! Estimated noise=0.372262 mJy/beam.
! Added new window around map position (47.6, 1.2).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00392579 Jy
! Component: 100  -  total flux cleaned = 0.00578109 Jy
! Total flux subtracted in 100 components = 0.00578109 Jy
! Clean residual min=-0.005327 max=0.005017 Jy/beam
! Clean residual mean=0.000277 rms=0.001062 Jy/beam
! Combined flux in latest and established models = 1.49499 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 29 model components to the UV plane model.
! The established model now contains 187 components and 1.49499 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.268052Jy  sigma=1.790702
! Fit after  self-cal, rms=0.271420Jy  sigma=1.791195
! Inverting map and beam 
! Estimated beam: bmin=2.797 mas, bmaj=4.066 mas, bpa=11.56 degrees
! Estimated noise=0.3718 mJy/beam.
! Added new window around map position (-5.6, -19.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00328127 Jy
! Component: 100  -  total flux cleaned = 0.00510821 Jy
! Total flux subtracted in 100 components = 0.00510821 Jy
! Clean residual min=-0.005218 max=0.004275 Jy/beam
! Clean residual mean=0.000235 rms=0.001006 Jy/beam
! Combined flux in latest and established models = 1.50009 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 28 model components to the UV plane model.
! The established model now contains 195 components and 1.50009 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.271257Jy  sigma=1.788907
! Fit after  self-cal, rms=0.274295Jy  sigma=1.789122
! Inverting map and beam 
! Estimated beam: bmin=2.797 mas, bmaj=4.065 mas, bpa=11.54 degrees
! Estimated noise=0.371363 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.00165094 Jy
! Component: 100  -  total flux cleaned = 0.00327491 Jy
! Total flux subtracted in 100 components = 0.00327491 Jy
! Clean residual min=-0.005015 max=0.004185 Jy/beam
! Clean residual mean=0.000214 rms=0.000971 Jy/beam
! Combined flux in latest and established models = 1.50337 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 200 components and 1.50337 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.274168Jy  sigma=1.787986
! Fit after  self-cal, rms=0.273492Jy  sigma=1.771291
! Inverting map 
print "***********************************"
! *********************************** 
print "*** Finished Amplitude Self-Cal ***"
! *** Finished Amplitude Self-Cal *** 
print "***********************************"
! *********************************** 
save 0710+439_S
! Writing UV FITS file: 0710+439_S.uvf
! Writing 200 model components to file: 0710+439_S.mod
! wwins: Wrote 12 windows to 0710+439_S.win
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.797 x 4.065 at 11.54 degrees (North through East)
! Clean map  min=-0.0045973  max=0.5802 Jy/beam
! Writing clean map to FITS file: 0710+439_S.fits
! Writing difmap environment to: 0710+439_S.par
print "****************************"
! **************************** 
print "*** Now, make some plots ***"
! *** Now, make some plots *** 
print "****************************"
! **************************** 
![@gonzo_plot 0710+439_S]
device 0710+439_Suv.ps/vps,1,2
! Attempting to open device: '0710+439_Suv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device 0710+439_Sim.ps/vps
! Attempting to open device: '0710+439_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.797 mas, bmaj=4.065 mas, bpa=11.54 degrees
! Estimated noise=0.371363 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.797 x 4.065 at 11.54 degrees (North through East)
! Clean map  min=-0.0045973  max=0.5802 Jy/beam
![Exited script file: gonzo_plot]
![Exited script file: gonzo7-0.4]
quit
! Quitting program
! Log file 0710+439_S.log closed on Wed Aug 20 11:30:34 2008