! Started logfile: /image_prod/ug002_uva//J0706-4122_S_b_dfm.log on Wed Mar  4 21:39:17 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//J0706-4122_S_b_uva.fits
! Reading UV FITS file: /image_prod/ug002_uva//J0706-4122_S_b_uva.fits
! AN table 1: 50 integrations on 6 of 6 possible baselines.
! AN table 2: 50 integrations on 10 of 10 possible baselines.
! Apparent sampling: 0.995 visibilities/baseline/integration-bin.
! Found source: J0706-4122
! 
! 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 2388 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 J0706-4122_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 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=1.018614Jy  sigma=6.282181
! Fit after  self-cal, rms=0.904359Jy  sigma=5.538746
! 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=8.751 mas, bmaj=70.52 mas, bpa=5.053 degrees
! Estimated noise=3.90785 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.074864 Jy
! Component: 100  -  total flux cleaned = 0.0948201 Jy
! Total flux subtracted in 100 components = 0.0948201 Jy
! Clean residual min=-0.009804 max=0.009921 Jy/beam
! Clean residual mean=-0.000047 rms=0.003692 Jy/beam
! Combined flux in latest and established models = 0.0948201 Jy
! Performing phase self-cal
! Adding 28 model components to the UV plane model.
! The established model now contains 28 components and 0.0948201 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.087798Jy  sigma=0.480316
! Fit after  self-cal, rms=0.087378Jy  sigma=0.477427
! 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=10.54 mas, bmaj=71.87 mas, bpa=3.09 degrees
! Estimated noise=2.39827 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.0114113 Jy
! Component: 100  -  total flux cleaned = 0.0122665 Jy
! Component: 150  -  total flux cleaned = 0.012752 Jy
! Component: 200  -  total flux cleaned = 0.0126902 Jy
! Total flux subtracted in 200 components = 0.0126902 Jy
! Clean residual min=-0.002640 max=0.002664 Jy/beam
! Clean residual mean=-0.000035 rms=0.001252 Jy/beam
! Combined flux in latest and established models = 0.10751 Jy
! Performing phase self-cal
! Adding 112 model components to the UV plane model.
! The established model now contains 139 components and 0.10751 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.084665Jy  sigma=0.455871
! Fit after  self-cal, rms=0.084741Jy  sigma=0.454752
! Inverting map 
print "============== Finished natural weighting clean =============="
! ============== Finished natural weighting clean ============== 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 10.54 x 71.87 at 3.09 degrees (North through East)
! Clean map  min=-0.0077461  max=0.1005 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   KP        0.98     LA        0.92     OV        1.05     PT        1.03     
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        1.00     KP        0.97     LA        1.00     NL        1.05     
!   OV        1.03     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   KP        0.95     LA        1.00     OV        1.03     PT        0.99     
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.99     KP        0.94     LA        1.05     NL        0.94     
!   OV        1.09     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   KP        0.93     LA        1.06     OV        1.03     PT        0.89     
! 
!  Telescope amplitude corrections in sub-array 2:
!   FD        0.97     KP        1.01     LA        0.93     NL        0.91     
!   OV        0.95     
! 
! 
! Fit before self-cal, rms=0.084741Jy  sigma=0.454752
! Fit after  self-cal, rms=0.081798Jy  sigma=0.452605
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=10.56 mas, bmaj=72.1 mas, bpa=3.06 degrees
! Estimated noise=2.33142 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.00290256 Jy
! Component: 100  -  total flux cleaned = 0.00404569 Jy
! Component: 150  -  total flux cleaned = 0.00467271 Jy
! Component: 200  -  total flux cleaned = 0.00441389 Jy
! Total flux subtracted in 200 components = 0.00441389 Jy
! Clean residual min=-0.001426 max=0.001412 Jy/beam
! Clean residual mean=-0.000018 rms=0.000580 Jy/beam
! Combined flux in latest and established models = 0.111924 Jy
! Performing phase self-cal
! Adding 94 model components to the UV plane model.
! The established model now contains 233 components and 0.111924 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.081452Jy  sigma=0.450455
! Fit after  self-cal, rms=0.081496Jy  sigma=0.450139
! 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.081496Jy  sigma=0.450139
! Fit after  self-cal, rms=0.080927Jy  sigma=0.450048
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.55 mas, bmaj=72.11 mas, bpa=3.075 degrees
! Estimated noise=2.32017 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00166581 Jy
! Component: 100  -  total flux cleaned = 0.0013725 Jy
! Component: 150  -  total flux cleaned = 0.000923729 Jy
! Component: 200  -  total flux cleaned = 0.000404303 Jy
! Total flux subtracted in 200 components = 0.000404303 Jy
! Clean residual min=-0.000935 max=0.000933 Jy/beam
! Clean residual mean=-0.000017 rms=0.000393 Jy/beam
! Combined flux in latest and established models = 0.112329 Jy
selfcal
! Performing phase self-cal
! Adding 101 model components to the UV plane model.
! The established model now contains 330 components and 0.112328 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080784Jy  sigma=0.449270
! Fit after  self-cal, rms=0.080789Jy  sigma=0.449112
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.000555701 Jy
! Component: 100  -  total flux cleaned = 0.000442783 Jy
! Component: 150  -  total flux cleaned = 0.00016208 Jy
! Component: 200  -  total flux cleaned = -4.8421e-05 Jy
! Total flux subtracted in 200 components = -4.8421e-05 Jy
! Clean residual min=-0.000854 max=0.000853 Jy/beam
! Clean residual mean=-0.000014 rms=0.000356 Jy/beam
! Combined flux in latest and established models = 0.11228 Jy
! Performing phase self-cal
! Adding 95 model components to the UV plane model.
! The established model now contains 416 components and 0.11228 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080697Jy  sigma=0.448663
! Fit after  self-cal, rms=0.080719Jy  sigma=0.448608
! 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.080719Jy  sigma=0.448608
! Fit after  self-cal, rms=0.080584Jy  sigma=0.448553
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=10.54 mas, bmaj=72.03 mas, bpa=3.106 degrees
! Estimated noise=2.31452 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00101972 Jy
! Component: 100  -  total flux cleaned = 0.000918074 Jy
! Component: 150  -  total flux cleaned = 0.00077213 Jy
! Component: 200  -  total flux cleaned = 0.000502587 Jy
! Total flux subtracted in 200 components = 0.000502587 Jy
! Clean residual min=-0.000734 max=0.000737 Jy/beam
! Clean residual mean=-0.000007 rms=0.000308 Jy/beam
! Combined flux in latest and established models = 0.112783 Jy
selfcal
! Performing phase self-cal
! Adding 91 model components to the UV plane model.
! The established model now contains 488 components and 0.112783 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080512Jy  sigma=0.448207
! Fit after  self-cal, rms=0.080543Jy  sigma=0.448167
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.000487668 Jy
! Component: 100  -  total flux cleaned = 0.000490546 Jy
! Component: 150  -  total flux cleaned = 0.000358349 Jy
! Component: 200  -  total flux cleaned = 0.000231461 Jy
! Total flux subtracted in 200 components = 0.000231461 Jy
! Clean residual min=-0.000692 max=0.000690 Jy/beam
! Clean residual mean=-0.000006 rms=0.000289 Jy/beam
! Combined flux in latest and established models = 0.113014 Jy
! Performing phase self-cal
! Adding 96 model components to the UV plane model.
! The established model now contains 557 components and 0.113014 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080491Jy  sigma=0.447914
! Fit after  self-cal, rms=0.080512Jy  sigma=0.447892
! 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.080512Jy  sigma=0.447892
! Fit after  self-cal, rms=0.080347Jy  sigma=0.447867
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=11.64 mas, bmaj=71.8 mas, bpa=2.477 degrees
! Estimated noise=2.29382 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000436565 Jy
! Component: 100  -  total flux cleaned = 0.000471882 Jy
! Component: 150  -  total flux cleaned = 0.000407921 Jy
! Component: 200  -  total flux cleaned = 0.000377124 Jy
! Total flux subtracted in 200 components = 0.000377124 Jy
! Clean residual min=-0.000518 max=0.000516 Jy/beam
! Clean residual mean=-0.000002 rms=0.000218 Jy/beam
! Combined flux in latest and established models = 0.113391 Jy
selfcal
! Performing phase self-cal
! Adding 80 model components to the UV plane model.
! The established model now contains 600 components and 0.113391 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080308Jy  sigma=0.447685
! Fit after  self-cal, rms=0.080325Jy  sigma=0.447671
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.000276372 Jy
! Component: 100  -  total flux cleaned = 0.000276804 Jy
! Component: 150  -  total flux cleaned = 0.000188575 Jy
! Component: 200  -  total flux cleaned = 0.000160052 Jy
! Total flux subtracted in 200 components = 0.000160052 Jy
! Clean residual min=-0.000479 max=0.000479 Jy/beam
! Clean residual mean=-0.000002 rms=0.000201 Jy/beam
! Combined flux in latest and established models = 0.113551 Jy
! Performing phase self-cal
! Adding 85 model components to the UV plane model.
! The established model now contains 637 components and 0.113551 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080295Jy  sigma=0.447526
! Fit after  self-cal, rms=0.080310Jy  sigma=0.447517
! 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=10.53 mas, bmaj=71.98 mas, bpa=3.12 degrees
! Estimated noise=2.30894 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 10.53 x 71.98 at 3.12 degrees (North through East)
! Clean map  min=-0.0069428  max=0.098484 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=8.657 mas, bmaj=70.99 mas, bpa=5.029 degrees
! Estimated noise=3.75545 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.068539 Jy
! Component: 100  -  total flux cleaned = 0.0862548 Jy
! Component: 150  -  total flux cleaned = 0.0893201 Jy
! Component: 200  -  total flux cleaned = 0.0894251 Jy
! Total flux subtracted in 200 components = 0.0894251 Jy
! Clean residual min=-0.005501 max=0.005574 Jy/beam
! Clean residual mean=-0.000007 rms=0.002365 Jy/beam
! Combined flux in latest and established models = 0.0894251 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 102 model components to the UV plane model.
! The established model now contains 102 components and 0.0894251 Jy
! Inverting map and beam 
! Estimated beam: bmin=10.53 mas, bmaj=71.98 mas, bpa=3.12 degrees
! Estimated noise=2.30894 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.00282015 Jy
! Component: 100  -  total flux cleaned = 0.00253546 Jy
! Component: 150  -  total flux cleaned = 0.00200996 Jy
! Component: 200  -  total flux cleaned = 0.00188457 Jy
! Total flux subtracted in 200 components = 0.00188457 Jy
! Clean residual min=-0.002089 max=0.002085 Jy/beam
! Clean residual mean=0.000009 rms=0.000932 Jy/beam
! Combined flux in latest and established models = 0.0913097 Jy
! Adding 102 model components to the UV plane model.
! The established model now contains 204 components and 0.0913096 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=11.64 mas, bmaj=71.8 mas, bpa=2.477 degrees
! Estimated noise=2.29382 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000465508 Jy
! Component: 100  -  total flux cleaned = 6.83128e-05 Jy
! Component: 150  -  total flux cleaned = 0.000770201 Jy
! Component: 200  -  total flux cleaned = 0.000767176 Jy
! Total flux subtracted in 200 components = 0.000767176 Jy
! Clean residual min=-0.001121 max=0.001108 Jy/beam
! Clean residual mean=0.000007 rms=0.000495 Jy/beam
! Combined flux in latest and established models = 0.0920768 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 = -6.87025e-05 Jy
! Component: 100  -  total flux cleaned = -0.000131775 Jy
! Component: 150  -  total flux cleaned = -0.000413451 Jy
! Component: 200  -  total flux cleaned = -0.00079491 Jy
! Total flux subtracted in 200 components = -0.00079491 Jy
! Clean residual min=-0.000892 max=0.000894 Jy/beam
! Clean residual mean=0.000008 rms=0.000368 Jy/beam
! Combined flux in latest and established models = 0.0912819 Jy
! Adding 157 model components to the UV plane model.
! The established model now contains 351 components and 0.0912819 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=10.53 mas, bmaj=71.98 mas, bpa=3.12 degrees
! Estimated noise=2.30894 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -7.96973e-05 Jy
! Component: 100  -  total flux cleaned = -0.000396949 Jy
! Component: 150  -  total flux cleaned = -0.000649733 Jy
! Component: 200  -  total flux cleaned = -0.000929749 Jy
! Total flux subtracted in 200 components = -0.000929749 Jy
! Clean residual min=-0.000650 max=0.000647 Jy/beam
! Clean residual mean=0.000002 rms=0.000281 Jy/beam
! Combined flux in latest and established models = 0.0903522 Jy
selfcal
! Performing phase self-cal
! Adding 99 model components to the UV plane model.
! The established model now contains 450 components and 0.0903521 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.080351Jy  sigma=0.447888
! Fit after  self-cal, rms=0.080343Jy  sigma=0.447807
wmodel J0706-4122_S_b_map.mod
! Writing 450 model components to file: J0706-4122_S_b_map.mod
wobs   J0706-4122_S_b_uvs.fits
! Writing UV FITS file: J0706-4122_S_b_uvs.fits
wwins  J0706-4122_S_b_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
! 
wmap   /image_prod/ug002_uva//J0706-4122_S_b_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 10.53 x 71.98 at 3.12 degrees (North through East)
! Clean map  min=-0.0093612  max=0.09064 Jy/beam
! Writing clean map to FITS file: /image_prod/ug002_uva//J0706-4122_S_b_map.fits
wmodel /image_prod/ug002_uva//J0706-4122_S_b_map.mod
! Writing 450 model components to file: /image_prod/ug002_uva//J0706-4122_S_b_map.mod
wobs   /image_prod/ug002_uva//J0706-4122_S_b_uvs.fits
! Writing UV FITS file: /image_prod/ug002_uva//J0706-4122_S_b_uvs.fits
wwins  /image_prod/ug002_uva//J0706-4122_S_b_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.0906403 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000292838 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.96923 0.96923 1.93846 3.87692 7.75384 15.5077 31.0154 62.0307
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0008 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       0.09 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  309.524 
quit
! Quitting program
! Log file /image_prod/ug002_uva//J0706-4122_S_b_dfm.log closed on Wed Mar  4 21:39:21 2020