! Started logfile: /scr/pima/rdv62_uvs/J0339-0133_X_uvs.log on Fri Nov 22 22:32:21 2013
obs /scr/pima/rdv62_uvs/J0339-0133_X_uva.fits
! Reading UV FITS file: /scr/pima/rdv62_uvs/J0339-0133_X_uva.fits
! AN table 1: 23 integrations on 28 of 28 possible baselines.
! Apparent sampling: 0.835404 visibilities/baseline/integration-bin.
! Found source: J0339-0133
! 
! 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 2152 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/J0339-0133_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:BR       self-cal  status=correctable  weight=1
! 1:FD       self-cal  status=correctable  weight=1
! 1:HN       self-cal  status=correctable  weight=1
! 1:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:MK       self-cal  status=correctable  weight=1
! 1:NL       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
system "rm /tmp/ff.ff"
mapsize field_size, field_cell
! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize.
startmod
! Applying default point source starting model.
! Performing phase self-cal
! Adding 1 model components to the UV plane model.
! The established model now contains 1 components and 1 Jy
! 
! Correcting IF 1.
!  A total of 8 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 8 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 8 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 8 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.028995Jy  sigma=2.578197
! Fit after  self-cal, rms=0.805160Jy  sigma=1.809257
! 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=1.344 mas, bmaj=2.675 mas, bpa=-14.83 degrees
! Estimated noise=15.5086 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (0, 0).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.177203 Jy
! Component: 100  -  total flux cleaned = 0.236623 Jy
! Total flux subtracted in 100 components = 0.236623 Jy
! Clean residual min=-0.021841 max=0.032534 Jy/beam
! Clean residual mean=-0.000026 rms=0.007592 Jy/beam
! Combined flux in latest and established models = 0.236624 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 16 components and 0.236624 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.090583Jy  sigma=0.220685
! Fit after  self-cal, rms=0.091049Jy  sigma=0.220472
! 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=1.848 mas, bmaj=3.036 mas, bpa=-7.228 degrees
! Estimated noise=6.94862 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.0366963 Jy
! Component: 100  -  total flux cleaned = 0.0545182 Jy
! Component: 150  -  total flux cleaned = 0.0643008 Jy
! Component: 200  -  total flux cleaned = 0.0684583 Jy
! Total flux subtracted in 200 components = 0.0684583 Jy
! Clean residual min=-0.013770 max=0.014364 Jy/beam
! Clean residual mean=-0.000017 rms=0.003488 Jy/beam
! Combined flux in latest and established models = 0.305082 Jy
! Performing phase self-cal
! Adding 34 model components to the UV plane model.
! The established model now contains 46 components and 0.305082 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.082122Jy  sigma=0.170963
! Fit after  self-cal, rms=0.082413Jy  sigma=0.170874
! 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:
!   BR        1.05     FD        0.97     HN        1.17     KP        0.98     
!   LA        0.93     MK        0.84     NL        1.01     PT        0.97     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.02     FD        1.00     HN        1.14     KP        0.91     
!   LA        0.94     MK        0.98     NL        0.98     PT        1.02     
! 
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.05     FD        1.05     HN        1.26     KP        0.85     
!   LA        1.02     MK        1.05     NL        1.00     PT        1.07     
! 
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.07     FD        1.10     HN        1.28     KP        0.90     
!   LA        0.98     MK        1.14     NL        1.05     PT        1.06     
! 
! 
! Fit before self-cal, rms=0.082413Jy  sigma=0.170874
! Fit after  self-cal, rms=0.082991Jy  sigma=0.158457
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=1.846 mas, bmaj=3.07 mas, bpa=-7.82 degrees
! Estimated noise=6.98664 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 (-7.4, -6.2).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00581594 Jy
! Component: 100  -  total flux cleaned = 0.00369449 Jy
! Component: 150  -  total flux cleaned = 0.00296753 Jy
! Component: 200  -  total flux cleaned = 0.00266123 Jy
! Total flux subtracted in 200 components = 0.00266123 Jy
! Clean residual min=-0.005828 max=0.007728 Jy/beam
! Clean residual mean=0.000000 rms=0.001622 Jy/beam
! Combined flux in latest and established models = 0.307743 Jy
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 65 components and 0.307743 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.081225Jy  sigma=0.155496
! Fit after  self-cal, rms=0.081481Jy  sigma=0.155382
! 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.081481Jy  sigma=0.155382
! Fit after  self-cal, rms=0.080744Jy  sigma=0.155236
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.845 mas, bmaj=3.067 mas, bpa=-7.814 degrees
! Estimated noise=6.95802 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000780912 Jy
! Component: 100  -  total flux cleaned = 0.000252089 Jy
! Component: 150  -  total flux cleaned = 3.26701e-06 Jy
! Component: 200  -  total flux cleaned = -3.57224e-07 Jy
! Total flux subtracted in 200 components = -3.57224e-07 Jy
! Clean residual min=-0.005583 max=0.006081 Jy/beam
! Clean residual mean=0.000000 rms=0.001456 Jy/beam
! Combined flux in latest and established models = 0.307743 Jy
selfcal
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 73 components and 0.307743 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.079899Jy  sigma=0.154695
! Fit after  self-cal, rms=0.080049Jy  sigma=0.154660
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.000358657 Jy
! Component: 100  -  total flux cleaned = 0.000355163 Jy
! Component: 150  -  total flux cleaned = 0.000637044 Jy
! Component: 200  -  total flux cleaned = 0.000905417 Jy
! Total flux subtracted in 200 components = 0.000905417 Jy
! Clean residual min=-0.005438 max=0.005510 Jy/beam
! Clean residual mean=0.000000 rms=0.001382 Jy/beam
! Combined flux in latest and established models = 0.308648 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 77 components and 0.308648 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.079500Jy  sigma=0.154402
! Fit after  self-cal, rms=0.079613Jy  sigma=0.154369
! 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.079613Jy  sigma=0.154369
! Fit after  self-cal, rms=0.079173Jy  sigma=0.154303
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.854 mas, bmaj=3.074 mas, bpa=-7.775 degrees
! Estimated noise=6.92984 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000585973 Jy
! Component: 100  -  total flux cleaned = -2.95093e-05 Jy
! Component: 150  -  total flux cleaned = 0.000393137 Jy
! Component: 200  -  total flux cleaned = 0.000553674 Jy
! Total flux subtracted in 200 components = 0.000553674 Jy
! Clean residual min=-0.005038 max=0.005137 Jy/beam
! Clean residual mean=-0.000000 rms=0.001289 Jy/beam
! Combined flux in latest and established models = 0.309202 Jy
selfcal
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 82 components and 0.309202 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.078608Jy  sigma=0.154121
! Fit after  self-cal, rms=0.078726Jy  sigma=0.154087
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.000236776 Jy
! Component: 100  -  total flux cleaned = 0.000466346 Jy
! Component: 150  -  total flux cleaned = 0.000537313 Jy
! Component: 200  -  total flux cleaned = 0.000604936 Jy
! Total flux subtracted in 200 components = 0.000604936 Jy
! Clean residual min=-0.004927 max=0.004968 Jy/beam
! Clean residual mean=-0.000000 rms=0.001247 Jy/beam
! Combined flux in latest and established models = 0.309807 Jy
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 87 components and 0.309807 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.078355Jy  sigma=0.153961
! Fit after  self-cal, rms=0.078448Jy  sigma=0.153930
! 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.078448Jy  sigma=0.153930
! Fit after  self-cal, rms=0.078108Jy  sigma=0.153878
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.267 mas, bmaj=4.297 mas, bpa=30.46 degrees
! Estimated noise=6.53155 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000301977 Jy
! Component: 100  -  total flux cleaned = 0.000753414 Jy
! Component: 150  -  total flux cleaned = 0.00127545 Jy
! Component: 200  -  total flux cleaned = 0.00164479 Jy
! Total flux subtracted in 200 components = 0.00164479 Jy
! Clean residual min=-0.003674 max=0.004099 Jy/beam
! Clean residual mean=0.000000 rms=0.001110 Jy/beam
! Combined flux in latest and established models = 0.311451 Jy
selfcal
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 94 components and 0.311451 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077813Jy  sigma=0.153789
! Fit after  self-cal, rms=0.077886Jy  sigma=0.153780
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.00022779 Jy
! Component: 100  -  total flux cleaned = 0.000468876 Jy
! Component: 150  -  total flux cleaned = 0.000651743 Jy
! Component: 200  -  total flux cleaned = 0.00074152 Jy
! Total flux subtracted in 200 components = 0.00074152 Jy
! Clean residual min=-0.003698 max=0.004003 Jy/beam
! Clean residual mean=0.000000 rms=0.001091 Jy/beam
! Combined flux in latest and established models = 0.312193 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 98 components and 0.312193 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077708Jy  sigma=0.153718
! Fit after  self-cal, rms=0.077758Jy  sigma=0.153708
! 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:BR,true
! 1:BR       self-cal  status=fixed  weight=1
selfant 1:FD,true
! 1:FD       self-cal  status=fixed  weight=1
selfant 1:HN,true
! 1:HN       self-cal  status=fixed  weight=1
selfant 1:KP,true
! 1:KP       self-cal  status=fixed  weight=1
selfant 1:LA,true
! 1:LA       self-cal  status=fixed  weight=1
selfant 1:MK,true
! 1:MK       self-cal  status=fixed  weight=1
selfant 1:NL,true
! 1:NL       self-cal  status=fixed  weight=1
selfant 1:PT,true
! 1:PT       self-cal  status=fixed  weight=1
![Exited script file: /tmp/ff.ff]
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.078165Jy  sigma=0.154812
! Fit after  self-cal, rms=0.078165Jy  sigma=0.154812
selfant "",false
! 1:BR       self-cal  status=correctable  weight=1
! 1:FD       self-cal  status=correctable  weight=1
! 1:HN       self-cal  status=correctable  weight=1
! 1:KP       self-cal  status=correctable  weight=1
! 1:LA       self-cal  status=correctable  weight=1
! 1:MK       self-cal  status=correctable  weight=1
! 1:NL       self-cal  status=correctable  weight=1
! 1:PT       self-cal  status=correctable  weight=1
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.078165Jy  sigma=0.154812
! Fit after  self-cal, rms=0.077763Jy  sigma=0.153665
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=1.363 mas, bmaj=2.717 mas, bpa=-14.57 degrees
! Estimated noise=16.004 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.00165532 Jy
! Component: 100  -  total flux cleaned = -0.00148751 Jy
! Component: 150  -  total flux cleaned = -0.000757785 Jy
! Component: 200  -  total flux cleaned = -0.00022757 Jy
! Total flux subtracted in 200 components = -0.00022757 Jy
! Clean residual min=-0.006386 max=0.006369 Jy/beam
! Clean residual mean=0.000002 rms=0.001735 Jy/beam
! Combined flux in latest and established models = 0.311965 Jy
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 110 components and 0.311965 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077363Jy  sigma=0.153877
! Fit after  self-cal, rms=0.077490Jy  sigma=0.153747
! 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=1.864 mas, bmaj=3.083 mas, bpa=-7.771 degrees
! Estimated noise=6.89685 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 (-1.4, 29.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000608467 Jy
! Component: 100  -  total flux cleaned = 0.000577386 Jy
! Component: 150  -  total flux cleaned = 0.00109271 Jy
! Component: 200  -  total flux cleaned = 0.00109227 Jy
! Component: 250  -  total flux cleaned = 0.000927068 Jy
! Component: 300  -  total flux cleaned = 0.000701162 Jy
! Component: 350  -  total flux cleaned = 0.000772857 Jy
! Component: 400  -  total flux cleaned = 0.000640325 Jy
! Total flux subtracted in 400 components = 0.000640325 Jy
! Clean residual min=-0.003379 max=0.003533 Jy/beam
! Clean residual mean=-0.000001 rms=0.000823 Jy/beam
! Combined flux in latest and established models = 0.312606 Jy
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 153 components and 0.312605 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.076607Jy  sigma=0.152723
! Fit after  self-cal, rms=0.076620Jy  sigma=0.152689
! 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.000365074 Jy
! Component: 100  -  total flux cleaned = -0.000618025 Jy
! Component: 150  -  total flux cleaned = -0.000916992 Jy
! Component: 200  -  total flux cleaned = -0.00114031 Jy
! Component: 250  -  total flux cleaned = -0.00129971 Jy
! Component: 300  -  total flux cleaned = -0.00145266 Jy
! Component: 350  -  total flux cleaned = -0.00169311 Jy
! Component: 400  -  total flux cleaned = -0.00183003 Jy
! Total flux subtracted in 400 components = -0.00183003 Jy
! Clean residual min=-0.003176 max=0.003073 Jy/beam
! Clean residual mean=-0.000001 rms=0.000730 Jy/beam
! Combined flux in latest and established models = 0.310775 Jy
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 162 components and 0.310775 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.076301Jy  sigma=0.152480
! Fit after  self-cal, rms=0.076329Jy  sigma=0.152468
! Inverting map 
print "=== Finished natural weighting clean for al sites ==="
! === Finished natural weighting clean for al sites === 
unflag *
! Unflagging all baselines in the currently selected channels.
uncalib false, false, true
! uncal: All telescope correction flags have been un-done.
selfcal true, true, soltime1
! Performing amp+phase self-cal over 800000 minute time intervals
! 
! Correcting IF 1.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 11 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.076791Jy  sigma=0.153631
! Fit after  self-cal, rms=0.075981Jy  sigma=0.152432
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.849 mas, bmaj=3.08 mas, bpa=-8.132 degrees
! Estimated noise=6.87709 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000297567 Jy
! Component: 100  -  total flux cleaned = -0.000255288 Jy
! Component: 150  -  total flux cleaned = -0.000341681 Jy
! Component: 200  -  total flux cleaned = -0.000381966 Jy
! Component: 250  -  total flux cleaned = -0.00046026 Jy
! Component: 300  -  total flux cleaned = -0.000646288 Jy
! Component: 350  -  total flux cleaned = -0.000788297 Jy
! Component: 400  -  total flux cleaned = -0.000856894 Jy
! Total flux subtracted in 400 components = -0.000856894 Jy
! Clean residual min=-0.002888 max=0.002598 Jy/beam
! Clean residual mean=-0.000001 rms=0.000665 Jy/beam
! Combined flux in latest and established models = 0.309919 Jy
selfcal
! Performing phase self-cal
! Adding 39 model components to the UV plane model.
! The established model now contains 185 components and 0.309919 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.075779Jy  sigma=0.152339
! Fit after  self-cal, rms=0.075799Jy  sigma=0.152327
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.000122241 Jy
! Component: 100  -  total flux cleaned = -0.000227272 Jy
! Component: 150  -  total flux cleaned = -0.000359597 Jy
! Component: 200  -  total flux cleaned = -0.000455977 Jy
! Component: 250  -  total flux cleaned = -0.000518908 Jy
! Component: 300  -  total flux cleaned = -0.000639621 Jy
! Component: 350  -  total flux cleaned = -0.000727452 Jy
! Component: 400  -  total flux cleaned = -0.000756578 Jy
! Total flux subtracted in 400 components = -0.000756578 Jy
! Clean residual min=-0.002794 max=0.002390 Jy/beam
! Clean residual mean=-0.000001 rms=0.000630 Jy/beam
! Combined flux in latest and established models = 0.309162 Jy
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 198 components and 0.309162 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.075679Jy  sigma=0.152270
! Fit after  self-cal, rms=0.075695Jy  sigma=0.152264
! 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.075695Jy  sigma=0.152264
! Fit after  self-cal, rms=0.075584Jy  sigma=0.152255
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.842 mas, bmaj=3.079 mas, bpa=-8.29 degrees
! Estimated noise=6.86457 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000271783 Jy
! Component: 100  -  total flux cleaned = -0.000272078 Jy
! Component: 150  -  total flux cleaned = -0.000328288 Jy
! Component: 200  -  total flux cleaned = -0.000409003 Jy
! Component: 250  -  total flux cleaned = -0.000460295 Jy
! Component: 300  -  total flux cleaned = -0.000535398 Jy
! Component: 350  -  total flux cleaned = -0.000559427 Jy
! Component: 400  -  total flux cleaned = -0.000607198 Jy
! Total flux subtracted in 400 components = -0.000607198 Jy
! Clean residual min=-0.002694 max=0.002328 Jy/beam
! Clean residual mean=-0.000001 rms=0.000600 Jy/beam
! Combined flux in latest and established models = 0.308555 Jy
selfcal
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 217 components and 0.308555 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.075487Jy  sigma=0.152217
! Fit after  self-cal, rms=0.075501Jy  sigma=0.152212
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 = -5.5485e-05 Jy
! Component: 100  -  total flux cleaned = -8.01196e-05 Jy
! Component: 150  -  total flux cleaned = -0.000175179 Jy
! Component: 200  -  total flux cleaned = -0.000198051 Jy
! Component: 250  -  total flux cleaned = -0.000176132 Jy
! Component: 300  -  total flux cleaned = -0.000220049 Jy
! Component: 350  -  total flux cleaned = -0.00024222 Jy
! Component: 400  -  total flux cleaned = -0.000263284 Jy
! Total flux subtracted in 400 components = -0.000263284 Jy
! Clean residual min=-0.002643 max=0.002373 Jy/beam
! Clean residual mean=-0.000001 rms=0.000583 Jy/beam
! Combined flux in latest and established models = 0.308292 Jy
! Performing phase self-cal
! Adding 49 model components to the UV plane model.
! The established model now contains 224 components and 0.308292 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.075439Jy  sigma=0.152184
! Fit after  self-cal, rms=0.075449Jy  sigma=0.152181
! 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.075449Jy  sigma=0.152181
! Fit after  self-cal, rms=0.075049Jy  sigma=0.147472
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.273 mas, bmaj=4.285 mas, bpa=30.19 degrees
! Estimated noise=6.51229 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000178188 Jy
! Component: 100  -  total flux cleaned = 0.000264125 Jy
! Component: 150  -  total flux cleaned = 0.000302988 Jy
! Component: 200  -  total flux cleaned = 0.000339473 Jy
! Component: 250  -  total flux cleaned = 0.00033877 Jy
! Component: 300  -  total flux cleaned = 0.000338532 Jy
! Component: 350  -  total flux cleaned = 0.000338471 Jy
! Component: 400  -  total flux cleaned = 0.000338686 Jy
! Total flux subtracted in 400 components = 0.000338686 Jy
! Clean residual min=-0.002161 max=0.002571 Jy/beam
! Clean residual mean=-0.000001 rms=0.000554 Jy/beam
! Combined flux in latest and established models = 0.30863 Jy
selfcal
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 247 components and 0.30863 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.075010Jy  sigma=0.147447
! Fit after  self-cal, rms=0.075013Jy  sigma=0.147434
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 (-46.2, -45.6).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00112863 Jy
! Component: 100  -  total flux cleaned = 0.00120044 Jy
! Component: 150  -  total flux cleaned = 0.00115411 Jy
! Component: 200  -  total flux cleaned = 0.00100116 Jy
! Component: 250  -  total flux cleaned = 0.000870494 Jy
! Component: 300  -  total flux cleaned = 0.000692926 Jy
! Component: 350  -  total flux cleaned = 0.00058312 Jy
! Component: 400  -  total flux cleaned = 0.000583543 Jy
! Total flux subtracted in 400 components = 0.000583543 Jy
! Clean residual min=-0.002301 max=0.002119 Jy/beam
! Clean residual mean=-0.000000 rms=0.000487 Jy/beam
! Combined flux in latest and established models = 0.309214 Jy
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 292 components and 0.309214 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.074892Jy  sigma=0.147150
! Fit after  self-cal, rms=0.074895Jy  sigma=0.147140
! 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=1.848 mas, bmaj=3.081 mas, bpa=-8.185 degrees
! Estimated noise=6.87669 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 1.848 x 3.081 at -8.185 degrees (North through East)
! Clean map  min=-0.0057562  max=0.18576 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=1.365 mas, bmaj=2.729 mas, bpa=-14.3 degrees
! Estimated noise=16.0105 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.176113 Jy
! Component: 100  -  total flux cleaned = 0.234013 Jy
! Component: 150  -  total flux cleaned = 0.260793 Jy
! Component: 200  -  total flux cleaned = 0.271666 Jy
! Component: 250  -  total flux cleaned = 0.278603 Jy
! Component: 300  -  total flux cleaned = 0.28396 Jy
! Component: 350  -  total flux cleaned = 0.287124 Jy
! Component: 400  -  total flux cleaned = 0.288982 Jy
! Total flux subtracted in 400 components = 0.288982 Jy
! Clean residual min=-0.007654 max=0.012729 Jy/beam
! Clean residual mean=-0.000003 rms=0.002332 Jy/beam
! Combined flux in latest and established models = 0.288982 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 61 model components to the UV plane model.
! The established model now contains 61 components and 0.288982 Jy
! Inverting map and beam 
! Estimated beam: bmin=1.848 mas, bmaj=3.081 mas, bpa=-8.185 degrees
! Estimated noise=6.87669 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)
! Added new window around map position (-2.2, -0.8).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0118169 Jy
! Component: 100  -  total flux cleaned = 0.0165339 Jy
! Component: 150  -  total flux cleaned = 0.019055 Jy
! Component: 200  -  total flux cleaned = 0.0200049 Jy
! Component: 250  -  total flux cleaned = 0.020482 Jy
! Component: 300  -  total flux cleaned = 0.0210063 Jy
! Component: 350  -  total flux cleaned = 0.021291 Jy
! Component: 400  -  total flux cleaned = 0.0216422 Jy
! Total flux subtracted in 400 components = 0.0216422 Jy
! Clean residual min=-0.003246 max=0.002942 Jy/beam
! Clean residual mean=0.000000 rms=0.000752 Jy/beam
! Combined flux in latest and established models = 0.310624 Jy
! Adding 46 model components to the UV plane model.
! The established model now contains 106 components and 0.310624 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=3.273 mas, bmaj=4.285 mas, bpa=30.19 degrees
! Estimated noise=6.51229 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000232361 Jy
! Component: 100  -  total flux cleaned = -0.000432708 Jy
! Component: 150  -  total flux cleaned = -0.000679449 Jy
! Component: 200  -  total flux cleaned = -0.000796216 Jy
! Component: 250  -  total flux cleaned = -0.000961637 Jy
! Component: 300  -  total flux cleaned = -0.00101255 Jy
! Component: 350  -  total flux cleaned = -0.00126112 Jy
! Component: 400  -  total flux cleaned = -0.00144952 Jy
! Total flux subtracted in 400 components = -0.00144952 Jy
! Clean residual min=-0.002117 max=0.002210 Jy/beam
! Clean residual mean=0.000001 rms=0.000568 Jy/beam
! Combined flux in latest and established models = 0.309175 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 (48.6, -17).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00142616 Jy
! Component: 100  -  total flux cleaned = 0.00207939 Jy
! Component: 150  -  total flux cleaned = 0.00245664 Jy
! Component: 200  -  total flux cleaned = 0.00270286 Jy
! Component: 250  -  total flux cleaned = 0.00279612 Jy
! Component: 300  -  total flux cleaned = 0.0027963 Jy
! Component: 350  -  total flux cleaned = 0.0027546 Jy
! Component: 400  -  total flux cleaned = 0.0026736 Jy
! Total flux subtracted in 400 components = 0.0026736 Jy
! Clean residual min=-0.001665 max=0.001798 Jy/beam
! Clean residual mean=-0.000000 rms=0.000436 Jy/beam
! Combined flux in latest and established models = 0.311848 Jy
! Adding 106 model components to the UV plane model.
! The established model now contains 200 components and 0.311848 Jy
! Inverting map 
! Added new window around map position (-49.8, 47.2).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00100046 Jy
! Component: 100  -  total flux cleaned = 0.00150274 Jy
! Component: 150  -  total flux cleaned = 0.00161957 Jy
! Component: 200  -  total flux cleaned = 0.00191987 Jy
! Component: 250  -  total flux cleaned = 0.0022462 Jy
! Component: 300  -  total flux cleaned = 0.0025986 Jy
! Component: 350  -  total flux cleaned = 0.0028034 Jy
! Component: 400  -  total flux cleaned = 0.00300489 Jy
! Total flux subtracted in 400 components = 0.00300489 Jy
! Clean residual min=-0.001357 max=0.001484 Jy/beam
! Clean residual mean=0.000000 rms=0.000351 Jy/beam
! Combined flux in latest and established models = 0.314853 Jy
! Adding 49 model components to the UV plane model.
! The established model now contains 232 components and 0.314853 Jy
! Inverting map 
! Added new window around map position (50.8, 45).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000151915 Jy
! Component: 100  -  total flux cleaned = 0.000153278 Jy
! Component: 150  -  total flux cleaned = 0.000190386 Jy
! Component: 200  -  total flux cleaned = 0.000304489 Jy
! Component: 250  -  total flux cleaned = 0.000452071 Jy
! Component: 300  -  total flux cleaned = 0.000560617 Jy
! Component: 350  -  total flux cleaned = 0.000525419 Jy
! Component: 400  -  total flux cleaned = 0.000419321 Jy
! Total flux subtracted in 400 components = 0.000419321 Jy
! Clean residual min=-0.001427 max=0.001529 Jy/beam
! Clean residual mean=-0.000002 rms=0.000345 Jy/beam
! Combined flux in latest and established models = 0.315273 Jy
! Adding 55 model components to the UV plane model.
! The established model now contains 281 components and 0.315273 Jy
! Inverting map 
! Added new window around map position (-47, -23).
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000211181 Jy
! Component: 100  -  total flux cleaned = 0.000124317 Jy
! Component: 150  -  total flux cleaned = 8.44152e-05 Jy
! Component: 200  -  total flux cleaned = 0.000120155 Jy
! Component: 250  -  total flux cleaned = 0.000223489 Jy
! Component: 300  -  total flux cleaned = 0.00039142 Jy
! Component: 350  -  total flux cleaned = 0.000621625 Jy
! Component: 400  -  total flux cleaned = 0.000718885 Jy
! Total flux subtracted in 400 components = 0.000718885 Jy
! Clean residual min=-0.001150 max=0.001038 Jy/beam
! Clean residual mean=-0.000001 rms=0.000276 Jy/beam
! Combined flux in latest and established models = 0.315992 Jy
! Adding 56 model components to the UV plane model.
! The established model now contains 328 components and 0.315992 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=1.848 mas, bmaj=3.081 mas, bpa=-8.185 degrees
! Estimated noise=6.87669 mJy/beam.
! clean: niter=400  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000173514 Jy
! Component: 100  -  total flux cleaned = 0.000115687 Jy
! Component: 150  -  total flux cleaned = 6.12033e-06 Jy
! Component: 200  -  total flux cleaned = 5.40587e-06 Jy
! Component: 250  -  total flux cleaned = 0.000104305 Jy
! Component: 300  -  total flux cleaned = 0.000246499 Jy
! Component: 350  -  total flux cleaned = 0.000383214 Jy
! Component: 400  -  total flux cleaned = 0.00051669 Jy
! Total flux subtracted in 400 components = 0.00051669 Jy
! Clean residual min=-0.001784 max=0.001309 Jy/beam
! Clean residual mean=-0.000000 rms=0.000348 Jy/beam
! Combined flux in latest and established models = 0.316508 Jy
selfcal
! Performing phase self-cal
! Adding 52 model components to the UV plane model.
! The established model now contains 379 components and 0.316508 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.074405Jy  sigma=0.146580
! Fit after  self-cal, rms=0.074370Jy  sigma=0.146545
wmodel /scr/pima/rdv62_uvs/J0339-0133_X_map.mod
! Writing 379 model components to file: /scr/pima/rdv62_uvs/J0339-0133_X_map.mod
wobs   /scr/pima/rdv62_uvs/J0339-0133_X_uvs.fits
! Applying 1288 buffered edits.
! Writing UV FITS file: /scr/pima/rdv62_uvs/J0339-0133_X_uvs.fits
wwins  /scr/pima/rdv62_uvs/J0339-0133_X_map.win
! wwins: Wrote 9 windows to /scr/pima/rdv62_uvs/J0339-0133_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.000482287 Jy
! Component: 100  -  total flux cleaned = -0.000506827 Jy
! Component: 150  -  total flux cleaned = -0.000462243 Jy
! Component: 200  -  total flux cleaned = -0.000419387 Jy
! Component: 250  -  total flux cleaned = -0.000356983 Jy
! Component: 300  -  total flux cleaned = -0.000296273 Jy
! Component: 350  -  total flux cleaned = -0.000276718 Jy
! Component: 400  -  total flux cleaned = -0.000200291 Jy
! Component: 450  -  total flux cleaned = -0.000181127 Jy
! Component: 500  -  total flux cleaned = -0.000143812 Jy
! Component: 550  -  total flux cleaned = -0.00012533 Jy
! Component: 600  -  total flux cleaned = -0.00014301 Jy
! Component: 650  -  total flux cleaned = -0.000160555 Jy
! Component: 700  -  total flux cleaned = -0.000177633 Jy
! Component: 750  -  total flux cleaned = -0.000177735 Jy
! Component: 800  -  total flux cleaned = -0.00022793 Jy
! Component: 850  -  total flux cleaned = -0.000260551 Jy
! Component: 900  -  total flux cleaned = -0.000260549 Jy
! Component: 950  -  total flux cleaned = -0.00029305 Jy
! Component: 1000  -  total flux cleaned = -0.00032503 Jy
! Component: 1050  -  total flux cleaned = -0.000341064 Jy
! Component: 1100  -  total flux cleaned = -0.000404222 Jy
! Component: 1150  -  total flux cleaned = -0.000497921 Jy
! Component: 1200  -  total flux cleaned = -0.000575345 Jy
! Component: 1250  -  total flux cleaned = -0.000606195 Jy
! Component: 1300  -  total flux cleaned = -0.000667302 Jy
! Component: 1350  -  total flux cleaned = -0.000667254 Jy
! Component: 1400  -  total flux cleaned = -0.000697531 Jy
! Component: 1450  -  total flux cleaned = -0.000712522 Jy
! Component: 1500  -  total flux cleaned = -0.000712529 Jy
! Component: 1550  -  total flux cleaned = -0.000712546 Jy
! Component: 1600  -  total flux cleaned = -0.000712525 Jy
! Component: 1650  -  total flux cleaned = -0.000712482 Jy
! Component: 1700  -  total flux cleaned = -0.000668202 Jy
! Component: 1750  -  total flux cleaned = -0.000624117 Jy
! Component: 1800  -  total flux cleaned = -0.000609537 Jy
! Component: 1850  -  total flux cleaned = -0.000594942 Jy
! Component: 1900  -  total flux cleaned = -0.000580392 Jy
! Component: 1950  -  total flux cleaned = -0.000536524 Jy
! Component: 2000  -  total flux cleaned = -0.000507364 Jy
! Component: 2050  -  total flux cleaned = -0.000463545 Jy
! Component: 2100  -  total flux cleaned = -0.000434417 Jy
! Component: 2150  -  total flux cleaned = -0.000434433 Jy
! Component: 2200  -  total flux cleaned = -0.000434536 Jy
! Component: 2250  -  total flux cleaned = -0.000405239 Jy
! Component: 2300  -  total flux cleaned = -0.00039058 Jy
! Component: 2350  -  total flux cleaned = -0.000375867 Jy
! Component: 2400  -  total flux cleaned = -0.000346555 Jy
! Component: 2450  -  total flux cleaned = -0.00031697 Jy
! Component: 2500  -  total flux cleaned = -0.000287479 Jy
! Component: 2550  -  total flux cleaned = -0.000243159 Jy
! Component: 2600  -  total flux cleaned = -0.000228396 Jy
! Component: 2650  -  total flux cleaned = -0.000213633 Jy
! Component: 2700  -  total flux cleaned = -0.000198787 Jy
! Component: 2750  -  total flux cleaned = -0.000183931 Jy
! Component: 2800  -  total flux cleaned = -0.000183936 Jy
! Component: 2850  -  total flux cleaned = -0.000153674 Jy
! Component: 2900  -  total flux cleaned = -0.000107952 Jy
! Component: 2950  -  total flux cleaned = -0.000123167 Jy
! Component: 3000  -  total flux cleaned = -7.72159e-05 Jy
! Component: 3050  -  total flux cleaned = -9.26622e-05 Jy
! Component: 3100  -  total flux cleaned = -7.71874e-05 Jy
! Component: 3150  -  total flux cleaned = -9.29645e-05 Jy
! Component: 3200  -  total flux cleaned = -9.30169e-05 Jy
! Component: 3250  -  total flux cleaned = -0.000108773 Jy
! Component: 3300  -  total flux cleaned = -0.000124988 Jy
! Component: 3350  -  total flux cleaned = -0.000141044 Jy
! Component: 3400  -  total flux cleaned = -0.000189841 Jy
! Component: 3450  -  total flux cleaned = -0.000189904 Jy
! Component: 3500  -  total flux cleaned = -0.000222967 Jy
! Component: 3550  -  total flux cleaned = -0.000239514 Jy
! Component: 3600  -  total flux cleaned = -0.000273107 Jy
! Component: 3650  -  total flux cleaned = -0.000290233 Jy
! Component: 3700  -  total flux cleaned = -0.000324571 Jy
! Component: 3750  -  total flux cleaned = -0.000324509 Jy
! Component: 3800  -  total flux cleaned = -0.000359741 Jy
! Component: 3850  -  total flux cleaned = -0.00041303 Jy
! Component: 3900  -  total flux cleaned = -0.000412858 Jy
! Component: 3950  -  total flux cleaned = -0.000430976 Jy
! Component: 4000  -  total flux cleaned = -0.000449349 Jy
! Component: 4050  -  total flux cleaned = -0.000467969 Jy
! Total flux subtracted in 4096 components = -0.000486917 Jy
! Clean residual min=-0.001040 max=0.001188 Jy/beam
! Clean residual mean=-0.000003 rms=0.000306 Jy/beam
! Combined flux in latest and established models = 0.316021 Jy
keep
! Adding 176 model components to the UV plane model.
! The established model now contains 552 components and 0.316021 Jy
![Exited script file: /home/lpetrov/bin/rdv_muppet]
wmap "/scr/pima/rdv62_uvs/J0339-0133_X\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 1.848 x 3.081 at -8.185 degrees (North through East)
! Clean map  min=-0.0045118  max=0.19578 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J0339-0133_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=  0.1957  rms=  0.00031  Jy 
quit
! Quitting program
! Log file /scr/pima/rdv62_uvs/J0339-0133_X_uvs.log closed on Fri Nov 22 22:32:28 2013