! Started logfile: J1625+4134_S_2001_07_05_pus_uvs.log on Sun Aug 12 23:55:23 2007
obs J1625+4134_S_2001_07_05_pus_uva.fits
! Reading UV FITS file: J1625+4134_S_2001_07_05_pus_uva.fits
! AN table 1: 214 integrations on 171 of 171 possible baselines.
! AN table 2: 130 integrations on 171 of 171 possible baselines.
! Apparent sampling: 0.201533 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J1625+4134
! Warning: Unknown AN-table MNTSTA value (3).
!          Will assume that it is ground-based.
! Warning: Unknown AN-table MNTSTA value (3).
!          Will assume that it is ground-based.
! 
! 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.22997e+09        8e+06          1        8e+06
!  02        2  2.25997e+09        8e+06          1        8e+06
!  03        3  2.34997e+09        8e+06          1        8e+06
!  04        4  2.36997e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1778 lines of history.
! 
! Reading 47420 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 = 1024 
print "field_cell =",field_cell
! field_cell = 0.5 
print "taper_size =",taper_size
! taper_size = 20 
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 = 1024x1024 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 2075 telescope corrections were flagged in sub-array 1.
!  A total of 2075 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 2075 telescope corrections were flagged in sub-array 1.
!  A total of 2075 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 2075 telescope corrections were flagged in sub-array 1.
!  A total of 2075 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 2075 telescope corrections were flagged in sub-array 1.
!  A total of 2081 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.890898Jy  sigma=0.008208
! Fit after  self-cal, rms=0.492322Jy  sigma=0.004178
! 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.199 mas, bmaj=3.16 mas, bpa=-4.31 degrees
! Estimated noise=908.456 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.424344 Jy
! Component: 100  -  total flux cleaned = 0.589538 Jy
! Total flux subtracted in 100 components = 0.589538 Jy
! Clean residual min=-0.022869 max=0.079932 Jy/beam
! Clean residual mean=0.000286 rms=0.006823 Jy/beam
! Combined flux in latest and established models = 0.589538 Jy
! Performing phase self-cal
! Adding 13 model components to the UV plane model.
! The established model now contains 13 components and 0.589538 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267449Jy  sigma=0.002374
! Fit after  self-cal, rms=0.267223Jy  sigma=0.002373
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.097584 Jy
! Component: 100  -  total flux cleaned = 0.161679 Jy
! Total flux subtracted in 100 components = 0.161679 Jy
! Clean residual min=-0.015817 max=0.036347 Jy/beam
! Clean residual mean=0.000179 rms=0.004136 Jy/beam
! Combined flux in latest and established models = 0.751218 Jy
! Performing phase self-cal
! Adding 22 model components to the UV plane model.
! The established model now contains 32 components and 0.751218 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.214023Jy  sigma=0.001715
! Fit after  self-cal, rms=0.213978Jy  sigma=0.001715
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0441526 Jy
! Component: 100  -  total flux cleaned = 0.0747705 Jy
! Total flux subtracted in 100 components = 0.0747705 Jy
! Clean residual min=-0.013770 max=0.027408 Jy/beam
! Clean residual mean=0.000130 rms=0.003283 Jy/beam
! Combined flux in latest and established models = 0.825988 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 39 components and 0.825988 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.201558Jy  sigma=0.001535
! Fit after  self-cal, rms=0.201563Jy  sigma=0.001534
! Inverting map 
! Added new window around map position (-0.5, 5.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0267553 Jy
! Component: 100  -  total flux cleaned = 0.0457238 Jy
! Total flux subtracted in 100 components = 0.0457238 Jy
! Clean residual min=-0.011714 max=0.021486 Jy/beam
! Clean residual mean=0.000102 rms=0.002904 Jy/beam
! Combined flux in latest and established models = 0.871712 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 54 components and 0.871712 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.196364Jy  sigma=0.001461
! Fit after  self-cal, rms=0.196241Jy  sigma=0.001460
! Inverting map 
! Added new window around map position (-3, -0.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0199942 Jy
! Component: 100  -  total flux cleaned = 0.0340467 Jy
! Total flux subtracted in 100 components = 0.0340467 Jy
! Clean residual min=-0.011465 max=0.013654 Jy/beam
! Clean residual mean=0.000080 rms=0.002663 Jy/beam
! Combined flux in latest and established models = 0.905759 Jy
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 65 components and 0.905759 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.194015Jy  sigma=0.001427
! Fit after  self-cal, rms=0.193942Jy  sigma=0.001426
! 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.849 mas, bmaj=3.891 mas, bpa=-16.86 degrees
! Estimated noise=372.753 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 (-3.5, 7.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0189753 Jy
! Component: 100  -  total flux cleaned = 0.0313846 Jy
! Component: 150  -  total flux cleaned = 0.0407645 Jy
! Component: 200  -  total flux cleaned = 0.0480139 Jy
! Total flux subtracted in 200 components = 0.0480139 Jy
! Clean residual min=-0.009346 max=0.010074 Jy/beam
! Clean residual mean=0.000333 rms=0.002442 Jy/beam
! Combined flux in latest and established models = 0.953773 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 99 components and 0.953773 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.191886Jy  sigma=0.001398
! Fit after  self-cal, rms=0.191754Jy  sigma=0.001396
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 60 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.99     FD        1.00     GC        1.02     GG        1.00*    
!   HH        1.00*    HN        0.99     KK        0.98     KP        1.08     
!   LA        1.05     MC        1.00*    MK        0.98     NL        1.01     
!   NY        0.99     OV        0.97     PT        0.96     SC        1.05     
!   TS        1.00     WF        0.98     WZ        1.00*    
! 
!  A total of 18 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.96     FD        0.94     GC        1.00*    GG        1.00*    
!   HH        1.00*    HN        0.92     KK        1.00*    KP        0.88     
!   LA        0.91     MC        1.00*    MK        1.00*    NL        0.90     
!   NY        1.00*    OV        0.97     PT        0.87     SC        1.03     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
! 
! Correcting IF 2.
!  A total of 60 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.98     FD        1.02     GC        0.99     GG        1.00*    
!   HH        1.00*    HN        1.00     KK        0.99     KP        0.94     
!   LA        1.02     MC        1.00*    MK        0.97     NL        0.99     
!   NY        1.00     OV        1.00     PT        1.00     SC        1.08     
!   TS        0.99     WF        0.98     WZ        1.00*    
! 
!  A total of 18 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.95     FD        0.94     GC        1.00*    GG        1.00*    
!   HH        1.00*    HN        0.92     KK        1.00*    KP        0.90     
!   LA        0.86     MC        1.00*    MK        1.00*    NL        0.88     
!   NY        1.00*    OV        0.98     PT        0.89     SC        1.05     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
! 
! Correcting IF 3.
!  A total of 60 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.03     FD        0.98     GC        1.04     GG        1.00*    
!   HH        1.00*    HN        1.07     KK        1.01     KP        1.08     
!   LA        1.07     MC        1.00*    MK        1.04     NL        1.04     
!   NY        1.01     OV        1.03     PT        0.99     SC        1.10     
!   TS        1.04     WF        1.07     WZ        1.00*    
! 
!  A total of 18 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        0.93     GC        1.00*    GG        1.00*    
!   HH        1.00*    HN        1.08     KK        1.00*    KP        0.91     
!   LA        0.90     MC        1.00*    MK        1.00*    NL        0.94     
!   NY        1.00*    OV        1.05     PT        0.90     SC        1.05     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
! 
! Correcting IF 4.
!  A total of 60 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.98     FD        1.03     GC        1.04     GG        1.00*    
!   HH        1.00*    HN        0.97     KK        0.99     KP        1.08     
!   LA        1.04     MC        1.00*    MK        1.03     NL        1.03     
!   NY        1.02     OV        1.01     PT        0.98     SC        1.10     
!   TS        1.04     WF        1.02     WZ        1.00*    
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.97     FD        0.95     GC        1.00*    GG        1.00*    
!   HH        1.00*    HN        0.92     KK        1.00*    KP        0.97     
!   LA        0.94     MC        1.00*    MK        1.00*    NL        0.92     
!   NY        1.00*    OV        0.97     PT        1.04     SC        1.05     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
! 
! Fit before self-cal, rms=0.191754Jy  sigma=0.001396
! Fit after  self-cal, rms=0.187120Jy  sigma=0.001289
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.879 mas, bmaj=3.927 mas, bpa=-17 degrees
! Estimated noise=374.649 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 (-4, 27).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00882898 Jy
! Component: 100  -  total flux cleaned = 0.0141647 Jy
! Component: 150  -  total flux cleaned = 0.0181159 Jy
! Component: 200  -  total flux cleaned = 0.0204886 Jy
! Total flux subtracted in 200 components = 0.0204886 Jy
! Clean residual min=-0.002819 max=0.004361 Jy/beam
! Clean residual mean=0.000306 rms=0.000829 Jy/beam
! Combined flux in latest and established models = 0.974261 Jy
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 141 components and 0.974261 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186649Jy  sigma=0.001283
! Fit after  self-cal, rms=0.186687Jy  sigma=0.001283
! Inverting map 
! Added new window around map position (-25, -50).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0039565 Jy
! Component: 100  -  total flux cleaned = 0.00606448 Jy
! Component: 150  -  total flux cleaned = 0.00730073 Jy
! Component: 200  -  total flux cleaned = 0.00814467 Jy
! Total flux subtracted in 200 components = 0.00814467 Jy
! Clean residual min=-0.002694 max=0.003683 Jy/beam
! Clean residual mean=0.000279 rms=0.000760 Jy/beam
! Combined flux in latest and established models = 0.982406 Jy
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 158 components and 0.982406 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186534Jy  sigma=0.001281
! Fit after  self-cal, rms=0.186525Jy  sigma=0.001281
! 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.186525Jy  sigma=0.001281
! Fit after  self-cal, rms=0.186515Jy  sigma=0.001281
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.878 mas, bmaj=3.924 mas, bpa=-16.89 degrees
! Estimated noise=374.627 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00242069 Jy
! Component: 100  -  total flux cleaned = 0.00315789 Jy
! Component: 150  -  total flux cleaned = 0.00373953 Jy
! Component: 200  -  total flux cleaned = 0.00433228 Jy
! Total flux subtracted in 200 components = 0.00433228 Jy
! Clean residual min=-0.002608 max=0.003464 Jy/beam
! Clean residual mean=0.000298 rms=0.000706 Jy/beam
! Combined flux in latest and established models = 0.986738 Jy
selfcal
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 173 components and 0.986738 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186446Jy  sigma=0.001280
! Fit after  self-cal, rms=0.186442Jy  sigma=0.001280
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)
! Added new window around map position (4, -28).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00290085 Jy
! Component: 100  -  total flux cleaned = 0.00421801 Jy
! Component: 150  -  total flux cleaned = 0.00471395 Jy
! Component: 200  -  total flux cleaned = 0.00529352 Jy
! Total flux subtracted in 200 components = 0.00529352 Jy
! Clean residual min=-0.002569 max=0.003195 Jy/beam
! Clean residual mean=0.000278 rms=0.000685 Jy/beam
! Combined flux in latest and established models = 0.992032 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 193 components and 0.992032 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186368Jy  sigma=0.001280
! Fit after  self-cal, rms=0.186365Jy  sigma=0.001280
! Inverting map 
! Added new window around map position (-2.5, 37.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00293001 Jy
! Component: 100  -  total flux cleaned = 0.00466275 Jy
! Component: 150  -  total flux cleaned = 0.00604427 Jy
! Component: 200  -  total flux cleaned = 0.00693382 Jy
! Total flux subtracted in 200 components = 0.00693382 Jy
! Clean residual min=-0.002451 max=0.002681 Jy/beam
! Clean residual mean=0.000255 rms=0.000649 Jy/beam
! Combined flux in latest and established models = 0.998966 Jy
! Performing phase self-cal
! Adding 36 model components to the UV plane model.
! The established model now contains 214 components and 0.998966 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186287Jy  sigma=0.001279
! Fit after  self-cal, rms=0.186283Jy  sigma=0.001279
! 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.186283Jy  sigma=0.001279
! Fit after  self-cal, rms=0.186345Jy  sigma=0.001277
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.878 mas, bmaj=3.926 mas, bpa=-16.88 degrees
! Estimated noise=374.839 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00153837 Jy
! Component: 100  -  total flux cleaned = 0.00265406 Jy
! Component: 150  -  total flux cleaned = 0.00343073 Jy
! Component: 200  -  total flux cleaned = 0.0038851 Jy
! Total flux subtracted in 200 components = 0.0038851 Jy
! Clean residual min=-0.002473 max=0.002604 Jy/beam
! Clean residual mean=0.000255 rms=0.000634 Jy/beam
! Combined flux in latest and established models = 1.00285 Jy
selfcal
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 226 components and 1.00285 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186310Jy  sigma=0.001277
! Fit after  self-cal, rms=0.186308Jy  sigma=0.001277
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.000525137 Jy
! Component: 100  -  total flux cleaned = 0.000952696 Jy
! Component: 150  -  total flux cleaned = 0.00127649 Jy
! Component: 200  -  total flux cleaned = 0.00147115 Jy
! Total flux subtracted in 200 components = 0.00147115 Jy
! Clean residual min=-0.002499 max=0.002561 Jy/beam
! Clean residual mean=0.000250 rms=0.000628 Jy/beam
! Combined flux in latest and established models = 1.00432 Jy
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 237 components and 1.00432 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186288Jy  sigma=0.001277
! Fit after  self-cal, rms=0.186286Jy  sigma=0.001277
! 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.186286Jy  sigma=0.001277
! Fit after  self-cal, rms=0.186244Jy  sigma=0.001276
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.363 mas, bmaj=8.963 mas, bpa=-28.68 degrees
! Estimated noise=495.542 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00280898 Jy
! Component: 100  -  total flux cleaned = 0.00470197 Jy
! Component: 150  -  total flux cleaned = 0.00619084 Jy
! Component: 200  -  total flux cleaned = 0.00726705 Jy
! Total flux subtracted in 200 components = 0.00726705 Jy
! Clean residual min=-0.002364 max=0.004815 Jy/beam
! Clean residual mean=0.000667 rms=0.001014 Jy/beam
! Combined flux in latest and established models = 1.01159 Jy
selfcal
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 250 components and 1.01159 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186228Jy  sigma=0.001275
! Fit after  self-cal, rms=0.186224Jy  sigma=0.001275
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)
! Added new window around map position (8.5, 54.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0048663 Jy
! Component: 100  -  total flux cleaned = 0.0079833 Jy
! Component: 150  -  total flux cleaned = 0.0103027 Jy
! Component: 200  -  total flux cleaned = 0.0122285 Jy
! Total flux subtracted in 200 components = 0.0122285 Jy
! Clean residual min=-0.002085 max=0.003333 Jy/beam
! Clean residual mean=0.000551 rms=0.000808 Jy/beam
! Combined flux in latest and established models = 1.02382 Jy
! Performing phase self-cal
! Adding 73 model components to the UV plane model.
! The established model now contains 319 components and 1.02382 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.186140Jy  sigma=0.001275
! Fit after  self-cal, rms=0.186140Jy  sigma=0.001275
! 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.878 mas, bmaj=3.927 mas, bpa=-16.88 degrees
! Estimated noise=374.936 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.878 x 3.927 at -16.88 degrees (North through East)
! Clean map  min=-0.0023025  max=0.59658 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.186140Jy  sigma=0.001275
! Fit after  self-cal, rms=0.276068Jy  sigma=0.001172
! Inverting map and beam 
! Estimated beam: bmin=2.876 mas, bmaj=3.919 mas, bpa=-16.94 degrees
! Estimated noise=374.53 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00118567 Jy
! Component: 100  -  total flux cleaned = 0.00201454 Jy
! Component: 150  -  total flux cleaned = 0.00244017 Jy
! Component: 200  -  total flux cleaned = 0.00278661 Jy
! Total flux subtracted in 200 components = 0.00278661 Jy
! Clean residual min=-0.001884 max=0.002368 Jy/beam
! Clean residual mean=0.000191 rms=0.000485 Jy/beam
! Combined flux in latest and established models = 1.0266 Jy
! Performing phase self-cal
! Adding 35 model components to the UV plane model.
! The established model now contains 350 components and 1.0266 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276041Jy  sigma=0.001172
! Fit after  self-cal, rms=0.276049Jy  sigma=0.001172
! Inverting map 
! Added new window around map position (127.5, 95.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00200631 Jy
! Component: 100  -  total flux cleaned = 0.00277639 Jy
! Component: 150  -  total flux cleaned = 0.003155 Jy
! Component: 200  -  total flux cleaned = 0.00355784 Jy
! Total flux subtracted in 200 components = 0.00355784 Jy
! Clean residual min=-0.001899 max=0.002218 Jy/beam
! Clean residual mean=0.000183 rms=0.000469 Jy/beam
! Combined flux in latest and established models = 1.03016 Jy
! Performing phase self-cal
! Adding 43 model components to the UV plane model.
! The established model now contains 368 components and 1.03016 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276019Jy  sigma=0.001171
! Fit after  self-cal, rms=0.276024Jy  sigma=0.001171
! Inverting map 
! Added new window around map position (1, 16).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00166585 Jy
! Component: 100  -  total flux cleaned = 0.00181991 Jy
! Component: 150  -  total flux cleaned = 0.00223698 Jy
! Component: 200  -  total flux cleaned = 0.00280276 Jy
! Total flux subtracted in 200 components = 0.00280276 Jy
! Clean residual min=-0.001823 max=0.002073 Jy/beam
! Clean residual mean=0.000173 rms=0.000449 Jy/beam
! Combined flux in latest and established models = 1.03296 Jy
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 389 components and 1.03296 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276004Jy  sigma=0.001171
! Fit after  self-cal, rms=0.276008Jy  sigma=0.001171
! Inverting map 
! Added new window around map position (-72, -75).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0018305 Jy
! Component: 100  -  total flux cleaned = 0.0027964 Jy
! Component: 150  -  total flux cleaned = 0.00323534 Jy
! Component: 200  -  total flux cleaned = 0.00344212 Jy
! Total flux subtracted in 200 components = 0.00344212 Jy
! Clean residual min=-0.001772 max=0.002019 Jy/beam
! Clean residual mean=0.000164 rms=0.000435 Jy/beam
! Combined flux in latest and established models = 1.03641 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 408 components and 1.03641 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.275985Jy  sigma=0.001171
! Fit after  self-cal, rms=0.275992Jy  sigma=0.001171
! Inverting map 
! Added new window around map position (-64, -9.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00124669 Jy
! Component: 100  -  total flux cleaned = 0.00194073 Jy
! Component: 150  -  total flux cleaned = 0.00246733 Jy
! Component: 200  -  total flux cleaned = 0.00289363 Jy
! Total flux subtracted in 200 components = 0.00289363 Jy
! Clean residual min=-0.001763 max=0.001961 Jy/beam
! Clean residual mean=0.000155 rms=0.000421 Jy/beam
! Combined flux in latest and established models = 1.0393 Jy
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 428 components and 1.0393 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.275975Jy  sigma=0.001170
! Fit after  self-cal, rms=0.275979Jy  sigma=0.001170
! Inverting map 
! Added new window around map position (-47.5, 102.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00193993 Jy
! Component: 100  -  total flux cleaned = 0.00286797 Jy
! Component: 150  -  total flux cleaned = 0.00355951 Jy
! Component: 200  -  total flux cleaned = 0.00409953 Jy
! Total flux subtracted in 200 components = 0.00409953 Jy
! Clean residual min=-0.001634 max=0.001789 Jy/beam
! Clean residual mean=0.000142 rms=0.000407 Jy/beam
! Combined flux in latest and established models = 1.0434 Jy
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 451 components and 1.0434 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.275956Jy  sigma=0.001170
! Fit after  self-cal, rms=0.275956Jy  sigma=0.001170
! Inverting map 
! Added new window around map position (-17.5, -48.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00126614 Jy
! Component: 100  -  total flux cleaned = 0.00171888 Jy
! Component: 150  -  total flux cleaned = 0.0022243 Jy
! Component: 200  -  total flux cleaned = 0.00278326 Jy
! Total flux subtracted in 200 components = 0.00278326 Jy
! Clean residual min=-0.001592 max=0.001640 Jy/beam
! Clean residual mean=0.000134 rms=0.000396 Jy/beam
! Combined flux in latest and established models = 1.04618 Jy
! Performing phase self-cal
! Adding 47 model components to the UV plane model.
! The established model now contains 474 components and 1.04618 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.275937Jy  sigma=0.001170
! Fit after  self-cal, rms=0.275937Jy  sigma=0.001170
! 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.275937Jy  sigma=0.001170
! Fit after  self-cal, rms=0.283633Jy  sigma=0.001170
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.218 mas, bmaj=3.172 mas, bpa=-4.179 degrees
! Estimated noise=818.683 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.42341 Jy
! Component: 100  -  total flux cleaned = 0.586659 Jy
! Component: 150  -  total flux cleaned = 0.680287 Jy
! Component: 200  -  total flux cleaned = 0.741076 Jy
! Total flux subtracted in 200 components = 0.741076 Jy
! Clean residual min=-0.010784 max=0.033028 Jy/beam
! Clean residual mean=0.000186 rms=0.003400 Jy/beam
! Combined flux in latest and established models = 0.741076 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 29 model components to the UV plane model.
! The established model now contains 29 components and 0.741076 Jy
! Inverting map and beam 
! Estimated beam: bmin=2.876 mas, bmaj=3.919 mas, bpa=-16.96 degrees
! Estimated noise=374.533 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.0679842 Jy
! Component: 100  -  total flux cleaned = 0.108903 Jy
! Component: 150  -  total flux cleaned = 0.135919 Jy
! Component: 200  -  total flux cleaned = 0.155594 Jy
! Total flux subtracted in 200 components = 0.155594 Jy
! Clean residual min=-0.005593 max=0.011771 Jy/beam
! Clean residual mean=0.000679 rms=0.002211 Jy/beam
! Combined flux in latest and established models = 0.896671 Jy
! Adding 75 model components to the UV plane model.
! The established model now contains 96 components and 0.896671 Jy
! Inverting map 
! Added new window around map position (-0.5, -4).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0153996 Jy
! Component: 100  -  total flux cleaned = 0.0279384 Jy
! Component: 150  -  total flux cleaned = 0.0382819 Jy
! Component: 200  -  total flux cleaned = 0.0469336 Jy
! Total flux subtracted in 200 components = 0.0469336 Jy
! Clean residual min=-0.003615 max=0.006102 Jy/beam
! Clean residual mean=0.000496 rms=0.001220 Jy/beam
! Combined flux in latest and established models = 0.943604 Jy
! Adding 73 model components to the UV plane model.
! The established model now contains 155 components and 0.943604 Jy
! Inverting map 
! Added new window around map position (-4.5, 16).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00751371 Jy
! Component: 100  -  total flux cleaned = 0.0138719 Jy
! Component: 150  -  total flux cleaned = 0.0193654 Jy
! Component: 200  -  total flux cleaned = 0.0241501 Jy
! Total flux subtracted in 200 components = 0.0241501 Jy
! Clean residual min=-0.002692 max=0.003378 Jy/beam
! Clean residual mean=0.000405 rms=0.000821 Jy/beam
! Combined flux in latest and established models = 0.967754 Jy
! Adding 60 model components to the UV plane model.
! The established model now contains 195 components and 0.967755 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.36 mas, bmaj=8.958 mas, bpa=-28.58 degrees
! Estimated noise=496.057 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00934706 Jy
! Component: 100  -  total flux cleaned = 0.0161815 Jy
! Component: 150  -  total flux cleaned = 0.0217422 Jy
! Component: 200  -  total flux cleaned = 0.0264179 Jy
! Total flux subtracted in 200 components = 0.0264179 Jy
! Clean residual min=-0.002116 max=0.003785 Jy/beam
! Clean residual mean=0.000897 rms=0.000902 Jy/beam
! Combined flux in latest and established models = 0.994173 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.00405781 Jy
! Component: 100  -  total flux cleaned = 0.00766891 Jy
! Component: 150  -  total flux cleaned = 0.0109553 Jy
! Component: 200  -  total flux cleaned = 0.0139974 Jy
! Total flux subtracted in 200 components = 0.0139974 Jy
! Clean residual min=-0.001701 max=0.003060 Jy/beam
! Clean residual mean=0.000759 rms=0.000661 Jy/beam
! Combined flux in latest and established models = 1.00817 Jy
! Adding 199 model components to the UV plane model.
! The established model now contains 371 components and 1.00817 Jy
! Inverting map 
! Added new window around map position (39.5, -15).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0031772 Jy
! Component: 100  -  total flux cleaned = 0.00589127 Jy
! Component: 150  -  total flux cleaned = 0.00843592 Jy
! Component: 200  -  total flux cleaned = 0.0108317 Jy
! Total flux subtracted in 200 components = 0.0108317 Jy
! Clean residual min=-0.001387 max=0.002608 Jy/beam
! Clean residual mean=0.000657 rms=0.000523 Jy/beam
! Combined flux in latest and established models = 1.019 Jy
! Adding 79 model components to the UV plane model.
! The established model now contains 436 components and 1.019 Jy
! Inverting map 
! Added new window around map position (2.5, -124).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00266639 Jy
! Component: 100  -  total flux cleaned = 0.00486262 Jy
! Component: 150  -  total flux cleaned = 0.00691824 Jy
! Component: 200  -  total flux cleaned = 0.00886378 Jy
! Total flux subtracted in 200 components = 0.00886378 Jy
! Clean residual min=-0.001162 max=0.002085 Jy/beam
! Clean residual mean=0.000574 rms=0.000441 Jy/beam
! Combined flux in latest and established models = 1.02786 Jy
! Adding 79 model components to the UV plane model.
! The established model now contains 496 components and 1.02787 Jy
! Inverting map 
! Added new window around map position (-21.5, 90.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00213544 Jy
! Component: 100  -  total flux cleaned = 0.00391627 Jy
! Component: 150  -  total flux cleaned = 0.00559562 Jy
! Component: 200  -  total flux cleaned = 0.00719065 Jy
! Total flux subtracted in 200 components = 0.00719065 Jy
! Clean residual min=-0.000944 max=0.001837 Jy/beam
! Clean residual mean=0.000506 rms=0.000389 Jy/beam
! Combined flux in latest and established models = 1.03506 Jy
! Adding 61 model components to the UV plane model.
! The established model now contains 542 components and 1.03506 Jy
! Inverting map 
! Added new window around map position (49, -86).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00182446 Jy
! Component: 100  -  total flux cleaned = 0.00333305 Jy
! Component: 150  -  total flux cleaned = 0.00476975 Jy
! Component: 200  -  total flux cleaned = 0.00614313 Jy
! Total flux subtracted in 200 components = 0.00614313 Jy
! Clean residual min=-0.000870 max=0.001761 Jy/beam
! Clean residual mean=0.000448 rms=0.000365 Jy/beam
! Combined flux in latest and established models = 1.0412 Jy
! Adding 64 model components to the UV plane model.
! The established model now contains 583 components and 1.0412 Jy
! Inverting map 
! Added new window around map position (36.5, 108).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00178936 Jy
! Component: 100  -  total flux cleaned = 0.00316674 Jy
! Component: 150  -  total flux cleaned = 0.0044628 Jy
! Component: 200  -  total flux cleaned = 0.0056989 Jy
! Total flux subtracted in 200 components = 0.0056989 Jy
! Clean residual min=-0.000922 max=0.001642 Jy/beam
! Clean residual mean=0.000393 rms=0.000341 Jy/beam
! Combined flux in latest and established models = 1.0469 Jy
! Adding 58 model components to the UV plane model.
! The established model now contains 619 components and 1.0469 Jy
! Inverting map 
! Added new window around map position (88.5, 127).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00151181 Jy
! Component: 100  -  total flux cleaned = 0.00270598 Jy
! Component: 150  -  total flux cleaned = 0.00384004 Jy
! Component: 200  -  total flux cleaned = 0.00492207 Jy
! Total flux subtracted in 200 components = 0.00492207 Jy
! Clean residual min=-0.000837 max=0.001427 Jy/beam
! Clean residual mean=0.000348 rms=0.000314 Jy/beam
! Combined flux in latest and established models = 1.05182 Jy
! Adding 68 model components to the UV plane model.
! The established model now contains 670 components and 1.05182 Jy
! Inverting map 
! Added new window around map position (5, 104.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00133426 Jy
! Component: 100  -  total flux cleaned = 0.00239802 Jy
! Component: 150  -  total flux cleaned = 0.00341146 Jy
! Component: 200  -  total flux cleaned = 0.00438672 Jy
! Total flux subtracted in 200 components = 0.00438672 Jy
! Clean residual min=-0.000877 max=0.001297 Jy/beam
! Clean residual mean=0.000308 rms=0.000297 Jy/beam
! Combined flux in latest and established models = 1.05621 Jy
! Adding 62 model components to the UV plane model.
! The established model now contains 710 components and 1.05621 Jy
! Inverting map 
! Added new window around map position (-112, -59.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0012347 Jy
! Component: 100  -  total flux cleaned = 0.00219405 Jy
! Component: 150  -  total flux cleaned = 0.00310908 Jy
! Component: 200  -  total flux cleaned = 0.00399018 Jy
! Total flux subtracted in 200 components = 0.00399018 Jy
! Clean residual min=-0.000827 max=0.001216 Jy/beam
! Clean residual mean=0.000276 rms=0.000279 Jy/beam
! Combined flux in latest and established models = 1.0602 Jy
! Adding 58 model components to the UV plane model.
! The established model now contains 745 components and 1.0602 Jy
! Inverting map 
! Added new window around map position (35.5, -54).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00111674 Jy
! Component: 100  -  total flux cleaned = 0.00199129 Jy
! Component: 150  -  total flux cleaned = 0.00282712 Jy
! Component: 200  -  total flux cleaned = 0.00363246 Jy
! Total flux subtracted in 200 components = 0.00363246 Jy
! Clean residual min=-0.000900 max=0.001147 Jy/beam
! Clean residual mean=0.000243 rms=0.000267 Jy/beam
! Combined flux in latest and established models = 1.06383 Jy
! Adding 63 model components to the UV plane model.
! The established model now contains 788 components and 1.06383 Jy
! Inverting map 
! Added new window around map position (89, -122.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00127724 Jy
! Component: 100  -  total flux cleaned = 0.00215922 Jy
! Component: 150  -  total flux cleaned = 0.0029713 Jy
! Component: 200  -  total flux cleaned = 0.00374524 Jy
! Total flux subtracted in 200 components = 0.00374524 Jy
! Clean residual min=-0.000956 max=0.001000 Jy/beam
! Clean residual mean=0.000210 rms=0.000251 Jy/beam
! Combined flux in latest and established models = 1.06757 Jy
! Adding 58 model components to the UV plane model.
! The established model now contains 825 components and 1.06757 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=2.876 mas, bmaj=3.919 mas, bpa=-16.96 degrees
! Estimated noise=374.533 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000646793 Jy
! Component: 100  -  total flux cleaned = 0.00146774 Jy
! Component: 150  -  total flux cleaned = 0.00198001 Jy
! Component: 200  -  total flux cleaned = 0.00264973 Jy
! Total flux subtracted in 200 components = 0.00264973 Jy
! Clean residual min=-0.001377 max=0.001719 Jy/beam
! Clean residual mean=0.000064 rms=0.000364 Jy/beam
! Combined flux in latest and established models = 1.07022 Jy
selfcal
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 868 components and 1.07022 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.283639Jy  sigma=0.001170
! Fit after  self-cal, rms=0.283607Jy  sigma=0.001170
wmodel J1625+4134_S_2001_07_05_pus_map.mod
! Writing 868 model components to file: J1625+4134_S_2001_07_05_pus_map.mod
wobs J1625+4134_S_2001_07_05_pus_uvs.fits
! Writing UV FITS file: J1625+4134_S_2001_07_05_pus_uvs.fits
wwins J1625+4134_S_2001_07_05_pus_map.win
! wwins: Wrote 27 windows to J1625+4134_S_2001_07_05_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=4096  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.000678458 Jy
! Component: 100  -  total flux cleaned = 0.000971781 Jy
! Component: 150  -  total flux cleaned = 0.00110033 Jy
! Component: 200  -  total flux cleaned = 0.00117568 Jy
! Component: 250  -  total flux cleaned = 0.00129895 Jy
! Component: 300  -  total flux cleaned = 0.0013958 Jy
! Component: 350  -  total flux cleaned = 0.0015629 Jy
! Component: 400  -  total flux cleaned = 0.00175156 Jy
! Component: 450  -  total flux cleaned = 0.00191499 Jy
! Component: 500  -  total flux cleaned = 0.00214591 Jy
! Component: 550  -  total flux cleaned = 0.00237389 Jy
! Component: 600  -  total flux cleaned = 0.00248683 Jy
! Component: 650  -  total flux cleaned = 0.00264392 Jy
! Component: 700  -  total flux cleaned = 0.00279917 Jy
! Component: 750  -  total flux cleaned = 0.00297508 Jy
! Component: 800  -  total flux cleaned = 0.00308421 Jy
! Component: 850  -  total flux cleaned = 0.00325715 Jy
! Component: 900  -  total flux cleaned = 0.00342877 Jy
! Component: 950  -  total flux cleaned = 0.00353517 Jy
! Component: 1000  -  total flux cleaned = 0.00372535 Jy
! Component: 1050  -  total flux cleaned = 0.00391373 Jy
! Component: 1100  -  total flux cleaned = 0.00401777 Jy
! Component: 1150  -  total flux cleaned = 0.00424507 Jy
! Component: 1200  -  total flux cleaned = 0.00444989 Jy
! Component: 1250  -  total flux cleaned = 0.00457195 Jy
! Component: 1300  -  total flux cleaned = 0.00473389 Jy
! Component: 1350  -  total flux cleaned = 0.00495513 Jy
! Component: 1400  -  total flux cleaned = 0.00517467 Jy
! Component: 1450  -  total flux cleaned = 0.0053334 Jy
! Component: 1500  -  total flux cleaned = 0.0054716 Jy
! Component: 1550  -  total flux cleaned = 0.00570695 Jy
! Component: 1600  -  total flux cleaned = 0.00586291 Jy
! Component: 1650  -  total flux cleaned = 0.00597936 Jy
! Component: 1700  -  total flux cleaned = 0.0061336 Jy
! Component: 1750  -  total flux cleaned = 0.0063637 Jy
! Component: 1800  -  total flux cleaned = 0.00647816 Jy
! Component: 1850  -  total flux cleaned = 0.0066679 Jy
! Component: 1900  -  total flux cleaned = 0.00679981 Jy
! Component: 1950  -  total flux cleaned = 0.00691228 Jy
! Component: 2000  -  total flux cleaned = 0.00708027 Jy
! Component: 2050  -  total flux cleaned = 0.00715437 Jy
! Component: 2100  -  total flux cleaned = 0.00730219 Jy
! Component: 2150  -  total flux cleaned = 0.00743081 Jy
! Component: 2200  -  total flux cleaned = 0.0075771 Jy
! Component: 2250  -  total flux cleaned = 0.00774074 Jy
! Component: 2300  -  total flux cleaned = 0.00779502 Jy
! Component: 2350  -  total flux cleaned = 0.00788523 Jy
! Component: 2400  -  total flux cleaned = 0.00802876 Jy
! Component: 2450  -  total flux cleaned = 0.00813592 Jy
! Component: 2500  -  total flux cleaned = 0.00827812 Jy
! Component: 2550  -  total flux cleaned = 0.00838421 Jy
! Component: 2600  -  total flux cleaned = 0.00847214 Jy
! Component: 2650  -  total flux cleaned = 0.00864759 Jy
! Component: 2700  -  total flux cleaned = 0.00875229 Jy
! Component: 2750  -  total flux cleaned = 0.00889134 Jy
! Component: 2800  -  total flux cleaned = 0.00899523 Jy
! Component: 2850  -  total flux cleaned = 0.00911593 Jy
! Component: 2900  -  total flux cleaned = 0.00920179 Jy
! Component: 2950  -  total flux cleaned = 0.00932153 Jy
! Component: 3000  -  total flux cleaned = 0.00949171 Jy
! Component: 3050  -  total flux cleaned = 0.00952555 Jy
! Component: 3100  -  total flux cleaned = 0.0096438 Jy
! Component: 3150  -  total flux cleaned = 0.00972795 Jy
! Component: 3200  -  total flux cleaned = 0.009862 Jy
! Component: 3250  -  total flux cleaned = 0.00996215 Jy
! Component: 3300  -  total flux cleaned = 0.0100453 Jy
! Component: 3350  -  total flux cleaned = 0.0101778 Jy
! Component: 3400  -  total flux cleaned = 0.0102109 Jy
! Component: 3450  -  total flux cleaned = 0.0103095 Jy
! Component: 3500  -  total flux cleaned = 0.010457 Jy
! Component: 3550  -  total flux cleaned = 0.0105385 Jy
! Component: 3600  -  total flux cleaned = 0.0106199 Jy
! Component: 3650  -  total flux cleaned = 0.0107171 Jy
! Component: 3700  -  total flux cleaned = 0.0107656 Jy
! Component: 3750  -  total flux cleaned = 0.0109264 Jy
! Component: 3800  -  total flux cleaned = 0.0110228 Jy
! Component: 3850  -  total flux cleaned = 0.0110868 Jy
! Component: 3900  -  total flux cleaned = 0.0111664 Jy
! Component: 3950  -  total flux cleaned = 0.0113092 Jy
! Component: 4000  -  total flux cleaned = 0.0113882 Jy
! Component: 4050  -  total flux cleaned = 0.0114512 Jy
! Total flux subtracted in 4096 components = 0.0115454 Jy
! Clean residual min=-0.001293 max=0.000905 Jy/beam
! Clean residual mean=0.000032 rms=0.000270 Jy/beam
! Combined flux in latest and established models = 1.08177 Jy
keep
! Adding 294 model components to the UV plane model.
! The established model now contains 1152 components and 1.08177 Jy
wmap "J1625+4134_S_2001_07_05_pus\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.876 x 3.919 at -16.96 degrees (North through East)
! Clean map  min=-0.0012746  max=0.59679 Jy/beam
! Writing clean map to FITS file: J1625+4134_S_2001_07_05_pus_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 0.596794 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000269591 
device "J1625+4134_S_2001_07_05_pus\_map.eps/VPS"
! Attempting to open device: 'J1625+4134_S_2001_07_05_pus_map.eps/VPS'
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.13552 0.13552 0.271039 0.542079 1.08416 2.16831 4.33663 8.67326 17.3465 34.693 69.3861
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0008 
print "peak=",int(1000*peak_flux)/1000
! peak= 0.596 
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.876 x 3.919 at -16.96 degrees (North through East)
! Clean map  min=-0.0012746  max=0.59679 Jy/beam
mapcolor color, 0.9
! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5.
device "J1625+4134_S_2001_07_05_pus\_map.gif/vgif"
! Attempting to open device: 'J1625+4134_S_2001_07_05_pus_map.gif/vgif'
mapplot cln
quit
! Quitting program
! Log file J1625+4134_S_2001_07_05_pus_uvs.log closed on Mon Aug 13 00:02:55 2007