! Started logfile: /scr/pima/rdv62_uvs/J0735-1735_S_uvs.log on Fri Nov 22 22:33:46 2013
obs /scr/pima/rdv62_uvs/J0735-1735_S_uva.fits
! Reading UV FITS file: /scr/pima/rdv62_uvs/J0735-1735_S_uva.fits
! AN table 1: 53 integrations on 10 of 10 possible baselines.
! AN table 2: 95 integrations on 55 of 55 possible baselines.
! Apparent sampling: 0.96073 visibilities/baseline/integration-bin.
! Found source: J0735-1735
! 
! 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 22116 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/rdv62_uvs/J0735-1735_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:OV       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
! 2:BR       self-cal  status=correctable  weight=1
! 2:FD       self-cal  status=correctable  weight=1
! 2:HN       self-cal  status=correctable  weight=1
! 2:KK       self-cal  status=correctable  weight=1
! 2:KP       self-cal  status=correctable  weight=1
! 2:LA       self-cal  status=correctable  weight=1
! 2:MK       self-cal  status=correctable  weight=1
! 2:NL       self-cal  status=correctable  weight=1
! 2:OV       self-cal  status=correctable  weight=1
! 2:PT       self-cal  status=correctable  weight=1
! 2:SC       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 10 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 10 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 10 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 10 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=1.782482Jy  sigma=31.233451
! Fit after  self-cal, rms=0.824709Jy  sigma=15.462833
! 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 2:KK
! Flagging baselines of antenna 2: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=3.244 mas, bmaj=8.874 mas, bpa=13.75 degrees
! Estimated noise=1.44348 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.09694 Jy
! Component: 100  -  total flux cleaned = 1.52248 Jy
! Total flux subtracted in 100 components = 1.52248 Jy
! Clean residual min=-0.102962 max=0.233040 Jy/beam
! Clean residual mean=0.001089 rms=0.038725 Jy/beam
! Combined flux in latest and established models = 1.52248 Jy
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 42 components and 1.52248 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.513640Jy  sigma=9.935321
! Fit after  self-cal, rms=0.512362Jy  sigma=9.928595
! Inverting map 
! Added new window around map position (-5.2, -7.6).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.25157 Jy
! Component: 100  -  total flux cleaned = 0.405362 Jy
! Total flux subtracted in 100 components = 0.405362 Jy
! Clean residual min=-0.064235 max=0.154589 Jy/beam
! Clean residual mean=0.000755 rms=0.023912 Jy/beam
! Combined flux in latest and established models = 1.92784 Jy
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 87 components and 1.92784 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.368109Jy  sigma=6.740356
! Fit after  self-cal, rms=0.364631Jy  sigma=6.619981
! 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.829 mas, bmaj=12.57 mas, bpa=12.93 degrees
! Estimated noise=0.454653 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.200409 Jy
! Component: 100  -  total flux cleaned = 0.306388 Jy
! Component: 150  -  total flux cleaned = 0.365835 Jy
! Component: 200  -  total flux cleaned = 0.391981 Jy
! Total flux subtracted in 200 components = 0.391981 Jy
! Clean residual min=-0.069141 max=0.159149 Jy/beam
! Clean residual mean=0.000149 rms=0.027778 Jy/beam
! Combined flux in latest and established models = 2.31982 Jy
! Performing phase self-cal
! Adding 64 model components to the UV plane model.
! The established model now contains 140 components and 2.31983 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.305800Jy  sigma=4.905191
! Fit after  self-cal, rms=0.304746Jy  sigma=4.826538
! 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.92     KP        0.98     LA        1.12     OV        0.99     
!   PT        0.98     
! 
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        0.97     HN        0.98     KK        1.00*    
!   KP        0.99     LA        1.09     MK        1.20     NL        1.25     
!   OV        1.02     PT        1.01     SC        0.95     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.01     KP        1.02     LA        1.04     OV        1.05     
!   PT        0.99     
! 
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.04     FD        1.03     HN        1.11     KK        1.00*    
!   KP        1.01     LA        1.02     MK        1.24     NL        1.31     
!   OV        1.07     PT        1.00     SC        1.03     
! 
! 
! 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     KP        1.04     LA        0.84     OV        1.11     
!   PT        1.01     
! 
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.07     FD        1.04     HN        1.05     KK        1.00*    
!   KP        1.02     LA        0.81     MK        1.19     NL        1.40     
!   OV        1.13     PT        1.02     SC        1.05     
! 
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        1.00     KP        0.96     LA        0.75     OV        0.90     
!   PT        0.97     
! 
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.04     FD        1.01     HN        0.98     KK        1.00*    
!   KP        0.97     LA        0.73     MK        1.13     NL        1.29     
!   OV        0.94     PT        0.98     SC        1.15     
! 
! 
! Fit before self-cal, rms=0.304746Jy  sigma=4.826538
! Fit after  self-cal, rms=0.236287Jy  sigma=3.035956
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=4.98 mas, bmaj=13.26 mas, bpa=13.81 degrees
! Estimated noise=0.465608 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (14.4, 3.6).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.100546 Jy
! Component: 100  -  total flux cleaned = 0.111206 Jy
! Component: 150  -  total flux cleaned = 0.121333 Jy
! Component: 200  -  total flux cleaned = 0.131283 Jy
! Total flux subtracted in 200 components = 0.131283 Jy
! Clean residual min=-0.033483 max=0.045753 Jy/beam
! Clean residual mean=0.000240 rms=0.011004 Jy/beam
! Combined flux in latest and established models = 2.45111 Jy
! Performing phase self-cal
! Adding 40 model components to the UV plane model.
! The established model now contains 178 components and 2.45111 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.192209Jy  sigma=2.107777
! Fit after  self-cal, rms=0.187089Jy  sigma=2.020570
! 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.187089Jy  sigma=2.020570
! Fit after  self-cal, rms=0.178309Jy  sigma=1.939568
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=4.978 mas, bmaj=13.29 mas, bpa=13.94 degrees
! Estimated noise=0.460959 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0214541 Jy
! Component: 100  -  total flux cleaned = 0.0288023 Jy
! Component: 150  -  total flux cleaned = 0.0359315 Jy
! Component: 200  -  total flux cleaned = 0.0420988 Jy
! Total flux subtracted in 200 components = 0.0420988 Jy
! Clean residual min=-0.024353 max=0.028958 Jy/beam
! Clean residual mean=0.000085 rms=0.007002 Jy/beam
! Combined flux in latest and established models = 2.49321 Jy
selfcal
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 204 components and 2.49321 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.161326Jy  sigma=1.691743
! Fit after  self-cal, rms=0.160342Jy  sigma=1.675461
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.00162426 Jy
! Component: 100  -  total flux cleaned = 0.00463603 Jy
! Component: 150  -  total flux cleaned = 0.00673943 Jy
! Component: 200  -  total flux cleaned = 0.00867536 Jy
! Total flux subtracted in 200 components = 0.00867536 Jy
! Clean residual min=-0.022467 max=0.024008 Jy/beam
! Clean residual mean=0.000086 rms=0.005937 Jy/beam
! Combined flux in latest and established models = 2.50188 Jy
! Performing phase self-cal
! Adding 27 model components to the UV plane model.
! The established model now contains 225 components and 2.50188 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.151840Jy  sigma=1.595122
! Fit after  self-cal, rms=0.151021Jy  sigma=1.586057
! 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.151021Jy  sigma=1.586057
! Fit after  self-cal, rms=0.146023Jy  sigma=1.516766
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=5.003 mas, bmaj=13.45 mas, bpa=14.24 degrees
! Estimated noise=0.457757 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00368651 Jy
! Component: 100  -  total flux cleaned = 0.00651749 Jy
! Component: 150  -  total flux cleaned = 0.00903713 Jy
! Component: 200  -  total flux cleaned = 0.0112625 Jy
! Total flux subtracted in 200 components = 0.0112625 Jy
! Clean residual min=-0.017085 max=0.021124 Jy/beam
! Clean residual mean=0.000031 rms=0.004592 Jy/beam
! Combined flux in latest and established models = 2.51314 Jy
selfcal
! Performing phase self-cal
! Adding 24 model components to the UV plane model.
! The established model now contains 234 components and 2.51315 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.138255Jy  sigma=1.441260
! Fit after  self-cal, rms=0.137845Jy  sigma=1.434999
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.00274383 Jy
! Component: 100  -  total flux cleaned = 0.00471401 Jy
! Component: 150  -  total flux cleaned = 0.00649262 Jy
! Component: 200  -  total flux cleaned = 0.00689744 Jy
! Total flux subtracted in 200 components = 0.00689744 Jy
! Clean residual min=-0.016090 max=0.018521 Jy/beam
! Clean residual mean=0.000037 rms=0.004040 Jy/beam
! Combined flux in latest and established models = 2.52004 Jy
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 248 components and 2.52004 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.133410Jy  sigma=1.398980
! Fit after  self-cal, rms=0.133059Jy  sigma=1.394285
! 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.133059Jy  sigma=1.394285
! Fit after  self-cal, rms=0.132299Jy  sigma=1.378493
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.53 mas, bmaj=22.84 mas, bpa=13.75 degrees
! Estimated noise=0.394298 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00420442 Jy
! Component: 100  -  total flux cleaned = 0.00514203 Jy
! Component: 150  -  total flux cleaned = 0.0067719 Jy
! Component: 200  -  total flux cleaned = 0.0078229 Jy
! Total flux subtracted in 200 components = 0.0078229 Jy
! Clean residual min=-0.012140 max=0.012714 Jy/beam
! Clean residual mean=0.000013 rms=0.003472 Jy/beam
! Combined flux in latest and established models = 2.52787 Jy
selfcal
! Performing phase self-cal
! Adding 13 model components to the UV plane model.
! The established model now contains 258 components and 2.52787 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.130370Jy  sigma=1.343554
! Fit after  self-cal, rms=0.130165Jy  sigma=1.340275
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.000568623 Jy
! Component: 100  -  total flux cleaned = 0.00147083 Jy
! Component: 150  -  total flux cleaned = 0.00146433 Jy
! Component: 200  -  total flux cleaned = 0.00146639 Jy
! Total flux subtracted in 200 components = 0.00146639 Jy
! Clean residual min=-0.010597 max=0.010678 Jy/beam
! Clean residual mean=0.000025 rms=0.003069 Jy/beam
! Combined flux in latest and established models = 2.52933 Jy
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 273 components and 2.52933 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.129511Jy  sigma=1.325512
! Fit after  self-cal, rms=0.129351Jy  sigma=1.322727
! 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.
! Applying 2560 buffered edits.
! Applying 2560 buffered edits.
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:OV,true
! 1:OV       self-cal  status=fixed  weight=1
selfant 1:PT,true
! 1:PT       self-cal  status=fixed  weight=1
selfant 2:BR,true
! 2:BR       self-cal  status=fixed  weight=1
selfant 2:FD,true
! 2:FD       self-cal  status=fixed  weight=1
selfant 2:HN,true
! 2:HN       self-cal  status=fixed  weight=1
selfant 2:KP,true
! 2:KP       self-cal  status=fixed  weight=1
selfant 2:LA,true
! 2:LA       self-cal  status=fixed  weight=1
selfant 2:MK,true
! 2:MK       self-cal  status=fixed  weight=1
selfant 2:NL,true
! 2:NL       self-cal  status=fixed  weight=1
selfant 2:OV,true
! 2:OV       self-cal  status=fixed  weight=1
selfant 2:PT,true
! 2:PT       self-cal  status=fixed  weight=1
selfant 2:SC,true
! 2:SC       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.
!  A total of 95 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 95 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 95 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 95 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.129439Jy  sigma=1.326058
! Fit after  self-cal, rms=0.129439Jy  sigma=1.326058
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:OV       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
! 2:BR       self-cal  status=correctable  weight=1
! 2:FD       self-cal  status=correctable  weight=1
! 2:HN       self-cal  status=correctable  weight=1
! 2:KK       self-cal  status=correctable  weight=1
! 2:KP       self-cal  status=correctable  weight=1
! 2:LA       self-cal  status=correctable  weight=1
! 2:MK       self-cal  status=correctable  weight=1
! 2:NL       self-cal  status=correctable  weight=1
! 2:OV       self-cal  status=correctable  weight=1
! 2:PT       self-cal  status=correctable  weight=1
! 2:SC       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 13 telescope corrections were flagged in sub-array 1.
!  A total of 16 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 16 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 16 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 16 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.129439Jy  sigma=1.326058
! Fit after  self-cal, rms=0.128736Jy  sigma=1.318960
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=3.353 mas, bmaj=9.59 mas, bpa=15.39 degrees
! Estimated noise=1.48459 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.00297563 Jy
! Component: 100  -  total flux cleaned = 0.00452163 Jy
! Component: 150  -  total flux cleaned = 0.00592702 Jy
! Component: 200  -  total flux cleaned = 0.005947 Jy
! Total flux subtracted in 200 components = 0.005947 Jy
! Clean residual min=-0.014091 max=0.012938 Jy/beam
! Clean residual mean=0.000006 rms=0.003072 Jy/beam
! Combined flux in latest and established models = 2.53528 Jy
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 301 components and 2.53528 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.123692Jy  sigma=1.302651
! Fit after  self-cal, rms=0.123661Jy  sigma=1.302018
! 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=5.026 mas, bmaj=13.55 mas, bpa=14.36 degrees
! Estimated noise=0.455066 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.000485357 Jy
! Component: 100  -  total flux cleaned = -0.0011513 Jy
! Component: 150  -  total flux cleaned = -0.0028454 Jy
! Component: 200  -  total flux cleaned = -0.00328898 Jy
! Component: 250  -  total flux cleaned = -0.00393614 Jy
! Component: 300  -  total flux cleaned = -0.0041461 Jy
! Component: 350  -  total flux cleaned = -0.00454005 Jy
! Component: 400  -  total flux cleaned = -0.00473106 Jy
! Total flux subtracted in 400 components = -0.00473106 Jy
! Clean residual min=-0.009027 max=0.011115 Jy/beam
! Clean residual mean=-0.000003 rms=0.002259 Jy/beam
! Combined flux in latest and established models = 2.53055 Jy
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 325 components and 2.53055 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.120658Jy  sigma=1.286301
! Fit after  self-cal, rms=0.120435Jy  sigma=1.284271
! 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.000227295 Jy
! Component: 100  -  total flux cleaned = -0.000165291 Jy
! Component: 150  -  total flux cleaned = -0.000520025 Jy
! Component: 200  -  total flux cleaned = -0.000865274 Jy
! Component: 250  -  total flux cleaned = -0.00103293 Jy
! Component: 300  -  total flux cleaned = -0.00119928 Jy
! Component: 350  -  total flux cleaned = -0.00151838 Jy
! Component: 400  -  total flux cleaned = -0.00182483 Jy
! Total flux subtracted in 400 components = -0.00182483 Jy
! Clean residual min=-0.007976 max=0.010011 Jy/beam
! Clean residual mean=-0.000002 rms=0.002049 Jy/beam
! Combined flux in latest and established models = 2.52872 Jy
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 342 components and 2.52872 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.118358Jy  sigma=1.276036
! Fit after  self-cal, rms=0.118211Jy  sigma=1.274564
! 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.
! Applying 2560 buffered edits.
! Applying 2560 buffered edits.
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 13 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 13 telescope corrections were flagged in sub-array 1.
!  A total of 111 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.118322Jy  sigma=1.278269
! Fit after  self-cal, rms=0.118754Jy  sigma=1.272047
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=5.042 mas, bmaj=13.66 mas, bpa=14.55 degrees
! Estimated noise=0.454329 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00187489 Jy
! Component: 100  -  total flux cleaned = 0.00166413 Jy
! Component: 150  -  total flux cleaned = 0.00183688 Jy
! Component: 200  -  total flux cleaned = 0.00167454 Jy
! Component: 250  -  total flux cleaned = 0.00138005 Jy
! Component: 300  -  total flux cleaned = 0.000956399 Jy
! Component: 350  -  total flux cleaned = 0.000691775 Jy
! Component: 400  -  total flux cleaned = 0.000435212 Jy
! Total flux subtracted in 400 components = 0.000435212 Jy
! Clean residual min=-0.006946 max=0.009005 Jy/beam
! Clean residual mean=-0.000003 rms=0.001756 Jy/beam
! Combined flux in latest and established models = 2.52916 Jy
selfcal
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 359 components and 2.52916 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.117049Jy  sigma=1.263992
! Fit after  self-cal, rms=0.116908Jy  sigma=1.262961
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.000325834 Jy
! Component: 100  -  total flux cleaned = 6.65784e-05 Jy
! Component: 150  -  total flux cleaned = -0.000178349 Jy
! Component: 200  -  total flux cleaned = -0.000294964 Jy
! Component: 250  -  total flux cleaned = -0.000409359 Jy
! Component: 300  -  total flux cleaned = -0.000634199 Jy
! Component: 350  -  total flux cleaned = -0.000852162 Jy
! Component: 400  -  total flux cleaned = -0.000957813 Jy
! Total flux subtracted in 400 components = -0.000957813 Jy
! Clean residual min=-0.006683 max=0.008135 Jy/beam
! Clean residual mean=-0.000003 rms=0.001644 Jy/beam
! Combined flux in latest and established models = 2.5282 Jy
! Performing phase self-cal
! Adding 49 model components to the UV plane model.
! The established model now contains 377 components and 2.5282 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.115649Jy  sigma=1.259603
! Fit after  self-cal, rms=0.115537Jy  sigma=1.258810
! 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.115537Jy  sigma=1.258810
! Fit after  self-cal, rms=0.115781Jy  sigma=1.250309
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=5.051 mas, bmaj=13.73 mas, bpa=14.64 degrees
! Estimated noise=0.453909 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0010228 Jy
! Component: 100  -  total flux cleaned = 0.000889052 Jy
! Component: 150  -  total flux cleaned = 0.000524422 Jy
! Component: 200  -  total flux cleaned = 0.000186761 Jy
! Component: 250  -  total flux cleaned = -3.36832e-05 Jy
! Component: 300  -  total flux cleaned = -0.00034157 Jy
! Component: 350  -  total flux cleaned = -0.000635251 Jy
! Component: 400  -  total flux cleaned = -0.000727271 Jy
! Total flux subtracted in 400 components = -0.000727271 Jy
! Clean residual min=-0.006450 max=0.007619 Jy/beam
! Clean residual mean=-0.000005 rms=0.001476 Jy/beam
! Combined flux in latest and established models = 2.52747 Jy
selfcal
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 390 components and 2.52747 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.114679Jy  sigma=1.246657
! Fit after  self-cal, rms=0.114606Jy  sigma=1.246039
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.000208558 Jy
! Component: 100  -  total flux cleaned = -8.56223e-05 Jy
! Component: 150  -  total flux cleaned = -0.0002675 Jy
! Component: 200  -  total flux cleaned = -0.000353807 Jy
! Component: 250  -  total flux cleaned = -0.00043998 Jy
! Component: 300  -  total flux cleaned = -0.000524391 Jy
! Component: 350  -  total flux cleaned = -0.000608927 Jy
! Component: 400  -  total flux cleaned = -0.000692935 Jy
! Total flux subtracted in 400 components = -0.000692935 Jy
! Clean residual min=-0.006159 max=0.006997 Jy/beam
! Clean residual mean=-0.000007 rms=0.001415 Jy/beam
! Combined flux in latest and established models = 2.52678 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 400 components and 2.52678 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.113687Jy  sigma=1.244346
! Fit after  self-cal, rms=0.113611Jy  sigma=1.243864
! 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.113611Jy  sigma=1.243864
! Fit after  self-cal, rms=0.109211Jy  sigma=1.037206
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.48 mas, bmaj=22.76 mas, bpa=13.45 degrees
! Estimated noise=0.391817 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.42777e-05 Jy
! Component: 100  -  total flux cleaned = -7.63533e-05 Jy
! Component: 150  -  total flux cleaned = -0.000149378 Jy
! Component: 200  -  total flux cleaned = -0.000215959 Jy
! Component: 250  -  total flux cleaned = -0.000217339 Jy
! Component: 300  -  total flux cleaned = -0.000218653 Jy
! Component: 350  -  total flux cleaned = -0.000166294 Jy
! Component: 400  -  total flux cleaned = -0.000116461 Jy
! Total flux subtracted in 400 components = -0.000116461 Jy
! Clean residual min=-0.003918 max=0.004548 Jy/beam
! Clean residual mean=-0.000013 rms=0.001235 Jy/beam
! Combined flux in latest and established models = 2.52667 Jy
selfcal
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 421 components and 2.52667 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.108963Jy  sigma=1.035454
! Fit after  self-cal, rms=0.108931Jy  sigma=1.034947
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 (-53.6, 38).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00246519 Jy
! Component: 100  -  total flux cleaned = 0.00260327 Jy
! Component: 150  -  total flux cleaned = 0.00285173 Jy
! Component: 200  -  total flux cleaned = 0.00318309 Jy
! Component: 250  -  total flux cleaned = 0.00338562 Jy
! Component: 300  -  total flux cleaned = 0.00356855 Jy
! Component: 350  -  total flux cleaned = 0.00373427 Jy
! Component: 400  -  total flux cleaned = 0.00397537 Jy
! Total flux subtracted in 400 components = 0.00397537 Jy
! Clean residual min=-0.002974 max=0.002466 Jy/beam
! Clean residual mean=0.000001 rms=0.000784 Jy/beam
! Combined flux in latest and established models = 2.53064 Jy
! Performing phase self-cal
! Adding 63 model components to the UV plane model.
! The established model now contains 481 components and 2.53064 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.108526Jy  sigma=1.028355
! Fit after  self-cal, rms=0.108494Jy  sigma=1.028083
! 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=5.057 mas, bmaj=13.76 mas, bpa=14.69 degrees
! Estimated noise=0.453699 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 5.057 x 13.76 at 14.69 degrees (North through East)
! Clean map  min=-0.013176  max=1.3435 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.108494Jy  sigma=1.028083
! Fit after  self-cal, rms=0.104506Jy  sigma=0.870179
! Inverting map and beam 
! Estimated beam: bmin=5.059 mas, bmaj=13.78 mas, bpa=14.73 degrees
! Estimated noise=0.453775 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00165449 Jy
! Component: 100  -  total flux cleaned = 0.00240705 Jy
! Component: 150  -  total flux cleaned = 0.00254997 Jy
! Component: 200  -  total flux cleaned = 0.00281382 Jy
! Component: 250  -  total flux cleaned = 0.00268941 Jy
! Component: 300  -  total flux cleaned = 0.00281023 Jy
! Component: 350  -  total flux cleaned = 0.00292738 Jy
! Component: 400  -  total flux cleaned = 0.00315255 Jy
! Total flux subtracted in 400 components = 0.00315255 Jy
! Clean residual min=-0.004802 max=0.004641 Jy/beam
! Clean residual mean=0.000011 rms=0.000899 Jy/beam
! Combined flux in latest and established models = 2.53379 Jy
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 531 components and 2.53379 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.102917Jy  sigma=0.865060
! Fit after  self-cal, rms=0.102860Jy  sigma=0.864891
! Inverting map 
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000117081 Jy
! Component: 100  -  total flux cleaned = 0.000437571 Jy
! Component: 150  -  total flux cleaned = 0.00054381 Jy
! Component: 200  -  total flux cleaned = 0.000641023 Jy
! Component: 250  -  total flux cleaned = 0.000740301 Jy
! Component: 300  -  total flux cleaned = 0.00102483 Jy
! Component: 350  -  total flux cleaned = 0.000929792 Jy
! Component: 400  -  total flux cleaned = 0.00102126 Jy
! Total flux subtracted in 400 components = 0.00102126 Jy
! Clean residual min=-0.004454 max=0.004152 Jy/beam
! Clean residual mean=0.000014 rms=0.000834 Jy/beam
! Combined flux in latest and established models = 2.53482 Jy
! Performing phase self-cal
! Adding 67 model components to the UV plane model.
! The established model now contains 568 components and 2.53481 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101936Jy  sigma=0.861990
! Fit after  self-cal, rms=0.101908Jy  sigma=0.861894
! 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.101908Jy  sigma=0.861894
! Fit after  self-cal, rms=0.101998Jy  sigma=0.861588
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=3.38 mas, bmaj=9.837 mas, bpa=15.92 degrees
! Estimated noise=1.48558 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.11287 Jy
! Component: 100  -  total flux cleaned = 1.53589 Jy
! Component: 150  -  total flux cleaned = 1.77333 Jy
! Component: 200  -  total flux cleaned = 1.93508 Jy
! Component: 250  -  total flux cleaned = 2.05761 Jy
! Component: 300  -  total flux cleaned = 2.15281 Jy
! Component: 350  -  total flux cleaned = 2.22748 Jy
! Component: 400  -  total flux cleaned = 2.27989 Jy
! Total flux subtracted in 400 components = 2.27989 Jy
! Clean residual min=-0.034995 max=0.043390 Jy/beam
! Clean residual mean=0.000174 rms=0.008924 Jy/beam
! Combined flux in latest and established models = 2.27989 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 166 model components to the UV plane model.
! The established model now contains 166 components and 2.27989 Jy
! Inverting map and beam 
! Estimated beam: bmin=5.058 mas, bmaj=13.79 mas, bpa=14.74 degrees
! Estimated noise=0.453763 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.0812707 Jy
! Component: 100  -  total flux cleaned = 0.132846 Jy
! Component: 150  -  total flux cleaned = 0.16824 Jy
! Component: 200  -  total flux cleaned = 0.188964 Jy
! Component: 250  -  total flux cleaned = 0.198293 Jy
! Component: 300  -  total flux cleaned = 0.205692 Jy
! Component: 350  -  total flux cleaned = 0.212209 Jy
! Component: 400  -  total flux cleaned = 0.218562 Jy
! Total flux subtracted in 400 components = 0.218562 Jy
! Clean residual min=-0.010393 max=0.015142 Jy/beam
! Clean residual mean=-0.000009 rms=0.003270 Jy/beam
! Combined flux in latest and established models = 2.49845 Jy
! Adding 114 model components to the UV plane model.
! The established model now contains 268 components and 2.49845 Jy
! Inverting map 
! Added new window around map position (-12, -10.4).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00768001 Jy
! Component: 100  -  total flux cleaned = 0.0102802 Jy
! Component: 150  -  total flux cleaned = 0.0125804 Jy
! Component: 200  -  total flux cleaned = 0.0151133 Jy
! Component: 250  -  total flux cleaned = 0.0174489 Jy
! Component: 300  -  total flux cleaned = 0.0188932 Jy
! Component: 350  -  total flux cleaned = 0.0202598 Jy
! Component: 400  -  total flux cleaned = 0.0215294 Jy
! Total flux subtracted in 400 components = 0.0215294 Jy
! Clean residual min=-0.006484 max=0.007205 Jy/beam
! Clean residual mean=-0.000001 rms=0.001870 Jy/beam
! Combined flux in latest and established models = 2.51998 Jy
! Adding 73 model components to the UV plane model.
! The established model now contains 332 components and 2.51998 Jy
! Inverting map 
! Added new window around map position (4.8, 11.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000854134 Jy
! Component: 100  -  total flux cleaned = 3.9254e-05 Jy
! Component: 150  -  total flux cleaned = 3.57985e-05 Jy
! Component: 200  -  total flux cleaned = 0.000557988 Jy
! Component: 250  -  total flux cleaned = 0.000560911 Jy
! Component: 300  -  total flux cleaned = 0.000804538 Jy
! Component: 350  -  total flux cleaned = 0.00103027 Jy
! Component: 400  -  total flux cleaned = 0.00103365 Jy
! Total flux subtracted in 400 components = 0.00103365 Jy
! Clean residual min=-0.005149 max=0.005170 Jy/beam
! Clean residual mean=-0.000026 rms=0.001398 Jy/beam
! Combined flux in latest and established models = 2.52102 Jy
! Adding 63 model components to the UV plane model.
! The established model now contains 377 components and 2.52102 Jy
! Inverting map 
! Added new window around map position (-4.8, 62.4).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00350723 Jy
! Component: 100  -  total flux cleaned = 0.00370616 Jy
! Component: 150  -  total flux cleaned = 0.00371562 Jy
! Component: 200  -  total flux cleaned = 0.00353744 Jy
! Component: 250  -  total flux cleaned = 0.00371443 Jy
! Component: 300  -  total flux cleaned = 0.00371475 Jy
! Component: 350  -  total flux cleaned = 0.00371237 Jy
! Component: 400  -  total flux cleaned = 0.00338991 Jy
! Total flux subtracted in 400 components = 0.00338991 Jy
! Clean residual min=-0.004325 max=0.004101 Jy/beam
! Clean residual mean=-0.000016 rms=0.001107 Jy/beam
! Combined flux in latest and established models = 2.52441 Jy
! Adding 66 model components to the UV plane model.
! The established model now contains 417 components and 2.52441 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=10.48 mas, bmaj=22.77 mas, bpa=13.46 degrees
! Estimated noise=0.391838 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00185063 Jy
! Component: 100  -  total flux cleaned = 0.00229514 Jy
! Component: 150  -  total flux cleaned = 0.00258656 Jy
! Component: 200  -  total flux cleaned = 0.0028615 Jy
! Component: 250  -  total flux cleaned = 0.00321056 Jy
! Component: 300  -  total flux cleaned = 0.00329149 Jy
! Component: 350  -  total flux cleaned = 0.00344635 Jy
! Component: 400  -  total flux cleaned = 0.00344471 Jy
! Total flux subtracted in 400 components = 0.00344471 Jy
! Clean residual min=-0.002400 max=0.002504 Jy/beam
! Clean residual mean=-0.000007 rms=0.000700 Jy/beam
! Combined flux in latest and established models = 2.52785 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 (95.2, 60).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000272909 Jy
! Component: 100  -  total flux cleaned = -9.73462e-05 Jy
! Component: 150  -  total flux cleaned = -9.91682e-05 Jy
! Component: 200  -  total flux cleaned = -3.91319e-05 Jy
! Component: 250  -  total flux cleaned = 2.46094e-05 Jy
! Component: 300  -  total flux cleaned = 0.000331394 Jy
! Component: 350  -  total flux cleaned = 0.000448959 Jy
! Component: 400  -  total flux cleaned = 0.000676974 Jy
! Total flux subtracted in 400 components = 0.000676974 Jy
! Clean residual min=-0.001893 max=0.001807 Jy/beam
! Clean residual mean=-0.000009 rms=0.000541 Jy/beam
! Combined flux in latest and established models = 2.52853 Jy
! Adding 131 model components to the UV plane model.
! The established model now contains 541 components and 2.52853 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=5.058 mas, bmaj=13.79 mas, bpa=14.74 degrees
! Estimated noise=0.453763 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00103426 Jy
! Component: 100  -  total flux cleaned = -0.00076544 Jy
! Component: 150  -  total flux cleaned = -0.000501549 Jy
! Component: 200  -  total flux cleaned = -0.000250391 Jy
! Component: 250  -  total flux cleaned = -0.000248743 Jy
! Component: 300  -  total flux cleaned = -0.000244468 Jy
! Component: 350  -  total flux cleaned = -0.000132591 Jy
! Component: 400  -  total flux cleaned = -0.000243177 Jy
! Total flux subtracted in 400 components = -0.000243177 Jy
! Clean residual min=-0.002746 max=0.002896 Jy/beam
! Clean residual mean=-0.000001 rms=0.000719 Jy/beam
! Combined flux in latest and established models = 2.52828 Jy
selfcal
! Performing phase self-cal
! Adding 50 model components to the UV plane model.
! The established model now contains 580 components and 2.52828 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.098263Jy  sigma=0.857179
! Fit after  self-cal, rms=0.097748Jy  sigma=0.856032
wmodel /scr/pima/rdv62_uvs/J0735-1735_S_map.mod
! Writing 580 model components to file: /scr/pima/rdv62_uvs/J0735-1735_S_map.mod
wobs   /scr/pima/rdv62_uvs/J0735-1735_S_uvs.fits
! Applying 1365 buffered edits.
! Writing UV FITS file: /scr/pima/rdv62_uvs/J0735-1735_S_uvs.fits
wwins  /scr/pima/rdv62_uvs/J0735-1735_S_map.win
! wwins: Wrote 8 windows to /scr/pima/rdv62_uvs/J0735-1735_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.000772951 Jy
! Component: 100  -  total flux cleaned = 0.00110561 Jy
! Component: 150  -  total flux cleaned = 0.00124235 Jy
! Component: 200  -  total flux cleaned = 0.00141676 Jy
! Component: 250  -  total flux cleaned = 0.00158674 Jy
! Component: 300  -  total flux cleaned = 0.00191674 Jy
! Component: 350  -  total flux cleaned = 0.00223784 Jy
! Component: 400  -  total flux cleaned = 0.00235599 Jy
! Component: 450  -  total flux cleaned = 0.00239553 Jy
! Component: 500  -  total flux cleaned = 0.00235782 Jy
! Component: 550  -  total flux cleaned = 0.00235807 Jy
! Component: 600  -  total flux cleaned = 0.00228454 Jy
! Component: 650  -  total flux cleaned = 0.00217581 Jy
! Component: 700  -  total flux cleaned = 0.00206812 Jy
! Component: 750  -  total flux cleaned = 0.00196149 Jy
! Component: 800  -  total flux cleaned = 0.00171614 Jy
! Component: 850  -  total flux cleaned = 0.00161269 Jy
! Component: 900  -  total flux cleaned = 0.00140713 Jy
! Component: 950  -  total flux cleaned = 0.00137319 Jy
! Component: 1000  -  total flux cleaned = 0.00113786 Jy
! Component: 1050  -  total flux cleaned = 0.00107114 Jy
! Component: 1100  -  total flux cleaned = 0.000906693 Jy
! Component: 1150  -  total flux cleaned = 0.000874137 Jy
! Component: 1200  -  total flux cleaned = 0.000777293 Jy
! Component: 1250  -  total flux cleaned = 0.000745081 Jy
! Component: 1300  -  total flux cleaned = 0.000554851 Jy
! Component: 1350  -  total flux cleaned = 0.000429503 Jy
! Component: 1400  -  total flux cleaned = 0.000336343 Jy
! Component: 1450  -  total flux cleaned = 0.000398669 Jy
! Component: 1500  -  total flux cleaned = 0.00030658 Jy
! Component: 1550  -  total flux cleaned = 0.000427653 Jy
! Component: 1600  -  total flux cleaned = 0.000277121 Jy
! Component: 1650  -  total flux cleaned = 0.000277015 Jy
! Component: 1700  -  total flux cleaned = 0.000277195 Jy
! Component: 1750  -  total flux cleaned = 0.00033586 Jy
! Component: 1800  -  total flux cleaned = 0.000365275 Jy
! Component: 1850  -  total flux cleaned = 0.000423205 Jy
! Component: 1900  -  total flux cleaned = 0.000480893 Jy
! Component: 1950  -  total flux cleaned = 0.000538062 Jy
! Component: 2000  -  total flux cleaned = 0.000594999 Jy
! Component: 2050  -  total flux cleaned = 0.000708413 Jy
! Component: 2100  -  total flux cleaned = 0.000792507 Jy
! Component: 2150  -  total flux cleaned = 0.000903872 Jy
! Component: 2200  -  total flux cleaned = 0.000876376 Jy
! Component: 2250  -  total flux cleaned = 0.00101395 Jy
! Component: 2300  -  total flux cleaned = 0.00101385 Jy
! Component: 2350  -  total flux cleaned = 0.00112224 Jy
! Component: 2400  -  total flux cleaned = 0.00106823 Jy
! Component: 2450  -  total flux cleaned = 0.00120218 Jy
! Component: 2500  -  total flux cleaned = 0.00122914 Jy
! Component: 2550  -  total flux cleaned = 0.00130835 Jy
! Component: 2600  -  total flux cleaned = 0.00130821 Jy
! Component: 2650  -  total flux cleaned = 0.00133438 Jy
! Component: 2700  -  total flux cleaned = 0.00123078 Jy
! Component: 2750  -  total flux cleaned = 0.00133446 Jy
! Component: 2800  -  total flux cleaned = 0.00130874 Jy
! Component: 2850  -  total flux cleaned = 0.00130905 Jy
! Component: 2900  -  total flux cleaned = 0.00128402 Jy
! Component: 2950  -  total flux cleaned = 0.00138487 Jy
! Component: 3000  -  total flux cleaned = 0.00140985 Jy
! Component: 3050  -  total flux cleaned = 0.00133499 Jy
! Component: 3100  -  total flux cleaned = 0.00143408 Jy
! Component: 3150  -  total flux cleaned = 0.00143398 Jy
! Component: 3200  -  total flux cleaned = 0.00140955 Jy
! Component: 3250  -  total flux cleaned = 0.00148281 Jy
! Component: 3300  -  total flux cleaned = 0.00153179 Jy
! Component: 3350  -  total flux cleaned = 0.0016283 Jy
! Component: 3400  -  total flux cleaned = 0.00158004 Jy
! Component: 3450  -  total flux cleaned = 0.00160382 Jy
! Component: 3500  -  total flux cleaned = 0.00160401 Jy
! Component: 3550  -  total flux cleaned = 0.0016511 Jy
! Component: 3600  -  total flux cleaned = 0.00160408 Jy
! Component: 3650  -  total flux cleaned = 0.00162751 Jy
! Component: 3700  -  total flux cleaned = 0.00167386 Jy
! Component: 3750  -  total flux cleaned = 0.00167382 Jy
! Component: 3800  -  total flux cleaned = 0.00167378 Jy
! Component: 3850  -  total flux cleaned = 0.00176532 Jy
! Component: 3900  -  total flux cleaned = 0.00169673 Jy
! Component: 3950  -  total flux cleaned = 0.00181004 Jy
! Component: 4000  -  total flux cleaned = 0.0018776 Jy
! Component: 4050  -  total flux cleaned = 0.00185516 Jy
! Total flux subtracted in 4096 components = 0.00187757 Jy
! Clean residual min=-0.001121 max=0.001119 Jy/beam
! Clean residual mean=0.000005 rms=0.000422 Jy/beam
! Combined flux in latest and established models = 2.53016 Jy
keep
! Adding 179 model components to the UV plane model.
! The established model now contains 740 components and 2.53016 Jy
![Exited script file: /home/lpetrov/bin/rdv_muppet]
wmap "/scr/pima/rdv62_uvs/J0735-1735_S\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 5.058 x 13.79 at 14.74 degrees (North through East)
! Clean map  min=-0.010309  max=1.3331 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J0735-1735_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.3331  rms=  0.00042  Jy 
quit
! Quitting program
! Log file /scr/pima/rdv62_uvs/J0735-1735_S_uvs.log closed on Fri Nov 22 22:33:56 2013