! Started logfile: /scr/pima/rdv62_uvs/J2206-1835_X_uvs.log on Fri Nov 22 22:42:40 2013
obs /scr/pima/rdv62_uvs/J2206-1835_X_uva.fits
! Reading UV FITS file: /scr/pima/rdv62_uvs/J2206-1835_X_uva.fits
! AN table 1: 37 integrations on 10 of 10 possible baselines.
! AN table 2: 12 integrations on 3 of 3 possible baselines.
! AN table 3: 12 integrations on 1 of 1 possible baselines.
! AN table 4: 12 integrations on 1 of 1 possible baselines.
! AN table 5: 32 integrations on 36 of 36 possible baselines.
! Apparent sampling: 0.3622 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J2206-1835
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  8.40599e+09        8e+06          1        8e+06
!  02        2  8.47599e+09        8e+06          1        8e+06
!  03        3  8.79099e+09        8e+06          1        8e+06
!  04        4  8.89599e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 0 lines of history.
! 
! Reading 2292 visibilities.
select RR
! Selecting polarization: RR,  channels: 1..4
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
print "field_size =",field_size
! field_size = 1024 
print "field_cell =",field_cell
! field_cell = 0.2 
print "taper_size =",taper_size
! taper_size = 70 
![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J2206-1835_X]
integer clean_niter
clean_niter = 100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float soltime1
soltime1 = 800000
float soltime2
soltime2 = 15
float soltime3
soltime3 = 1
float thresh
thresh = 0.5
float win_mult
win_mult = 1.8
float old_peak
float new_peak
float flux_cutoff
selfant "",false
! 1:FD       self-cal  status=correctable  weight=1
! 1:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:NL       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
! 2:KP       self-cal  status=correctable  weight=1
! 2:LA       self-cal  status=correctable  weight=1
! 2:PT       self-cal  status=correctable  weight=1
! 3:NL       self-cal  status=correctable  weight=1
! 3:SC       self-cal  status=correctable  weight=1
! 4:MK       self-cal  status=correctable  weight=1
! 4:OV       self-cal  status=correctable  weight=1
! 5:BR       self-cal  status=correctable  weight=1
! 5:FD       self-cal  status=correctable  weight=1
! 5:HN       self-cal  status=correctable  weight=1
! 5:KP       self-cal  status=correctable  weight=1
! 5:LA       self-cal  status=correctable  weight=1
! 5:MK       self-cal  status=correctable  weight=1
! 5:NL       self-cal  status=correctable  weight=1
! 5:OV       self-cal  status=correctable  weight=1
! 5:PT       self-cal  status=correctable  weight=1
system "rm /tmp/ff.ff"
mapsize field_size, field_cell
! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize.
startmod
! Applying default point source starting model.
! Performing phase self-cal
! Adding 1 model components to the UV plane model.
! The established model now contains 1 components and 1 Jy
! 
! Correcting IF 1.
!  A total of 74 telescope corrections were flagged in sub-array 1.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 2.
!  A total of 74 telescope corrections were flagged in sub-array 1.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 3.
!  A total of 74 telescope corrections were flagged in sub-array 1.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 4.
!  A total of 74 telescope corrections were flagged in sub-array 1.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Fit before self-cal, rms=0.997713Jy  sigma=3.080485
! Fit after  self-cal, rms=0.823686Jy  sigma=1.759064
! 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=9.898 mas, bmaj=19.74 mas, bpa=-14.07 degrees
! Estimated noise=12.4785 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.352234 Jy
! Component: 100  -  total flux cleaned = 0.405297 Jy
! Total flux subtracted in 100 components = 0.405297 Jy
! Clean residual min=-0.066136 max=0.066511 Jy/beam
! Clean residual mean=0.000488 rms=0.041257 Jy/beam
! Combined flux in latest and established models = 0.405297 Jy
! Performing phase self-cal
! Adding 82 model components to the UV plane model.
! The established model now contains 82 components and 0.405297 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.126140Jy  sigma=0.390619
! Fit after  self-cal, rms=0.123917Jy  sigma=0.384042
! 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=10.11 mas, bmaj=22.32 mas, bpa=-20.53 degrees
! Estimated noise=7.23397 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.0484289 Jy
! Component: 100  -  total flux cleaned = -0.0559724 Jy
! Component: 150  -  total flux cleaned = -0.0793693 Jy
! Component: 200  -  total flux cleaned = -0.0918981 Jy
! Total flux subtracted in 200 components = -0.0918981 Jy
! Clean residual min=-0.023429 max=0.023019 Jy/beam
! Clean residual mean=0.000194 rms=0.008524 Jy/beam
! Combined flux in latest and established models = 0.313399 Jy
! Performing phase self-cal
! Adding 72 model components to the UV plane model.
! The established model now contains 154 components and 0.313399 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.087490Jy  sigma=0.273672
! Fit after  self-cal, rms=0.084230Jy  sigma=0.272280
! 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 111 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    KP        0.00*    LA        0.00*    NL        0.00*    
!   PT        0.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    LA        0.00*    PT        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   NL        0.00*    SC        0.00*    
! 
!  Telescope amplitude corrections in sub-array 4:
!   MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 5:
!   BR        1.00*    FD        0.85     HN        1.00*    KP        1.44     
!   LA        0.97     MK        1.00*    NL        1.00*    OV        0.73     
!   PT        0.91     
! 
! 
! Correcting IF 2.
!  A total of 111 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    KP        0.00*    LA        0.00*    NL        0.00*    
!   PT        0.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    LA        0.00*    PT        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   NL        0.00*    SC        0.00*    
! 
!  Telescope amplitude corrections in sub-array 4:
!   MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 5:
!   BR        1.00*    FD        0.91     HN        1.00*    KP        1.59     
!   LA        1.00     MK        1.00*    NL        1.00*    OV        0.45     
!   PT        0.93     
! 
! 
! Correcting IF 3.
!  A total of 111 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    KP        0.00*    LA        0.00*    NL        0.00*    
!   PT        0.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    LA        0.00*    PT        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   NL        0.00*    SC        0.00*    
! 
!  Telescope amplitude corrections in sub-array 4:
!   MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 5:
!   BR        1.00*    FD        0.92     HN        1.00*    KP        0.95     
!   LA        1.08     MK        1.00*    NL        1.00*    OV        0.99     
!   PT        0.97     
! 
! 
! Correcting IF 4.
!  A total of 111 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    KP        0.00*    LA        0.00*    NL        0.00*    
!   PT        0.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    LA        0.00*    PT        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   NL        0.00*    SC        0.00*    
! 
!  Telescope amplitude corrections in sub-array 4:
!   MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 5:
!   BR        1.00*    FD        1.01     HN        1.00*    KP        0.90     
!   LA        1.06     MK        1.00*    NL        1.00*    OV        1.04     
!   PT        0.95     
! 
! 
! Fit before self-cal, rms=0.084230Jy  sigma=0.272280
! Fit after  self-cal, rms=0.070962Jy  sigma=0.107404
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=7.794 mas, bmaj=13.7 mas, bpa=-1.77 degrees
! Estimated noise=16.8067 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.00314827 Jy
! Component: 100  -  total flux cleaned = 0.00236688 Jy
! Component: 150  -  total flux cleaned = 0.00799303 Jy
! Component: 200  -  total flux cleaned = 0.0135513 Jy
! Total flux subtracted in 200 components = 0.0135513 Jy
! Clean residual min=-0.004728 max=0.004697 Jy/beam
! Clean residual mean=-0.000010 rms=0.002095 Jy/beam
! Combined flux in latest and established models = 0.326951 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 200 components and 0.326951 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062140Jy  sigma=0.105952
! Fit after  self-cal, rms=0.062773Jy  sigma=0.105851
! 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.062773Jy  sigma=0.105851
! Fit after  self-cal, rms=0.063995Jy  sigma=0.105705
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.808 mas, bmaj=13.84 mas, bpa=-2.256 degrees
! Estimated noise=17.0533 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00155254 Jy
! Component: 100  -  total flux cleaned = 0.00139753 Jy
! Component: 150  -  total flux cleaned = 0.00514294 Jy
! Component: 200  -  total flux cleaned = 0.00795219 Jy
! Total flux subtracted in 200 components = 0.00795219 Jy
! Clean residual min=-0.002654 max=0.002679 Jy/beam
! Clean residual mean=-0.000003 rms=0.001226 Jy/beam
! Combined flux in latest and established models = 0.334903 Jy
selfcal
! Performing phase self-cal
! Adding 36 model components to the UV plane model.
! The established model now contains 225 components and 0.334903 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060798Jy  sigma=0.105334
! Fit after  self-cal, rms=0.060996Jy  sigma=0.105317
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.00229409 Jy
! Component: 100  -  total flux cleaned = 0.00544213 Jy
! Component: 150  -  total flux cleaned = 0.00831718 Jy
! Component: 200  -  total flux cleaned = 0.0105519 Jy
! Total flux subtracted in 200 components = 0.0105519 Jy
! Clean residual min=-0.001930 max=0.001930 Jy/beam
! Clean residual mean=0.000001 rms=0.000850 Jy/beam
! Combined flux in latest and established models = 0.345455 Jy
! Performing phase self-cal
! Adding 49 model components to the UV plane model.
! The established model now contains 272 components and 0.345455 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.059808Jy  sigma=0.105171
! Fit after  self-cal, rms=0.059861Jy  sigma=0.105154
! 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.059861Jy  sigma=0.105154
! Fit after  self-cal, rms=0.060668Jy  sigma=0.105113
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.885 mas, bmaj=13.91 mas, bpa=-2.98 degrees
! Estimated noise=17.2634 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000925514 Jy
! Component: 100  -  total flux cleaned = -0.00101114 Jy
! Component: 150  -  total flux cleaned = -0.0016542 Jy
! Component: 200  -  total flux cleaned = -0.00259102 Jy
! Total flux subtracted in 200 components = -0.00259102 Jy
! Clean residual min=-0.001035 max=0.001027 Jy/beam
! Clean residual mean=0.000003 rms=0.000545 Jy/beam
! Combined flux in latest and established models = 0.342864 Jy
selfcal
! Performing phase self-cal
! Adding 53 model components to the UV plane model.
! The established model now contains 321 components and 0.342864 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060226Jy  sigma=0.105066
! Fit after  self-cal, rms=0.060250Jy  sigma=0.105057
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.000112108 Jy
! Component: 100  -  total flux cleaned = -8.74852e-05 Jy
! Component: 150  -  total flux cleaned = -0.000864099 Jy
! Component: 200  -  total flux cleaned = -0.000471997 Jy
! Total flux subtracted in 200 components = -0.000471997 Jy
! Clean residual min=-0.000820 max=0.000838 Jy/beam
! Clean residual mean=-0.000001 rms=0.000395 Jy/beam
! Combined flux in latest and established models = 0.342392 Jy
! Performing phase self-cal
! Adding 67 model components to the UV plane model.
! The established model now contains 388 components and 0.342392 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060058Jy  sigma=0.105036
! Fit after  self-cal, rms=0.060070Jy  sigma=0.105035
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime2
! Performing amp+phase self-cal over 15 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060070Jy  sigma=0.105035
! Fit after  self-cal, rms=0.060549Jy  sigma=0.104997
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.063 mas, bmaj=14.03 mas, bpa=-3.896 degrees
! Estimated noise=16.9166 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00054184 Jy
! Component: 100  -  total flux cleaned = -0.00128614 Jy
! Component: 150  -  total flux cleaned = -0.00184336 Jy
! Component: 200  -  total flux cleaned = -0.00218274 Jy
! Total flux subtracted in 200 components = -0.00218274 Jy
! Clean residual min=-0.000565 max=0.000561 Jy/beam
! Clean residual mean=-0.000002 rms=0.000268 Jy/beam
! Combined flux in latest and established models = 0.340209 Jy
selfcal
! Performing phase self-cal
! Adding 54 model components to the UV plane model.
! The established model now contains 437 components and 0.340209 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060440Jy  sigma=0.104979
! Fit after  self-cal, rms=0.060440Jy  sigma=0.104978
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.000250536 Jy
! Component: 100  -  total flux cleaned = 0.000791773 Jy
! Component: 150  -  total flux cleaned = 0.00109875 Jy
! Component: 200  -  total flux cleaned = 0.00127494 Jy
! Total flux subtracted in 200 components = 0.00127494 Jy
! Clean residual min=-0.000630 max=0.000627 Jy/beam
! Clean residual mean=0.000003 rms=0.000228 Jy/beam
! Combined flux in latest and established models = 0.341484 Jy
! Performing phase self-cal
! Adding 47 model components to the UV plane model.
! The established model now contains 478 components and 0.341484 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060375Jy  sigma=0.104974
! Fit after  self-cal, rms=0.060378Jy  sigma=0.104970
! 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:NL,true
! 1:NL       self-cal  status=fixed  weight=1
selfant 1:PT,true
! 1:PT       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:PT,true
! 2:PT       self-cal  status=fixed  weight=1
selfant 3:NL,true
! 3:NL       self-cal  status=fixed  weight=1
selfant 3:SC,true
! 3:SC       self-cal  status=fixed  weight=1
selfant 4:MK,true
! 4:MK       self-cal  status=fixed  weight=1
selfant 4:OV,true
! 4:OV       self-cal  status=fixed  weight=1
selfant 5:BR,true
! 5:BR       self-cal  status=fixed  weight=1
selfant 5:FD,true
! 5:FD       self-cal  status=fixed  weight=1
selfant 5:HN,true
! 5:HN       self-cal  status=fixed  weight=1
selfant 5:KP,true
! 5:KP       self-cal  status=fixed  weight=1
selfant 5:LA,true
! 5:LA       self-cal  status=fixed  weight=1
selfant 5:MK,true
! 5:MK       self-cal  status=fixed  weight=1
selfant 5:NL,true
! 5:NL       self-cal  status=fixed  weight=1
selfant 5:OV,true
! 5:OV       self-cal  status=fixed  weight=1
selfant 5:PT,true
! 5: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.108539Jy  sigma=0.233035
! Fit after  self-cal, rms=0.108539Jy  sigma=0.233035
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:NL       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
! 2:KP       self-cal  status=correctable  weight=1
! 2:LA       self-cal  status=correctable  weight=1
! 2:PT       self-cal  status=correctable  weight=1
! 3:NL       self-cal  status=correctable  weight=1
! 3:SC       self-cal  status=correctable  weight=1
! 4:MK       self-cal  status=correctable  weight=1
! 4:OV       self-cal  status=correctable  weight=1
! 5:BR       self-cal  status=correctable  weight=1
! 5:FD       self-cal  status=correctable  weight=1
! 5:HN       self-cal  status=correctable  weight=1
! 5:KP       self-cal  status=correctable  weight=1
! 5:LA       self-cal  status=correctable  weight=1
! 5:MK       self-cal  status=correctable  weight=1
! 5:NL       self-cal  status=correctable  weight=1
! 5:OV       self-cal  status=correctable  weight=1
! 5: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 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 2.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 3.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 4.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Fit before self-cal, rms=0.108539Jy  sigma=0.233035
! Fit after  self-cal, rms=0.060288Jy  sigma=0.104961
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=6.882 mas, bmaj=13.05 mas, bpa=-1.647 degrees
! Estimated noise=30.799 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 = -3.47806e-05 Jy
! Component: 100  -  total flux cleaned = 0.000908263 Jy
! Component: 150  -  total flux cleaned = 0.00104118 Jy
! Component: 200  -  total flux cleaned = 0.00158879 Jy
! Total flux subtracted in 200 components = 0.00158879 Jy
! Clean residual min=-0.002310 max=0.002316 Jy/beam
! Clean residual mean=-0.000011 rms=0.001054 Jy/beam
! Combined flux in latest and established models = 0.343073 Jy
! Performing phase self-cal
! Adding 70 model components to the UV plane model.
! The established model now contains 544 components and 0.343073 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060359Jy  sigma=0.105032
! Fit after  self-cal, rms=0.060342Jy  sigma=0.104972
! 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=7.923 mas, bmaj=13.94 mas, bpa=-3.325 degrees
! Estimated noise=17.3497 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.000815443 Jy
! Component: 100  -  total flux cleaned = 0.00159756 Jy
! Component: 150  -  total flux cleaned = 0.00230348 Jy
! Component: 200  -  total flux cleaned = 0.00303002 Jy
! Component: 250  -  total flux cleaned = 0.00371599 Jy
! Component: 300  -  total flux cleaned = 0.00429247 Jy
! Component: 350  -  total flux cleaned = 0.00495753 Jy
! Component: 400  -  total flux cleaned = 0.00524461 Jy
! Total flux subtracted in 400 components = 0.00524461 Jy
! Clean residual min=-0.001049 max=0.001040 Jy/beam
! Clean residual mean=0.000020 rms=0.000464 Jy/beam
! Combined flux in latest and established models = 0.348318 Jy
! Performing phase self-cal
! Adding 117 model components to the UV plane model.
! The established model now contains 660 components and 0.348318 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060204Jy  sigma=0.104980
! Fit after  self-cal, rms=0.060202Jy  sigma=0.104965
! 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.000376456 Jy
! Component: 100  -  total flux cleaned = -0.000431157 Jy
! Component: 150  -  total flux cleaned = -0.000404828 Jy
! Component: 200  -  total flux cleaned = -0.000453409 Jy
! Component: 250  -  total flux cleaned = -0.000316834 Jy
! Component: 300  -  total flux cleaned = -0.000495342 Jy
! Component: 350  -  total flux cleaned = -0.000495821 Jy
! Component: 400  -  total flux cleaned = -0.000489499 Jy
! Total flux subtracted in 400 components = -0.000489499 Jy
! Clean residual min=-0.001294 max=0.001291 Jy/beam
! Clean residual mean=-0.000012 rms=0.000507 Jy/beam
! Combined flux in latest and established models = 0.347828 Jy
! Performing phase self-cal
! Adding 91 model components to the UV plane model.
! The established model now contains 747 components and 0.347828 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060182Jy  sigma=0.104993
! Fit after  self-cal, rms=0.060176Jy  sigma=0.104968
! 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.
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 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 2.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 3.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 4.
!  A total of 185 telescope corrections were flagged in sub-array 1.
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  A total of 24 telescope corrections were flagged in sub-array 3.
!  A total of 24 telescope corrections were flagged in sub-array 4.
!  A total of 129 telescope corrections were flagged in sub-array 5.
! 
! Fit before self-cal, rms=0.109073Jy  sigma=0.237170
! Fit after  self-cal, rms=0.060317Jy  sigma=0.104916
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.903 mas, bmaj=13.96 mas, bpa=-3.17 degrees
! Estimated noise=17.3432 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000245813 Jy
! Component: 100  -  total flux cleaned = -0.000232191 Jy
! Component: 150  -  total flux cleaned = -0.000247549 Jy
! Component: 200  -  total flux cleaned = -0.000274069 Jy
! Component: 250  -  total flux cleaned = -0.000220482 Jy
! Component: 300  -  total flux cleaned = -7.69533e-05 Jy
! Component: 350  -  total flux cleaned = 2.61818e-05 Jy
! Component: 400  -  total flux cleaned = 0.000155366 Jy
! Total flux subtracted in 400 components = 0.000155366 Jy
! Clean residual min=-0.000401 max=0.000401 Jy/beam
! Clean residual mean=0.000001 rms=0.000159 Jy/beam
! Combined flux in latest and established models = 0.347984 Jy
selfcal
! Performing phase self-cal
! Adding 100 model components to the UV plane model.
! The established model now contains 845 components and 0.347984 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060286Jy  sigma=0.104914
! Fit after  self-cal, rms=0.060282Jy  sigma=0.104908
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 (-50.2, -46.6).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 5.72732e-05 Jy
! Component: 100  -  total flux cleaned = 0.000199028 Jy
! Component: 150  -  total flux cleaned = 0.000351355 Jy
! Component: 200  -  total flux cleaned = 0.000469408 Jy
! Component: 250  -  total flux cleaned = 0.000515356 Jy
! Component: 300  -  total flux cleaned = 0.000535998 Jy
! Component: 350  -  total flux cleaned = 0.000550329 Jy
! Component: 400  -  total flux cleaned = 0.000528767 Jy
! Total flux subtracted in 400 components = 0.000528767 Jy
! Clean residual min=-0.000235 max=0.000270 Jy/beam
! Clean residual mean=0.000001 rms=0.000078 Jy/beam
! Combined flux in latest and established models = 0.348512 Jy
! Performing phase self-cal
! Adding 62 model components to the UV plane model.
! The established model now contains 904 components and 0.348513 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060269Jy  sigma=0.104906
! Fit after  self-cal, rms=0.060268Jy  sigma=0.104906
! Inverting map 
! Added new window around map position (23.8, -1).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -1.401e-06 Jy
! Component: 100  -  total flux cleaned = -0.000140302 Jy
! Component: 150  -  total flux cleaned = -0.000182956 Jy
! Component: 200  -  total flux cleaned = -0.000136762 Jy
! Component: 250  -  total flux cleaned = -0.000136526 Jy
! Component: 300  -  total flux cleaned = -0.000118035 Jy
! Component: 350  -  total flux cleaned = -0.000111962 Jy
! Component: 400  -  total flux cleaned = -0.000123483 Jy
! Total flux subtracted in 400 components = -0.000123483 Jy
! Clean residual min=-0.000160 max=0.000182 Jy/beam
! Clean residual mean=0.000000 rms=0.000052 Jy/beam
! Combined flux in latest and established models = 0.348389 Jy
! Performing phase self-cal
! Adding 77 model components to the UV plane model.
! The established model now contains 978 components and 0.348389 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060262Jy  sigma=0.104903
! Fit after  self-cal, rms=0.060259Jy  sigma=0.104903
! 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.060259Jy  sigma=0.104903
! Fit after  self-cal, rms=0.060597Jy  sigma=0.104890
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.926 mas, bmaj=13.98 mas, bpa=-3.403 degrees
! Estimated noise=17.4021 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 8.55502e-05 Jy
! Component: 100  -  total flux cleaned = -0.000136941 Jy
! Component: 150  -  total flux cleaned = -0.000324204 Jy
! Component: 200  -  total flux cleaned = -0.000416545 Jy
! Component: 250  -  total flux cleaned = -0.000470405 Jy
! Component: 300  -  total flux cleaned = -0.000522618 Jy
! Component: 350  -  total flux cleaned = -0.00057322 Jy
! Component: 400  -  total flux cleaned = -0.000598179 Jy
! Total flux subtracted in 400 components = -0.000598179 Jy
! Clean residual min=-0.000237 max=0.000221 Jy/beam
! Clean residual mean=0.000002 rms=0.000056 Jy/beam
! Combined flux in latest and established models = 0.347791 Jy
selfcal
! Performing phase self-cal
! Adding 74 model components to the UV plane model.
! The established model now contains 1050 components and 0.347791 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060587Jy  sigma=0.104886
! Fit after  self-cal, rms=0.060586Jy  sigma=0.104886
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 = -2.68271e-05 Jy
! Component: 100  -  total flux cleaned = -0.000136566 Jy
! Component: 150  -  total flux cleaned = -0.000197496 Jy
! Component: 200  -  total flux cleaned = -0.000224171 Jy
! Component: 250  -  total flux cleaned = -0.000245114 Jy
! Component: 300  -  total flux cleaned = -0.000271232 Jy
! Component: 350  -  total flux cleaned = -0.00029742 Jy
! Component: 400  -  total flux cleaned = -0.000318687 Jy
! Total flux subtracted in 400 components = -0.000318687 Jy
! Clean residual min=-0.000198 max=0.000189 Jy/beam
! Clean residual mean=0.000001 rms=0.000049 Jy/beam
! Combined flux in latest and established models = 0.347472 Jy
! Performing phase self-cal
! Adding 55 model components to the UV plane model.
! The established model now contains 1096 components and 0.347472 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060581Jy  sigma=0.104884
! Fit after  self-cal, rms=0.060579Jy  sigma=0.104884
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime3
! Performing amp+phase self-cal over 1 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.060579Jy  sigma=0.104884
! Fit after  self-cal, rms=0.063014Jy  sigma=0.092297
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.131 mas, bmaj=14.15 mas, bpa=-4.509 degrees
! Estimated noise=17.179 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 8.19546e-05 Jy
! Component: 100  -  total flux cleaned = -0.000218318 Jy
! Component: 150  -  total flux cleaned = -0.000351002 Jy
! Component: 200  -  total flux cleaned = -0.000456863 Jy
! Component: 250  -  total flux cleaned = -0.000553387 Jy
! Component: 300  -  total flux cleaned = -0.000598356 Jy
! Component: 350  -  total flux cleaned = -0.000618702 Jy
! Component: 400  -  total flux cleaned = -0.000628577 Jy
! Total flux subtracted in 400 components = -0.000628577 Jy
! Clean residual min=-0.000316 max=0.000285 Jy/beam
! Clean residual mean=0.000001 rms=0.000090 Jy/beam
! Combined flux in latest and established models = 0.346844 Jy
selfcal
! Performing phase self-cal
! Adding 93 model components to the UV plane model.
! The established model now contains 1177 components and 0.346844 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062986Jy  sigma=0.092259
! Fit after  self-cal, rms=0.063024Jy  sigma=0.092257
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 (-44.2, 48.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 3.15733e-05 Jy
! Component: 100  -  total flux cleaned = 1.90069e-05 Jy
! Component: 150  -  total flux cleaned = -8.35689e-06 Jy
! Component: 200  -  total flux cleaned = -7.74019e-05 Jy
! Component: 250  -  total flux cleaned = -0.000179582 Jy
! Component: 300  -  total flux cleaned = -0.00028935 Jy
! Component: 350  -  total flux cleaned = -0.000392219 Jy
! Component: 400  -  total flux cleaned = -0.000451028 Jy
! Total flux subtracted in 400 components = -0.000451028 Jy
! Clean residual min=-0.000606 max=0.000459 Jy/beam
! Clean residual mean=0.000002 rms=0.000153 Jy/beam
! Combined flux in latest and established models = 0.346393 Jy
! Performing phase self-cal
! Adding 88 model components to the UV plane model.
! The established model now contains 1263 components and 0.346393 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.062996Jy  sigma=0.092258
! Fit after  self-cal, rms=0.062985Jy  sigma=0.092256
! 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=7.996 mas, bmaj=14.06 mas, bpa=-4.026 degrees
! Estimated noise=17.6278 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 7.996 x 14.06 at -4.026 degrees (North through East)
! Clean map  min=-0.059016  max=0.20231 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
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=6.969 mas, bmaj=13.2 mas, bpa=-1.922 degrees
! Estimated noise=31.5756 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0252339 Jy
! Component: 100  -  total flux cleaned = 0.0440032 Jy
! Component: 150  -  total flux cleaned = 0.0606924 Jy
! Component: 200  -  total flux cleaned = 0.0781783 Jy
! Component: 250  -  total flux cleaned = 0.102497 Jy
! Component: 300  -  total flux cleaned = 0.125142 Jy
! Component: 350  -  total flux cleaned = 0.154717 Jy
! Component: 400  -  total flux cleaned = 0.187846 Jy
! Total flux subtracted in 400 components = 0.187846 Jy
! Clean residual min=-0.076489 max=0.084162 Jy/beam
! Clean residual mean=0.000387 rms=0.022400 Jy/beam
! Combined flux in latest and established models = 0.187846 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 170 model components to the UV plane model.
! The established model now contains 170 components and 0.187846 Jy
! Inverting map and beam 
! Estimated beam: bmin=7.996 mas, bmaj=14.06 mas, bpa=-4.026 degrees
! Estimated noise=17.6278 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.0592902 Jy
! Component: 100  -  total flux cleaned = -0.0810361 Jy
! Component: 150  -  total flux cleaned = -0.10511 Jy
! Component: 200  -  total flux cleaned = -0.133987 Jy
! Component: 250  -  total flux cleaned = -0.166453 Jy
! Component: 300  -  total flux cleaned = -0.203757 Jy
! Component: 350  -  total flux cleaned = -0.251003 Jy
! Component: 400  -  total flux cleaned = -0.297666 Jy
! Total flux subtracted in 400 components = -0.297666 Jy
! Clean residual min=-0.121118 max=0.115418 Jy/beam
! Clean residual mean=-0.000877 rms=0.038201 Jy/beam
! Combined flux in latest and established models = -0.10982 Jy
! Adding 79 model components to the UV plane model.
! The established model now contains 221 components and -0.10982 Jy
! Inverting map 
! Added new window around map position (15.6, 33.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0682755 Jy
! Component: 100  -  total flux cleaned = 0.0907688 Jy
! Component: 150  -  total flux cleaned = 0.124003 Jy
! Component: 200  -  total flux cleaned = 0.167131 Jy
! Component: 250  -  total flux cleaned = 0.214973 Jy
! Component: 300  -  total flux cleaned = 0.26913 Jy
! Component: 350  -  total flux cleaned = 0.333579 Jy
! Component: 400  -  total flux cleaned = 0.39711 Jy
! Total flux subtracted in 400 components = 0.39711 Jy
! Clean residual min=-0.182583 max=0.208958 Jy/beam
! Clean residual mean=0.000665 rms=0.056949 Jy/beam
! Combined flux in latest and established models = 0.28729 Jy
! Adding 75 model components to the UV plane model.
! The established model now contains 257 components and 0.28729 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.131 mas, bmaj=14.15 mas, bpa=-4.509 degrees
! Estimated noise=17.179 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0457269 Jy
! Component: 100  -  total flux cleaned = -0.0889552 Jy
! Component: 150  -  total flux cleaned = -0.151906 Jy
! Component: 200  -  total flux cleaned = -0.220733 Jy
! Component: 250  -  total flux cleaned = -0.30435 Jy
! Component: 300  -  total flux cleaned = -0.409217 Jy
! Component: 350  -  total flux cleaned = -0.52351 Jy
! Component: 400  -  total flux cleaned = -0.659101 Jy
! Total flux subtracted in 400 components = -0.659101 Jy
! Clean residual min=-0.391873 max=0.346889 Jy/beam
! Clean residual mean=-0.001228 rms=0.106183 Jy/beam
! Combined flux in latest and established models = -0.371811 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.162423 Jy
! Component: 100  -  total flux cleaned = -0.322631 Jy
! Component: 150  -  total flux cleaned = -0.531683 Jy
! Component: 200  -  total flux cleaned = -0.781664 Jy
! Component: 250  -  total flux cleaned = -1.05108 Jy
! Component: 300  -  total flux cleaned = -1.36856 Jy
! Component: 350  -  total flux cleaned = -1.78558 Jy
! Component: 400  -  total flux cleaned = -2.23471 Jy
! Total flux subtracted in 400 components = -2.23471 Jy
! Clean residual min=-1.688902 max=1.431152 Jy/beam
! Clean residual mean=-0.003376 rms=0.425161 Jy/beam
! Combined flux in latest and established models = -2.60652 Jy
! Adding 101 model components to the UV plane model.
! The established model now contains 307 components and -2.60652 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=7.996 mas, bmaj=14.06 mas, bpa=-4.026 degrees
! Estimated noise=17.6278 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.318328 Jy
! Component: 100  -  total flux cleaned = 0.809906 Jy
! Component: 150  -  total flux cleaned = 1.33608 Jy
! Component: 200  -  total flux cleaned = 1.91078 Jy
! Component: 250  -  total flux cleaned = 2.6122 Jy
! Component: 300  -  total flux cleaned = 3.63026 Jy
! Component: 350  -  total flux cleaned = 4.6708 Jy
! Component: 400  -  total flux cleaned = 5.61041 Jy
! Total flux subtracted in 400 components = 5.61041 Jy
! Clean residual min=-2.836576 max=3.188854 Jy/beam
! Clean residual mean=0.009073 rms=0.863025 Jy/beam
! Combined flux in latest and established models = 3.00389 Jy
selfcal
! Performing phase self-cal
! Adding 70 model components to the UV plane model.
! The established model now contains 322 components and 3.00388 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=3.387178Jy  sigma=5.396346
! Fit after  self-cal, rms=3.349257Jy  sigma=5.314090
wmodel /scr/pima/rdv62_uvs/J2206-1835_X_map.mod
! Writing 322 model components to file: /scr/pima/rdv62_uvs/J2206-1835_X_map.mod
wobs   /scr/pima/rdv62_uvs/J2206-1835_X_uvs.fits
! Applying 604 buffered edits.
! Writing UV FITS file: /scr/pima/rdv62_uvs/J2206-1835_X_uvs.fits
wwins  /scr/pima/rdv62_uvs/J2206-1835_X_map.win
! wwins: Wrote 4 windows to /scr/pima/rdv62_uvs/J2206-1835_X_map.win
x = (field_size-8) * field_cell / 4
addwin -x,x,-x,x
clean (field_size*4),0.01
! Inverting map 
! clean: niter=4096  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.523519 Jy
! Component: 100  -  total flux cleaned = -0.55358 Jy
! Component: 150  -  total flux cleaned = -0.55447 Jy
! Component: 200  -  total flux cleaned = -0.604919 Jy
! Component: 250  -  total flux cleaned = -0.802444 Jy
! Component: 300  -  total flux cleaned = -1.07597 Jy
! Component: 350  -  total flux cleaned = -1.35441 Jy
! Component: 400  -  total flux cleaned = -1.63997 Jy
! Component: 450  -  total flux cleaned = -1.93377 Jy
! Component: 500  -  total flux cleaned = -2.20815 Jy
! Component: 550  -  total flux cleaned = -2.35198 Jy
! Component: 600  -  total flux cleaned = -2.53655 Jy
! Component: 650  -  total flux cleaned = -2.73527 Jy
! Component: 700  -  total flux cleaned = -2.94849 Jy
! Component: 750  -  total flux cleaned = -3.216 Jy
! Component: 800  -  total flux cleaned = -3.46423 Jy
! Component: 850  -  total flux cleaned = -3.77539 Jy
! Component: 900  -  total flux cleaned = -4.01338 Jy
! Component: 950  -  total flux cleaned = -4.32209 Jy
! Component: 1000  -  total flux cleaned = -4.546 Jy
! Component: 1050  -  total flux cleaned = -4.84716 Jy
! Component: 1100  -  total flux cleaned = -5.10894 Jy
! Component: 1150  -  total flux cleaned = -5.4603 Jy
! Component: 1200  -  total flux cleaned = -5.7667 Jy
! Component: 1250  -  total flux cleaned = -6.17669 Jy
! Component: 1300  -  total flux cleaned = -6.70839 Jy
! Component: 1350  -  total flux cleaned = -7.18465 Jy
! Component: 1400  -  total flux cleaned = -7.80096 Jy
! Component: 1450  -  total flux cleaned = -8.47195 Jy
! Component: 1500  -  total flux cleaned = -9.19997 Jy
! Component: 1550  -  total flux cleaned = -9.85502 Jy
! Component: 1600  -  total flux cleaned = -10.7027 Jy
! Component: 1650  -  total flux cleaned = -11.4631 Jy
! Component: 1700  -  total flux cleaned = -12.4485 Jy
! Component: 1750  -  total flux cleaned = -13.5173 Jy
! Component: 1800  -  total flux cleaned = -14.2787 Jy
! Component: 1850  -  total flux cleaned = -15.509 Jy
! Component: 1900  -  total flux cleaned = -17.079 Jy
! Component: 1950  -  total flux cleaned = -18.2863 Jy
! Component: 2000  -  total flux cleaned = -20.1254 Jy
! Component: 2050  -  total flux cleaned = -21.5617 Jy
! Component: 2100  -  total flux cleaned = -23.0956 Jy
! Component: 2150  -  total flux cleaned = -25.4255 Jy
! Component: 2200  -  total flux cleaned = -27.2334 Jy
! Component: 2250  -  total flux cleaned = -29.5779 Jy
! Component: 2300  -  total flux cleaned = -31.6901 Jy
! Component: 2350  -  total flux cleaned = -33.9372 Jy
! Component: 2400  -  total flux cleaned = -37.3558 Jy
! Component: 2450  -  total flux cleaned = -40.0156 Jy
! Component: 2500  -  total flux cleaned = -42.8435 Jy
! Component: 2550  -  total flux cleaned = -47.1648 Jy
! Component: 2600  -  total flux cleaned = -50.508 Jy
! Component: 2650  -  total flux cleaned = -54.8429 Jy
! Component: 2700  -  total flux cleaned = -58.7315 Jy
! Component: 2750  -  total flux cleaned = -63.7443 Jy
! Component: 2800  -  total flux cleaned = -69.1769 Jy
! Component: 2850  -  total flux cleaned = -76.1249 Jy
! Component: 2900  -  total flux cleaned = -80.3945 Jy
! Component: 2950  -  total flux cleaned = -87.2655 Jy
! Component: 3000  -  total flux cleaned = -96.021 Jy
! Component: 3050  -  total flux cleaned = -102.82 Jy
! Component: 3100  -  total flux cleaned = -111.63 Jy
! Component: 3150  -  total flux cleaned = -119.486 Jy
! Component: 3200  -  total flux cleaned = -131.471 Jy
! Component: 3250  -  total flux cleaned = -140.812 Jy
! Component: 3300  -  total flux cleaned = -150.745 Jy
! Component: 3350  -  total flux cleaned = -163.628 Jy
! Component: 3400  -  total flux cleaned = -177.59 Jy
! Component: 3450  -  total flux cleaned = -192.738 Jy
! Component: 3500  -  total flux cleaned = -206.281 Jy
! Component: 3550  -  total flux cleaned = -223.882 Jy
! Component: 3600  -  total flux cleaned = -243.028 Jy
! Component: 3650  -  total flux cleaned = -260.172 Jy
! Component: 3700  -  total flux cleaned = -282.382 Jy
! Component: 3750  -  total flux cleaned = -306.42 Jy
! Component: 3800  -  total flux cleaned = -337.086 Jy
! Component: 3850  -  total flux cleaned = -356.027 Jy
! Component: 3900  -  total flux cleaned = -386.389 Jy
! Component: 3950  -  total flux cleaned = -419.368 Jy
! Component: 4000  -  total flux cleaned = -455.186 Jy
! Component: 4050  -  total flux cleaned = -487.302 Jy
! Total flux subtracted in 4096 components = -528.841 Jy
! Clean residual min=-370.092896 max=385.031891 Jy/beam
! Clean residual mean=-5.416779 rms=165.757141 Jy/beam
! Combined flux in latest and established models = -525.837 Jy
keep
! Adding 67 model components to the UV plane model.
! The established model now contains 341 components and -525.837 Jy
![Exited script file: /home/lpetrov/bin/rdv_muppet]
wmap "/scr/pima/rdv62_uvs/J2206-1835_X\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 7.996 x 14.06 at -4.026 degrees (North through East)
! Clean map  min=-484.86  max=332.87 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J2206-1835_X_map.fits
float peak_flux
peak_flux = peak(flux,max)
invert
! Inverting map 
float image_peak
image_peak = int(10000*peak_flux)/10000
float image_rms
image_rms  = int(100000*imstat(rms))/100000
print "========================================="
! ========================================= 
print "The image peak= ", image_peak, " rms= ", image_rms, " Jy"
! The image peak=  332.874  rms=  201.096  Jy 
quit
! Quitting program
! Log file /scr/pima/rdv62_uvs/J2206-1835_X_uvs.log closed on Fri Nov 22 22:42:50 2013