! Started logfile: /image_prod/ug002_uva//J1941-3007_S_dfm.log on Thu Mar  5 04:22:27 2020
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     /image_prod/ug002_uva//J1941-3007_S_uva.fits
! Reading UV FITS file: /image_prod/ug002_uva//J1941-3007_S_uva.fits
! AN table 1: 60 integrations on 45 of 45 possible baselines.
! Apparent sampling: 0.886667 visibilities/baseline/integration-bin.
! Found source: J1941-3007
! 
! 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 44 lines of history.
! 
! Reading 7182 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 J1941-3007_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 34 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 34 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 34 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.031425Jy  sigma=6.291833
! Fit after  self-cal, rms=0.874201Jy  sigma=5.146994
! 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.482 mas, bmaj=6.556 mas, bpa=-4.227 degrees
! Estimated noise=2.35861 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.0974429 Jy
! Component: 100  -  total flux cleaned = 0.145187 Jy
! Total flux subtracted in 100 components = 0.145187 Jy
! Clean residual min=-0.014529 max=0.042176 Jy/beam
! Clean residual mean=-0.000129 rms=0.004924 Jy/beam
! Combined flux in latest and established models = 0.145187 Jy
! Performing phase self-cal
! Adding 22 model components to the UV plane model.
! The established model now contains 22 components and 0.145187 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.147004Jy  sigma=0.743453
! Fit after  self-cal, rms=0.145842Jy  sigma=0.735687
! Inverting map 
! Added new window around map position (4, -0.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.040195 Jy
! Component: 100  -  total flux cleaned = 0.0630347 Jy
! Total flux subtracted in 100 components = 0.0630347 Jy
! Clean residual min=-0.011311 max=0.022597 Jy/beam
! Clean residual mean=-0.000065 rms=0.003475 Jy/beam
! Combined flux in latest and established models = 0.208222 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 37 components and 0.208222 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.132695Jy  sigma=0.593599
! Fit after  self-cal, rms=0.131304Jy  sigma=0.583668
! 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.071 mas, bmaj=7.955 mas, bpa=-3.692 degrees
! Estimated noise=1.45157 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.0273635 Jy
! Component: 100  -  total flux cleaned = 0.0420387 Jy
! Component: 150  -  total flux cleaned = 0.0506306 Jy
! Component: 200  -  total flux cleaned = 0.0556788 Jy
! Total flux subtracted in 200 components = 0.0556788 Jy
! Clean residual min=-0.009118 max=0.014552 Jy/beam
! Clean residual mean=-0.000028 rms=0.002652 Jy/beam
! Combined flux in latest and established models = 0.263901 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 45 components and 0.263901 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.127238Jy  sigma=0.523010
! Fit after  self-cal, rms=0.126804Jy  sigma=0.518457
! 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.071 x 7.955 at -3.692 degrees (North through East)
! Clean map  min=-0.0089632  max=0.11363 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        2.21     FD        0.93     HN        1.37     KP        0.95     
!   LA        0.97     MK        0.88     NL        0.98     OV        0.92     
!   PT        0.99     SC        0.79     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        2.16     FD        0.92     HN        1.21     KP        0.94     
!   LA        0.97     MK        0.91     NL        1.09     OV        0.95     
!   PT        0.99     SC        0.77     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        2.56     FD        0.92     HN        1.40     KP        0.98     
!   LA        0.98     MK        0.86     NL        0.97     OV        0.94     
!   PT        0.97     SC        0.85     
! 
! 
! Fit before self-cal, rms=0.126804Jy  sigma=0.518457
! Fit after  self-cal, rms=0.231029Jy  sigma=0.479427
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.919 mas, bmaj=7.62 mas, bpa=-2.684 degrees
! Estimated noise=1.38291 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.00335284 Jy
! Component: 100  -  total flux cleaned = 0.00321315 Jy
! Component: 150  -  total flux cleaned = 0.00289557 Jy
! Component: 200  -  total flux cleaned = 0.00262921 Jy
! Total flux subtracted in 200 components = 0.00262921 Jy
! Clean residual min=-0.005420 max=0.004702 Jy/beam
! Clean residual mean=-0.000005 rms=0.001198 Jy/beam
! Combined flux in latest and established models = 0.26653 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 63 components and 0.26653 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.230673Jy  sigma=0.477278
! Fit after  self-cal, rms=0.231065Jy  sigma=0.476282
! 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.231065Jy  sigma=0.476282
! Fit after  self-cal, rms=0.208695Jy  sigma=0.472357
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.895 mas, bmaj=7.604 mas, bpa=-1.582 degrees
! Estimated noise=1.3782 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00286196 Jy
! Component: 100  -  total flux cleaned = 0.00305385 Jy
! Component: 150  -  total flux cleaned = 0.00327347 Jy
! Component: 200  -  total flux cleaned = 0.00326994 Jy
! Total flux subtracted in 200 components = 0.00326994 Jy
! Clean residual min=-0.004195 max=0.004663 Jy/beam
! Clean residual mean=0.000000 rms=0.000961 Jy/beam
! Combined flux in latest and established models = 0.2698 Jy
selfcal
! Performing phase self-cal
! Adding 22 model components to the UV plane model.
! The established model now contains 77 components and 0.2698 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.208417Jy  sigma=0.471555
! Fit after  self-cal, rms=0.208429Jy  sigma=0.471187
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 (-19.5, 16).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00150614 Jy
! Component: 100  -  total flux cleaned = 0.00142207 Jy
! Component: 150  -  total flux cleaned = 0.00107681 Jy
! Component: 200  -  total flux cleaned = 0.000769901 Jy
! Total flux subtracted in 200 components = 0.000769901 Jy
! Clean residual min=-0.003721 max=0.003912 Jy/beam
! Clean residual mean=0.000002 rms=0.000882 Jy/beam
! Combined flux in latest and established models = 0.27057 Jy
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 90 components and 0.27057 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.208241Jy  sigma=0.469826
! Fit after  self-cal, rms=0.208237Jy  sigma=0.469694
! Inverting map 
! Added new window around map position (-35, 37.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00189093 Jy
! Component: 100  -  total flux cleaned = 0.00248001 Jy
! Component: 150  -  total flux cleaned = 0.00273391 Jy
! Component: 200  -  total flux cleaned = 0.00273627 Jy
! Total flux subtracted in 200 components = 0.00273627 Jy
! Clean residual min=-0.003323 max=0.002947 Jy/beam
! Clean residual mean=0.000008 rms=0.000812 Jy/beam
! Combined flux in latest and established models = 0.273306 Jy
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 104 components and 0.273306 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.208071Jy  sigma=0.468462
! Fit after  self-cal, rms=0.208081Jy  sigma=0.468390
! 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.208081Jy  sigma=0.468390
! Fit after  self-cal, rms=0.183677Jy  sigma=0.465920
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.893 mas, bmaj=7.648 mas, bpa=-1.531 degrees
! Estimated noise=1.37143 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0007664 Jy
! Component: 100  -  total flux cleaned = -0.000827029 Jy
! Component: 150  -  total flux cleaned = -0.000998964 Jy
! Component: 200  -  total flux cleaned = -0.00110688 Jy
! Total flux subtracted in 200 components = -0.00110688 Jy
! Clean residual min=-0.003105 max=0.002718 Jy/beam
! Clean residual mean=-0.000000 rms=0.000674 Jy/beam
! Combined flux in latest and established models = 0.272199 Jy
selfcal
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 115 components and 0.272199 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.183600Jy  sigma=0.465477
! Fit after  self-cal, rms=0.183582Jy  sigma=0.465441
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.97796e-05 Jy
! Component: 100  -  total flux cleaned = -0.000250938 Jy
! Component: 150  -  total flux cleaned = -0.000346706 Jy
! Component: 200  -  total flux cleaned = -0.000347908 Jy
! Total flux subtracted in 200 components = -0.000347908 Jy
! Clean residual min=-0.003008 max=0.002628 Jy/beam
! Clean residual mean=-0.000001 rms=0.000659 Jy/beam
! Combined flux in latest and established models = 0.271851 Jy
! Performing phase self-cal
! Adding 28 model components to the UV plane model.
! The established model now contains 124 components and 0.271851 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.183522Jy  sigma=0.465230
! Fit after  self-cal, rms=0.183527Jy  sigma=0.465215
! 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.183527Jy  sigma=0.465215
! Fit after  self-cal, rms=0.183470Jy  sigma=0.465192
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.09 mas, bmaj=20.89 mas, bpa=1.777 degrees
! Estimated noise=1.86928 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 8.42225e-05 Jy
! Component: 100  -  total flux cleaned = 0.000113587 Jy
! Component: 150  -  total flux cleaned = 5.72658e-05 Jy
! Component: 200  -  total flux cleaned = 3.42452e-06 Jy
! Total flux subtracted in 200 components = 3.42452e-06 Jy
! Clean residual min=-0.002235 max=0.002430 Jy/beam
! Clean residual mean=-0.000007 rms=0.000693 Jy/beam
! Combined flux in latest and established models = 0.271855 Jy
selfcal
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 135 components and 0.271855 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.183454Jy  sigma=0.465172
! Fit after  self-cal, rms=0.183458Jy  sigma=0.465168
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 = 8.32101e-05 Jy
! Component: 100  -  total flux cleaned = 5.83353e-05 Jy
! Component: 150  -  total flux cleaned = -1.28319e-05 Jy
! Component: 200  -  total flux cleaned = -3.62191e-05 Jy
! Total flux subtracted in 200 components = -3.62191e-05 Jy
! Clean residual min=-0.002186 max=0.002412 Jy/beam
! Clean residual mean=-0.000007 rms=0.000688 Jy/beam
! Combined flux in latest and established models = 0.271819 Jy
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 137 components and 0.271819 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.183446Jy  sigma=0.465158
! Fit after  self-cal, rms=0.183446Jy  sigma=0.465155
! 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=2.892 mas, bmaj=7.669 mas, bpa=-1.669 degrees
! Estimated noise=1.37026 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.892 x 7.669 at -1.669 degrees (North through East)
! Clean map  min=-0.0028126  max=0.10722 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=2.359 mas, bmaj=6.405 mas, bpa=-2.425 degrees
! Estimated noise=2.15095 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.097086 Jy
! Component: 100  -  total flux cleaned = 0.150492 Jy
! Component: 150  -  total flux cleaned = 0.18541 Jy
! Component: 200  -  total flux cleaned = 0.209126 Jy
! Total flux subtracted in 200 components = 0.209126 Jy
! Clean residual min=-0.006455 max=0.014604 Jy/beam
! Clean residual mean=-0.000057 rms=0.001997 Jy/beam
! Combined flux in latest and established models = 0.209126 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 45 model components to the UV plane model.
! The established model now contains 45 components and 0.209126 Jy
! Inverting map and beam 
! Estimated beam: bmin=2.892 mas, bmaj=7.669 mas, bpa=-1.669 degrees
! Estimated noise=1.37026 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.0247483 Jy
! Component: 100  -  total flux cleaned = 0.0381922 Jy
! Component: 150  -  total flux cleaned = 0.0459887 Jy
! Component: 200  -  total flux cleaned = 0.0499802 Jy
! Total flux subtracted in 200 components = 0.0499802 Jy
! Clean residual min=-0.003324 max=0.005095 Jy/beam
! Clean residual mean=-0.000024 rms=0.000881 Jy/beam
! Combined flux in latest and established models = 0.259106 Jy
! Adding 62 model components to the UV plane model.
! The established model now contains 84 components and 0.259106 Jy
! Inverting map 
! Added new window around map position (-4, -4).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00437207 Jy
! Component: 100  -  total flux cleaned = 0.00657807 Jy
! Component: 150  -  total flux cleaned = 0.00785008 Jy
! Component: 200  -  total flux cleaned = 0.00848081 Jy
! Total flux subtracted in 200 components = 0.00848081 Jy
! Clean residual min=-0.002699 max=0.002507 Jy/beam
! Clean residual mean=-0.000009 rms=0.000666 Jy/beam
! Combined flux in latest and established models = 0.267587 Jy
! Adding 29 model components to the UV plane model.
! The established model now contains 109 components and 0.267587 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=8.09 mas, bmaj=20.89 mas, bpa=1.777 degrees
! Estimated noise=1.86928 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00175491 Jy
! Component: 100  -  total flux cleaned = 0.00249613 Jy
! Component: 150  -  total flux cleaned = 0.00314212 Jy
! Component: 200  -  total flux cleaned = 0.00359453 Jy
! Total flux subtracted in 200 components = 0.00359453 Jy
! Clean residual min=-0.002342 max=0.002361 Jy/beam
! Clean residual mean=-0.000009 rms=0.000673 Jy/beam
! Combined flux in latest and established models = 0.271181 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.000214728 Jy
! Component: 100  -  total flux cleaned = 0.00040772 Jy
! Component: 150  -  total flux cleaned = 0.000496197 Jy
! Component: 200  -  total flux cleaned = 0.000607215 Jy
! Total flux subtracted in 200 components = 0.000607215 Jy
! Clean residual min=-0.002311 max=0.002426 Jy/beam
! Clean residual mean=-0.000006 rms=0.000668 Jy/beam
! Combined flux in latest and established models = 0.271789 Jy
! Adding 39 model components to the UV plane model.
! The established model now contains 146 components and 0.271789 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=2.892 mas, bmaj=7.669 mas, bpa=-1.669 degrees
! Estimated noise=1.37026 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000148883 Jy
! Component: 100  -  total flux cleaned = -0.00021745 Jy
! Component: 150  -  total flux cleaned = -0.000218326 Jy
! Component: 200  -  total flux cleaned = -0.000112668 Jy
! Total flux subtracted in 200 components = -0.000112668 Jy
! Clean residual min=-0.002589 max=0.002534 Jy/beam
! Clean residual mean=-0.000001 rms=0.000606 Jy/beam
! Combined flux in latest and established models = 0.271676 Jy
selfcal
! Performing phase self-cal
! Adding 24 model components to the UV plane model.
! The established model now contains 168 components and 0.271676 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.183341Jy  sigma=0.464537
! Fit after  self-cal, rms=0.183357Jy  sigma=0.464260
wmodel J1941-3007_S_map.mod
! Writing 168 model components to file: J1941-3007_S_map.mod
wobs   J1941-3007_S_uvs.fits
! Writing UV FITS file: J1941-3007_S_uvs.fits
wwins  J1941-3007_S_map.win
! wwins: Wrote 5 windows to J1941-3007_S_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 = 9.77413e-05 Jy
! Component: 100  -  total flux cleaned = 0.000101882 Jy
! Component: 150  -  total flux cleaned = -4.94e-05 Jy
! Component: 200  -  total flux cleaned = -0.000122316 Jy
! Component: 250  -  total flux cleaned = -0.000297095 Jy
! Component: 300  -  total flux cleaned = -0.00050008 Jy
! Component: 350  -  total flux cleaned = -0.000729763 Jy
! Component: 400  -  total flux cleaned = -0.000858448 Jy
! Component: 450  -  total flux cleaned = -0.00107689 Jy
! Component: 500  -  total flux cleaned = -0.00119843 Jy
! Component: 550  -  total flux cleaned = -0.00131747 Jy
! Component: 600  -  total flux cleaned = -0.00137621 Jy
! Component: 650  -  total flux cleaned = -0.00146189 Jy
! Component: 700  -  total flux cleaned = -0.00146172 Jy
! Component: 750  -  total flux cleaned = -0.00137934 Jy
! Component: 800  -  total flux cleaned = -0.00140645 Jy
! Component: 850  -  total flux cleaned = -0.00135287 Jy
! Component: 900  -  total flux cleaned = -0.00135258 Jy
! Component: 950  -  total flux cleaned = -0.00140415 Jy
! Component: 1000  -  total flux cleaned = -0.00132767 Jy
! Component: 1050  -  total flux cleaned = -0.0013273 Jy
! Component: 1100  -  total flux cleaned = -0.00125297 Jy
! Component: 1150  -  total flux cleaned = -0.00117965 Jy
! Component: 1200  -  total flux cleaned = -0.00115549 Jy
! Component: 1250  -  total flux cleaned = -0.00110749 Jy
! Component: 1300  -  total flux cleaned = -0.0010836 Jy
! Component: 1350  -  total flux cleaned = -0.00101309 Jy
! Component: 1400  -  total flux cleaned = -0.000943365 Jy
! Component: 1450  -  total flux cleaned = -0.00092033 Jy
! Component: 1500  -  total flux cleaned = -0.00092049 Jy
! Component: 1550  -  total flux cleaned = -0.000965908 Jy
! Component: 1600  -  total flux cleaned = -0.000920876 Jy
! Component: 1650  -  total flux cleaned = -0.000876398 Jy
! Component: 1700  -  total flux cleaned = -0.000942311 Jy
! Component: 1750  -  total flux cleaned = -0.000876676 Jy
! Component: 1800  -  total flux cleaned = -0.000833134 Jy
! Component: 1850  -  total flux cleaned = -0.000790118 Jy
! Component: 1900  -  total flux cleaned = -0.00074745 Jy
! Component: 1950  -  total flux cleaned = -0.000768592 Jy
! Component: 2000  -  total flux cleaned = -0.000684351 Jy
! Component: 2050  -  total flux cleaned = -0.000663526 Jy
! Component: 2100  -  total flux cleaned = -0.000684102 Jy
! Component: 2150  -  total flux cleaned = -0.000622268 Jy
! Component: 2200  -  total flux cleaned = -0.000622346 Jy
! Component: 2250  -  total flux cleaned = -0.000520841 Jy
! Component: 2300  -  total flux cleaned = -0.000521022 Jy
! Component: 2350  -  total flux cleaned = -0.000460892 Jy
! Component: 2400  -  total flux cleaned = -0.000440851 Jy
! Component: 2450  -  total flux cleaned = -0.000420951 Jy
! Component: 2500  -  total flux cleaned = -0.00036189 Jy
! Component: 2550  -  total flux cleaned = -0.000401096 Jy
! Component: 2600  -  total flux cleaned = -0.000381779 Jy
! Component: 2650  -  total flux cleaned = -0.00036256 Jy
! Component: 2700  -  total flux cleaned = -0.000343356 Jy
! Component: 2750  -  total flux cleaned = -0.000343235 Jy
! Component: 2800  -  total flux cleaned = -0.00041911 Jy
! Component: 2850  -  total flux cleaned = -0.000381293 Jy
! Component: 2900  -  total flux cleaned = -0.000437551 Jy
! Component: 2950  -  total flux cleaned = -0.000418763 Jy
! Component: 3000  -  total flux cleaned = -0.0004931 Jy
! Component: 3050  -  total flux cleaned = -0.000400551 Jy
! Component: 3100  -  total flux cleaned = -0.000474406 Jy
! Component: 3150  -  total flux cleaned = -0.00051102 Jy
! Component: 3200  -  total flux cleaned = -0.000474623 Jy
! Component: 3250  -  total flux cleaned = -0.000583653 Jy
! Component: 3300  -  total flux cleaned = -0.000637832 Jy
! Component: 3350  -  total flux cleaned = -0.000637791 Jy
! Component: 3400  -  total flux cleaned = -0.000673555 Jy
! Component: 3450  -  total flux cleaned = -0.000691417 Jy
! Component: 3500  -  total flux cleaned = -0.000726848 Jy
! Component: 3550  -  total flux cleaned = -0.000762097 Jy
! Component: 3600  -  total flux cleaned = -0.00077966 Jy
! Component: 3650  -  total flux cleaned = -0.00079713 Jy
! Component: 3700  -  total flux cleaned = -0.0008319 Jy
! Component: 3750  -  total flux cleaned = -0.000849231 Jy
! Component: 3800  -  total flux cleaned = -0.000900694 Jy
! Component: 3850  -  total flux cleaned = -0.000934943 Jy
! Component: 3900  -  total flux cleaned = -0.000952057 Jy
! Component: 3950  -  total flux cleaned = -0.00100299 Jy
! Component: 4000  -  total flux cleaned = -0.00107065 Jy
! Component: 4050  -  total flux cleaned = -0.00103703 Jy
! Total flux subtracted in 4096 components = -0.00110406 Jy
! Clean residual min=-0.000904 max=0.001037 Jy/beam
! Clean residual mean=0.000001 rms=0.000307 Jy/beam
! Combined flux in latest and established models = 0.270572 Jy
keep
! Adding 260 model components to the UV plane model.
! The established model now contains 428 components and 0.270572 Jy
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
wmap   /image_prod/ug002_uva//J1941-3007_S_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.892 x 7.669 at -1.669 degrees (North through East)
! Clean map  min=-0.0029891  max=0.11673 Jy/beam
! Writing clean map to FITS file: /image_prod/ug002_uva//J1941-3007_S_map.fits
wmodel /image_prod/ug002_uva//J1941-3007_S_map.mod
! Writing 428 model components to file: /image_prod/ug002_uva//J1941-3007_S_map.mod
wobs   /image_prod/ug002_uva//J1941-3007_S_uvs.fits
! Writing UV FITS file: /image_prod/ug002_uva//J1941-3007_S_uvs.fits
wwins  /image_prod/ug002_uva//J1941-3007_S_map.win
! wwins: Wrote 6 windows to /image_prod/ug002_uva//J1941-3007_S_map.win
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 0.116725 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000301693 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.775395 0.775395 1.55079 3.10158 6.20316 12.4063 24.8126 49.6253 99.2506
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0009 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       0.116 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  386.899 
quit
! Quitting program
! Log file /image_prod/ug002_uva//J1941-3007_S_dfm.log closed on Thu Mar  5 04:22:31 2020