! Started logfile: /scr/pima/rdv62_uvs/J2056-4714_X_uvs.log on Fri Nov 22 22:41:57 2013
obs /scr/pima/rdv62_uvs/J2056-4714_X_uva.fits
! Reading UV FITS file: /scr/pima/rdv62_uvs/J2056-4714_X_uva.fits
! AN table 1: 7 integrations on 15 of 15 possible baselines.
! Apparent sampling: 0.866667 visibilities/baseline/integration-bin.
! Found source: J2056-4714
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  8.40599e+09        8e+06          1        8e+06
!  02        2  8.47599e+09        8e+06          1        8e+06
!  03        3  8.79099e+09        8e+06          1        8e+06
!  04        4  8.89599e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 0 lines of history.
! 
! Reading 364 visibilities.
select RR
! 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.2 
print "taper_size =",taper_size
! taper_size = 70 
![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J2056-4714_X]
integer clean_niter
clean_niter = 100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float soltime1
soltime1 = 800000
float soltime2
soltime2 = 15
float soltime3
soltime3 = 1
float thresh
thresh = 0.5
float win_mult
win_mult = 1.8
float old_peak
float new_peak
float flux_cutoff
selfant "",false
! 1:FD       self-cal  status=correctable  weight=1
! 1:KK       self-cal  status=correctable  weight=1
! 1:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:MK       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
system "rm /tmp/ff.ff"
mapsize field_size, field_cell
! Map grid = 1024x1024 pixels with 0.200x0.200 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 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=2.373533Jy  sigma=46.815540
! Fit after  self-cal, rms=1.190470Jy  sigma=23.947434
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
print "=== Flag out IVS stations"
! === Flag out IVS stations 
integer lun
string  cfg_file
lun = outfile ( "/tmp/ff.ff" ) 
do i=1,nsub
  do j=1,ntel(i)
    fprintf lun, "%s\n", telname(j,i)
    end do
  end do
close ( lun )
system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff"
![@/tmp/ff.ff \]
flag 1:KK
! Flagging baselines of antenna 1:KK in the currently selected channels.
![Exited script file: /tmp/ff.ff]
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=0.7053 mas, bmaj=5.607 mas, bpa=-16.03 degrees
! Estimated noise=3.1048 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 = 1.73313 Jy
! Component: 100  -  total flux cleaned = 2.11224 Jy
! Total flux subtracted in 100 components = 2.11224 Jy
! Clean residual min=-0.103313 max=0.118902 Jy/beam
! Clean residual mean=0.000488 rms=0.033709 Jy/beam
! Combined flux in latest and established models = 2.11224 Jy
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 6 components and 2.11224 Jy
! 
! Correcting IF 1.
!  A total of 7 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 7 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 7 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 7 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.311739Jy  sigma=6.857709
! Fit after  self-cal, rms=0.311502Jy  sigma=6.852530
! Inverting map and beam 
! Estimated beam: bmin=0.7053 mas, bmaj=5.607 mas, bpa=-16.03 degrees
! Estimated noise=3.1048 mJy/beam.
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=0.9275 mas, bmaj=7.184 mas, bpa=-16.33 degrees
! Estimated noise=2.20753 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.186511 Jy
! Component: 100  -  total flux cleaned = 0.256674 Jy
! Component: 150  -  total flux cleaned = 0.286195 Jy
! Component: 200  -  total flux cleaned = 0.289737 Jy
! Total flux subtracted in 200 components = 0.289737 Jy
! Clean residual min=-0.086926 max=0.071743 Jy/beam
! Clean residual mean=-0.000500 rms=0.023990 Jy/beam
! Combined flux in latest and established models = 2.40198 Jy
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 37 components and 2.40198 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.208806Jy  sigma=4.400938
! Fit after  self-cal, rms=0.208766Jy  sigma=4.398764
! Inverting map 
print "=== Finished natural weighting clean for VLBI sites ==="
! === Finished natural weighting clean for VLBI sites === 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.99     KK        1.00*    KP        0.97     LA        0.96     
!   MK        0.97     PT        1.01     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.98     KK        1.00*    KP        0.92     LA        0.96     
!   MK        1.00     PT        1.07     
! 
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.03     KK        1.00*    KP        0.87     LA        1.04     
!   MK        0.99     PT        1.11     
! 
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.05     KK        1.00*    KP        0.93     LA        1.02     
!   MK        1.03     PT        1.10     
! 
! 
! Fit before self-cal, rms=0.208766Jy  sigma=4.398764
! Fit after  self-cal, rms=0.104572Jy  sigma=2.098081
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.9496 mas, bmaj=7.3 mas, bpa=-16.97 degrees
! Estimated noise=2.24005 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.0053178 Jy
! Component: 100  -  total flux cleaned = 0.00592058 Jy
! Component: 150  -  total flux cleaned = 0.00619093 Jy
! Component: 200  -  total flux cleaned = 0.00645807 Jy
! Total flux subtracted in 200 components = 0.00645807 Jy
! Clean residual min=-0.029257 max=0.021339 Jy/beam
! Clean residual mean=0.000370 rms=0.006838 Jy/beam
! Combined flux in latest and established models = 2.40843 Jy
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 45 components and 2.40843 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.102875Jy  sigma=2.060966
! Fit after  self-cal, rms=0.102960Jy  sigma=2.058726
! Inverting map 
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.102960Jy  sigma=2.058726
! Fit after  self-cal, rms=0.102686Jy  sigma=2.055273
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.948 mas, bmaj=7.29 mas, bpa=-16.97 degrees
! Estimated noise=2.23938 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00469465 Jy
! Component: 100  -  total flux cleaned = 0.00494573 Jy
! Component: 150  -  total flux cleaned = 0.00517682 Jy
! Component: 200  -  total flux cleaned = 0.00540101 Jy
! Total flux subtracted in 200 components = 0.00540101 Jy
! Clean residual min=-0.027541 max=0.020880 Jy/beam
! Clean residual mean=0.000364 rms=0.006507 Jy/beam
! Combined flux in latest and established models = 2.41384 Jy
selfcal
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 48 components and 2.41384 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101489Jy  sigma=2.029333
! Fit after  self-cal, rms=0.101473Jy  sigma=2.028193
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.000336304 Jy
! Component: 100  -  total flux cleaned = 0.000542834 Jy
! Component: 150  -  total flux cleaned = 0.000745914 Jy
! Component: 200  -  total flux cleaned = 0.00112835 Jy
! Total flux subtracted in 200 components = 0.00112835 Jy
! Clean residual min=-0.026725 max=0.020398 Jy/beam
! Clean residual mean=0.000356 rms=0.006330 Jy/beam
! Combined flux in latest and established models = 2.41496 Jy
! Performing phase self-cal
! Adding 11 model components to the UV plane model.
! The established model now contains 51 components and 2.41496 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.100766Jy  sigma=2.012311
! Fit after  self-cal, rms=0.100770Jy  sigma=2.011759
! Inverting map 
selfcal true, true, soltime2
! Performing amp+phase self-cal over 15 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.100770Jy  sigma=2.011759
! Fit after  self-cal, rms=0.093194Jy  sigma=1.825674
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.9469 mas, bmaj=7.296 mas, bpa=-16.98 degrees
! Estimated noise=2.23859 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000143351 Jy
! Component: 100  -  total flux cleaned = -0.000114581 Jy
! Component: 150  -  total flux cleaned = 6.19847e-06 Jy
! Component: 200  -  total flux cleaned = -0.000100703 Jy
! Total flux subtracted in 200 components = -0.000100703 Jy
! Clean residual min=-0.013414 max=0.015690 Jy/beam
! Clean residual mean=-0.000029 rms=0.003626 Jy/beam
! Combined flux in latest and established models = 2.41486 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 57 components and 2.41486 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.092827Jy  sigma=1.819770
! Fit after  self-cal, rms=0.092832Jy  sigma=1.819378
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 (-7.2, 0.4).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00877838 Jy
! Component: 100  -  total flux cleaned = 0.010611 Jy
! Component: 150  -  total flux cleaned = 0.0108877 Jy
! Component: 200  -  total flux cleaned = 0.0114465 Jy
! Total flux subtracted in 200 components = 0.0114465 Jy
! Clean residual min=-0.011106 max=0.010577 Jy/beam
! Clean residual mean=-0.000042 rms=0.002835 Jy/beam
! Combined flux in latest and established models = 2.42631 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 77 components and 2.42631 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.089677Jy  sigma=1.763597
! Fit after  self-cal, rms=0.089585Jy  sigma=1.762115
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime2
! Performing amp+phase self-cal over 15 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.089585Jy  sigma=1.762115
! Fit after  self-cal, rms=0.089425Jy  sigma=1.759496
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.702 mas, bmaj=18.33 mas, bpa=-2.41 degrees
! Estimated noise=2.5738 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00278619 Jy
! Component: 100  -  total flux cleaned = 0.00279068 Jy
! Component: 150  -  total flux cleaned = 0.00289083 Jy
! Component: 200  -  total flux cleaned = 0.00322477 Jy
! Total flux subtracted in 200 components = 0.00322477 Jy
! Clean residual min=-0.008241 max=0.008347 Jy/beam
! Clean residual mean=-0.000076 rms=0.002870 Jy/beam
! Combined flux in latest and established models = 2.42954 Jy
selfcal
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 90 components and 2.42954 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.089221Jy  sigma=1.755382
! Fit after  self-cal, rms=0.089160Jy  sigma=1.754721
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.000785092 Jy
! Component: 100  -  total flux cleaned = 0.00110113 Jy
! Component: 150  -  total flux cleaned = 0.00109991 Jy
! Component: 200  -  total flux cleaned = 0.00109748 Jy
! Total flux subtracted in 200 components = 0.00109748 Jy
! Clean residual min=-0.007868 max=0.008258 Jy/beam
! Clean residual mean=-0.000069 rms=0.002727 Jy/beam
! Combined flux in latest and established models = 2.43063 Jy
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 94 components and 2.43063 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.089144Jy  sigma=1.753918
! Fit after  self-cal, rms=0.089114Jy  sigma=1.753675
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
win_mult = win_mult / 1.2
print "=== Finished amplitude self-cal for VLBA stations ==="
! === Finished amplitude self-cal for VLBA stations === 
unflag *
! Unflagging all baselines in the currently selected channels.
uncalib false, false, true
! uncal: All telescope correction flags have been un-done.
integer lun
lun = outfile ( "/tmp/ff.ff" )
do i=1,nsub
  do j=1,ntel(i)
    fprintf lun, "%s\n", telname(j,i)
    end do
  end do
close ( lun )
system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff"
![@/tmp/ff.ff]
selfant 1:FD,true
! 1:FD       self-cal  status=fixed  weight=1
selfant 1:KP,true
! 1:KP       self-cal  status=fixed  weight=1
selfant 1:LA,true
! 1:LA       self-cal  status=fixed  weight=1
selfant 1:MK,true
! 1:MK       self-cal  status=fixed  weight=1
selfant 1:PT,true
! 1:PT       self-cal  status=fixed  weight=1
![Exited script file: /tmp/ff.ff]
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.435537Jy  sigma=4.649099
! Fit after  self-cal, rms=0.422054Jy  sigma=3.890535
selfant "",false
! 1:FD       self-cal  status=correctable  weight=1
! 1:KK       self-cal  status=correctable  weight=1
! 1:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:MK       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.422054Jy  sigma=3.890535
! Fit after  self-cal, rms=0.421946Jy  sigma=3.873805
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=0.6312 mas, bmaj=5.095 mas, bpa=-15.64 degrees
! Estimated noise=3.10911 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.0143526 Jy
! Component: 100  -  total flux cleaned = 0.015744 Jy
! Component: 150  -  total flux cleaned = 0.0171468 Jy
! Component: 200  -  total flux cleaned = 0.0133334 Jy
! Total flux subtracted in 200 components = 0.0133334 Jy
! Clean residual min=-0.042801 max=0.040831 Jy/beam
! Clean residual mean=-0.000093 rms=0.010617 Jy/beam
! Combined flux in latest and established models = 2.44397 Jy
! Performing phase self-cal
! Adding 26 model components to the UV plane model.
! The established model now contains 120 components and 2.44397 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.403159Jy  sigma=3.779091
! Fit after  self-cal, rms=0.402716Jy  sigma=3.773070
! Inverting map 
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=0.7905 mas, bmaj=6.314 mas, bpa=-16.76 degrees
! Estimated noise=2.12528 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 (-9.6, 34).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00804115 Jy
! Component: 100  -  total flux cleaned = 0.0137729 Jy
! Component: 150  -  total flux cleaned = 0.0163921 Jy
! Component: 200  -  total flux cleaned = 0.0181776 Jy
! Component: 250  -  total flux cleaned = 0.0192629 Jy
! Component: 300  -  total flux cleaned = 0.0207382 Jy
! Component: 350  -  total flux cleaned = 0.0216407 Jy
! Component: 400  -  total flux cleaned = 0.0224823 Jy
! Total flux subtracted in 400 components = 0.0224823 Jy
! Clean residual min=-0.025589 max=0.018933 Jy/beam
! Clean residual mean=-0.000074 rms=0.005491 Jy/beam
! Combined flux in latest and established models = 2.46645 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 162 components and 2.46645 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.394158Jy  sigma=3.709856
! Fit after  self-cal, rms=0.393877Jy  sigma=3.703634
! Inverting map 
dynam = dynam - 0.5
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)
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00294251 Jy
! Component: 100  -  total flux cleaned = 0.00290233 Jy
! Component: 150  -  total flux cleaned = 0.00326432 Jy
! Component: 200  -  total flux cleaned = 0.00326258 Jy
! Component: 250  -  total flux cleaned = 0.00358683 Jy
! Component: 300  -  total flux cleaned = 0.00328189 Jy
! Component: 350  -  total flux cleaned = 0.00358042 Jy
! Component: 400  -  total flux cleaned = 0.0033128 Jy
! Total flux subtracted in 400 components = 0.0033128 Jy
! Clean residual min=-0.023612 max=0.017783 Jy/beam
! Clean residual mean=-0.000082 rms=0.005019 Jy/beam
! Combined flux in latest and established models = 2.46976 Jy
! Performing phase self-cal
! Adding 39 model components to the UV plane model.
! The established model now contains 188 components and 2.46976 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.391171Jy  sigma=3.688749
! Fit after  self-cal, rms=0.391083Jy  sigma=3.686639
! Inverting map 
print "=== Finished natural weighting clean for al sites ==="
! === Finished natural weighting clean for al sites === 
unflag *
! Unflagging all baselines in the currently selected channels.
uncalib false, false, true
! uncal: All telescope correction flags have been un-done.
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.391083Jy  sigma=3.686639
! Fit after  self-cal, rms=0.386160Jy  sigma=3.669788
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.7878 mas, bmaj=6.297 mas, bpa=-16.76 degrees
! Estimated noise=2.11707 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00473695 Jy
! Component: 100  -  total flux cleaned = 0.00441925 Jy
! Component: 150  -  total flux cleaned = 0.00411893 Jy
! Component: 200  -  total flux cleaned = 0.00384373 Jy
! Component: 250  -  total flux cleaned = 0.00384367 Jy
! Component: 300  -  total flux cleaned = 0.00334244 Jy
! Component: 350  -  total flux cleaned = 0.0033499 Jy
! Component: 400  -  total flux cleaned = 0.0033487 Jy
! Total flux subtracted in 400 components = 0.0033487 Jy
! Clean residual min=-0.017602 max=0.017531 Jy/beam
! Clean residual mean=-0.000074 rms=0.004281 Jy/beam
! Combined flux in latest and established models = 2.47311 Jy
selfcal
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 218 components and 2.47311 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.383509Jy  sigma=3.661830
! Fit after  self-cal, rms=0.383485Jy  sigma=3.661123
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=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000648699 Jy
! Component: 100  -  total flux cleaned = 0.000421686 Jy
! Component: 150  -  total flux cleaned = -6.30545e-06 Jy
! Component: 200  -  total flux cleaned = -5.60794e-06 Jy
! Component: 250  -  total flux cleaned = -0.000208205 Jy
! Component: 300  -  total flux cleaned = -0.000398309 Jy
! Component: 350  -  total flux cleaned = -0.000399329 Jy
! Component: 400  -  total flux cleaned = -0.000583172 Jy
! Total flux subtracted in 400 components = -0.000583172 Jy
! Clean residual min=-0.017028 max=0.016897 Jy/beam
! Clean residual mean=-0.000084 rms=0.004175 Jy/beam
! Combined flux in latest and established models = 2.47253 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 237 components and 2.47253 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.382446Jy  sigma=3.656875
! Fit after  self-cal, rms=0.382429Jy  sigma=3.656482
! Inverting map 
selfcal true, true, soltime2
! Performing amp+phase self-cal over 15 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.382429Jy  sigma=3.656482
! Fit after  self-cal, rms=0.381545Jy  sigma=3.634830
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.7868 mas, bmaj=6.29 mas, bpa=-16.77 degrees
! Estimated noise=2.11221 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000407903 Jy
! Component: 100  -  total flux cleaned = 0.000183751 Jy
! Component: 150  -  total flux cleaned = -0.000438235 Jy
! Component: 200  -  total flux cleaned = -0.000816334 Jy
! Component: 250  -  total flux cleaned = -0.00135886 Jy
! Component: 300  -  total flux cleaned = -0.00188153 Jy
! Component: 350  -  total flux cleaned = -0.00222038 Jy
! Component: 400  -  total flux cleaned = -0.0025545 Jy
! Total flux subtracted in 400 components = -0.0025545 Jy
! Clean residual min=-0.014992 max=0.017073 Jy/beam
! Clean residual mean=-0.000079 rms=0.003656 Jy/beam
! Combined flux in latest and established models = 2.46997 Jy
selfcal
! Performing phase self-cal
! Adding 48 model components to the UV plane model.
! The established model now contains 270 components and 2.46997 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.380406Jy  sigma=3.631432
! Fit after  self-cal, rms=0.380387Jy  sigma=3.630963
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 (49.6, -45.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00257254 Jy
! Component: 100  -  total flux cleaned = 0.00309921 Jy
! Component: 150  -  total flux cleaned = 0.00307355 Jy
! Component: 200  -  total flux cleaned = 0.00308147 Jy
! Component: 250  -  total flux cleaned = 0.00341115 Jy
! Component: 300  -  total flux cleaned = 0.00341512 Jy
! Component: 350  -  total flux cleaned = 0.00312923 Jy
! Component: 400  -  total flux cleaned = 0.00285251 Jy
! Total flux subtracted in 400 components = 0.00285251 Jy
! Clean residual min=-0.014927 max=0.014878 Jy/beam
! Clean residual mean=-0.000098 rms=0.003599 Jy/beam
! Combined flux in latest and established models = 2.47283 Jy
! Performing phase self-cal
! Adding 57 model components to the UV plane model.
! The established model now contains 322 components and 2.47282 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.376631Jy  sigma=3.606678
! Fit after  self-cal, rms=0.376417Jy  sigma=3.602000
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime3
! Performing amp+phase self-cal over 1 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.376417Jy  sigma=3.602000
! Fit after  self-cal, rms=0.745955Jy  sigma=4.572256
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.511 mas, bmaj=17.86 mas, bpa=-4.642 degrees
! Estimated noise=2.54699 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00639301 Jy
! Component: 100  -  total flux cleaned = -0.00447165 Jy
! Component: 150  -  total flux cleaned = -0.00448525 Jy
! Component: 200  -  total flux cleaned = -0.00450301 Jy
! Component: 250  -  total flux cleaned = -0.00775692 Jy
! Component: 300  -  total flux cleaned = -0.00773615 Jy
! Component: 350  -  total flux cleaned = -0.00772801 Jy
! Component: 400  -  total flux cleaned = -0.0091167 Jy
! Total flux subtracted in 400 components = -0.0091167 Jy
! Clean residual min=-0.041444 max=0.047135 Jy/beam
! Clean residual mean=0.000207 rms=0.019363 Jy/beam
! Combined flux in latest and established models = 2.46371 Jy
selfcal
! Performing phase self-cal
! Adding 64 model components to the UV plane model.
! The established model now contains 384 components and 2.46371 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.709877Jy  sigma=4.285316
! Fit after  self-cal, rms=0.709829Jy  sigma=4.277442
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=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00152287 Jy
! Component: 100  -  total flux cleaned = -0.000318281 Jy
! Component: 150  -  total flux cleaned = 0.00202768 Jy
! Component: 200  -  total flux cleaned = 0.00200833 Jy
! Component: 250  -  total flux cleaned = 0.00414203 Jy
! Component: 300  -  total flux cleaned = 0.00413246 Jy
! Component: 350  -  total flux cleaned = 0.00508766 Jy
! Component: 400  -  total flux cleaned = 0.00508932 Jy
! Total flux subtracted in 400 components = 0.00508932 Jy
! Clean residual min=-0.037358 max=0.047656 Jy/beam
! Clean residual mean=0.000280 rms=0.016605 Jy/beam
! Combined flux in latest and established models = 2.4688 Jy
! Performing phase self-cal
! Adding 59 model components to the UV plane model.
! The established model now contains 423 components and 2.4688 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.688742Jy  sigma=4.137385
! Fit after  self-cal, rms=0.688682Jy  sigma=4.128964
! Inverting map 
! Added new window around map position (18.2, -31.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0142129 Jy
! Component: 100  -  total flux cleaned = 0.026957 Jy
! Component: 150  -  total flux cleaned = 0.0382518 Jy
! Component: 200  -  total flux cleaned = 0.0508967 Jy
! Component: 250  -  total flux cleaned = 0.0618993 Jy
! Component: 300  -  total flux cleaned = 0.0713797 Jy
! Component: 350  -  total flux cleaned = 0.0768093 Jy
! Component: 400  -  total flux cleaned = 0.0819939 Jy
! Total flux subtracted in 400 components = 0.0819939 Jy
! Clean residual min=-0.023743 max=0.023723 Jy/beam
! Clean residual mean=0.000161 rms=0.010127 Jy/beam
! Combined flux in latest and established models = 2.55079 Jy
! Performing phase self-cal
! Adding 88 model components to the UV plane model.
! The established model now contains 508 components and 2.55079 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.672403Jy  sigma=3.977625
! Fit after  self-cal, rms=0.672012Jy  sigma=3.968826
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
win_mult = win_mult / 1.2
print "=== Finished amplitude self-cal for all stations ==="
! === Finished amplitude self-cal for all stations === 
restore
! Inverting map and beam 
! Estimated beam: bmin=0.8094 mas, bmaj=6.425 mas, bpa=-16.07 degrees
! Estimated noise=2.21286 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.8094 x 6.425 at -16.07 degrees (North through East)
! Clean map  min=-0.068352  max=2.2831 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.672012Jy  sigma=3.968826
! Fit after  self-cal, rms=0.843751Jy  sigma=2.135137
! Inverting map and beam 
! Estimated beam: bmin=0.8065 mas, bmaj=6.398 mas, bpa=-16.22 degrees
! Estimated noise=2.17844 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0434703 Jy
! Component: 100  -  total flux cleaned = -0.0594929 Jy
! Component: 150  -  total flux cleaned = -0.0689681 Jy
! Component: 200  -  total flux cleaned = -0.0783311 Jy
! Component: 250  -  total flux cleaned = -0.0850498 Jy
! Component: 300  -  total flux cleaned = -0.0866473 Jy
! Component: 350  -  total flux cleaned = -0.085904 Jy
! Component: 400  -  total flux cleaned = -0.0858578 Jy
! Total flux subtracted in 400 components = -0.0858578 Jy
! Clean residual min=-0.015657 max=0.018391 Jy/beam
! Clean residual mean=-0.000162 rms=0.005046 Jy/beam
! Combined flux in latest and established models = 2.46493 Jy
! Performing phase self-cal
! Adding 61 model components to the UV plane model.
! The established model now contains 565 components and 2.46493 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.839181Jy  sigma=1.757955
! Fit after  self-cal, rms=0.839004Jy  sigma=1.751365
! Inverting map 
! Added new window around map position (-21, -47.4).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00466167 Jy
! Component: 100  -  total flux cleaned = 0.00401926 Jy
! Component: 150  -  total flux cleaned = 0.00288929 Jy
! Component: 200  -  total flux cleaned = 0.00286817 Jy
! Component: 250  -  total flux cleaned = 0.00283927 Jy
! Component: 300  -  total flux cleaned = 0.00282163 Jy
! Component: 350  -  total flux cleaned = 0.0019471 Jy
! Component: 400  -  total flux cleaned = 0.00154389 Jy
! Total flux subtracted in 400 components = 0.00154389 Jy
! Clean residual min=-0.011302 max=0.013578 Jy/beam
! Clean residual mean=-0.000154 rms=0.003406 Jy/beam
! Combined flux in latest and established models = 2.46648 Jy
! Performing phase self-cal
! Adding 76 model components to the UV plane model.
! The established model now contains 616 components and 2.46648 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.836830Jy  sigma=1.634937
! Fit after  self-cal, rms=0.836809Jy  sigma=1.631480
! Inverting map 
! Added new window around map position (-46.6, 48.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00605874 Jy
! Component: 100  -  total flux cleaned = 0.00812143 Jy
! Component: 150  -  total flux cleaned = 0.00814506 Jy
! Component: 200  -  total flux cleaned = 0.0077933 Jy
! Component: 250  -  total flux cleaned = 0.00813506 Jy
! Component: 300  -  total flux cleaned = 0.00812809 Jy
! Component: 350  -  total flux cleaned = 0.00811166 Jy
! Component: 400  -  total flux cleaned = 0.00842093 Jy
! Total flux subtracted in 400 components = 0.00842093 Jy
! Clean residual min=-0.008779 max=0.009542 Jy/beam
! Clean residual mean=-0.000142 rms=0.002461 Jy/beam
! Combined flux in latest and established models = 2.4749 Jy
! Performing phase self-cal
! Adding 79 model components to the UV plane model.
! The established model now contains 677 components and 2.4749 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.836384Jy  sigma=1.571112
! Fit after  self-cal, rms=0.836464Jy  sigma=1.569041
! Inverting map 
! Added new window around map position (50.2, 9.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00322521 Jy
! Component: 100  -  total flux cleaned = 0.00357619 Jy
! Component: 150  -  total flux cleaned = 0.00526576 Jy
! Component: 200  -  total flux cleaned = 0.00763605 Jy
! Component: 250  -  total flux cleaned = 0.00913575 Jy
! Component: 300  -  total flux cleaned = 0.0103346 Jy
! Component: 350  -  total flux cleaned = 0.0112614 Jy
! Component: 400  -  total flux cleaned = 0.0121573 Jy
! Total flux subtracted in 400 components = 0.0121573 Jy
! Clean residual min=-0.006853 max=0.006786 Jy/beam
! Clean residual mean=-0.000112 rms=0.001898 Jy/beam
! Combined flux in latest and established models = 2.48706 Jy
! Performing phase self-cal
! Adding 72 model components to the UV plane model.
! The established model now contains 726 components and 2.48706 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.835373Jy  sigma=1.535675
! Fit after  self-cal, rms=0.835375Jy  sigma=1.534831
! 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.835375Jy  sigma=1.534831
! Fit after  self-cal, rms=0.905397Jy  sigma=1.513591
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=0.6373 mas, bmaj=5.167 mas, bpa=-15.07 degrees
! Estimated noise=3.20634 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.72919 Jy
! Component: 100  -  total flux cleaned = 2.10927 Jy
! Component: 150  -  total flux cleaned = 2.23314 Jy
! Component: 200  -  total flux cleaned = 2.29527 Jy
! Component: 250  -  total flux cleaned = 2.30588 Jy
! Component: 300  -  total flux cleaned = 2.30135 Jy
! Component: 350  -  total flux cleaned = 2.3024 Jy
! Component: 400  -  total flux cleaned = 2.30862 Jy
! Total flux subtracted in 400 components = 2.30862 Jy
! Clean residual min=-0.024703 max=0.023025 Jy/beam
! Clean residual mean=-0.000004 rms=0.006497 Jy/beam
! Combined flux in latest and established models = 2.30862 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 81 model components to the UV plane model.
! The established model now contains 81 components and 2.30862 Jy
! Inverting map and beam 
! Estimated beam: bmin=0.8069 mas, bmaj=6.4 mas, bpa=-16.21 degrees
! Estimated noise=2.17795 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=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0374325 Jy
! Component: 100  -  total flux cleaned = 0.0442891 Jy
! Component: 150  -  total flux cleaned = 0.0515819 Jy
! Component: 200  -  total flux cleaned = 0.0548286 Jy
! Component: 250  -  total flux cleaned = 0.0578086 Jy
! Component: 300  -  total flux cleaned = 0.0592149 Jy
! Component: 350  -  total flux cleaned = 0.0598678 Jy
! Component: 400  -  total flux cleaned = 0.0593178 Jy
! Total flux subtracted in 400 components = 0.0593178 Jy
! Clean residual min=-0.014536 max=0.017228 Jy/beam
! Clean residual mean=-0.000098 rms=0.005022 Jy/beam
! Combined flux in latest and established models = 2.36793 Jy
! Adding 103 model components to the UV plane model.
! The established model now contains 181 components and 2.36793 Jy
! Inverting map 
! Added new window around map position (-21.2, -21.6).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00666277 Jy
! Component: 100  -  total flux cleaned = 0.0045394 Jy
! Component: 150  -  total flux cleaned = 0.00453112 Jy
! Component: 200  -  total flux cleaned = 0.004065 Jy
! Component: 250  -  total flux cleaned = 0.00448411 Jy
! Component: 300  -  total flux cleaned = 0.00487292 Jy
! Component: 350  -  total flux cleaned = 0.00486688 Jy
! Component: 400  -  total flux cleaned = 0.00450254 Jy
! Total flux subtracted in 400 components = 0.00450254 Jy
! Clean residual min=-0.010052 max=0.010743 Jy/beam
! Clean residual mean=-0.000062 rms=0.003292 Jy/beam
! Combined flux in latest and established models = 2.37244 Jy
! Adding 92 model components to the UV plane model.
! The established model now contains 257 components and 2.37244 Jy
! Inverting map 
! Added new window around map position (-49, 23.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00616234 Jy
! Component: 100  -  total flux cleaned = 0.0079407 Jy
! Component: 150  -  total flux cleaned = 0.00858915 Jy
! Component: 200  -  total flux cleaned = 0.00919735 Jy
! Component: 250  -  total flux cleaned = 0.0109017 Jy
! Component: 300  -  total flux cleaned = 0.0117014 Jy
! Component: 350  -  total flux cleaned = 0.0119612 Jy
! Component: 400  -  total flux cleaned = 0.0122287 Jy
! Total flux subtracted in 400 components = 0.0122287 Jy
! Clean residual min=-0.007925 max=0.007770 Jy/beam
! Clean residual mean=-0.000066 rms=0.002368 Jy/beam
! Combined flux in latest and established models = 2.38466 Jy
! Adding 85 model components to the UV plane model.
! The established model now contains 317 components and 2.38466 Jy
! Inverting map 
! Added new window around map position (42.2, 46.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00518232 Jy
! Component: 100  -  total flux cleaned = 0.00644148 Jy
! Component: 150  -  total flux cleaned = 0.00621637 Jy
! Component: 200  -  total flux cleaned = 0.00687875 Jy
! Component: 250  -  total flux cleaned = 0.00644303 Jy
! Component: 300  -  total flux cleaned = 0.00644759 Jy
! Component: 350  -  total flux cleaned = 0.00685611 Jy
! Component: 400  -  total flux cleaned = 0.00705888 Jy
! Total flux subtracted in 400 components = 0.00705888 Jy
! Clean residual min=-0.006834 max=0.006352 Jy/beam
! Clean residual mean=-0.000068 rms=0.001895 Jy/beam
! Combined flux in latest and established models = 2.39172 Jy
! Adding 88 model components to the UV plane model.
! The established model now contains 382 components and 2.39172 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.528 mas, bmaj=17.87 mas, bpa=-4.704 degrees
! Estimated noise=2.51081 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00178018 Jy
! Component: 100  -  total flux cleaned = -0.00244199 Jy
! Component: 150  -  total flux cleaned = -0.0032219 Jy
! Component: 200  -  total flux cleaned = -0.00394385 Jy
! Component: 250  -  total flux cleaned = -0.00444354 Jy
! Component: 300  -  total flux cleaned = -0.0050726 Jy
! Component: 350  -  total flux cleaned = -0.00566644 Jy
! Component: 400  -  total flux cleaned = -0.00650708 Jy
! Total flux subtracted in 400 components = -0.00650708 Jy
! Clean residual min=-0.005415 max=0.005428 Jy/beam
! Clean residual mean=-0.000058 rms=0.001757 Jy/beam
! Combined flux in latest and established models = 2.38522 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)
! Added new window around map position (-51, -23.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00311616 Jy
! Component: 100  -  total flux cleaned = 0.00517903 Jy
! Component: 150  -  total flux cleaned = 0.0060092 Jy
! Component: 200  -  total flux cleaned = 0.00679891 Jy
! Component: 250  -  total flux cleaned = 0.00743211 Jy
! Component: 300  -  total flux cleaned = 0.00829669 Jy
! Component: 350  -  total flux cleaned = 0.00901577 Jy
! Component: 400  -  total flux cleaned = 0.00996439 Jy
! Total flux subtracted in 400 components = 0.00996439 Jy
! Clean residual min=-0.003918 max=0.004238 Jy/beam
! Clean residual mean=-0.000005 rms=0.001313 Jy/beam
! Combined flux in latest and established models = 2.39518 Jy
! Adding 234 model components to the UV plane model.
! The established model now contains 606 components and 2.39518 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=0.8069 mas, bmaj=6.4 mas, bpa=-16.21 degrees
! Estimated noise=2.17795 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.8443e-05 Jy
! Component: 100  -  total flux cleaned = -0.000708945 Jy
! Component: 150  -  total flux cleaned = -0.000894281 Jy
! Component: 200  -  total flux cleaned = -0.000896398 Jy
! Component: 250  -  total flux cleaned = -0.000295996 Jy
! Component: 300  -  total flux cleaned = 0.00012722 Jy
! Component: 350  -  total flux cleaned = 0.000262648 Jy
! Component: 400  -  total flux cleaned = 0.000925618 Jy
! Total flux subtracted in 400 components = 0.000925618 Jy
! Clean residual min=-0.004383 max=0.004809 Jy/beam
! Clean residual mean=-0.000026 rms=0.001151 Jy/beam
! Combined flux in latest and established models = 2.39611 Jy
selfcal
! Performing phase self-cal
! Adding 70 model components to the UV plane model.
! The established model now contains 670 components and 2.39611 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.905331Jy  sigma=1.484289
! Fit after  self-cal, rms=0.905201Jy  sigma=1.475708
wmodel /scr/pima/rdv62_uvs/J2056-4714_X_map.mod
! Writing 670 model components to file: /scr/pima/rdv62_uvs/J2056-4714_X_map.mod
wobs   /scr/pima/rdv62_uvs/J2056-4714_X_uvs.fits
! Applying 217 buffered edits.
! Writing UV FITS file: /scr/pima/rdv62_uvs/J2056-4714_X_uvs.fits
wwins  /scr/pima/rdv62_uvs/J2056-4714_X_map.win
! wwins: Wrote 12 windows to /scr/pima/rdv62_uvs/J2056-4714_X_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.000221513 Jy
! Component: 100  -  total flux cleaned = -0.00114212 Jy
! Component: 150  -  total flux cleaned = -0.00157106 Jy
! Component: 200  -  total flux cleaned = -0.00190612 Jy
! Component: 250  -  total flux cleaned = -0.00196904 Jy
! Component: 300  -  total flux cleaned = -0.00221366 Jy
! Component: 350  -  total flux cleaned = -0.00191821 Jy
! Component: 400  -  total flux cleaned = -0.00162883 Jy
! Component: 450  -  total flux cleaned = -0.00129102 Jy
! Component: 500  -  total flux cleaned = -0.00107094 Jy
! Component: 550  -  total flux cleaned = -0.000747446 Jy
! Component: 600  -  total flux cleaned = -0.000535829 Jy
! Component: 650  -  total flux cleaned = -0.00032791 Jy
! Component: 700  -  total flux cleaned = -0.000226378 Jy
! Component: 750  -  total flux cleaned = -0.000126928 Jy
! Component: 800  -  total flux cleaned = -0.000176693 Jy
! Component: 850  -  total flux cleaned = -0.000273916 Jy
! Component: 900  -  total flux cleaned = -0.000130274 Jy
! Component: 950  -  total flux cleaned = -0.000177982 Jy
! Component: 1000  -  total flux cleaned = -0.000224327 Jy
! Component: 1050  -  total flux cleaned = -8.58671e-05 Jy
! Component: 1100  -  total flux cleaned = -4.04902e-05 Jy
! Component: 1150  -  total flux cleaned = 4.47544e-06 Jy
! Component: 1200  -  total flux cleaned = 0.00022782 Jy
! Component: 1250  -  total flux cleaned = 0.000492903 Jy
! Component: 1300  -  total flux cleaned = 0.000842566 Jy
! Component: 1350  -  total flux cleaned = 0.000928739 Jy
! Component: 1400  -  total flux cleaned = 0.00114324 Jy
! Component: 1450  -  total flux cleaned = 0.00156816 Jy
! Component: 1500  -  total flux cleaned = 0.00190395 Jy
! Component: 1550  -  total flux cleaned = 0.00202896 Jy
! Component: 1600  -  total flux cleaned = 0.00219402 Jy
! Component: 1650  -  total flux cleaned = 0.00223535 Jy
! Component: 1700  -  total flux cleaned = 0.00268163 Jy
! Component: 1750  -  total flux cleaned = 0.00280218 Jy
! Component: 1800  -  total flux cleaned = 0.00308159 Jy
! Component: 1850  -  total flux cleaned = 0.00339855 Jy
! Component: 1900  -  total flux cleaned = 0.00343817 Jy
! Component: 1950  -  total flux cleaned = 0.00375056 Jy
! Component: 2000  -  total flux cleaned = 0.00398331 Jy
! Component: 2050  -  total flux cleaned = 0.00421472 Jy
! Component: 2100  -  total flux cleaned = 0.00440668 Jy
! Component: 2150  -  total flux cleaned = 0.00459776 Jy
! Component: 2200  -  total flux cleaned = 0.00475011 Jy
! Component: 2250  -  total flux cleaned = 0.00490177 Jy
! Component: 2300  -  total flux cleaned = 0.00505306 Jy
! Component: 2350  -  total flux cleaned = 0.0050906 Jy
! Component: 2400  -  total flux cleaned = 0.00531616 Jy
! Component: 2450  -  total flux cleaned = 0.00535378 Jy
! Component: 2500  -  total flux cleaned = 0.00557827 Jy
! Component: 2550  -  total flux cleaned = 0.00569002 Jy
! Component: 2600  -  total flux cleaned = 0.00583926 Jy
! Component: 2650  -  total flux cleaned = 0.00598755 Jy
! Component: 2700  -  total flux cleaned = 0.00621005 Jy
! Component: 2750  -  total flux cleaned = 0.00635811 Jy
! Component: 2800  -  total flux cleaned = 0.0065058 Jy
! Component: 2850  -  total flux cleaned = 0.0065796 Jy
! Component: 2900  -  total flux cleaned = 0.00661655 Jy
! Component: 2950  -  total flux cleaned = 0.00676394 Jy
! Component: 3000  -  total flux cleaned = 0.00680076 Jy
! Component: 3050  -  total flux cleaned = 0.00698462 Jy
! Component: 3100  -  total flux cleaned = 0.0071682 Jy
! Component: 3150  -  total flux cleaned = 0.007205 Jy
! Component: 3200  -  total flux cleaned = 0.00724171 Jy
! Component: 3250  -  total flux cleaned = 0.00731525 Jy
! Component: 3300  -  total flux cleaned = 0.00735173 Jy
! Component: 3350  -  total flux cleaned = 0.00749888 Jy
! Component: 3400  -  total flux cleaned = 0.00746192 Jy
! Component: 3450  -  total flux cleaned = 0.00760952 Jy
! Component: 3500  -  total flux cleaned = 0.0076464 Jy
! Component: 3550  -  total flux cleaned = 0.00775744 Jy
! Component: 3600  -  total flux cleaned = 0.00779469 Jy
! Component: 3650  -  total flux cleaned = 0.0079062 Jy
! Component: 3700  -  total flux cleaned = 0.00786887 Jy
! Component: 3750  -  total flux cleaned = 0.00801813 Jy
! Component: 3800  -  total flux cleaned = 0.00805568 Jy
! Component: 3850  -  total flux cleaned = 0.00816796 Jy
! Component: 3900  -  total flux cleaned = 0.00824298 Jy
! Component: 3950  -  total flux cleaned = 0.00835571 Jy
! Component: 4000  -  total flux cleaned = 0.008431 Jy
! Component: 4050  -  total flux cleaned = 0.00846895 Jy
! Total flux subtracted in 4096 components = 0.00858286 Jy
! Clean residual min=-0.002139 max=0.001925 Jy/beam
! Clean residual mean=0.000025 rms=0.000679 Jy/beam
! Combined flux in latest and established models = 2.40469 Jy
keep
! Adding 181 model components to the UV plane model.
! The established model now contains 847 components and 2.40469 Jy
![Exited script file: /home/lpetrov/bin/rdv_muppet]
wmap "/scr/pima/rdv62_uvs/J2056-4714_X\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.8069 x 6.4 at -16.21 degrees (North through East)
! Clean map  min=-0.013283  max=2.1536 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J2056-4714_X_map.fits
float peak_flux
peak_flux = peak(flux,max)
invert
! Inverting map 
float image_peak
image_peak = int(10000*peak_flux)/10000
float image_rms
image_rms  = int(100000*imstat(rms))/100000
print "========================================="
! ========================================= 
print "The image peak= ", image_peak, " rms= ", image_rms, " Jy"
! The image peak=  2.1535  rms=  0.00065  Jy 
quit
! Quitting program
! Log file /scr/pima/rdv62_uvs/J2056-4714_X_uvs.log closed on Fri Nov 22 22:42:05 2013