! Started logfile: /scr/pima/rdv95_uvs/J0922-3959_S_uvs.log on Wed Nov  6 20:45:09 2013
obs /scr/pima/rdv95_uvs/J0922-3959_S_uva.fits
! Reading UV FITS file: /scr/pima/rdv95_uvs/J0922-3959_S_uva.fits
! AN table 1: 15 integrations on 15 of 15 possible baselines.
! Apparent sampling: 0.831111 visibilities/baseline/integration-bin.
! Found source: J0922-3959
! 
! 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.23299e+09        8e+06          1        8e+06
!  02        2  2.26299e+09        8e+06          1        8e+06
!  03        3  2.35299e+09        8e+06          1        8e+06
!  04        4  2.37299e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 0 lines of history.
! 
! Reading 748 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.4 
print "taper_size =",taper_size
! taper_size = 20 
![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv95_uvs/J0922-3959_S]
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:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:MK       self-cal  status=correctable  weight=1
! 1:OV       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.400x0.400 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 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.465728Jy  sigma=18.645961
! Fit after  self-cal, rms=0.276820Jy  sigma=3.703606
! 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 \]
![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=2.985 mas, bmaj=28.63 mas, bpa=-11.86 degrees
! Estimated noise=3.59388 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=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.88468 Jy
! Component: 100  -  total flux cleaned = 1.09947 Jy
! Total flux subtracted in 100 components = 1.09947 Jy
! Clean residual min=-0.055017 max=0.087354 Jy/beam
! Clean residual mean=0.000032 rms=0.023118 Jy/beam
! Combined flux in latest and established models = 1.09947 Jy
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 17 components and 1.09947 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.215487Jy  sigma=2.856143
! Fit after  self-cal, rms=0.215464Jy  sigma=2.855867
! 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=4.185 mas, bmaj=29.32 mas, bpa=-10.95 degrees
! Estimated noise=2.0558 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.114356 Jy
! Component: 100  -  total flux cleaned = 0.149577 Jy
! Component: 150  -  total flux cleaned = 0.149669 Jy
! Component: 200  -  total flux cleaned = 0.148219 Jy
! Total flux subtracted in 200 components = 0.148219 Jy
! Clean residual min=-0.021715 max=0.021495 Jy/beam
! Clean residual mean=-0.000387 rms=0.009136 Jy/beam
! Combined flux in latest and established models = 1.24769 Jy
! Performing phase self-cal
! Adding 85 model components to the UV plane model.
! The established model now contains 100 components and 1.24769 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.156289Jy  sigma=2.008162
! Fit after  self-cal, rms=0.155956Jy  sigma=2.004589
! 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.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.97     KP        0.94     LA        1.19     MK        0.99     
!   OV        0.98     PT        1.03     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.04     KP        1.00     LA        1.07     MK        1.01     
!   OV        1.03     PT        1.03     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.06     KP        1.01     LA        1.00     MK        0.98     
!   OV        1.04     PT        1.06     
! 
! 
! Correcting IF 4.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.04     KP        0.94     LA        0.84     MK        1.08     
!   OV        0.88     PT        0.94     
! 
! 
! Fit before self-cal, rms=0.155956Jy  sigma=2.004589
! Fit after  self-cal, rms=0.072480Jy  sigma=0.904524
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=4.216 mas, bmaj=29.25 mas, bpa=-10.94 degrees
! Estimated noise=2.05728 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.0112717 Jy
! Component: 100  -  total flux cleaned = 0.0156066 Jy
! Component: 150  -  total flux cleaned = 0.018957 Jy
! Component: 200  -  total flux cleaned = 0.0216267 Jy
! Total flux subtracted in 200 components = 0.0216267 Jy
! Clean residual min=-0.005127 max=0.005147 Jy/beam
! Clean residual mean=-0.000009 rms=0.002161 Jy/beam
! Combined flux in latest and established models = 1.26931 Jy
! Performing phase self-cal
! Adding 56 model components to the UV plane model.
! The established model now contains 155 components and 1.26931 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.065298Jy  sigma=0.813319
! Fit after  self-cal, rms=0.065080Jy  sigma=0.811093
! 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.065080Jy  sigma=0.811093
! Fit after  self-cal, rms=0.064900Jy  sigma=0.809213
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.215 mas, bmaj=29.23 mas, bpa=-10.94 degrees
! Estimated noise=2.0546 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00102053 Jy
! Component: 100  -  total flux cleaned = 0.00142743 Jy
! Component: 150  -  total flux cleaned = 0.00236007 Jy
! Component: 200  -  total flux cleaned = 0.0022103 Jy
! Total flux subtracted in 200 components = 0.0022103 Jy
! Clean residual min=-0.002632 max=0.002632 Jy/beam
! Clean residual mean=-0.000035 rms=0.001095 Jy/beam
! Combined flux in latest and established models = 1.27152 Jy
selfcal
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 195 components and 1.27152 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.063509Jy  sigma=0.793022
! Fit after  self-cal, rms=0.063440Jy  sigma=0.792582
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 = 3.94013e-05 Jy
! Component: 100  -  total flux cleaned = -0.00035293 Jy
! Component: 150  -  total flux cleaned = -0.000107801 Jy
! Component: 200  -  total flux cleaned = -0.000111664 Jy
! Total flux subtracted in 200 components = -0.000111664 Jy
! Clean residual min=-0.001778 max=0.001771 Jy/beam
! Clean residual mean=-0.000038 rms=0.000730 Jy/beam
! Combined flux in latest and established models = 1.27141 Jy
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 233 components and 1.27141 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062930Jy  sigma=0.787473
! Fit after  self-cal, rms=0.062894Jy  sigma=0.787185
! 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.062894Jy  sigma=0.787185
! Fit after  self-cal, rms=0.062821Jy  sigma=0.786025
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.214 mas, bmaj=29.23 mas, bpa=-10.94 degrees
! Estimated noise=2.05504 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000498897 Jy
! Component: 100  -  total flux cleaned = -0.00102805 Jy
! Component: 150  -  total flux cleaned = -0.00121757 Jy
! Component: 200  -  total flux cleaned = -0.00130821 Jy
! Total flux subtracted in 200 components = -0.00130821 Jy
! Clean residual min=-0.001439 max=0.001447 Jy/beam
! Clean residual mean=-0.000026 rms=0.000586 Jy/beam
! Combined flux in latest and established models = 1.2701 Jy
selfcal
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 266 components and 1.2701 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062498Jy  sigma=0.782772
! Fit after  self-cal, rms=0.062473Jy  sigma=0.782626
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.000801629 Jy
! Component: 100  -  total flux cleaned = -0.000573101 Jy
! Component: 150  -  total flux cleaned = -0.00019254 Jy
! Component: 200  -  total flux cleaned = 0.000322679 Jy
! Total flux subtracted in 200 components = 0.000322679 Jy
! Clean residual min=-0.001205 max=0.001195 Jy/beam
! Clean residual mean=-0.000027 rms=0.000484 Jy/beam
! Combined flux in latest and established models = 1.27043 Jy
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 300 components and 1.27043 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062268Jy  sigma=0.780773
! Fit after  self-cal, rms=0.062260Jy  sigma=0.780694
! 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 15 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062260Jy  sigma=0.780694
! Fit after  self-cal, rms=0.062214Jy  sigma=0.779997
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=11.61 mas, bmaj=47.07 mas, bpa=7.121 degrees
! Estimated noise=2.2426 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000589751 Jy
! Component: 100  -  total flux cleaned = -0.00124692 Jy
! Component: 150  -  total flux cleaned = -0.00166772 Jy
! Component: 200  -  total flux cleaned = -0.00150122 Jy
! Total flux subtracted in 200 components = -0.00150122 Jy
! Clean residual min=-0.000912 max=0.000910 Jy/beam
! Clean residual mean=-0.000034 rms=0.000392 Jy/beam
! Combined flux in latest and established models = 1.26893 Jy
selfcal
! Performing phase self-cal
! Adding 59 model components to the UV plane model.
! The established model now contains 357 components and 1.26892 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062120Jy  sigma=0.778843
! Fit after  self-cal, rms=0.062100Jy  sigma=0.778710
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.000206756 Jy
! Component: 100  -  total flux cleaned = 9.83726e-05 Jy
! Component: 150  -  total flux cleaned = -0.000205407 Jy
! Component: 200  -  total flux cleaned = -0.000449318 Jy
! Total flux subtracted in 200 components = -0.000449318 Jy
! Clean residual min=-0.000793 max=0.000792 Jy/beam
! Clean residual mean=-0.000031 rms=0.000351 Jy/beam
! Combined flux in latest and established models = 1.26848 Jy
! Performing phase self-cal
! Adding 55 model components to the UV plane model.
! The established model now contains 400 components and 1.26848 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062033Jy  sigma=0.777902
! Fit after  self-cal, rms=0.062027Jy  sigma=0.777847
! 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:OV,true
! 1:OV       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.132920Jy  sigma=1.589390
! Fit after  self-cal, rms=0.132920Jy  sigma=1.589390
selfant "",false
! 1:FD       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:OV       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 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.132920Jy  sigma=1.589390
! Fit after  self-cal, rms=0.062028Jy  sigma=0.777770
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.996 mas, bmaj=28.57 mas, bpa=-11.9 degrees
! Estimated noise=3.58484 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.00114959 Jy
! Component: 100  -  total flux cleaned = 0.00156806 Jy
! Component: 150  -  total flux cleaned = 0.00173975 Jy
! Component: 200  -  total flux cleaned = 0.00187441 Jy
! Total flux subtracted in 200 components = 0.00187441 Jy
! Clean residual min=-0.002174 max=0.002169 Jy/beam
! Clean residual mean=-0.000009 rms=0.000916 Jy/beam
! Combined flux in latest and established models = 1.27035 Jy
! Performing phase self-cal
! Adding 75 model components to the UV plane model.
! The established model now contains 471 components and 1.27035 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061869Jy  sigma=0.776563
! Fit after  self-cal, rms=0.061848Jy  sigma=0.776398
! 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=4.212 mas, bmaj=29.24 mas, bpa=-10.94 degrees
! Estimated noise=2.05574 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=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000956991 Jy
! Component: 100  -  total flux cleaned = 0.00150097 Jy
! Component: 150  -  total flux cleaned = 0.00205966 Jy
! Component: 200  -  total flux cleaned = 0.00252675 Jy
! Component: 250  -  total flux cleaned = 0.00291839 Jy
! Component: 300  -  total flux cleaned = 0.00313301 Jy
! Component: 350  -  total flux cleaned = 0.00318442 Jy
! Component: 400  -  total flux cleaned = 0.00323693 Jy
! Total flux subtracted in 400 components = 0.00323693 Jy
! Clean residual min=-0.000866 max=0.000861 Jy/beam
! Clean residual mean=-0.000031 rms=0.000355 Jy/beam
! Combined flux in latest and established models = 1.27359 Jy
! Performing phase self-cal
! Adding 81 model components to the UV plane model.
! The established model now contains 548 components and 1.27359 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061636Jy  sigma=0.774204
! Fit after  self-cal, rms=0.061629Jy  sigma=0.774133
! 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.000219805 Jy
! Component: 100  -  total flux cleaned = -0.000319341 Jy
! Component: 150  -  total flux cleaned = -0.000366604 Jy
! Component: 200  -  total flux cleaned = -0.000503553 Jy
! Component: 250  -  total flux cleaned = -0.000457897 Jy
! Component: 300  -  total flux cleaned = -0.000371463 Jy
! Component: 350  -  total flux cleaned = -0.000542122 Jy
! Component: 400  -  total flux cleaned = -0.000542699 Jy
! Total flux subtracted in 400 components = -0.000542699 Jy
! Clean residual min=-0.000710 max=0.000707 Jy/beam
! Clean residual mean=-0.000028 rms=0.000301 Jy/beam
! Combined flux in latest and established models = 1.27304 Jy
! Performing phase self-cal
! Adding 80 model components to the UV plane model.
! The established model now contains 610 components and 1.27304 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061498Jy  sigma=0.772803
! Fit after  self-cal, rms=0.061497Jy  sigma=0.772763
! 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 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.132824Jy  sigma=1.588370
! Fit after  self-cal, rms=0.061449Jy  sigma=0.772492
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.209 mas, bmaj=29.24 mas, bpa=-10.94 degrees
! Estimated noise=2.05517 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000294848 Jy
! Component: 100  -  total flux cleaned = 0.000201392 Jy
! Component: 150  -  total flux cleaned = 0.000201116 Jy
! Component: 200  -  total flux cleaned = 0.000244737 Jy
! Component: 250  -  total flux cleaned = 0.000203081 Jy
! Component: 300  -  total flux cleaned = 3.39655e-05 Jy
! Component: 350  -  total flux cleaned = 3.43074e-05 Jy
! Component: 400  -  total flux cleaned = -0.000189012 Jy
! Total flux subtracted in 400 components = -0.000189012 Jy
! Clean residual min=-0.000764 max=0.000761 Jy/beam
! Clean residual mean=-0.000025 rms=0.000289 Jy/beam
! Combined flux in latest and established models = 1.27285 Jy
selfcal
! Performing phase self-cal
! Adding 61 model components to the UV plane model.
! The established model now contains 662 components and 1.27285 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061328Jy  sigma=0.771395
! Fit after  self-cal, rms=0.061328Jy  sigma=0.771368
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 (-9.6, 21.6).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -3.03301e-05 Jy
! Component: 100  -  total flux cleaned = -0.000346921 Jy
! Component: 150  -  total flux cleaned = -0.000501723 Jy
! Component: 200  -  total flux cleaned = -0.000647089 Jy
! Component: 250  -  total flux cleaned = -0.000793542 Jy
! Component: 300  -  total flux cleaned = -0.000887454 Jy
! Component: 350  -  total flux cleaned = -0.0009874 Jy
! Component: 400  -  total flux cleaned = -0.00103468 Jy
! Total flux subtracted in 400 components = -0.00103468 Jy
! Clean residual min=-0.001088 max=0.001077 Jy/beam
! Clean residual mean=-0.000023 rms=0.000261 Jy/beam
! Combined flux in latest and established models = 1.27182 Jy
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 707 components and 1.27182 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061304Jy  sigma=0.771130
! Fit after  self-cal, rms=0.061304Jy  sigma=0.771116
! 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.061304Jy  sigma=0.771116
! Fit after  self-cal, rms=0.061278Jy  sigma=0.770598
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.21 mas, bmaj=29.24 mas, bpa=-10.94 degrees
! Estimated noise=2.05545 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -6.2362e-05 Jy
! Component: 100  -  total flux cleaned = -0.000135046 Jy
! Component: 150  -  total flux cleaned = -0.000202405 Jy
! Component: 200  -  total flux cleaned = -0.000251597 Jy
! Component: 250  -  total flux cleaned = -0.000275071 Jy
! Component: 300  -  total flux cleaned = -0.00031898 Jy
! Component: 350  -  total flux cleaned = -0.000350082 Jy
! Component: 400  -  total flux cleaned = -0.000359881 Jy
! Total flux subtracted in 400 components = -0.000359881 Jy
! Clean residual min=-0.001157 max=0.001139 Jy/beam
! Clean residual mean=-0.000022 rms=0.000274 Jy/beam
! Combined flux in latest and established models = 1.27146 Jy
selfcal
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 744 components and 1.27146 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061265Jy  sigma=0.770516
! Fit after  self-cal, rms=0.061267Jy  sigma=0.770505
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 = 9.0072e-06 Jy
! Component: 100  -  total flux cleaned = -3.53318e-05 Jy
! Component: 150  -  total flux cleaned = -5.56282e-05 Jy
! Component: 200  -  total flux cleaned = -6.53638e-05 Jy
! Component: 250  -  total flux cleaned = -6.53421e-05 Jy
! Component: 300  -  total flux cleaned = -6.50288e-05 Jy
! Component: 350  -  total flux cleaned = -6.46477e-05 Jy
! Component: 400  -  total flux cleaned = -5.5968e-05 Jy
! Total flux subtracted in 400 components = -5.5968e-05 Jy
! Clean residual min=-0.001142 max=0.001111 Jy/beam
! Clean residual mean=-0.000023 rms=0.000280 Jy/beam
! Combined flux in latest and established models = 1.2714 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 768 components and 1.2714 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061259Jy  sigma=0.770452
! Fit after  self-cal, rms=0.061260Jy  sigma=0.770444
! Inverting map 
! Added new window around map position (-102.4, -71.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000355918 Jy
! Component: 100  -  total flux cleaned = 0.000185634 Jy
! Component: 150  -  total flux cleaned = 2.20232e-05 Jy
! Component: 200  -  total flux cleaned = 4.28052e-05 Jy
! Component: 250  -  total flux cleaned = 6.17981e-05 Jy
! Component: 300  -  total flux cleaned = 2.58686e-05 Jy
! Component: 350  -  total flux cleaned = -2.60612e-05 Jy
! Component: 400  -  total flux cleaned = -4.34623e-05 Jy
! Total flux subtracted in 400 components = -4.34623e-05 Jy
! Clean residual min=-0.001144 max=0.000978 Jy/beam
! Clean residual mean=-0.000028 rms=0.000284 Jy/beam
! Combined flux in latest and established models = 1.27136 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 809 components and 1.27136 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061219Jy  sigma=0.770076
! Fit after  self-cal, rms=0.061218Jy  sigma=0.770059
! Inverting map 
! Added new window around map position (102, -99.6).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00021344 Jy
! Component: 100  -  total flux cleaned = -0.000574127 Jy
! Component: 150  -  total flux cleaned = -0.000670798 Jy
! Component: 200  -  total flux cleaned = -0.000671522 Jy
! Component: 250  -  total flux cleaned = -0.000671523 Jy
! Component: 300  -  total flux cleaned = -0.000725504 Jy
! Component: 350  -  total flux cleaned = -0.000778788 Jy
! Component: 400  -  total flux cleaned = -0.000831031 Jy
! Total flux subtracted in 400 components = -0.000831031 Jy
! Clean residual min=-0.000870 max=0.000869 Jy/beam
! Clean residual mean=-0.000021 rms=0.000258 Jy/beam
! Combined flux in latest and established models = 1.27053 Jy
! Performing phase self-cal
! Adding 47 model components to the UV plane model.
! The established model now contains 841 components and 1.27053 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061149Jy  sigma=0.769481
! Fit after  self-cal, rms=0.061148Jy  sigma=0.769458
! Inverting map 
! Added new window around map position (-88, 74.4).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000590636 Jy
! Component: 100  -  total flux cleaned = 0.00072229 Jy
! Component: 150  -  total flux cleaned = 0.000783408 Jy
! Component: 200  -  total flux cleaned = 0.000754236 Jy
! Component: 250  -  total flux cleaned = 0.000781458 Jy
! Component: 300  -  total flux cleaned = 0.000753837 Jy
! Component: 350  -  total flux cleaned = 0.000754407 Jy
! Component: 400  -  total flux cleaned = 0.00078192 Jy
! Total flux subtracted in 400 components = 0.00078192 Jy
! Clean residual min=-0.000833 max=0.000917 Jy/beam
! Clean residual mean=-0.000020 rms=0.000281 Jy/beam
! Combined flux in latest and established models = 1.27131 Jy
! Performing phase self-cal
! Adding 60 model components to the UV plane model.
! The established model now contains 880 components and 1.27131 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061092Jy  sigma=0.768845
! Fit after  self-cal, rms=0.061091Jy  sigma=0.768819
! Inverting map 
! Added new window around map position (102, 12.4).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 7.81854e-05 Jy
! Component: 100  -  total flux cleaned = 1.90219e-05 Jy
! Component: 150  -  total flux cleaned = 0.00012944 Jy
! Component: 200  -  total flux cleaned = 0.000182765 Jy
! Component: 250  -  total flux cleaned = 0.000235265 Jy
! Component: 300  -  total flux cleaned = 0.000235957 Jy
! Component: 350  -  total flux cleaned = 0.000210848 Jy
! Component: 400  -  total flux cleaned = 0.000185094 Jy
! Total flux subtracted in 400 components = 0.000185094 Jy
! Clean residual min=-0.000872 max=0.000800 Jy/beam
! Clean residual mean=-0.000018 rms=0.000289 Jy/beam
! Combined flux in latest and established models = 1.2715 Jy
! Performing phase self-cal
! Adding 53 model components to the UV plane model.
! The established model now contains 909 components and 1.2715 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.061043Jy  sigma=0.768335
! Fit after  self-cal, rms=0.061043Jy  sigma=0.768308
! 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, 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.061043Jy  sigma=0.768308
! Fit after  self-cal, rms=0.059190Jy  sigma=0.742550
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=11.6 mas, bmaj=47 mas, bpa=7.09 degrees
! Estimated noise=2.24429 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -3.24666e-05 Jy
! Component: 100  -  total flux cleaned = -9.51973e-05 Jy
! Component: 150  -  total flux cleaned = -7.00932e-05 Jy
! Component: 200  -  total flux cleaned = -9.37564e-05 Jy
! Component: 250  -  total flux cleaned = -9.3679e-05 Jy
! Component: 300  -  total flux cleaned = -7.07079e-05 Jy
! Component: 350  -  total flux cleaned = -0.000121266 Jy
! Component: 400  -  total flux cleaned = -0.000230089 Jy
! Total flux subtracted in 400 components = -0.000230089 Jy
! Clean residual min=-0.000970 max=0.000869 Jy/beam
! Clean residual mean=-0.000014 rms=0.000330 Jy/beam
! Combined flux in latest and established models = 1.27127 Jy
selfcal
! Performing phase self-cal
! Adding 86 model components to the UV plane model.
! The established model now contains 979 components and 1.27127 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.059166Jy  sigma=0.742215
! Fit after  self-cal, rms=0.059164Jy  sigma=0.742173
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.000283394 Jy
! Component: 100  -  total flux cleaned = 0.000375222 Jy
! Component: 150  -  total flux cleaned = 0.000464358 Jy
! Component: 200  -  total flux cleaned = 0.00057293 Jy
! Component: 250  -  total flux cleaned = 0.000658875 Jy
! Component: 300  -  total flux cleaned = 0.000743237 Jy
! Component: 350  -  total flux cleaned = 0.000742985 Jy
! Component: 400  -  total flux cleaned = 0.000721453 Jy
! Total flux subtracted in 400 components = 0.000721453 Jy
! Clean residual min=-0.001005 max=0.000934 Jy/beam
! Clean residual mean=-0.000022 rms=0.000325 Jy/beam
! Combined flux in latest and established models = 1.27199 Jy
! Performing phase self-cal
! Adding 76 model components to the UV plane model.
! The established model now contains 1029 components and 1.27199 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.059133Jy  sigma=0.741863
! Fit after  self-cal, rms=0.059132Jy  sigma=0.741848
! Inverting map 
! Added new window around map position (42.8, -80).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000146266 Jy
! Component: 100  -  total flux cleaned = 3.65046e-05 Jy
! Component: 150  -  total flux cleaned = 0.000187759 Jy
! Component: 200  -  total flux cleaned = 0.000186972 Jy
! Component: 250  -  total flux cleaned = 0.00016379 Jy
! Component: 300  -  total flux cleaned = 0.00014131 Jy
! Component: 350  -  total flux cleaned = 9.63157e-05 Jy
! Component: 400  -  total flux cleaned = 5.17033e-05 Jy
! Total flux subtracted in 400 components = 5.17033e-05 Jy
! Clean residual min=-0.000387 max=0.000366 Jy/beam
! Clean residual mean=-0.000024 rms=0.000184 Jy/beam
! Combined flux in latest and established models = 1.27204 Jy
! Performing phase self-cal
! Adding 109 model components to the UV plane model.
! The established model now contains 1114 components and 1.27204 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.059098Jy  sigma=0.741322
! Fit after  self-cal, rms=0.059095Jy  sigma=0.741293
! 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=4.208 mas, bmaj=29.24 mas, bpa=-10.95 degrees
! Estimated noise=2.0555 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.208 x 29.24 at -10.95 degrees (North through East)
! Clean map  min=-0.014416  max=1.117 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.059095Jy  sigma=0.741293
! Fit after  self-cal, rms=0.053031Jy  sigma=0.648044
! Inverting map and beam 
! Estimated beam: bmin=4.208 mas, bmaj=29.24 mas, bpa=-10.95 degrees
! Estimated noise=2.05561 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000688003 Jy
! Component: 100  -  total flux cleaned = 0.000955548 Jy
! Component: 150  -  total flux cleaned = 0.00101538 Jy
! Component: 200  -  total flux cleaned = 0.00112956 Jy
! Component: 250  -  total flux cleaned = 0.00129438 Jy
! Component: 300  -  total flux cleaned = 0.00140155 Jy
! Component: 350  -  total flux cleaned = 0.00148103 Jy
! Component: 400  -  total flux cleaned = 0.00158611 Jy
! Total flux subtracted in 400 components = 0.00158611 Jy
! Clean residual min=-0.000587 max=0.000660 Jy/beam
! Clean residual mean=-0.000014 rms=0.000192 Jy/beam
! Combined flux in latest and established models = 1.27363 Jy
! Performing phase self-cal
! Adding 68 model components to the UV plane model.
! The established model now contains 1168 components and 1.27363 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.052968Jy  sigma=0.647458
! Fit after  self-cal, rms=0.052965Jy  sigma=0.647412
! Inverting map 
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -4.62323e-05 Jy
! Component: 100  -  total flux cleaned = -0.000149639 Jy
! Component: 150  -  total flux cleaned = -0.000198254 Jy
! Component: 200  -  total flux cleaned = -0.000174491 Jy
! Component: 250  -  total flux cleaned = -0.000151683 Jy
! Component: 300  -  total flux cleaned = -0.000129358 Jy
! Component: 350  -  total flux cleaned = -0.000129413 Jy
! Component: 400  -  total flux cleaned = -0.000194549 Jy
! Total flux subtracted in 400 components = -0.000194549 Jy
! Clean residual min=-0.000485 max=0.000523 Jy/beam
! Clean residual mean=-0.000019 rms=0.000169 Jy/beam
! Combined flux in latest and established models = 1.27343 Jy
! Performing phase self-cal
! Adding 65 model components to the UV plane model.
! The established model now contains 1209 components and 1.27343 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.052923Jy  sigma=0.647001
! Fit after  self-cal, rms=0.052921Jy  sigma=0.646980
! 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.052921Jy  sigma=0.646980
! Fit after  self-cal, rms=0.052925Jy  sigma=0.646935
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.994 mas, bmaj=28.57 mas, bpa=-11.9 degrees
! Estimated noise=3.58403 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.88225 Jy
! Component: 100  -  total flux cleaned = 1.09487 Jy
! Component: 150  -  total flux cleaned = 1.18358 Jy
! Component: 200  -  total flux cleaned = 1.22687 Jy
! Component: 250  -  total flux cleaned = 1.2469 Jy
! Component: 300  -  total flux cleaned = 1.25989 Jy
! Component: 350  -  total flux cleaned = 1.27046 Jy
! Component: 400  -  total flux cleaned = 1.27702 Jy
! Total flux subtracted in 400 components = 1.27702 Jy
! Clean residual min=-0.018295 max=0.013235 Jy/beam
! Clean residual mean=0.000081 rms=0.004008 Jy/beam
! Combined flux in latest and established models = 1.27702 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 124 model components to the UV plane model.
! The established model now contains 124 components and 1.27702 Jy
! Inverting map and beam 
! Estimated beam: bmin=4.208 mas, bmaj=29.24 mas, bpa=-10.95 degrees
! Estimated noise=2.0556 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.00311935 Jy
! Component: 100  -  total flux cleaned = 0.00799429 Jy
! Component: 150  -  total flux cleaned = 0.0130346 Jy
! Component: 200  -  total flux cleaned = 0.0173278 Jy
! Component: 250  -  total flux cleaned = 0.0212722 Jy
! Component: 300  -  total flux cleaned = 0.024213 Jy
! Component: 350  -  total flux cleaned = 0.0268248 Jy
! Component: 400  -  total flux cleaned = 0.028417 Jy
! Total flux subtracted in 400 components = 0.028417 Jy
! Clean residual min=-0.008005 max=0.003461 Jy/beam
! Clean residual mean=0.000024 rms=0.001829 Jy/beam
! Combined flux in latest and established models = 1.30543 Jy
! Adding 148 model components to the UV plane model.
! The established model now contains 267 components and 1.30543 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=11.6 mas, bmaj=47.01 mas, bpa=7.09 degrees
! Estimated noise=2.24428 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00113998 Jy
! Component: 100  -  total flux cleaned = 0.00233463 Jy
! Component: 150  -  total flux cleaned = 0.00309723 Jy
! Component: 200  -  total flux cleaned = 0.00344151 Jy
! Component: 250  -  total flux cleaned = 0.00376886 Jy
! Component: 300  -  total flux cleaned = 0.0040734 Jy
! Component: 350  -  total flux cleaned = 0.00436727 Jy
! Component: 400  -  total flux cleaned = 0.00456138 Jy
! Total flux subtracted in 400 components = 0.00456138 Jy
! Clean residual min=-0.004790 max=0.002552 Jy/beam
! Clean residual mean=-0.000029 rms=0.001231 Jy/beam
! Combined flux in latest and established models = 1.31 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=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000198678 Jy
! Component: 100  -  total flux cleaned = 0.000299944 Jy
! Component: 150  -  total flux cleaned = 0.000639773 Jy
! Component: 200  -  total flux cleaned = 0.000760352 Jy
! Component: 250  -  total flux cleaned = 0.000886485 Jy
! Component: 300  -  total flux cleaned = 0.00145795 Jy
! Component: 350  -  total flux cleaned = 0.00162052 Jy
! Component: 400  -  total flux cleaned = 0.00196187 Jy
! Total flux subtracted in 400 components = 0.00196187 Jy
! Clean residual min=-0.006067 max=0.004859 Jy/beam
! Clean residual mean=-0.000002 rms=0.001516 Jy/beam
! Combined flux in latest and established models = 1.31196 Jy
! Adding 221 model components to the UV plane model.
! The established model now contains 480 components and 1.31196 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=4.208 mas, bmaj=29.24 mas, bpa=-10.95 degrees
! Estimated noise=2.0556 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00341788 Jy
! Component: 100  -  total flux cleaned = 0.00344054 Jy
! Component: 150  -  total flux cleaned = 0.00305569 Jy
! Component: 200  -  total flux cleaned = 0.0029405 Jy
! Component: 250  -  total flux cleaned = 0.00282798 Jy
! Component: 300  -  total flux cleaned = 0.00323845 Jy
! Component: 350  -  total flux cleaned = 0.00363626 Jy
! Component: 400  -  total flux cleaned = 0.00402356 Jy
! Total flux subtracted in 400 components = 0.00402356 Jy
! Clean residual min=-0.003786 max=0.002485 Jy/beam
! Clean residual mean=-0.000058 rms=0.000919 Jy/beam
! Combined flux in latest and established models = 1.31598 Jy
selfcal
! Performing phase self-cal
! Adding 69 model components to the UV plane model.
! The established model now contains 545 components and 1.31598 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.054212Jy  sigma=0.662976
! Fit after  self-cal, rms=0.054089Jy  sigma=0.661493
wmodel /scr/pima/rdv95_uvs/J0922-3959_S_map.mod
! Writing 545 model components to file: /scr/pima/rdv95_uvs/J0922-3959_S_map.mod
wobs   /scr/pima/rdv95_uvs/J0922-3959_S_uvs.fits
! Applying 450 buffered edits.
! Writing UV FITS file: /scr/pima/rdv95_uvs/J0922-3959_S_uvs.fits
wwins  /scr/pima/rdv95_uvs/J0922-3959_S_map.win
! wwins: Wrote 6 windows to /scr/pima/rdv95_uvs/J0922-3959_S_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.00142134 Jy
! Component: 100  -  total flux cleaned = -0.00248376 Jy
! Component: 150  -  total flux cleaned = -0.00334711 Jy
! Component: 200  -  total flux cleaned = -0.00399355 Jy
! Component: 250  -  total flux cleaned = -0.00446416 Jy
! Component: 300  -  total flux cleaned = -0.00484422 Jy
! Component: 350  -  total flux cleaned = -0.00520432 Jy
! Component: 400  -  total flux cleaned = -0.00552229 Jy
! Component: 450  -  total flux cleaned = -0.00572652 Jy
! Component: 500  -  total flux cleaned = -0.00594876 Jy
! Component: 550  -  total flux cleaned = -0.00609286 Jy
! Component: 600  -  total flux cleaned = -0.00613989 Jy
! Component: 650  -  total flux cleaned = -0.00616289 Jy
! Component: 700  -  total flux cleaned = -0.00611708 Jy
! Component: 750  -  total flux cleaned = -0.0060944 Jy
! Component: 800  -  total flux cleaned = -0.00600468 Jy
! Component: 850  -  total flux cleaned = -0.00596032 Jy
! Component: 900  -  total flux cleaned = -0.00587226 Jy
! Component: 950  -  total flux cleaned = -0.0057626 Jy
! Component: 1000  -  total flux cleaned = -0.00574082 Jy
! Component: 1050  -  total flux cleaned = -0.00571868 Jy
! Component: 1100  -  total flux cleaned = -0.00571853 Jy
! Component: 1150  -  total flux cleaned = -0.00571855 Jy
! Component: 1200  -  total flux cleaned = -0.00571853 Jy
! Component: 1250  -  total flux cleaned = -0.0056957 Jy
! Component: 1300  -  total flux cleaned = -0.00569551 Jy
! Component: 1350  -  total flux cleaned = -0.00569608 Jy
! Component: 1400  -  total flux cleaned = -0.00569643 Jy
! Component: 1450  -  total flux cleaned = -0.00569653 Jy
! Component: 1500  -  total flux cleaned = -0.00569701 Jy
! Component: 1550  -  total flux cleaned = -0.00569765 Jy
! Component: 1600  -  total flux cleaned = -0.00569762 Jy
! Component: 1650  -  total flux cleaned = -0.00569732 Jy
! Component: 1700  -  total flux cleaned = -0.00572968 Jy
! Component: 1750  -  total flux cleaned = -0.00573041 Jy
! Component: 1800  -  total flux cleaned = -0.00573084 Jy
! Component: 1850  -  total flux cleaned = -0.00572933 Jy
! Component: 1900  -  total flux cleaned = -0.00572967 Jy
! Component: 1950  -  total flux cleaned = -0.00576882 Jy
! Component: 2000  -  total flux cleaned = -0.00576987 Jy
! Component: 2050  -  total flux cleaned = -0.00577062 Jy
! Component: 2100  -  total flux cleaned = -0.00576993 Jy
! Component: 2150  -  total flux cleaned = -0.00576903 Jy
! Component: 2200  -  total flux cleaned = -0.00576882 Jy
! Component: 2250  -  total flux cleaned = -0.00576974 Jy
! Component: 2300  -  total flux cleaned = -0.00577106 Jy
! Component: 2350  -  total flux cleaned = -0.0057183 Jy
! Component: 2400  -  total flux cleaned = -0.00571689 Jy
! Component: 2450  -  total flux cleaned = -0.00571719 Jy
! Component: 2500  -  total flux cleaned = -0.00571722 Jy
! Component: 2550  -  total flux cleaned = -0.00565745 Jy
! Component: 2600  -  total flux cleaned = -0.00559354 Jy
! Component: 2650  -  total flux cleaned = -0.00559355 Jy
! Component: 2700  -  total flux cleaned = -0.00552404 Jy
! Component: 2750  -  total flux cleaned = -0.00552307 Jy
! Component: 2800  -  total flux cleaned = -0.00552153 Jy
! Component: 2850  -  total flux cleaned = -0.00552101 Jy
! Component: 2900  -  total flux cleaned = -0.00552201 Jy
! Component: 2950  -  total flux cleaned = -0.0055225 Jy
! Component: 3000  -  total flux cleaned = -0.00543759 Jy
! Component: 3050  -  total flux cleaned = -0.00543763 Jy
! Component: 3100  -  total flux cleaned = -0.00543765 Jy
! Component: 3150  -  total flux cleaned = -0.00543804 Jy
! Component: 3200  -  total flux cleaned = -0.00543608 Jy
! Component: 3250  -  total flux cleaned = -0.00543686 Jy
! Component: 3300  -  total flux cleaned = -0.00543553 Jy
! Component: 3350  -  total flux cleaned = -0.00532198 Jy
! Component: 3400  -  total flux cleaned = -0.00532057 Jy
! Component: 3450  -  total flux cleaned = -0.00531981 Jy
! Component: 3500  -  total flux cleaned = -0.00531852 Jy
! Component: 3550  -  total flux cleaned = -0.00531549 Jy
! Component: 3600  -  total flux cleaned = -0.00545175 Jy
! Component: 3650  -  total flux cleaned = -0.00531187 Jy
! Component: 3700  -  total flux cleaned = -0.00531084 Jy
! Component: 3750  -  total flux cleaned = -0.00530982 Jy
! Component: 3800  -  total flux cleaned = -0.00530847 Jy
! Component: 3850  -  total flux cleaned = -0.00530726 Jy
! Component: 3900  -  total flux cleaned = -0.00530663 Jy
! Component: 3950  -  total flux cleaned = -0.00530515 Jy
! Component: 4000  -  total flux cleaned = -0.00530532 Jy
! Component: 4050  -  total flux cleaned = -0.00530369 Jy
! Total flux subtracted in 4096 components = -0.00529944 Jy
! Clean residual min=-0.010901 max=0.011301 Jy/beam
! Clean residual mean=-0.000025 rms=0.002792 Jy/beam
! Combined flux in latest and established models = 1.31068 Jy
keep
! Adding 169 model components to the UV plane model.
! The established model now contains 714 components and 1.31068 Jy
![Exited script file: /home/lpetrov/bin/rdv_muppet]
wmap "/scr/pima/rdv95_uvs/J0922-3959_S\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.208 x 29.24 at -10.95 degrees (North through East)
! Clean map  min=-0.033672  max=1.1219 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv95_uvs/J0922-3959_S_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=  1.1219  rms=  0.00341  Jy 
quit
! Quitting program
! Log file /scr/pima/rdv95_uvs/J0922-3959_S_uvs.log closed on Wed Nov  6 20:45:19 2013