! Started logfile: /vlbi/uf001/uva//J0423-0120_S_b_dfm.log on Sun Jan 13 22:10:00 2019
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     /vlbi/uf001/uva//J0423-0120_S_b_uva.fits
! Reading UV FITS file: /vlbi/uf001/uva//J0423-0120_S_b_uva.fits
! AN table 1: 14 integrations on 28 of 28 possible baselines.
! Apparent sampling: 0.770408 visibilities/baseline/integration-bin.
! Found source: J0423-0120
! 
! There are 3 IFs, and a total of 3 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
! 
! Polarization(s): RR
! 
! Read 40 lines of history.
! 
! Reading 906 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 3 channels.
! Selecting polarization: RR,  channels: 1..3
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
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 J0423-0120_S_b]
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 61 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 61 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 61 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.795209Jy  sigma=20.260994
! Fit after  self-cal, rms=0.552803Jy  sigma=5.996533
! 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=3.145 mas, bmaj=14.05 mas, bpa=-17.42 degrees
! Estimated noise=3.28383 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 = 1.09126 Jy
! Component: 100  -  total flux cleaned = 1.4085 Jy
! Total flux subtracted in 100 components = 1.4085 Jy
! Clean residual min=-0.068336 max=0.131193 Jy/beam
! Clean residual mean=0.000194 rms=0.019562 Jy/beam
! Combined flux in latest and established models = 1.4085 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 15 components and 1.4085 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.247245Jy  sigma=2.728580
! Fit after  self-cal, rms=0.246654Jy  sigma=2.723503
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.132795 Jy
! Component: 100  -  total flux cleaned = 0.188963 Jy
! Total flux subtracted in 100 components = 0.188963 Jy
! Clean residual min=-0.047239 max=0.040589 Jy/beam
! Clean residual mean=0.000067 rms=0.010850 Jy/beam
! Combined flux in latest and established models = 1.59747 Jy
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 23 components and 1.59747 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.163397Jy  sigma=1.800908
! Fit after  self-cal, rms=0.162935Jy  sigma=1.796511
! 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.473 mas, bmaj=15.02 mas, bpa=-17.3 degrees
! Estimated noise=2.96643 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.0333684 Jy
! Component: 100  -  total flux cleaned = 0.0545698 Jy
! Component: 150  -  total flux cleaned = 0.0671441 Jy
! Component: 200  -  total flux cleaned = 0.0697104 Jy
! Total flux subtracted in 200 components = 0.0697104 Jy
! Clean residual min=-0.042497 max=0.040179 Jy/beam
! Clean residual mean=0.000073 rms=0.010301 Jy/beam
! Combined flux in latest and established models = 1.66718 Jy
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 50 components and 1.66718 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.154830Jy  sigma=1.700425
! Fit after  self-cal, rms=0.154771Jy  sigma=1.698867
! 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.473 x 15.02 at -17.3 degrees (North through East)
! Clean map  min=-0.041978  max=1.4136 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.07     HN        0.98     KP        1.04     LA        0.90     
!   MK        1.01     NL        1.04     OV        1.01     PT        1.03     
! 
! 
! Correcting IF 2.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.90     HN        0.99     KP        1.04     LA        0.91     
!   MK        1.01     NL        1.06     OV        0.88     PT        1.05     
! 
! 
! Correcting IF 3.
!  A total of 6 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.05     HN        1.02     KP        1.04     LA        0.89     
!   MK        1.02     NL        1.08     OV        1.01     PT        1.04     
! 
! 
! Fit before self-cal, rms=0.154771Jy  sigma=1.698867
! Fit after  self-cal, rms=0.071807Jy  sigma=0.782821
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=3.559 mas, bmaj=14.72 mas, bpa=-16.91 degrees
! Estimated noise=3.03327 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.0113076 Jy
! Component: 100  -  total flux cleaned = -0.0123198 Jy
! Component: 150  -  total flux cleaned = -0.0123097 Jy
! Component: 200  -  total flux cleaned = -0.0129006 Jy
! Total flux subtracted in 200 components = -0.0129006 Jy
! Clean residual min=-0.010988 max=0.014706 Jy/beam
! Clean residual mean=0.000012 rms=0.003163 Jy/beam
! Combined flux in latest and established models = 1.65428 Jy
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 73 components and 1.65428 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.068081Jy  sigma=0.744035
! Fit after  self-cal, rms=0.068072Jy  sigma=0.743091
! 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.
! 
! Fit before self-cal, rms=0.068072Jy  sigma=0.743091
! Fit after  self-cal, rms=0.064792Jy  sigma=0.699001
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.56 mas, bmaj=14.74 mas, bpa=-16.92 degrees
! Estimated noise=3.03597 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000673897 Jy
! Component: 100  -  total flux cleaned = -0.00135832 Jy
! Component: 150  -  total flux cleaned = -0.00153677 Jy
! Component: 200  -  total flux cleaned = -0.00161412 Jy
! Total flux subtracted in 200 components = -0.00161412 Jy
! Clean residual min=-0.008454 max=0.011027 Jy/beam
! Clean residual mean=0.000004 rms=0.002486 Jy/beam
! Combined flux in latest and established models = 1.65266 Jy
selfcal
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 83 components and 1.65266 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.064347Jy  sigma=0.694272
! Fit after  self-cal, rms=0.064362Jy  sigma=0.693962
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 (2, -32.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0119599 Jy
! Component: 100  -  total flux cleaned = 0.0144423 Jy
! Component: 150  -  total flux cleaned = 0.0151554 Jy
! Component: 200  -  total flux cleaned = 0.0153363 Jy
! Total flux subtracted in 200 components = 0.0153363 Jy
! Clean residual min=-0.006438 max=0.007124 Jy/beam
! Clean residual mean=-0.000006 rms=0.001792 Jy/beam
! Combined flux in latest and established models = 1.668 Jy
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 120 components and 1.668 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.061369Jy  sigma=0.665453
! Fit after  self-cal, rms=0.061164Jy  sigma=0.664016
! 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.
! 
! Fit before self-cal, rms=0.061164Jy  sigma=0.664016
! Fit after  self-cal, rms=0.060822Jy  sigma=0.660759
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.564 mas, bmaj=14.72 mas, bpa=-16.91 degrees
! Estimated noise=3.03574 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00300364 Jy
! Component: 100  -  total flux cleaned = 0.00318327 Jy
! Component: 150  -  total flux cleaned = 0.00234042 Jy
! Component: 200  -  total flux cleaned = 0.00159549 Jy
! Total flux subtracted in 200 components = 0.00159549 Jy
! Clean residual min=-0.006137 max=0.005930 Jy/beam
! Clean residual mean=-0.000008 rms=0.001541 Jy/beam
! Combined flux in latest and established models = 1.66959 Jy
selfcal
! Performing phase self-cal
! Adding 39 model components to the UV plane model.
! The established model now contains 157 components and 1.66959 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.060186Jy  sigma=0.654378
! Fit after  self-cal, rms=0.060152Jy  sigma=0.654087
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 = 4.54098e-05 Jy
! Component: 100  -  total flux cleaned = -0.000496438 Jy
! Component: 150  -  total flux cleaned = -0.00114495 Jy
! Component: 200  -  total flux cleaned = -0.00181171 Jy
! Total flux subtracted in 200 components = -0.00181171 Jy
! Clean residual min=-0.006100 max=0.005615 Jy/beam
! Clean residual mean=-0.000008 rms=0.001471 Jy/beam
! Combined flux in latest and established models = 1.66778 Jy
! Performing phase self-cal
! Adding 34 model components to the UV plane model.
! The established model now contains 172 components and 1.66778 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.059912Jy  sigma=0.651846
! Fit after  self-cal, rms=0.059908Jy  sigma=0.651773
! 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.
! 
! Fit before self-cal, rms=0.059908Jy  sigma=0.651773
! Fit after  self-cal, rms=0.059850Jy  sigma=0.650915
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.413 mas, bmaj=15.74 mas, bpa=-10.81 degrees
! Estimated noise=3.75338 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000344804 Jy
! Component: 100  -  total flux cleaned = -0.00118867 Jy
! Component: 150  -  total flux cleaned = -0.00186162 Jy
! Component: 200  -  total flux cleaned = -0.00243845 Jy
! Total flux subtracted in 200 components = -0.00243845 Jy
! Clean residual min=-0.005516 max=0.005249 Jy/beam
! Clean residual mean=-0.000009 rms=0.001521 Jy/beam
! Combined flux in latest and established models = 1.66534 Jy
selfcal
! Performing phase self-cal
! Adding 36 model components to the UV plane model.
! The established model now contains 204 components and 1.66534 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.059677Jy  sigma=0.649887
! Fit after  self-cal, rms=0.059660Jy  sigma=0.649762
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.000413903 Jy
! Component: 100  -  total flux cleaned = -0.000713186 Jy
! Component: 150  -  total flux cleaned = -0.000973631 Jy
! Component: 200  -  total flux cleaned = -0.00117045 Jy
! Total flux subtracted in 200 components = -0.00117045 Jy
! Clean residual min=-0.005521 max=0.005170 Jy/beam
! Clean residual mean=-0.000009 rms=0.001500 Jy/beam
! Combined flux in latest and established models = 1.66417 Jy
! Performing phase self-cal
! Adding 25 model components to the UV plane model.
! The established model now contains 220 components and 1.66417 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.059584Jy  sigma=0.649628
! Fit after  self-cal, rms=0.059589Jy  sigma=0.649589
! 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.564 mas, bmaj=14.72 mas, bpa=-16.91 degrees
! Estimated noise=3.03547 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.564 x 14.72 at -16.91 degrees (North through East)
! Clean map  min=-0.011581  max=1.4211 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.
! 
! Fit before self-cal, rms=0.059589Jy  sigma=0.649589
! Fit after  self-cal, rms=0.057186Jy  sigma=0.612285
! Inverting map and beam 
! Estimated beam: bmin=3.565 mas, bmaj=14.73 mas, bpa=-16.91 degrees
! Estimated noise=3.03482 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000168393 Jy
! Component: 100  -  total flux cleaned = -0.000771712 Jy
! Component: 150  -  total flux cleaned = -0.0014011 Jy
! Component: 200  -  total flux cleaned = -0.00177603 Jy
! Total flux subtracted in 200 components = -0.00177603 Jy
! Clean residual min=-0.005623 max=0.005428 Jy/beam
! Clean residual mean=-0.000007 rms=0.001362 Jy/beam
! Combined flux in latest and established models = 1.6624 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 241 components and 1.6624 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.057061Jy  sigma=0.610992
! Fit after  self-cal, rms=0.057074Jy  sigma=0.610899
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000405727 Jy
! Component: 100  -  total flux cleaned = -0.000746903 Jy
! Component: 150  -  total flux cleaned = -0.000941087 Jy
! Component: 200  -  total flux cleaned = -0.0011602 Jy
! Total flux subtracted in 200 components = -0.0011602 Jy
! Clean residual min=-0.005559 max=0.005333 Jy/beam
! Clean residual mean=-0.000007 rms=0.001348 Jy/beam
! Combined flux in latest and established models = 1.66124 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 251 components and 1.66124 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.057011Jy  sigma=0.610407
! Fit after  self-cal, rms=0.057017Jy  sigma=0.610366
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.057017Jy  sigma=0.610366
! Fit after  self-cal, rms=0.057010Jy  sigma=0.610105
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=3.174 mas, bmaj=14.06 mas, bpa=-17.31 degrees
! Estimated noise=3.3395 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.09355 Jy
! Component: 100  -  total flux cleaned = 1.40458 Jy
! Component: 150  -  total flux cleaned = 1.53301 Jy
! Component: 200  -  total flux cleaned = 1.5869 Jy
! Total flux subtracted in 200 components = 1.5869 Jy
! Clean residual min=-0.014038 max=0.022919 Jy/beam
! Clean residual mean=0.000038 rms=0.004224 Jy/beam
! Combined flux in latest and established models = 1.5869 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 19 model components to the UV plane model.
! The established model now contains 19 components and 1.5869 Jy
! Inverting map and beam 
! Estimated beam: bmin=3.566 mas, bmaj=14.73 mas, bpa=-16.91 degrees
! Estimated noise=3.03472 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.0295069 Jy
! Component: 100  -  total flux cleaned = 0.0475421 Jy
! Component: 150  -  total flux cleaned = 0.0545944 Jy
! Component: 200  -  total flux cleaned = 0.0589336 Jy
! Total flux subtracted in 200 components = 0.0589336 Jy
! Clean residual min=-0.009478 max=0.009623 Jy/beam
! Clean residual mean=0.000009 rms=0.002370 Jy/beam
! Combined flux in latest and established models = 1.64583 Jy
! Adding 44 model components to the UV plane model.
! The established model now contains 60 components and 1.64583 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.416 mas, bmaj=15.75 mas, bpa=-10.8 degrees
! Estimated noise=3.75044 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00903262 Jy
! Component: 100  -  total flux cleaned = 0.0118465 Jy
! Component: 150  -  total flux cleaned = 0.0124001 Jy
! Component: 200  -  total flux cleaned = 0.0128442 Jy
! Total flux subtracted in 200 components = 0.0128442 Jy
! Clean residual min=-0.007863 max=0.005907 Jy/beam
! Clean residual mean=-0.000002 rms=0.001981 Jy/beam
! Combined flux in latest and established models = 1.65868 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.000374484 Jy
! Component: 100  -  total flux cleaned = 0.00070595 Jy
! Component: 150  -  total flux cleaned = 0.000997192 Jy
! Component: 200  -  total flux cleaned = 0.0012523 Jy
! Total flux subtracted in 200 components = 0.0012523 Jy
! Clean residual min=-0.006825 max=0.005389 Jy/beam
! Clean residual mean=-0.000006 rms=0.001698 Jy/beam
! Combined flux in latest and established models = 1.65993 Jy
! Adding 89 model components to the UV plane model.
! The established model now contains 144 components and 1.65993 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=3.566 mas, bmaj=14.73 mas, bpa=-16.91 degrees
! Estimated noise=3.03472 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000131166 Jy
! Component: 100  -  total flux cleaned = 0.000528817 Jy
! Component: 150  -  total flux cleaned = 0.000203991 Jy
! Component: 200  -  total flux cleaned = -9.04999e-05 Jy
! Total flux subtracted in 200 components = -9.04999e-05 Jy
! Clean residual min=-0.006459 max=0.005232 Jy/beam
! Clean residual mean=-0.000005 rms=0.001500 Jy/beam
! Combined flux in latest and established models = 1.65984 Jy
selfcal
! Performing phase self-cal
! Adding 22 model components to the UV plane model.
! The established model now contains 164 components and 1.65984 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.057478Jy  sigma=0.616561
! Fit after  self-cal, rms=0.057397Jy  sigma=0.615853
wmodel J0423-0120_S_b_map.mod
! Writing 164 model components to file: J0423-0120_S_b_map.mod
wobs   J0423-0120_S_b_uvs.fits
! Writing UV FITS file: J0423-0120_S_b_uvs.fits
wwins  J0423-0120_S_b_map.win
! wwins: Wrote 2 windows to J0423-0120_S_b_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.0027404 Jy
! Component: 100  -  total flux cleaned = -0.00418998 Jy
! Component: 150  -  total flux cleaned = -0.00490456 Jy
! Component: 200  -  total flux cleaned = -0.00540758 Jy
! Component: 250  -  total flux cleaned = -0.00620408 Jy
! Component: 300  -  total flux cleaned = -0.00681485 Jy
! Component: 350  -  total flux cleaned = -0.0074729 Jy
! Component: 400  -  total flux cleaned = -0.00782377 Jy
! Component: 450  -  total flux cleaned = -0.00816147 Jy
! Component: 500  -  total flux cleaned = -0.00855048 Jy
! Component: 550  -  total flux cleaned = -0.00886311 Jy
! Component: 600  -  total flux cleaned = -0.00916708 Jy
! Component: 650  -  total flux cleaned = -0.00940475 Jy
! Component: 700  -  total flux cleaned = -0.00992334 Jy
! Component: 750  -  total flux cleaned = -0.0102586 Jy
! Component: 800  -  total flux cleaned = -0.0105861 Jy
! Component: 850  -  total flux cleaned = -0.0109581 Jy
! Component: 900  -  total flux cleaned = -0.0111134 Jy
! Component: 950  -  total flux cleaned = -0.0113667 Jy
! Component: 1000  -  total flux cleaned = -0.0115156 Jy
! Component: 1050  -  total flux cleaned = -0.0117595 Jy
! Component: 1100  -  total flux cleaned = -0.0119019 Jy
! Component: 1150  -  total flux cleaned = -0.0119496 Jy
! Component: 1200  -  total flux cleaned = -0.0120862 Jy
! Component: 1250  -  total flux cleaned = -0.0121312 Jy
! Component: 1300  -  total flux cleaned = -0.0122188 Jy
! Component: 1350  -  total flux cleaned = -0.0122182 Jy
! Component: 1400  -  total flux cleaned = -0.0120898 Jy
! Component: 1450  -  total flux cleaned = -0.0120901 Jy
! Component: 1500  -  total flux cleaned = -0.0121735 Jy
! Component: 1550  -  total flux cleaned = -0.0121327 Jy
! Component: 1600  -  total flux cleaned = -0.0121726 Jy
! Component: 1650  -  total flux cleaned = -0.0121342 Jy
! Component: 1700  -  total flux cleaned = -0.012212 Jy
! Component: 1750  -  total flux cleaned = -0.0121361 Jy
! Component: 1800  -  total flux cleaned = -0.0120611 Jy
! Component: 1850  -  total flux cleaned = -0.0120243 Jy
! Component: 1900  -  total flux cleaned = -0.011769 Jy
! Component: 1950  -  total flux cleaned = -0.011661 Jy
! Component: 2000  -  total flux cleaned = -0.0115541 Jy
! Component: 2050  -  total flux cleaned = -0.0113784 Jy
! Component: 2100  -  total flux cleaned = -0.0112047 Jy
! Component: 2150  -  total flux cleaned = -0.0110674 Jy
! Component: 2200  -  total flux cleaned = -0.0108301 Jy
! Component: 2250  -  total flux cleaned = -0.0106962 Jy
! Component: 2300  -  total flux cleaned = -0.0105966 Jy
! Component: 2350  -  total flux cleaned = -0.0103997 Jy
! Component: 2400  -  total flux cleaned = -0.0102698 Jy
! Component: 2450  -  total flux cleaned = -0.0100775 Jy
! Component: 2500  -  total flux cleaned = -0.00998261 Jy
! Component: 2550  -  total flux cleaned = -0.00988828 Jy
! Component: 2600  -  total flux cleaned = -0.00976349 Jy
! Component: 2650  -  total flux cleaned = -0.00954773 Jy
! Component: 2700  -  total flux cleaned = -0.00945584 Jy
! Component: 2750  -  total flux cleaned = -0.00930475 Jy
! Component: 2800  -  total flux cleaned = -0.00930509 Jy
! Component: 2850  -  total flux cleaned = -0.00927542 Jy
! Component: 2900  -  total flux cleaned = -0.00924608 Jy
! Component: 2950  -  total flux cleaned = -0.0092459 Jy
! Component: 3000  -  total flux cleaned = -0.00915878 Jy
! Component: 3050  -  total flux cleaned = -0.00907265 Jy
! Component: 3100  -  total flux cleaned = -0.00907264 Jy
! Component: 3150  -  total flux cleaned = -0.00907256 Jy
! Component: 3200  -  total flux cleaned = -0.00904436 Jy
! Component: 3250  -  total flux cleaned = -0.00890488 Jy
! Component: 3300  -  total flux cleaned = -0.0088776 Jy
! Component: 3350  -  total flux cleaned = -0.00887764 Jy
! Component: 3400  -  total flux cleaned = -0.0088775 Jy
! Component: 3450  -  total flux cleaned = -0.00876909 Jy
! Component: 3500  -  total flux cleaned = -0.0087693 Jy
! Component: 3550  -  total flux cleaned = -0.00868905 Jy
! Component: 3600  -  total flux cleaned = -0.00863575 Jy
! Component: 3650  -  total flux cleaned = -0.00852993 Jy
! Component: 3700  -  total flux cleaned = -0.00845105 Jy
! Component: 3750  -  total flux cleaned = -0.00845108 Jy
! Component: 3800  -  total flux cleaned = -0.00832081 Jy
! Component: 3850  -  total flux cleaned = -0.00821728 Jy
! Component: 3900  -  total flux cleaned = -0.00819148 Jy
! Component: 3950  -  total flux cleaned = -0.008089 Jy
! Component: 4000  -  total flux cleaned = -0.0081401 Jy
! Component: 4050  -  total flux cleaned = -0.0080891 Jy
! Total flux subtracted in 4096 components = -0.00793743 Jy
! Clean residual min=-0.001346 max=0.001407 Jy/beam
! Clean residual mean=0.000001 rms=0.000479 Jy/beam
! Combined flux in latest and established models = 1.6519 Jy
keep
! Adding 276 model components to the UV plane model.
! The established model now contains 440 components and 1.6519 Jy
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
wmap   /vlbi/uf001/uva//J0423-0120_S_b_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.566 x 14.73 at -16.91 degrees (North through East)
! Clean map  min=-0.0083825  max=1.4253 Jy/beam
! Writing clean map to FITS file: /vlbi/uf001/uva//J0423-0120_S_b_map.fits
wmodel /vlbi/uf001/uva//J0423-0120_S_b_map.mod
! Writing 440 model components to file: /vlbi/uf001/uva//J0423-0120_S_b_map.mod
wobs   /vlbi/uf001/uva//J0423-0120_S_b_uvs.fits
! Writing UV FITS file: /vlbi/uf001/uva//J0423-0120_S_b_uvs.fits
wwins  /vlbi/uf001/uva//J0423-0120_S_b_map.win
! wwins: Wrote 3 windows to /vlbi/uf001/uva//J0423-0120_S_b_map.win
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 1.4253 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00046704 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.0983034 0.0983034 0.196607 0.393214 0.786427 1.57285 3.14571 6.29142 12.5828 25.1657 50.3313
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0014 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       1.425 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  3051.77 
quit
! Quitting program
! Log file /vlbi/uf001/uva//J0423-0120_S_b_dfm.log closed on Sun Jan 13 22:10:05 2019