! Started logfile: /scr/pima/rv123_uvs/J1709-3525_S_dfm.log on Sun Dec 23 05:05:13 2018
float field_size
field_size =  1024
float field_cell
field_cell =   0.5000
float taper_size
taper_size =   20.0
float freq
freq       =    2.23275
obs     /scr/pima/rv123_uvs/J1709-3525_S_uva.fits
! Reading UV FITS file: /scr/pima/rv123_uvs/J1709-3525_S_uva.fits
! AN table 1: 20 integrations on 1 of 1 possible baselines.
! AN table 2: 23 integrations on 6 of 6 possible baselines.
! AN table 3: 39 integrations on 10 of 10 possible baselines.
! Apparent sampling: 0.487226 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J1709-3525
! 
! 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.23275e+09        8e+06          1        8e+06
!  02        2  2.26275e+09        8e+06          1        8e+06
!  03        3  2.35275e+09        8e+06          1        8e+06
!  04        4  2.37275e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 89 lines of history.
! 
! Reading 1068 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 J1709-3525_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
print "==============================================================="
! =============================================================== 
! 
print "=== Difmap script: pima_mupet_01.dvm version of 2016.05.19 === "
! === Difmap script: pima_mupet_01.dvm version of 2016.05.19 ===  
! 
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 40 telescope corrections were flagged in sub-array 1.
!  A total of 19 telescope corrections were flagged in sub-array 2.
!  A total of 137 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 2.
!  A total of 40 telescope corrections were flagged in sub-array 1.
!  A total of 19 telescope corrections were flagged in sub-array 2.
!  A total of 137 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 3.
!  A total of 40 telescope corrections were flagged in sub-array 1.
!  A total of 92 telescope corrections were flagged in sub-array 2.
!  A total of 195 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 4.
!  A total of 40 telescope corrections were flagged in sub-array 1.
!  A total of 19 telescope corrections were flagged in sub-array 2.
!  A total of 137 telescope corrections were flagged in sub-array 3.
! 
! Fit before self-cal, rms=1.043627Jy  sigma=12.751271
! Fit after  self-cal, rms=0.812805Jy  sigma=9.967282
! 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=7.767 mas, bmaj=35.63 mas, bpa=10.56 degrees
! Estimated noise=4.06209 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.143585 Jy
! Component: 100  -  total flux cleaned = 0.203546 Jy
! Total flux subtracted in 100 components = 0.203546 Jy
! Clean residual min=-0.032044 max=0.032456 Jy/beam
! Clean residual mean=-0.000180 rms=0.012113 Jy/beam
! Combined flux in latest and established models = 0.203546 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 29 components and 0.203546 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.123940Jy  sigma=1.540698
! Fit after  self-cal, rms=0.123629Jy  sigma=1.535320
! 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=9.595 mas, bmaj=34.43 mas, bpa=9.144 degrees
! Estimated noise=2.35702 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.0297741 Jy
! Component: 100  -  total flux cleaned = 0.0429966 Jy
! Component: 150  -  total flux cleaned = 0.042906 Jy
! Component: 200  -  total flux cleaned = 0.0427679 Jy
! Total flux subtracted in 200 components = 0.0427679 Jy
! Clean residual min=-0.014693 max=0.014729 Jy/beam
! Clean residual mean=0.000095 rms=0.006584 Jy/beam
! Combined flux in latest and established models = 0.246314 Jy
! Performing phase self-cal
! Adding 130 model components to the UV plane model.
! The established model now contains 155 components and 0.246314 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.087280Jy  sigma=1.055273
! Fit after  self-cal, rms=0.086172Jy  sigma=1.043844
! Inverting map 
print "============== Finished natural weighting clean =============="
! ============== Finished natural weighting clean ============== 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 9.595 x 34.43 at 9.144 degrees (North through East)
! Clean map  min=-0.026711  max=0.18506 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   HN        0.00*    SC        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.90     KP        0.88     LA        1.13     OV        0.88     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.90     FD        0.98     KP        0.94     NL        1.27     
!   OV        0.98     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   HN        0.00*    SC        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        1.08     KP        0.87     LA        1.03     OV        0.74     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.64     FD        0.98     KP        0.88     NL        1.20     
!   OV        0.85     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   HN        0.00*    SC        0.00*    
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.00*    KP        0.00*    LA        0.00*    OV        0.00*    
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.00*    FD        0.00*    KP        0.00*    NL        0.00*    
!   OV        0.00*    
! 
! 
! Correcting IF 4.
!  Telescope amplitude corrections in sub-array 1:
!   HN        0.00*    SC        0.00*    
! 
!  A total of 33 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.90     KP        0.90     LA        1.04     OV        0.85     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.11     FD        1.00     KP        0.99     NL        1.39     
!   OV        0.99     
! 
! 
! Fit before self-cal, rms=0.086172Jy  sigma=1.043844
! Fit after  self-cal, rms=0.080377Jy  sigma=0.911543
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=9.541 mas, bmaj=36.39 mas, bpa=9.772 degrees
! Estimated noise=2.45929 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.0108457 Jy
! Component: 100  -  total flux cleaned = 0.0136115 Jy
! Component: 150  -  total flux cleaned = 0.0126451 Jy
! Component: 200  -  total flux cleaned = 0.0122397 Jy
! Total flux subtracted in 200 components = 0.0122397 Jy
! Clean residual min=-0.003138 max=0.003121 Jy/beam
! Clean residual mean=-0.000043 rms=0.001404 Jy/beam
! Combined flux in latest and established models = 0.258554 Jy
! Performing phase self-cal
! Adding 103 model components to the UV plane model.
! The established model now contains 258 components and 0.258554 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077822Jy  sigma=0.876812
! Fit after  self-cal, rms=0.077812Jy  sigma=0.874672
! 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.077812Jy  sigma=0.874672
! Fit after  self-cal, rms=0.076942Jy  sigma=0.873528
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=9.502 mas, bmaj=36.41 mas, bpa=9.724 degrees
! Estimated noise=2.44343 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00491161 Jy
! Component: 100  -  total flux cleaned = 0.00733083 Jy
! Component: 150  -  total flux cleaned = 0.00899416 Jy
! Component: 200  -  total flux cleaned = 0.00994786 Jy
! Total flux subtracted in 200 components = 0.00994786 Jy
! Clean residual min=-0.001664 max=0.001673 Jy/beam
! Clean residual mean=-0.000016 rms=0.000727 Jy/beam
! Combined flux in latest and established models = 0.268502 Jy
selfcal
! Performing phase self-cal
! Adding 103 model components to the UV plane model.
! The established model now contains 361 components and 0.268502 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.076395Jy  sigma=0.867536
! Fit after  self-cal, rms=0.076374Jy  sigma=0.867236
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.00155754 Jy
! Component: 100  -  total flux cleaned = 0.0020823 Jy
! Component: 150  -  total flux cleaned = 0.00279355 Jy
! Component: 200  -  total flux cleaned = 0.00329727 Jy
! Total flux subtracted in 200 components = 0.00329727 Jy
! Clean residual min=-0.001151 max=0.001147 Jy/beam
! Clean residual mean=-0.000014 rms=0.000502 Jy/beam
! Combined flux in latest and established models = 0.271799 Jy
! Performing phase self-cal
! Adding 83 model components to the UV plane model.
! The established model now contains 431 components and 0.271799 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.076158Jy  sigma=0.865307
! Fit after  self-cal, rms=0.076146Jy  sigma=0.865221
! 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.076146Jy  sigma=0.865221
! Fit after  self-cal, rms=0.077203Jy  sigma=0.865400
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=9.535 mas, bmaj=36.42 mas, bpa=9.702 degrees
! Estimated noise=2.45867 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000869414 Jy
! Component: 100  -  total flux cleaned = 0.00160328 Jy
! Component: 150  -  total flux cleaned = 0.00219922 Jy
! Component: 200  -  total flux cleaned = 0.00269796 Jy
! Total flux subtracted in 200 components = 0.00269796 Jy
! Clean residual min=-0.001001 max=0.001006 Jy/beam
! Clean residual mean=-0.000001 rms=0.000429 Jy/beam
! Combined flux in latest and established models = 0.274497 Jy
selfcal
! Performing phase self-cal
! Adding 92 model components to the UV plane model.
! The established model now contains 518 components and 0.274497 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077043Jy  sigma=0.864143
! Fit after  self-cal, rms=0.077033Jy  sigma=0.864074
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.000262223 Jy
! Component: 100  -  total flux cleaned = 0.000207288 Jy
! Component: 150  -  total flux cleaned = 0.000533396 Jy
! Component: 200  -  total flux cleaned = 0.000585322 Jy
! Total flux subtracted in 200 components = 0.000585322 Jy
! Clean residual min=-0.000853 max=0.000857 Jy/beam
! Clean residual mean=-0.000001 rms=0.000364 Jy/beam
! Combined flux in latest and established models = 0.275082 Jy
! Performing phase self-cal
! Adding 79 model components to the UV plane model.
! The established model now contains 582 components and 0.275082 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.076925Jy  sigma=0.863318
! Fit after  self-cal, rms=0.076922Jy  sigma=0.863296
! 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.076922Jy  sigma=0.863296
! Fit after  self-cal, rms=0.077404Jy  sigma=0.863335
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.86 mas, bmaj=36.42 mas, bpa=8.851 degrees
! Estimated noise=2.43439 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000152216 Jy
! Component: 100  -  total flux cleaned = -4.8883e-05 Jy
! Component: 150  -  total flux cleaned = 9.03039e-05 Jy
! Component: 200  -  total flux cleaned = 0.000223352 Jy
! Total flux subtracted in 200 components = 0.000223352 Jy
! Clean residual min=-0.000726 max=0.000728 Jy/beam
! Clean residual mean=0.000003 rms=0.000281 Jy/beam
! Combined flux in latest and established models = 0.275306 Jy
selfcal
! Performing phase self-cal
! Adding 72 model components to the UV plane model.
! The established model now contains 643 components and 0.275306 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077331Jy  sigma=0.862611
! Fit after  self-cal, rms=0.077330Jy  sigma=0.862578
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 = 9.43847e-05 Jy
! Component: 100  -  total flux cleaned = 0.00013096 Jy
! Component: 150  -  total flux cleaned = 0.000320134 Jy
! Component: 200  -  total flux cleaned = 0.000616737 Jy
! Total flux subtracted in 200 components = 0.000616737 Jy
! Clean residual min=-0.000610 max=0.000611 Jy/beam
! Clean residual mean=0.000003 rms=0.000246 Jy/beam
! Combined flux in latest and established models = 0.275922 Jy
! Performing phase self-cal
! Adding 68 model components to the UV plane model.
! The established model now contains 697 components and 0.275922 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077273Jy  sigma=0.862140
! Fit after  self-cal, rms=0.077268Jy  sigma=0.862125
! 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=9.544 mas, bmaj=36.45 mas, bpa=9.666 degrees
! Estimated noise=2.46537 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 9.544 x 36.45 at 9.666 degrees (North through East)
! Clean map  min=-0.012043  max=0.17245 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=7.977 mas, bmaj=37.06 mas, bpa=10.71 degrees
! Estimated noise=4.5523 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.131952 Jy
! Component: 100  -  total flux cleaned = 0.172896 Jy
! Component: 150  -  total flux cleaned = 0.188185 Jy
! Component: 200  -  total flux cleaned = 0.190843 Jy
! Total flux subtracted in 200 components = 0.190843 Jy
! Clean residual min=-0.017269 max=0.017209 Jy/beam
! Clean residual mean=0.000036 rms=0.007087 Jy/beam
! Combined flux in latest and established models = 0.190844 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 86 model components to the UV plane model.
! The established model now contains 86 components and 0.190844 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.544 mas, bmaj=36.45 mas, bpa=9.666 degrees
! Estimated noise=2.46537 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.0110379 Jy
! Component: 100  -  total flux cleaned = 0.015013 Jy
! Component: 150  -  total flux cleaned = 0.0186583 Jy
! Component: 200  -  total flux cleaned = 0.0245709 Jy
! Total flux subtracted in 200 components = 0.0245709 Jy
! Clean residual min=-0.007161 max=0.007220 Jy/beam
! Clean residual mean=0.000048 rms=0.003249 Jy/beam
! Combined flux in latest and established models = 0.215414 Jy
! Adding 100 model components to the UV plane model.
! The established model now contains 182 components and 0.215414 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.86 mas, bmaj=36.42 mas, bpa=8.851 degrees
! Estimated noise=2.43439 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00369294 Jy
! Component: 100  -  total flux cleaned = 0.0064367 Jy
! Component: 150  -  total flux cleaned = 0.00866631 Jy
! Component: 200  -  total flux cleaned = 0.00964006 Jy
! Total flux subtracted in 200 components = 0.00964006 Jy
! Clean residual min=-0.003026 max=0.003001 Jy/beam
! Clean residual mean=0.000018 rms=0.001485 Jy/beam
! Combined flux in latest and established models = 0.225054 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.000832534 Jy
! Component: 100  -  total flux cleaned = 0.00205289 Jy
! Component: 150  -  total flux cleaned = 0.00346187 Jy
! Component: 200  -  total flux cleaned = 0.00532733 Jy
! Total flux subtracted in 200 components = 0.00532733 Jy
! Clean residual min=-0.002168 max=0.002160 Jy/beam
! Clean residual mean=0.000027 rms=0.000851 Jy/beam
! Combined flux in latest and established models = 0.230382 Jy
! Adding 181 model components to the UV plane model.
! The established model now contains 356 components and 0.230382 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=9.544 mas, bmaj=36.45 mas, bpa=9.666 degrees
! Estimated noise=2.46537 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000397795 Jy
! Component: 100  -  total flux cleaned = 0.000707018 Jy
! Component: 150  -  total flux cleaned = 0.000803749 Jy
! Component: 200  -  total flux cleaned = 0.000648363 Jy
! Total flux subtracted in 200 components = 0.000648363 Jy
! Clean residual min=-0.001237 max=0.001236 Jy/beam
! Clean residual mean=-0.000012 rms=0.000497 Jy/beam
! Combined flux in latest and established models = 0.23103 Jy
selfcal
! Performing phase self-cal
! Adding 74 model components to the UV plane model.
! The established model now contains 430 components and 0.23103 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.077671Jy  sigma=0.865385
! Fit after  self-cal, rms=0.077609Jy  sigma=0.864887
wmodel J1709-3525_S_map.mod
! Writing 430 model components to file: J1709-3525_S_map.mod
wobs   J1709-3525_S_uvs.fits
! Writing UV FITS file: J1709-3525_S_uvs.fits
wwins  J1709-3525_S_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
! 
wmap   /scr/pima/rv123_uvs/J1709-3525_S_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 9.544 x 36.45 at 9.666 degrees (North through East)
! Clean map  min=-0.015036  max=0.16705 Jy/beam
! Writing clean map to FITS file: /scr/pima/rv123_uvs/J1709-3525_S_map.fits
wmodel /scr/pima/rv123_uvs/J1709-3525_S_map.mod
! Writing 430 model components to file: /scr/pima/rv123_uvs/J1709-3525_S_map.mod
wobs   /scr/pima/rv123_uvs/J1709-3525_S_uvs.fits
! Writing UV FITS file: /scr/pima/rv123_uvs/J1709-3525_S_uvs.fits
wwins  /scr/pima/rv123_uvs/J1709-3525_S_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
! 
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 0.167047 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000475834 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.854551 0.854551 1.7091 3.41821 6.83641 13.6728 27.3456 54.6913
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0014 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       0.167 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  351.061 
quit
! Quitting program
! Log file /scr/pima/rv123_uvs/J1709-3525_S_dfm.log closed on Sun Dec 23 05:05:16 2018