! Started logfile: /scr/pima/rdv93_uvs/J1107-4449_S_uvs.log on Tue Sep 24 17:06:57 2013
obs /scr/pima/rdv93_uvs/J1107-4449_S_uva.fits
! Reading UV FITS file: /scr/pima/rdv93_uvs/J1107-4449_S_uva.fits
! AN table 1: 23 integrations on 3 of 3 possible baselines.
! AN table 2: 11 integrations on 3 of 3 possible baselines.
! AN table 3: 11 integrations on 10 of 10 possible baselines.
! Apparent sampling: 0.778302 visibilities/baseline/integration-bin.
! Found source: J1107-4449
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  2.23693e+09      1.6e+07          1      1.6e+07
!  02        2  2.26693e+09      1.6e+07          1      1.6e+07
!  03        3  2.35693e+09      1.6e+07          1      1.6e+07
!  04        4  2.37693e+09      1.6e+07          1      1.6e+07
! 
! Polarization(s): RR
! 
! Read 0 lines of history.
! 
! Reading 660 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.5 
print "taper_size =",taper_size
! taper_size = 20 
![@/home/lpetrov/bin/01_muppet /scr/pima/rdv93_uvs/J1107-4449_S]
integer clean_niter
clean_niter = 100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float soltime1
soltime1 = 120
float soltime2
soltime2 = 20
float thresh
thresh = 0.5
float win_mult
win_mult = 1.8
float time_av
time_av = 16
float old_peak
float new_peak
float flux_cutoff
mapsize field_size, field_cell
! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize.
startmod
! Applying default point source starting model.
! Performing phase self-cal
! Adding 1 model components to the UV plane model.
! The established model now contains 1 components and 1 Jy
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  A total of 11 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  A total of 11 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  A total of 11 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  A total of 11 telescope corrections were flagged in sub-array 3.
! 
! Fit before self-cal, rms=1.827362Jy  sigma=18.144550
! Fit after  self-cal, rms=0.846233Jy  sigma=8.850061
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
uvw 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.862 mas, bmaj=28.42 mas, bpa=-15.38 degrees
! Estimated noise=6.46624 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 = 1.13318 Jy
! Component: 100  -  total flux cleaned = 1.6022 Jy
! Total flux subtracted in 100 components = 1.6022 Jy
! Clean residual min=-0.175812 max=0.205384 Jy/beam
! Clean residual mean=-0.006750 rms=0.069858 Jy/beam
! Combined flux in latest and established models = 1.60221 Jy
! Performing phase self-cal
! Adding 27 model components to the UV plane model.
! The established model now contains 27 components and 1.60221 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.402331Jy  sigma=4.101649
! Fit after  self-cal, rms=0.402140Jy  sigma=4.101054
! Inverting map 
print "************** FINISHED UNIFORM WEIGHTING CLEAN **************"
! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** 
! 
uvw 0,-1
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
win_mult = win_mult * 1.6
clean_niter = clean_niter * 2
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=4.711 mas, bmaj=39.93 mas, bpa=-15.8 degrees
! Estimated noise=3.32937 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.230245 Jy
! Component: 100  -  total flux cleaned = 0.187773 Jy
! Component: 150  -  total flux cleaned = 0.154589 Jy
! Component: 200  -  total flux cleaned = 0.14678 Jy
! Total flux subtracted in 200 components = 0.14678 Jy
! Clean residual min=-0.023264 max=0.023086 Jy/beam
! Clean residual mean=-0.001463 rms=0.011423 Jy/beam
! Combined flux in latest and established models = 1.74898 Jy
! Performing phase self-cal
! Adding 91 model components to the UV plane model.
! The established model now contains 106 components and 1.74898 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.262297Jy  sigma=2.415503
! Fit after  self-cal, rms=0.262204Jy  sigma=2.414959
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 66 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    LA        0.00*    PT        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.02     LA        1.07     MK        0.88     OV        1.00*    
!   PT        1.07     
! 
! 
! Correcting IF 2.
!  A total of 66 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    LA        0.00*    PT        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.04     LA        1.01     MK        0.96     OV        1.00*    
!   PT        1.06     
! 
! 
! Correcting IF 3.
!  A total of 66 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    LA        0.00*    PT        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   FD        1.03     LA        1.02     MK        1.01     OV        1.00*    
!   PT        1.01     
! 
! 
! Correcting IF 4.
!  A total of 66 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.00*    LA        0.00*    PT        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   KP        0.00*    MK        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   FD        0.94     LA        0.88     MK        1.13     OV        1.00*    
!   PT        0.86     
! 
! 
! Fit before self-cal, rms=0.262204Jy  sigma=2.414959
! Fit after  self-cal, rms=0.242061Jy  sigma=1.152553
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=3.17 mas, bmaj=62.32 mas, bpa=-11.84 degrees
! Estimated noise=6.85759 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.000682136 Jy
! Component: 100  -  total flux cleaned = 0.000613613 Jy
! Component: 150  -  total flux cleaned = 0.000537451 Jy
! Component: 200  -  total flux cleaned = -3.70782e-05 Jy
! Total flux subtracted in 200 components = -3.70782e-05 Jy
! Clean residual min=-0.009163 max=0.009065 Jy/beam
! Clean residual mean=0.000034 rms=0.002269 Jy/beam
! Combined flux in latest and established models = 1.74895 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 135 components and 1.74895 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.239178Jy  sigma=1.143659
! Fit after  self-cal, rms=0.239197Jy  sigma=1.143526
! Inverting map 
selfcal true, true, soltime1
! Performing amp+phase self-cal over 120 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.239197Jy  sigma=1.143526
! Fit after  self-cal, rms=0.236675Jy  sigma=1.142756
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.154 mas, bmaj=62.4 mas, bpa=-11.84 degrees
! Estimated noise=6.82663 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00185135 Jy
! Component: 100  -  total flux cleaned = 0.00185805 Jy
! Component: 150  -  total flux cleaned = 0.00255123 Jy
! Component: 200  -  total flux cleaned = 0.00450704 Jy
! Total flux subtracted in 200 components = 0.00450704 Jy
! Clean residual min=-0.006912 max=0.006970 Jy/beam
! Clean residual mean=-0.000050 rms=0.002102 Jy/beam
! Combined flux in latest and established models = 1.75345 Jy
selfcal
! Performing phase self-cal
! Adding 47 model components to the UV plane model.
! The established model now contains 180 components and 1.75345 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.234752Jy  sigma=1.139158
! Fit after  self-cal, rms=0.234794Jy  sigma=1.138844
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.00214609 Jy
! Component: 100  -  total flux cleaned = 0.002719 Jy
! Component: 150  -  total flux cleaned = 0.00391297 Jy
! Component: 200  -  total flux cleaned = 0.00391692 Jy
! Total flux subtracted in 200 components = 0.00391692 Jy
! Clean residual min=-0.006144 max=0.005956 Jy/beam
! Clean residual mean=-0.000177 rms=0.002096 Jy/beam
! Combined flux in latest and established models = 1.75737 Jy
! Performing phase self-cal
! Adding 50 model components to the UV plane model.
! The established model now contains 223 components and 1.75737 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.233684Jy  sigma=1.136382
! Fit after  self-cal, rms=0.233735Jy  sigma=1.136002
! Inverting map 
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.233735Jy  sigma=1.136002
! Fit after  self-cal, rms=0.229297Jy  sigma=1.134486
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.128 mas, bmaj=62.53 mas, bpa=-11.84 degrees
! Estimated noise=6.76511 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00180878 Jy
! Component: 100  -  total flux cleaned = 0.00212415 Jy
! Component: 150  -  total flux cleaned = 0.00119107 Jy
! Component: 200  -  total flux cleaned = -0.00129229 Jy
! Total flux subtracted in 200 components = -0.00129229 Jy
! Clean residual min=-0.006322 max=0.006518 Jy/beam
! Clean residual mean=0.000131 rms=0.001887 Jy/beam
! Combined flux in latest and established models = 1.75608 Jy
selfcal
! Performing phase self-cal
! Adding 54 model components to the UV plane model.
! The established model now contains 263 components and 1.75608 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.227558Jy  sigma=1.131676
! Fit after  self-cal, rms=0.227608Jy  sigma=1.131594
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.002924 Jy
! Component: 100  -  total flux cleaned = -0.0021434 Jy
! Component: 150  -  total flux cleaned = -0.000354833 Jy
! Component: 200  -  total flux cleaned = 0.00213322 Jy
! Total flux subtracted in 200 components = 0.00213322 Jy
! Clean residual min=-0.006549 max=0.006411 Jy/beam
! Clean residual mean=-0.000129 rms=0.002013 Jy/beam
! Combined flux in latest and established models = 1.75821 Jy
! Performing phase self-cal
! Adding 35 model components to the UV plane model.
! The established model now contains 281 components and 1.75821 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.226645Jy  sigma=1.128920
! Fit after  self-cal, rms=0.226675Jy  sigma=1.128583
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.226675Jy  sigma=1.128583
! Fit after  self-cal, rms=0.222419Jy  sigma=1.127095
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=18.38 mas, bmaj=61.24 mas, bpa=-4.541 degrees
! Estimated noise=6.78508 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00104351 Jy
! Component: 100  -  total flux cleaned = 0.00254262 Jy
! Component: 150  -  total flux cleaned = 0.00462723 Jy
! Component: 200  -  total flux cleaned = 0.00761264 Jy
! Total flux subtracted in 200 components = 0.00761264 Jy
! Clean residual min=-0.001495 max=0.002367 Jy/beam
! Clean residual mean=0.000031 rms=0.000867 Jy/beam
! Combined flux in latest and established models = 1.76583 Jy
selfcal
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 304 components and 1.76583 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.222121Jy  sigma=1.127050
! Fit after  self-cal, rms=0.222107Jy  sigma=1.127015
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.00079888 Jy
! Component: 100  -  total flux cleaned = 0.0027711 Jy
! Component: 150  -  total flux cleaned = 0.00567062 Jy
! Component: 200  -  total flux cleaned = 0.010266 Jy
! Total flux subtracted in 200 components = 0.010266 Jy
! Clean residual min=-0.002825 max=0.003890 Jy/beam
! Clean residual mean=0.000052 rms=0.001422 Jy/beam
! Combined flux in latest and established models = 1.77609 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 319 components and 1.77609 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.221946Jy  sigma=1.126721
! Fit after  self-cal, rms=0.221955Jy  sigma=1.126678
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
win_mult = win_mult / 1.2
print "************** FINISHED AMPLITUDE SELF-CAL **************"
! ************** FINISHED AMPLITUDE SELF-CAL ************** 
restore
! Inverting map and beam 
! Estimated beam: bmin=3.103 mas, bmaj=62.63 mas, bpa=-11.85 degrees
! Estimated noise=6.69878 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.103 x 62.63 at -11.85 degrees (North through East)
! Clean map  min=-0.03299  max=1.2249 Jy/beam
if(peak(flux) > thresh)
  selfcal true, true
  clean clean_niter,clean_gain
  selfcal
  flux_cutoff = imstat(rms) * dynam
  repeat
    if (peak(flux) > flux_cutoff) peakwin win_mult
    clean clean_niter,clean_gain
    flux_cutoff = imstat(rms) * dynam
    selfcal
    new_peak = peak(flux)
    until(new_peak<=flux_cutoff)
  selfcal true, true
  end if
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.221955Jy  sigma=1.126678
! Fit after  self-cal, rms=0.726571Jy  sigma=0.347212
! Inverting map and beam 
! Estimated beam: bmin=3.132 mas, bmaj=62.51 mas, bpa=-11.83 degrees
! Estimated noise=6.76191 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000850371 Jy
! Component: 100  -  total flux cleaned = -0.00201115 Jy
! Component: 150  -  total flux cleaned = -0.00177017 Jy
! Component: 200  -  total flux cleaned = -0.000626408 Jy
! Total flux subtracted in 200 components = -0.000626408 Jy
! Clean residual min=-0.001890 max=0.001872 Jy/beam
! Clean residual mean=-0.000014 rms=0.000771 Jy/beam
! Combined flux in latest and established models = 1.77547 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 360 components and 1.77547 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.725738Jy  sigma=0.344832
! Fit after  self-cal, rms=0.725967Jy  sigma=0.344723
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000750864 Jy
! Component: 100  -  total flux cleaned = 0.00127491 Jy
! Component: 150  -  total flux cleaned = 0.00127589 Jy
! Component: 200  -  total flux cleaned = 0.000176754 Jy
! Total flux subtracted in 200 components = 0.000176754 Jy
! Clean residual min=-0.002188 max=0.002197 Jy/beam
! Clean residual mean=0.000034 rms=0.000764 Jy/beam
! Combined flux in latest and established models = 1.77564 Jy
! Performing phase self-cal
! Adding 27 model components to the UV plane model.
! The established model now contains 385 components and 1.77564 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.725716Jy  sigma=0.344285
! Fit after  self-cal, rms=0.725915Jy  sigma=0.344199
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.725915Jy  sigma=0.344199
! Fit after  self-cal, rms=0.695647Jy  sigma=0.336456
print "************** CLEARING MODEL AND STARTING OVER **********"
! ************** CLEARING MODEL AND STARTING OVER ********** 
clrmod true
! clrmod: Cleared the established and tentative models.
uvw 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.326 mas, bmaj=69.64 mas, bpa=-11.98 degrees
! Estimated noise=11.8042 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.702532 Jy
! Component: 100  -  total flux cleaned = 1.08997 Jy
! Component: 150  -  total flux cleaned = 1.19824 Jy
! Component: 200  -  total flux cleaned = 1.18371 Jy
! Total flux subtracted in 200 components = 1.18371 Jy
! Clean residual min=-0.087438 max=0.086856 Jy/beam
! Clean residual mean=-0.009985 rms=0.043809 Jy/beam
! Combined flux in latest and established models = 1.18371 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 60 model components to the UV plane model.
! The established model now contains 60 components and 1.18371 Jy
! Inverting map and beam 
! Estimated beam: bmin=3.129 mas, bmaj=62.55 mas, bpa=-11.83 degrees
! Estimated noise=6.76103 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=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.120106 Jy
! Component: 100  -  total flux cleaned = 0.0978244 Jy
! Component: 150  -  total flux cleaned = 0.0622121 Jy
! Component: 200  -  total flux cleaned = 0.0222384 Jy
! Total flux subtracted in 200 components = 0.0222384 Jy
! Clean residual min=-0.028858 max=0.018743 Jy/beam
! Clean residual mean=-0.002766 rms=0.008826 Jy/beam
! Combined flux in latest and established models = 1.20595 Jy
! Adding 73 model components to the UV plane model.
! The established model now contains 121 components and 1.20595 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=18.58 mas, bmaj=61.41 mas, bpa=-4.339 degrees
! Estimated noise=6.78115 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00784728 Jy
! Component: 100  -  total flux cleaned = -0.0120843 Jy
! Component: 150  -  total flux cleaned = -0.0403478 Jy
! Component: 200  -  total flux cleaned = -0.0819612 Jy
! Total flux subtracted in 200 components = -0.0819612 Jy
! Clean residual min=-0.033384 max=0.022682 Jy/beam
! Clean residual mean=-0.000566 rms=0.010558 Jy/beam
! Combined flux in latest and established models = 1.12399 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=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0616273 Jy
! Component: 100  -  total flux cleaned = -0.153203 Jy
! Component: 150  -  total flux cleaned = -0.289323 Jy
! Component: 200  -  total flux cleaned = -0.491746 Jy
! Total flux subtracted in 200 components = -0.491746 Jy
! Clean residual min=-0.165149 max=0.136474 Jy/beam
! Clean residual mean=-0.001457 rms=0.071521 Jy/beam
! Combined flux in latest and established models = 0.632241 Jy
! Adding 81 model components to the UV plane model.
! The established model now contains 202 components and 0.632241 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=3.129 mas, bmaj=62.55 mas, bpa=-11.83 degrees
! Estimated noise=6.76103 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0241261 Jy
! Component: 100  -  total flux cleaned = 0.0723462 Jy
! Component: 150  -  total flux cleaned = 0.153869 Jy
! Component: 200  -  total flux cleaned = 0.266518 Jy
! Total flux subtracted in 200 components = 0.266518 Jy
! Clean residual min=-0.070145 max=0.084456 Jy/beam
! Clean residual mean=0.000136 rms=0.019294 Jy/beam
! Combined flux in latest and established models = 0.898759 Jy
selfcal
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 228 components and 0.89876 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.699681Jy  sigma=0.528553
! Fit after  self-cal, rms=0.696955Jy  sigma=0.514007
wmodel /scr/pima/rdv93_uvs/J1107-4449_S_map.mod
! Writing 228 model components to file: /scr/pima/rdv93_uvs/J1107-4449_S_map.mod
wobs /scr/pima/rdv93_uvs/J1107-4449_S_uvs.fits
! Writing UV FITS file: /scr/pima/rdv93_uvs/J1107-4449_S_uvs.fits
wwins /scr/pima/rdv93_uvs/J1107-4449_S_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
![Exited script file: /home/lpetrov/bin/01_muppet]
! 
wmap "/scr/pima/rdv93_uvs/J1107-4449_S\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.129 x 62.55 at -11.83 degrees (North through East)
! Clean map  min=-0.1371  max=1.14 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J1107-4449_S_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 1.13998 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.0153976 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -4.05209 4.05209 8.10418 16.2084 32.4167 64.8335
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0461 
print "peak=",int(1000*peak_flux)/1000
! peak= 1.139 
quit
! Quitting program
! Log file /scr/pima/rdv93_uvs/J1107-4449_S_uvs.log closed on Tue Sep 24 17:07:01 2013