! Started logfile: J0522-3627_S_1998_10_01_pus_uvs.log on Thu Feb  7 17:05:18 2008
obs J0522-3627_S_1998_10_01_pus_uva.fits
! Reading UV FITS file: J0522-3627_S_1998_10_01_pus_uva.fits
! AN table 1: 267 integrations on 136 of 136 possible baselines.
! Apparent sampling: 0.107237 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J0522-3627
! 
! 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.22298e+09        4e+06          1        4e+06
!  02        2  2.24298e+09        4e+06          1        4e+06
!  03        3  2.33298e+09        4e+06          1        4e+06
!  04        4  2.36298e+09        4e+06          1        4e+06
! 
! Polarization(s): RR
! 
! Read 1537 lines of history.
! 
! Reading 15576 visibilities.
select I
! Polarization I is unavailable.
! 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
print "field_size =",field_size
! field_size = 512 
print "field_cell =",field_cell
! field_cell = 0.5 
print "taper_size =",taper_size
! taper_size = 20 
![@muppet J0522-3627_S_1998_10_01_pus]
integer clean_niter
clean_niter = 100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float soltime1
soltime1 = 120
float soltime2
soltime2 = 20
float thresh
thresh = 0.5
float win_mult
win_mult = 1.8
float time_av
time_av = 16
float old_peak
float new_peak
float flux_cutoff
mapsize field_size, field_cell
! Map grid = 512x512 pixels with 0.500x0.500 milli-arcsec cellsize.
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 3286 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 3286 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 3286 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 3286 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.712955Jy  sigma=0.004502
! Fit after  self-cal, rms=0.554523Jy  sigma=0.003527
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
uvw 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.
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.352 mas, bmaj=8.81 mas, bpa=-9.677 degrees
! Estimated noise=1405.33 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (0, 0).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.607724 Jy
! Component: 100  -  total flux cleaned = 0.768293 Jy
! Total flux subtracted in 100 components = 0.768293 Jy
! Clean residual min=-0.051649 max=0.129411 Jy/beam
! Clean residual mean=0.003154 rms=0.019330 Jy/beam
! Combined flux in latest and established models = 0.768293 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 16 components and 0.768293 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.533969Jy  sigma=0.003474
! Fit after  self-cal, rms=0.533915Jy  sigma=0.003472
! Inverting map 
! Added new window around map position (-21, 19.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.122065 Jy
! Component: 100  -  total flux cleaned = 0.198914 Jy
! Total flux subtracted in 100 components = 0.198914 Jy
! Clean residual min=-0.039167 max=0.072588 Jy/beam
! Clean residual mean=0.002323 rms=0.014764 Jy/beam
! Combined flux in latest and established models = 0.967207 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 34 components and 0.967206 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.484848Jy  sigma=0.003093
! Fit after  self-cal, rms=0.482273Jy  sigma=0.003079
! Inverting map 
print "************** FINISHED UNIFORM WEIGHTING CLEAN **************"
! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** 
! 
uvw 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
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.882 mas, bmaj=9.985 mas, bpa=-12.04 degrees
! Estimated noise=877.508 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.132611 Jy
! Component: 100  -  total flux cleaned = 0.224452 Jy
! Component: 150  -  total flux cleaned = 0.291957 Jy
! Component: 200  -  total flux cleaned = 0.341899 Jy
! Total flux subtracted in 200 components = 0.341899 Jy
! Clean residual min=-0.034788 max=0.110604 Jy/beam
! Clean residual mean=0.005990 rms=0.021438 Jy/beam
! Combined flux in latest and established models = 1.30911 Jy
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 73 components and 1.30911 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.440748Jy  sigma=0.002705
! Fit after  self-cal, rms=0.437900Jy  sigma=0.002688
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 42 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        0.98     GC        1.00*    HN        1.18     
!   KK        0.98     KP        0.89     LA        0.84     MK        0.91     
!   NL        1.27     NY        1.00*    ON        1.00*    OV        1.15     
!   PT        0.92     SC        1.11     WF        1.14     MC        1.00*    
!   GN        1.04     
! 
! 
! Correcting IF 2.
!  A total of 42 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        0.96     GC        1.00*    HN        1.23     
!   KK        0.99     KP        0.88     LA        0.88     MK        0.99     
!   NL        1.23     NY        1.00*    ON        1.00*    OV        1.14     
!   PT        0.89     SC        1.11     WF        1.15     MC        1.00*    
!   GN        1.08     
! 
! 
! Correcting IF 3.
!  A total of 42 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00     GC        1.00*    HN        1.23     
!   KK        0.98     KP        0.92     LA        0.87     MK        0.91     
!   NL        1.24     NY        1.00*    ON        1.00*    OV        1.11     
!   PT        0.88     SC        1.02     WF        1.10     MC        1.00*    
!   GN        1.06     
! 
! 
! Correcting IF 4.
!  A total of 42 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        0.95     GC        1.00*    HN        1.09     
!   KK        0.99     KP        0.90     LA        0.90     MK        0.92     
!   NL        1.23     NY        1.00*    ON        1.00*    OV        1.11     
!   PT        0.86     SC        1.02     WF        1.07     MC        1.00*    
!   GN        1.02     
! 
! 
! Fit before self-cal, rms=0.437900Jy  sigma=0.002688
! Fit after  self-cal, rms=0.445464Jy  sigma=0.002521
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.885 mas, bmaj=10.1 mas, bpa=-12.36 degrees
! Estimated noise=860.079 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (-14.5, 21.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0800919 Jy
! Component: 100  -  total flux cleaned = 0.11343 Jy
! Component: 150  -  total flux cleaned = 0.118245 Jy
! Component: 200  -  total flux cleaned = 0.118267 Jy
! Total flux subtracted in 200 components = 0.118267 Jy
! Clean residual min=-0.042842 max=0.024736 Jy/beam
! Clean residual mean=0.000191 rms=0.007331 Jy/beam
! Combined flux in latest and established models = 1.42737 Jy
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 110 components and 1.42737 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.430987Jy  sigma=0.002404
! Fit after  self-cal, rms=0.431658Jy  sigma=0.002398
! Inverting map 
selfcal true, true, soltime1
! Performing amp+phase self-cal over 120 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.431658Jy  sigma=0.002398
! Fit after  self-cal, rms=0.420664Jy  sigma=0.002364
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.892 mas, bmaj=10.04 mas, bpa=-12.11 degrees
! Estimated noise=857.105 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0230377 Jy
! Component: 100  -  total flux cleaned = 0.0269029 Jy
! Component: 150  -  total flux cleaned = 0.0290355 Jy
! Component: 200  -  total flux cleaned = 0.0294528 Jy
! Total flux subtracted in 200 components = 0.0294528 Jy
! Clean residual min=-0.031820 max=0.022110 Jy/beam
! Clean residual mean=-0.000396 rms=0.005154 Jy/beam
! Combined flux in latest and established models = 1.45682 Jy
selfcal
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 142 components and 1.45683 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.418745Jy  sigma=0.002349
! Fit after  self-cal, rms=0.418291Jy  sigma=0.002347
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00141678 Jy
! Component: 100  -  total flux cleaned = 0.00142189 Jy
! Component: 150  -  total flux cleaned = 0.00107744 Jy
! Component: 200  -  total flux cleaned = 0.00140029 Jy
! Total flux subtracted in 200 components = 0.00140029 Jy
! Clean residual min=-0.030541 max=0.020822 Jy/beam
! Clean residual mean=-0.000432 rms=0.004936 Jy/beam
! Combined flux in latest and established models = 1.45823 Jy
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 151 components and 1.45823 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.417675Jy  sigma=0.002343
! Fit after  self-cal, rms=0.417649Jy  sigma=0.002343
! Inverting map 
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.417649Jy  sigma=0.002343
! Fit after  self-cal, rms=0.414884Jy  sigma=0.002327
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.894 mas, bmaj=10.03 mas, bpa=-12.07 degrees
! Estimated noise=855.122 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00976086 Jy
! Component: 100  -  total flux cleaned = 0.0107646 Jy
! Component: 150  -  total flux cleaned = 0.011361 Jy
! Component: 200  -  total flux cleaned = 0.0118959 Jy
! Total flux subtracted in 200 components = 0.0118959 Jy
! Clean residual min=-0.025804 max=0.021422 Jy/beam
! Clean residual mean=-0.000265 rms=0.004204 Jy/beam
! Combined flux in latest and established models = 1.47012 Jy
selfcal
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 171 components and 1.47012 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.414439Jy  sigma=0.002323
! Fit after  self-cal, rms=0.414364Jy  sigma=0.002323
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00115425 Jy
! Component: 100  -  total flux cleaned = 0.00165494 Jy
! Component: 150  -  total flux cleaned = 0.00188507 Jy
! Component: 200  -  total flux cleaned = 0.0018837 Jy
! Total flux subtracted in 200 components = 0.0018837 Jy
! Clean residual min=-0.025155 max=0.020700 Jy/beam
! Clean residual mean=-0.000320 rms=0.004136 Jy/beam
! Combined flux in latest and established models = 1.47201 Jy
! Performing phase self-cal
! Adding 26 model components to the UV plane model.
! The established model now contains 178 components and 1.47201 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.414127Jy  sigma=0.002321
! Fit after  self-cal, rms=0.414113Jy  sigma=0.002321
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.414113Jy  sigma=0.002321
! Fit after  self-cal, rms=0.413350Jy  sigma=0.002320
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.632 mas, bmaj=21.95 mas, bpa=-9.667 degrees
! Estimated noise=1109.06 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00602904 Jy
! Component: 100  -  total flux cleaned = 0.00635519 Jy
! Component: 150  -  total flux cleaned = 0.00663803 Jy
! Component: 200  -  total flux cleaned = 0.00711924 Jy
! Total flux subtracted in 200 components = 0.00711924 Jy
! Clean residual min=-0.021591 max=0.012989 Jy/beam
! Clean residual mean=-0.000520 rms=0.004953 Jy/beam
! Combined flux in latest and established models = 1.47912 Jy
selfcal
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 183 components and 1.47912 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.413733Jy  sigma=0.002322
! Fit after  self-cal, rms=0.413510Jy  sigma=0.002321
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000996647 Jy
! Component: 100  -  total flux cleaned = 0.000207652 Jy
! Component: 150  -  total flux cleaned = 0.000459616 Jy
! Component: 200  -  total flux cleaned = 0.000679853 Jy
! Total flux subtracted in 200 components = 0.000679853 Jy
! Clean residual min=-0.022403 max=0.012822 Jy/beam
! Clean residual mean=-0.000562 rms=0.004988 Jy/beam
! Combined flux in latest and established models = 1.4798 Jy
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 200 components and 1.4798 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.413801Jy  sigma=0.002323
! Fit after  self-cal, rms=0.413642Jy  sigma=0.002321
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
win_mult = win_mult / 1.2
print "************** FINISHED AMPLITUDE SELF-CAL **************"
! ************** FINISHED AMPLITUDE SELF-CAL ************** 
restore
! Inverting map and beam 
! Estimated beam: bmin=2.895 mas, bmaj=10.03 mas, bpa=-12.07 degrees
! Estimated noise=853.818 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.895 x 10.03 at -12.07 degrees (North through East)
! Clean map  min=-0.022774  max=0.8037 Jy/beam
if(peak(flux) > thresh)
  selfcal true, true
  clean clean_niter,clean_gain
  selfcal
  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
    new_peak = peak(flux)
    until(new_peak<=flux_cutoff)
  selfcal true, true
  end if
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.413642Jy  sigma=0.002321
! Fit after  self-cal, rms=0.450705Jy  sigma=0.002088
! Inverting map and beam 
! Estimated beam: bmin=2.904 mas, bmaj=10.07 mas, bpa=-12.13 degrees
! Estimated noise=853.379 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00153729 Jy
! Component: 100  -  total flux cleaned = 0.00156167 Jy
! Component: 150  -  total flux cleaned = 0.00207759 Jy
! Component: 200  -  total flux cleaned = 0.00254641 Jy
! Total flux subtracted in 200 components = 0.00254641 Jy
! Clean residual min=-0.020408 max=0.019911 Jy/beam
! Clean residual mean=-0.000113 rms=0.003669 Jy/beam
! Combined flux in latest and established models = 1.48235 Jy
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 216 components and 1.48235 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.450350Jy  sigma=0.002085
! Fit after  self-cal, rms=0.450386Jy  sigma=0.002085
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00123248 Jy
! Component: 100  -  total flux cleaned = 0.00144306 Jy
! Component: 150  -  total flux cleaned = 0.00164594 Jy
! Component: 200  -  total flux cleaned = 0.0018328 Jy
! Total flux subtracted in 200 components = 0.0018328 Jy
! Clean residual min=-0.019916 max=0.019328 Jy/beam
! Clean residual mean=-0.000161 rms=0.003611 Jy/beam
! Combined flux in latest and established models = 1.48418 Jy
! Performing phase self-cal
! Adding 24 model components to the UV plane model.
! The established model now contains 220 components and 1.48418 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.450174Jy  sigma=0.002083
! Fit after  self-cal, rms=0.450191Jy  sigma=0.002083
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.450191Jy  sigma=0.002083
! Fit after  self-cal, rms=0.455632Jy  sigma=0.002083
print "************** CLEARING MODEL AND STARTING OVER **********"
! ************** CLEARING MODEL AND STARTING OVER ********** 
clrmod true
! clrmod: Cleared the established and tentative models.
uvw 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.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.363 mas, bmaj=8.871 mas, bpa=-9.71 degrees
! Estimated noise=1384.77 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.596169 Jy
! Component: 100  -  total flux cleaned = 0.777942 Jy
! Component: 150  -  total flux cleaned = 0.891695 Jy
! Component: 200  -  total flux cleaned = 0.977105 Jy
! Total flux subtracted in 200 components = 0.977105 Jy
! Clean residual min=-0.030930 max=0.050497 Jy/beam
! Clean residual mean=0.002268 rms=0.011586 Jy/beam
! Combined flux in latest and established models = 0.977105 Jy
uvw 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.
flux_cutoff = imstat(rms) * dynam
! Adding 34 model components to the UV plane model.
! The established model now contains 34 components and 0.977105 Jy
! Inverting map and beam 
! Estimated beam: bmin=2.904 mas, bmaj=10.07 mas, bpa=-12.14 degrees
! Estimated noise=852.582 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  keep
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.133454 Jy
! Component: 100  -  total flux cleaned = 0.2265 Jy
! Component: 150  -  total flux cleaned = 0.295463 Jy
! Component: 200  -  total flux cleaned = 0.347366 Jy
! Total flux subtracted in 200 components = 0.347366 Jy
! Clean residual min=-0.015526 max=0.047479 Jy/beam
! Clean residual mean=0.004013 rms=0.009748 Jy/beam
! Combined flux in latest and established models = 1.32447 Jy
! Adding 67 model components to the UV plane model.
! The established model now contains 88 components and 1.32447 Jy
! Inverting map 
! Added new window around map position (-6.5, 7).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.046406 Jy
! Component: 100  -  total flux cleaned = 0.0775853 Jy
! Component: 150  -  total flux cleaned = 0.100842 Jy
! Component: 200  -  total flux cleaned = 0.118488 Jy
! Total flux subtracted in 200 components = 0.118487 Jy
! Clean residual min=-0.013937 max=0.012191 Jy/beam
! Clean residual mean=0.000985 rms=0.003972 Jy/beam
! Combined flux in latest and established models = 1.44296 Jy
! Adding 68 model components to the UV plane model.
! The established model now contains 145 components and 1.44296 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.622 mas, bmaj=22 mas, bpa=-9.748 degrees
! Estimated noise=1104.37 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0224306 Jy
! Component: 100  -  total flux cleaned = 0.0348034 Jy
! Component: 150  -  total flux cleaned = 0.041878 Jy
! Component: 200  -  total flux cleaned = 0.0447483 Jy
! Total flux subtracted in 200 components = 0.0447483 Jy
! Clean residual min=-0.014521 max=0.012426 Jy/beam
! Clean residual mean=-0.000368 rms=0.004220 Jy/beam
! Combined flux in latest and established models = 1.48771 Jy
flux_cutoff = imstat(rms) * dynam
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  keep
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000294875 Jy
! Component: 100  -  total flux cleaned = 0.000431844 Jy
! Component: 150  -  total flux cleaned = 0.000558832 Jy
! Component: 200  -  total flux cleaned = 0.00056134 Jy
! Total flux subtracted in 200 components = 0.00056134 Jy
! Clean residual min=-0.013833 max=0.012302 Jy/beam
! Clean residual mean=-0.000398 rms=0.004147 Jy/beam
! Combined flux in latest and established models = 1.48827 Jy
! Adding 70 model components to the UV plane model.
! The established model now contains 194 components and 1.48827 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=2.904 mas, bmaj=10.07 mas, bpa=-12.14 degrees
! Estimated noise=852.582 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00168778 Jy
! Component: 100  -  total flux cleaned = 0.00148342 Jy
! Component: 150  -  total flux cleaned = 0.00148745 Jy
! Component: 200  -  total flux cleaned = 0.00132771 Jy
! Total flux subtracted in 200 components = 0.00132771 Jy
! Clean residual min=-0.013199 max=0.010136 Jy/beam
! Clean residual mean=-0.000205 rms=0.002935 Jy/beam
! Combined flux in latest and established models = 1.4896 Jy
selfcal
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 212 components and 1.4896 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.453906Jy  sigma=0.002074
! Fit after  self-cal, rms=0.453567Jy  sigma=0.002071
wmodel J0522-3627_S_1998_10_01_pus_map.mod
! Writing 212 model components to file: J0522-3627_S_1998_10_01_pus_map.mod
wobs J0522-3627_S_1998_10_01_pus_uvs.fits
! Writing UV FITS file: J0522-3627_S_1998_10_01_pus_uvs.fits
wwins J0522-3627_S_1998_10_01_pus_map.win
! wwins: Wrote 4 windows to J0522-3627_S_1998_10_01_pus_map.win
x = (field_size-8) * field_cell / 4
addwin -x,x,-x,x
clean (field_size*4),0.01
! Inverting map 
! clean: niter=2048  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.00647649 Jy
! Component: 100  -  total flux cleaned = -0.0102675 Jy
! Component: 150  -  total flux cleaned = -0.012838 Jy
! Component: 200  -  total flux cleaned = -0.0123934 Jy
! Component: 250  -  total flux cleaned = -0.0115547 Jy
! Component: 300  -  total flux cleaned = -0.010886 Jy
! Component: 350  -  total flux cleaned = -0.00999393 Jy
! Component: 400  -  total flux cleaned = -0.0091478 Jy
! Component: 450  -  total flux cleaned = -0.00856414 Jy
! Component: 500  -  total flux cleaned = -0.00811451 Jy
! Component: 550  -  total flux cleaned = -0.00768109 Jy
! Component: 600  -  total flux cleaned = -0.00705254 Jy
! Component: 650  -  total flux cleaned = -0.00664732 Jy
! Component: 700  -  total flux cleaned = -0.00625442 Jy
! Component: 750  -  total flux cleaned = -0.00577497 Jy
! Component: 800  -  total flux cleaned = -0.00521119 Jy
! Component: 850  -  total flux cleaned = -0.00438684 Jy
! Component: 900  -  total flux cleaned = -0.00393855 Jy
! Component: 950  -  total flux cleaned = -0.00340907 Jy
! Component: 1000  -  total flux cleaned = -0.00280274 Jy
! Component: 1050  -  total flux cleaned = -0.00237737 Jy
! Component: 1100  -  total flux cleaned = -0.00204286 Jy
! Component: 1150  -  total flux cleaned = -0.00179656 Jy
! Component: 1200  -  total flux cleaned = -0.00155344 Jy
! Component: 1250  -  total flux cleaned = -0.00115579 Jy
! Component: 1300  -  total flux cleaned = -0.00115709 Jy
! Component: 1350  -  total flux cleaned = -0.00107939 Jy
! Component: 1400  -  total flux cleaned = -0.00100286 Jy
! Component: 1450  -  total flux cleaned = -0.00085233 Jy
! Component: 1500  -  total flux cleaned = -0.000852671 Jy
! Component: 1550  -  total flux cleaned = -0.000778734 Jy
! Component: 1600  -  total flux cleaned = -0.000633 Jy
! Component: 1650  -  total flux cleaned = -0.000490326 Jy
! Component: 1700  -  total flux cleaned = -0.000419794 Jy
! Component: 1750  -  total flux cleaned = -0.000279426 Jy
! Component: 1800  -  total flux cleaned = -3.88273e-06 Jy
! Component: 1850  -  total flux cleaned = 0.000132693 Jy
! Component: 1900  -  total flux cleaned = 0.000335285 Jy
! Component: 1950  -  total flux cleaned = 0.000468803 Jy
! Component: 2000  -  total flux cleaned = 0.000601051 Jy
! Total flux subtracted in 2048 components = 0.000797574 Jy
! Clean residual min=-0.003547 max=0.003933 Jy/beam
! Clean residual mean=-0.000072 rms=0.001402 Jy/beam
! Combined flux in latest and established models = 1.49039 Jy
keep
! Adding 153 model components to the UV plane model.
! The established model now contains 364 components and 1.49039 Jy
![Exited script file: muppet]
wmap "J0522-3627_S_1998_10_01_pus\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.904 x 10.07 at -12.14 degrees (North through East)
! Clean map  min=-0.0099782  max=0.80543 Jy/beam
! Writing clean map to FITS file: J0522-3627_S_1998_10_01_pus_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 0.805435 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00140448 
device "J0522-3627_S_1998_10_01_pus\_map.eps/VPS"
! Attempting to open device: 'J0522-3627_S_1998_10_01_pus_map.eps/VPS'
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.523128 0.523128 1.04626 2.09251 4.18502 8.37004 16.7401 33.4802 66.9603
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0042 
print "peak=",int(1000*peak_flux)/1000
! peak= 0.805 
mapcolor none
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
mapplot cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.904 x 10.07 at -12.14 degrees (North through East)
! Clean map  min=-0.0099782  max=0.80543 Jy/beam
mapcolor color, 0.9
! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5.
device "J0522-3627_S_1998_10_01_pus\_map.gif/vgif"
! Attempting to open device: 'J0522-3627_S_1998_10_01_pus_map.gif/vgif'
mapplot cln
quit
! Quitting program
! Log file J0522-3627_S_1998_10_01_pus_uvs.log closed on Thu Feb  7 17:05:29 2008