! Started logfile: /scr/pima/bg219h_uvs/J1642+6856_S_dfm.log on Thu May 19 22:33:32 2016
float field_size
field_size =  1024
float field_cell
field_cell =   0.5000
float taper_size
taper_size =   20.0
float freq
freq       =    2.22000
obs     /scr/pima/bg219h_uvs/J1642+6856_S_uva.fits
! Reading UV FITS file: /scr/pima/bg219h_uvs/J1642+6856_S_uva.fits
! AN table 1: 7 integrations on 36 of 36 possible baselines.
! Apparent sampling: 0.742063 visibilities/baseline/integration-bin.
! Found source: J1642+6856
! 
! 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.22e+09      3.2e+07          1      3.2e+07
!  02        2    2.252e+09      3.2e+07          1      3.2e+07
!  03        3    2.284e+09      3.2e+07          1      3.2e+07
!  04        4    2.348e+09      3.2e+07          1      3.2e+07
! 
! Polarization(s): RR
! 
! Read 54 lines of history.
! 
! Reading 748 visibilities.
select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 
! Restricting channel ranges to the available 4 channels.
! 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 
![@/opt64/bin/pima_mupet_01.dfm J1642+6856_S]
float lim_self
lim_self = 120.0
float lim_val
lim_va   = 1.2
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 15 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 15 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 15 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 15 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.705924Jy  sigma=24.380327
! Fit after  self-cal, rms=0.435022Jy  sigma=6.301713
! 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.483 mas, bmaj=9.626 mas, bpa=-69.07 degrees
! Estimated noise=3.58557 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.86102 Jy
! Component: 100  -  total flux cleaned = 1.05611 Jy
! Total flux subtracted in 100 components = 1.05611 Jy
! Clean residual min=-0.166840 max=0.251660 Jy/beam
! Clean residual mean=-0.000622 rms=0.047420 Jy/beam
! Combined flux in latest and established models = 1.05611 Jy
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 9 components and 1.05611 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.421338Jy  sigma=6.179606
! Fit after  self-cal, rms=0.421274Jy  sigma=6.178800
! 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=3.375 mas, bmaj=10.17 mas, bpa=-69.83 degrees
! Estimated noise=2.46719 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.0161012 Jy
! Component: 100  -  total flux cleaned = 7.7187e-05 Jy
! Component: 150  -  total flux cleaned = -0.00960756 Jy
! Component: 200  -  total flux cleaned = -0.0176842 Jy
! Total flux subtracted in 200 components = -0.0176842 Jy
! Clean residual min=-0.177887 max=0.267855 Jy/beam
! Clean residual mean=-0.001664 rms=0.055734 Jy/beam
! Combined flux in latest and established models = 1.03843 Jy
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 21 components and 1.03843 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.397427Jy  sigma=5.939388
! Fit after  self-cal, rms=0.397158Jy  sigma=5.937997
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.375 x 10.17 at -69.83 degrees (North through East)
! Clean map  min=-0.17588  max=1.1425 Jy/beam
if ( peak(flux,max)/imstat(rms) < lim_self ) selflim lim_val
! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2.
! Selfcal phase corrections will not be limited.
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00*    HN        1.00*    KP        1.00*    
!   LA        1.00*    MK        1.00*    NL        1.00*    OV        1.00*    
!   PT        1.00*    
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00*    HN        1.00*    KP        1.00*    
!   LA        1.00*    MK        1.00*    NL        1.00*    OV        1.00*    
!   PT        1.00*    
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00*    HN        1.00*    KP        1.00*    
!   LA        1.00*    MK        1.00*    NL        1.00*    OV        1.00*    
!   PT        1.00*    
! 
! 
! Correcting IF 4.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00*    HN        1.00*    KP        1.00*    
!   LA        1.00*    MK        1.00*    NL        1.00*    OV        1.00*    
!   PT        1.00*    
! 
! 
! Fit before self-cal, rms=0.397158Jy  sigma=5.937997
! Fit after  self-cal, rms=0.397158Jy  sigma=5.937997
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=3.375 mas, bmaj=10.17 mas, bpa=-69.83 degrees
! Estimated noise=2.46719 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.00809564 Jy
! Component: 100  -  total flux cleaned = -0.00987597 Jy
! Component: 150  -  total flux cleaned = -0.010691 Jy
! Component: 200  -  total flux cleaned = -0.011484 Jy
! Total flux subtracted in 200 components = -0.011484 Jy
! Clean residual min=-0.179947 max=0.264277 Jy/beam
! Clean residual mean=-0.001642 rms=0.055290 Jy/beam
! Combined flux in latest and established models = 1.02694 Jy
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 35 components and 1.02694 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.394260Jy  sigma=5.895186
! Fit after  self-cal, rms=0.394239Jy  sigma=5.894745
! 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.394239Jy  sigma=5.894745
! Fit after  self-cal, rms=0.394239Jy  sigma=5.894745
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.375 mas, bmaj=10.17 mas, bpa=-69.83 degrees
! Estimated noise=2.46719 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000775577 Jy
! Component: 100  -  total flux cleaned = -0.00150055 Jy
! Component: 150  -  total flux cleaned = -0.00219543 Jy
! Component: 200  -  total flux cleaned = -0.00288201 Jy
! Total flux subtracted in 200 components = -0.00288201 Jy
! Clean residual min=-0.179715 max=0.262173 Jy/beam
! Clean residual mean=-0.001638 rms=0.054928 Jy/beam
! Combined flux in latest and established models = 1.02406 Jy
selfcal
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 44 components and 1.02406 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.392263Jy  sigma=5.863863
! Fit after  self-cal, rms=0.392290Jy  sigma=5.863650
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 (61.5, -60.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.225932 Jy
! Component: 100  -  total flux cleaned = 0.326043 Jy
! Component: 150  -  total flux cleaned = 0.353408 Jy
! Component: 200  -  total flux cleaned = 0.378063 Jy
! Total flux subtracted in 200 components = 0.378063 Jy
! Clean residual min=-0.112994 max=0.148413 Jy/beam
! Clean residual mean=-0.000539 rms=0.029061 Jy/beam
! Combined flux in latest and established models = 1.40212 Jy
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 75 components and 1.40212 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.254220Jy  sigma=3.667928
! Fit after  self-cal, rms=0.238794Jy  sigma=3.480308
! 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.238794Jy  sigma=3.480308
! Fit after  self-cal, rms=0.238794Jy  sigma=3.480308
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.375 mas, bmaj=10.17 mas, bpa=-69.83 degrees
! Estimated noise=2.46719 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0505571 Jy
! Component: 100  -  total flux cleaned = 0.071276 Jy
! Component: 150  -  total flux cleaned = 0.0819736 Jy
! Component: 200  -  total flux cleaned = 0.089144 Jy
! Total flux subtracted in 200 components = 0.089144 Jy
! Clean residual min=-0.102973 max=0.101723 Jy/beam
! Clean residual mean=-0.000378 rms=0.023422 Jy/beam
! Combined flux in latest and established models = 1.49127 Jy
selfcal
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 102 components and 1.49127 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.220500Jy  sigma=3.158645
! Fit after  self-cal, rms=0.218195Jy  sigma=3.139039
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.00872731 Jy
! Component: 100  -  total flux cleaned = 0.0159208 Jy
! Component: 150  -  total flux cleaned = 0.02059 Jy
! Component: 200  -  total flux cleaned = 0.0248388 Jy
! Total flux subtracted in 200 components = 0.0248388 Jy
! Clean residual min=-0.096769 max=0.086587 Jy/beam
! Clean residual mean=-0.000339 rms=0.021698 Jy/beam
! Combined flux in latest and established models = 1.5161 Jy
! Performing phase self-cal
! Adding 25 model components to the UV plane model.
! The established model now contains 114 components and 1.5161 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211185Jy  sigma=3.006879
! Fit after  self-cal, rms=0.210230Jy  sigma=3.002412
! 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.210230Jy  sigma=3.002412
! Fit after  self-cal, rms=0.204546Jy  sigma=2.912805
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.486 mas, bmaj=14.09 mas, bpa=-87.57 degrees
! Estimated noise=3.06206 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00751139 Jy
! Component: 100  -  total flux cleaned = 0.0184025 Jy
! Component: 150  -  total flux cleaned = 0.0196635 Jy
! Component: 200  -  total flux cleaned = 0.0207474 Jy
! Total flux subtracted in 200 components = 0.0207474 Jy
! Clean residual min=-0.085290 max=0.082917 Jy/beam
! Clean residual mean=-0.000620 rms=0.026084 Jy/beam
! Combined flux in latest and established models = 1.53685 Jy
selfcal
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 136 components and 1.53685 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.202944Jy  sigma=2.891599
! Fit after  self-cal, rms=0.201357Jy  sigma=2.878520
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.00349438 Jy
! Component: 100  -  total flux cleaned = 0.00503835 Jy
! Component: 150  -  total flux cleaned = 0.00593361 Jy
! Component: 200  -  total flux cleaned = 0.00709546 Jy
! Total flux subtracted in 200 components = 0.00709546 Jy
! Clean residual min=-0.081663 max=0.083457 Jy/beam
! Clean residual mean=-0.000636 rms=0.025216 Jy/beam
! Combined flux in latest and established models = 1.54395 Jy
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 140 components and 1.54395 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.201244Jy  sigma=2.879505
! Fit after  self-cal, rms=0.200502Jy  sigma=2.875307
! 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.375 mas, bmaj=10.15 mas, bpa=-69.8 degrees
! Estimated noise=2.49145 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.375 x 10.15 at -69.8 degrees (North through East)
! Clean map  min=-0.077403  max=1.1616 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.200502Jy  sigma=2.875307
! Fit after  self-cal, rms=0.198994Jy  sigma=2.838917
! Inverting map and beam 
! Estimated beam: bmin=3.376 mas, bmaj=10.14 mas, bpa=-69.77 degrees
! Estimated noise=2.49857 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0067176 Jy
! Component: 100  -  total flux cleaned = -0.00870815 Jy
! Component: 150  -  total flux cleaned = -0.0092766 Jy
! Component: 200  -  total flux cleaned = -0.00932065 Jy
! Total flux subtracted in 200 components = -0.00932065 Jy
! Clean residual min=-0.074477 max=0.067474 Jy/beam
! Clean residual mean=-0.000297 rms=0.017879 Jy/beam
! Combined flux in latest and established models = 1.53463 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 159 components and 1.53463 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.195609Jy  sigma=2.783704
! Fit after  self-cal, rms=0.195307Jy  sigma=2.781596
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0015695 Jy
! Component: 100  -  total flux cleaned = -0.00204426 Jy
! Component: 150  -  total flux cleaned = -0.00247324 Jy
! Component: 200  -  total flux cleaned = -0.00172303 Jy
! Total flux subtracted in 200 components = -0.00172303 Jy
! Clean residual min=-0.074324 max=0.064597 Jy/beam
! Clean residual mean=-0.000286 rms=0.017783 Jy/beam
! Combined flux in latest and established models = 1.5329 Jy
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 166 components and 1.5329 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.193702Jy  sigma=2.756454
! Fit after  self-cal, rms=0.193669Jy  sigma=2.755276
! 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.193669Jy  sigma=2.755276
! Fit after  self-cal, rms=0.191065Jy  sigma=2.702390
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.504 mas, bmaj=9.581 mas, bpa=-69.05 degrees
! Estimated noise=3.66692 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.883373 Jy
! Component: 100  -  total flux cleaned = 1.1717 Jy
! Component: 150  -  total flux cleaned = 1.31011 Jy
! Component: 200  -  total flux cleaned = 1.40205 Jy
! Total flux subtracted in 200 components = 1.40205 Jy
! Clean residual min=-0.097153 max=0.115357 Jy/beam
! Clean residual mean=-0.000166 rms=0.024833 Jy/beam
! Combined flux in latest and established models = 1.40205 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 26 model components to the UV plane model.
! The established model now contains 26 components and 1.40205 Jy
! Inverting map and beam 
! Estimated beam: bmin=3.377 mas, bmaj=10.13 mas, bpa=-69.74 degrees
! Estimated noise=2.51305 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 (-61, 60.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0866829 Jy
! Component: 100  -  total flux cleaned = 0.112192 Jy
! Component: 150  -  total flux cleaned = 0.133615 Jy
! Component: 200  -  total flux cleaned = 0.142707 Jy
! Total flux subtracted in 200 components = 0.142707 Jy
! Clean residual min=-0.082472 max=0.072072 Jy/beam
! Clean residual mean=-0.000144 rms=0.017645 Jy/beam
! Combined flux in latest and established models = 1.54475 Jy
! Adding 37 model components to the UV plane model.
! The established model now contains 61 components and 1.54475 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=7.484 mas, bmaj=14.05 mas, bpa=-87.5 degrees
! Estimated noise=3.0946 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0558569 Jy
! Component: 100  -  total flux cleaned = 0.0775094 Jy
! Component: 150  -  total flux cleaned = 0.0836788 Jy
! Component: 200  -  total flux cleaned = 0.0877191 Jy
! Total flux subtracted in 200 components = 0.0877191 Jy
! Clean residual min=-0.069778 max=0.051763 Jy/beam
! Clean residual mean=-0.000120 rms=0.016982 Jy/beam
! Combined flux in latest and established models = 1.63247 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.00254457 Jy
! Component: 100  -  total flux cleaned = 0.00328658 Jy
! Component: 150  -  total flux cleaned = 0.00459476 Jy
! Component: 200  -  total flux cleaned = 0.00685199 Jy
! Total flux subtracted in 200 components = 0.00685199 Jy
! Clean residual min=-0.060046 max=0.046469 Jy/beam
! Clean residual mean=-0.000121 rms=0.014978 Jy/beam
! Combined flux in latest and established models = 1.63933 Jy
! Adding 33 model components to the UV plane model.
! The established model now contains 90 components and 1.63933 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=3.377 mas, bmaj=10.13 mas, bpa=-69.74 degrees
! Estimated noise=2.51305 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0156484 Jy
! Component: 100  -  total flux cleaned = -0.0204139 Jy
! Component: 150  -  total flux cleaned = -0.0234128 Jy
! Component: 200  -  total flux cleaned = -0.0261338 Jy
! Total flux subtracted in 200 components = -0.0261338 Jy
! Clean residual min=-0.058700 max=0.046987 Jy/beam
! Clean residual mean=-0.000088 rms=0.012966 Jy/beam
! Combined flux in latest and established models = 1.61319 Jy
selfcal
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 124 components and 1.61319 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.173248Jy  sigma=2.383486
! Fit after  self-cal, rms=0.167170Jy  sigma=2.298101
wmodel /J1642+6856_S_map.mod
! wmodel: Unable to open new model file: /J1642+6856_S_map.mod
! Error occured in command: wmodel
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
! 
wmap /scr/pima/bg219h_uvs/J1642+6856_S_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.377 x 10.13 at -69.74 degrees (North through East)
! Clean map  min=-0.062677  max=1.0839 Jy/beam
! Writing clean map to FITS file: /scr/pima/bg219h_uvs/J1642+6856_S_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 1.08393 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.0125845 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -3.48302 3.48302 6.96604 13.9321 27.8642 55.7283
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0377 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       1.083 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  86.132 
quit
! Quitting program
! Log file /scr/pima/bg219h_uvs/J1642+6856_S_dfm.log closed on Thu May 19 22:33:34 2016