! Started logfile: /image_prod/ug002_uva//J0854+2006_S_a_dfm.log on Wed Mar  4 22:38:11 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//J0854+2006_S_a_uva.fits
! Reading UV FITS file: /image_prod/ug002_uva//J0854+2006_S_a_uva.fits
! AN table 1: 11 integrations on 36 of 36 possible baselines.
! Apparent sampling: 0.911616 visibilities/baseline/integration-bin.
! Found source: J0854+2006
! 
! 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 1083 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 J0854+2006_S_a]
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 9 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 9 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 9 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=3.952481Jy  sigma=8.844193
! Fit after  self-cal, rms=2.876024Jy  sigma=6.129562
! 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.847 mas, bmaj=10.87 mas, bpa=-18.73 degrees
! Estimated noise=10.6569 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 = 2.69888 Jy
! Component: 100  -  total flux cleaned = 3.40077 Jy
! Total flux subtracted in 100 components = 3.40077 Jy
! Clean residual min=-0.138414 max=0.301793 Jy/beam
! Clean residual mean=-0.000239 rms=0.042184 Jy/beam
! Combined flux in latest and established models = 3.40077 Jy
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 7 components and 3.40077 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.922935Jy  sigma=1.473493
! Fit after  self-cal, rms=0.922096Jy  sigma=1.471136
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.32183 Jy
! Component: 100  -  total flux cleaned = 0.489468 Jy
! Total flux subtracted in 100 components = 0.489468 Jy
! Clean residual min=-0.128727 max=0.127452 Jy/beam
! Clean residual mean=-0.000239 rms=0.029278 Jy/beam
! Combined flux in latest and established models = 3.89024 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 17 components and 3.89024 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.735984Jy  sigma=1.129485
! Fit after  self-cal, rms=0.735175Jy  sigma=1.126885
! 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.313 mas, bmaj=11.06 mas, bpa=-18.67 degrees
! Estimated noise=9.74826 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.111384 Jy
! Component: 100  -  total flux cleaned = 0.116878 Jy
! Component: 150  -  total flux cleaned = 0.118983 Jy
! Component: 200  -  total flux cleaned = 0.118981 Jy
! Total flux subtracted in 200 components = 0.118981 Jy
! Clean residual min=-0.137408 max=0.101979 Jy/beam
! Clean residual mean=-0.000317 rms=0.028667 Jy/beam
! Combined flux in latest and established models = 4.00922 Jy
! Performing phase self-cal
! Adding 28 model components to the UV plane model.
! The established model now contains 39 components and 4.00922 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.692836Jy  sigma=1.095390
! Fit after  self-cal, rms=0.692340Jy  sigma=1.093528
! 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.313 x 11.06 at -18.67 degrees (North through East)
! Clean map  min=-0.13407  max=3.5259 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.04     FD        0.98     HN        1.02     KP        1.03     
!   LA        1.09     MK        1.02     NL        1.08     OV        1.06     
!   PT        0.98     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        1.00     HN        1.01     KP        1.06     
!   LA        1.08     MK        1.00     NL        0.89     OV        0.89     
!   PT        0.94     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.80     FD        1.00     HN        1.00     KP        1.03     
!   LA        1.08     MK        1.00     NL        0.99     OV        1.05     
!   PT        0.92     
! 
! 
! Fit before self-cal, rms=0.692340Jy  sigma=1.093528
! Fit after  self-cal, rms=0.440595Jy  sigma=0.817532
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=3.321 mas, bmaj=10.95 mas, bpa=-18.59 degrees
! Estimated noise=9.78377 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.0395088 Jy
! Component: 100  -  total flux cleaned = 0.0536431 Jy
! Component: 150  -  total flux cleaned = 0.0557979 Jy
! Component: 200  -  total flux cleaned = 0.0576021 Jy
! Total flux subtracted in 200 components = 0.0576021 Jy
! Clean residual min=-0.097108 max=0.097443 Jy/beam
! Clean residual mean=-0.000327 rms=0.022690 Jy/beam
! Combined flux in latest and established models = 4.06682 Jy
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 55 components and 4.06682 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.432271Jy  sigma=0.814979
! Fit after  self-cal, rms=0.432042Jy  sigma=0.814146
! 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.432042Jy  sigma=0.814146
! Fit after  self-cal, rms=0.095836Jy  sigma=0.200338
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.325 mas, bmaj=10.91 mas, bpa=-18.55 degrees
! Estimated noise=9.77303 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00502552 Jy
! Component: 100  -  total flux cleaned = 0.00715027 Jy
! Component: 150  -  total flux cleaned = 0.00799197 Jy
! Component: 200  -  total flux cleaned = 0.00905657 Jy
! Total flux subtracted in 200 components = 0.00905657 Jy
! Clean residual min=-0.016135 max=0.028048 Jy/beam
! Clean residual mean=0.000000 rms=0.004502 Jy/beam
! Combined flux in latest and established models = 4.07588 Jy
selfcal
! Performing phase self-cal
! Adding 26 model components to the UV plane model.
! The established model now contains 75 components and 4.07588 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.090407Jy  sigma=0.188143
! Fit after  self-cal, rms=0.089947Jy  sigma=0.186662
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 (-9, -1).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0207288 Jy
! Component: 100  -  total flux cleaned = 0.0202969 Jy
! Component: 150  -  total flux cleaned = 0.0183625 Jy
! Component: 200  -  total flux cleaned = 0.0173184 Jy
! Total flux subtracted in 200 components = 0.0173184 Jy
! Clean residual min=-0.016233 max=0.014764 Jy/beam
! Clean residual mean=-0.000002 rms=0.003611 Jy/beam
! Combined flux in latest and established models = 4.0932 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 98 components and 4.0932 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.081678Jy  sigma=0.168834
! Fit after  self-cal, rms=0.080787Jy  sigma=0.167038
! 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.080787Jy  sigma=0.167038
! Fit after  self-cal, rms=0.078010Jy  sigma=0.163964
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=3.325 mas, bmaj=10.9 mas, bpa=-18.53 degrees
! Estimated noise=9.78055 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00843768 Jy
! Component: 100  -  total flux cleaned = 0.0105843 Jy
! Component: 150  -  total flux cleaned = 0.0105895 Jy
! Component: 200  -  total flux cleaned = 0.0111417 Jy
! Total flux subtracted in 200 components = 0.0111417 Jy
! Clean residual min=-0.014911 max=0.013286 Jy/beam
! Clean residual mean=0.000002 rms=0.003098 Jy/beam
! Combined flux in latest and established models = 4.10434 Jy
selfcal
! Performing phase self-cal
! Adding 36 model components to the UV plane model.
! The established model now contains 119 components and 4.10434 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.074908Jy  sigma=0.157655
! Fit after  self-cal, rms=0.074658Jy  sigma=0.157115
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.00131027 Jy
! Component: 100  -  total flux cleaned = 0.000583877 Jy
! Component: 150  -  total flux cleaned = 0.000583496 Jy
! Component: 200  -  total flux cleaned = -3.86961e-05 Jy
! Total flux subtracted in 200 components = -3.86961e-05 Jy
! Clean residual min=-0.014171 max=0.012186 Jy/beam
! Clean residual mean=0.000002 rms=0.002918 Jy/beam
! Combined flux in latest and established models = 4.1043 Jy
! Performing phase self-cal
! Adding 43 model components to the UV plane model.
! The established model now contains 141 components and 4.1043 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.073117Jy  sigma=0.154101
! Fit after  self-cal, rms=0.073025Jy  sigma=0.153841
! 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.073025Jy  sigma=0.153841
! Fit after  self-cal, rms=0.072196Jy  sigma=0.152798
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.391 mas, bmaj=11.54 mas, bpa=-17.56 degrees
! Estimated noise=12.6093 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00099329 Jy
! Component: 100  -  total flux cleaned = 0.000465013 Jy
! Component: 150  -  total flux cleaned = 0.000237031 Jy
! Component: 200  -  total flux cleaned = -0.000143376 Jy
! Total flux subtracted in 200 components = -0.000143376 Jy
! Clean residual min=-0.010210 max=0.014320 Jy/beam
! Clean residual mean=0.000011 rms=0.003212 Jy/beam
! Combined flux in latest and established models = 4.10416 Jy
selfcal
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 165 components and 4.10416 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.070952Jy  sigma=0.150371
! Fit after  self-cal, rms=0.070884Jy  sigma=0.150030
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.5, -28.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0120132 Jy
! Component: 100  -  total flux cleaned = 0.0181839 Jy
! Component: 150  -  total flux cleaned = 0.0203587 Jy
! Component: 200  -  total flux cleaned = 0.0212393 Jy
! Total flux subtracted in 200 components = 0.0212393 Jy
! Clean residual min=-0.009789 max=0.010679 Jy/beam
! Clean residual mean=0.000014 rms=0.002735 Jy/beam
! Combined flux in latest and established models = 4.1254 Jy
! Performing phase self-cal
! Adding 41 model components to the UV plane model.
! The established model now contains 199 components and 4.12539 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.068385Jy  sigma=0.145799
! Fit after  self-cal, rms=0.068361Jy  sigma=0.145580
! 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.326 mas, bmaj=10.89 mas, bpa=-18.52 degrees
! Estimated noise=9.78939 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.326 x 10.89 at -18.52 degrees (North through East)
! Clean map  min=-0.014163  max=3.5333 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.068361Jy  sigma=0.145580
! Fit after  self-cal, rms=0.064156Jy  sigma=0.132651
! Inverting map and beam 
! Estimated beam: bmin=3.326 mas, bmaj=10.89 mas, bpa=-18.52 degrees
! Estimated noise=9.79103 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -4.38187e-05 Jy
! Component: 100  -  total flux cleaned = 0.000231044 Jy
! Component: 150  -  total flux cleaned = 0.000737037 Jy
! Component: 200  -  total flux cleaned = 0.00119752 Jy
! Total flux subtracted in 200 components = 0.00119752 Jy
! Clean residual min=-0.010920 max=0.008614 Jy/beam
! Clean residual mean=0.000008 rms=0.002112 Jy/beam
! Combined flux in latest and established models = 4.12659 Jy
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 226 components and 4.12659 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.062030Jy  sigma=0.128532
! Fit after  self-cal, rms=0.062035Jy  sigma=0.128415
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000709423 Jy
! Component: 100  -  total flux cleaned = 0.00132866 Jy
! Component: 150  -  total flux cleaned = 0.00190244 Jy
! Component: 200  -  total flux cleaned = 0.0024348 Jy
! Total flux subtracted in 200 components = 0.0024348 Jy
! Clean residual min=-0.009927 max=0.008025 Jy/beam
! Clean residual mean=0.000008 rms=0.001958 Jy/beam
! Combined flux in latest and established models = 4.12903 Jy
! Performing phase self-cal
! Adding 36 model components to the UV plane model.
! The established model now contains 240 components and 4.12902 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.060788Jy  sigma=0.126074
! Fit after  self-cal, rms=0.060803Jy  sigma=0.126009
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.060803Jy  sigma=0.126009
! Fit after  self-cal, rms=0.060517Jy  sigma=0.125724
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.861 mas, bmaj=10.68 mas, bpa=-18.61 degrees
! Estimated noise=10.7023 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 2.70867 Jy
! Component: 100  -  total flux cleaned = 3.40965 Jy
! Component: 150  -  total flux cleaned = 3.72879 Jy
! Component: 200  -  total flux cleaned = 3.89744 Jy
! Total flux subtracted in 200 components = 3.89744 Jy
! Clean residual min=-0.030642 max=0.082943 Jy/beam
! Clean residual mean=-0.000017 rms=0.010788 Jy/beam
! Combined flux in latest and established models = 3.89744 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 17 model components to the UV plane model.
! The established model now contains 17 components and 3.89744 Jy
! Inverting map and beam 
! Estimated beam: bmin=3.326 mas, bmaj=10.89 mas, bpa=-18.51 degrees
! Estimated noise=9.79193 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.111335 Jy
! Component: 100  -  total flux cleaned = 0.167821 Jy
! Component: 150  -  total flux cleaned = 0.200529 Jy
! Component: 200  -  total flux cleaned = 0.218368 Jy
! Total flux subtracted in 200 components = 0.218368 Jy
! Clean residual min=-0.013936 max=0.012014 Jy/beam
! Clean residual mean=0.000002 rms=0.003355 Jy/beam
! Combined flux in latest and established models = 4.1158 Jy
! Adding 37 model components to the UV plane model.
! The established model now contains 47 components and 4.1158 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.391 mas, bmaj=11.53 mas, bpa=-17.5 degrees
! Estimated noise=12.615 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0134562 Jy
! Component: 100  -  total flux cleaned = 0.0165203 Jy
! Component: 150  -  total flux cleaned = 0.0174533 Jy
! Component: 200  -  total flux cleaned = 0.0190112 Jy
! Total flux subtracted in 200 components = 0.0190112 Jy
! Clean residual min=-0.010905 max=0.010533 Jy/beam
! Clean residual mean=0.000010 rms=0.003034 Jy/beam
! Combined flux in latest and established models = 4.13481 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.00201486 Jy
! Component: 100  -  total flux cleaned = 0.00291279 Jy
! Component: 150  -  total flux cleaned = 0.00399753 Jy
! Component: 200  -  total flux cleaned = 0.00474379 Jy
! Total flux subtracted in 200 components = 0.00474379 Jy
! Clean residual min=-0.009184 max=0.008916 Jy/beam
! Clean residual mean=0.000011 rms=0.002632 Jy/beam
! Combined flux in latest and established models = 4.13956 Jy
! Adding 76 model components to the UV plane model.
! The established model now contains 111 components and 4.13956 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=3.326 mas, bmaj=10.89 mas, bpa=-18.51 degrees
! Estimated noise=9.79193 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00139118 Jy
! Component: 100  -  total flux cleaned = 0.00107838 Jy
! Component: 150  -  total flux cleaned = 0.00207666 Jy
! Component: 200  -  total flux cleaned = 0.0032201 Jy
! Total flux subtracted in 200 components = 0.0032201 Jy
! Clean residual min=-0.009594 max=0.008842 Jy/beam
! Clean residual mean=0.000006 rms=0.002168 Jy/beam
! Combined flux in latest and established models = 4.14278 Jy
selfcal
! Performing phase self-cal
! Adding 35 model components to the UV plane model.
! The established model now contains 144 components and 4.14278 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.062549Jy  sigma=0.129872
! Fit after  self-cal, rms=0.062157Jy  sigma=0.128703
wmodel J0854+2006_S_a_map.mod
! Writing 144 model components to file: J0854+2006_S_a_map.mod
wobs   J0854+2006_S_a_uvs.fits
! Writing UV FITS file: J0854+2006_S_a_uvs.fits
wwins  J0854+2006_S_a_map.win
! wwins: Wrote 3 windows to J0854+2006_S_a_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.0018949 Jy
! Component: 100  -  total flux cleaned = -0.00346593 Jy
! Component: 150  -  total flux cleaned = -0.0046751 Jy
! Component: 200  -  total flux cleaned = -0.00581619 Jy
! Component: 250  -  total flux cleaned = -0.00714529 Jy
! Component: 300  -  total flux cleaned = -0.00807109 Jy
! Component: 350  -  total flux cleaned = -0.00896421 Jy
! Component: 400  -  total flux cleaned = -0.0100353 Jy
! Component: 450  -  total flux cleaned = -0.0106548 Jy
! Component: 500  -  total flux cleaned = -0.0112527 Jy
! Component: 550  -  total flux cleaned = -0.0117308 Jy
! Component: 600  -  total flux cleaned = -0.0122861 Jy
! Component: 650  -  total flux cleaned = -0.012919 Jy
! Component: 700  -  total flux cleaned = -0.0133576 Jy
! Component: 750  -  total flux cleaned = -0.0139545 Jy
! Component: 800  -  total flux cleaned = -0.0146153 Jy
! Component: 850  -  total flux cleaned = -0.0150957 Jy
! Component: 900  -  total flux cleaned = -0.0155647 Jy
! Component: 950  -  total flux cleaned = -0.0160975 Jy
! Component: 1000  -  total flux cleaned = -0.0163209 Jy
! Component: 1050  -  total flux cleaned = -0.0165386 Jy
! Component: 1100  -  total flux cleaned = -0.0167521 Jy
! Component: 1150  -  total flux cleaned = -0.0168214 Jy
! Component: 1200  -  total flux cleaned = -0.0169562 Jy
! Component: 1250  -  total flux cleaned = -0.0171538 Jy
! Component: 1300  -  total flux cleaned = -0.0172822 Jy
! Component: 1350  -  total flux cleaned = -0.0173457 Jy
! Component: 1400  -  total flux cleaned = -0.0177797 Jy
! Component: 1450  -  total flux cleaned = -0.0180218 Jy
! Component: 1500  -  total flux cleaned = -0.0181415 Jy
! Component: 1550  -  total flux cleaned = -0.018377 Jy
! Component: 1600  -  total flux cleaned = -0.0186659 Jy
! Component: 1650  -  total flux cleaned = -0.0189499 Jy
! Component: 1700  -  total flux cleaned = -0.0192834 Jy
! Component: 1750  -  total flux cleaned = -0.0197735 Jy
! Component: 1800  -  total flux cleaned = -0.0200425 Jy
! Component: 1850  -  total flux cleaned = -0.0202002 Jy
! Component: 1900  -  total flux cleaned = -0.0203572 Jy
! Component: 1950  -  total flux cleaned = -0.0206142 Jy
! Component: 2000  -  total flux cleaned = -0.020766 Jy
! Component: 2050  -  total flux cleaned = -0.021162 Jy
! Component: 2100  -  total flux cleaned = -0.0213079 Jy
! Component: 2150  -  total flux cleaned = -0.0214525 Jy
! Component: 2200  -  total flux cleaned = -0.0216425 Jy
! Component: 2250  -  total flux cleaned = -0.0218755 Jy
! Component: 2300  -  total flux cleaned = -0.0218288 Jy
! Component: 2350  -  total flux cleaned = -0.0220561 Jy
! Component: 2400  -  total flux cleaned = -0.0223242 Jy
! Component: 2450  -  total flux cleaned = -0.0224115 Jy
! Component: 2500  -  total flux cleaned = -0.0224546 Jy
! Component: 2550  -  total flux cleaned = -0.0224534 Jy
! Component: 2600  -  total flux cleaned = -0.0226223 Jy
! Component: 2650  -  total flux cleaned = -0.0226637 Jy
! Component: 2700  -  total flux cleaned = -0.0227461 Jy
! Component: 2750  -  total flux cleaned = -0.0226648 Jy
! Component: 2800  -  total flux cleaned = -0.0227459 Jy
! Component: 2850  -  total flux cleaned = -0.0229044 Jy
! Component: 2900  -  total flux cleaned = -0.0230221 Jy
! Component: 2950  -  total flux cleaned = -0.0230608 Jy
! Component: 3000  -  total flux cleaned = -0.0230228 Jy
! Component: 3050  -  total flux cleaned = -0.0231365 Jy
! Component: 3100  -  total flux cleaned = -0.0231738 Jy
! Component: 3150  -  total flux cleaned = -0.0232108 Jy
! Component: 3200  -  total flux cleaned = -0.0232837 Jy
! Component: 3250  -  total flux cleaned = -0.0233913 Jy
! Component: 3300  -  total flux cleaned = -0.0233912 Jy
! Component: 3350  -  total flux cleaned = -0.0234621 Jy
! Component: 3400  -  total flux cleaned = -0.0234275 Jy
! Component: 3450  -  total flux cleaned = -0.0234613 Jy
! Component: 3500  -  total flux cleaned = -0.0233252 Jy
! Component: 3550  -  total flux cleaned = -0.0233257 Jy
! Component: 3600  -  total flux cleaned = -0.0232255 Jy
! Component: 3650  -  total flux cleaned = -0.0232584 Jy
! Component: 3700  -  total flux cleaned = -0.0231933 Jy
! Component: 3750  -  total flux cleaned = -0.0233232 Jy
! Component: 3800  -  total flux cleaned = -0.0233232 Jy
! Component: 3850  -  total flux cleaned = -0.0233868 Jy
! Component: 3900  -  total flux cleaned = -0.0234182 Jy
! Component: 3950  -  total flux cleaned = -0.0234495 Jy
! Component: 4000  -  total flux cleaned = -0.0235115 Jy
! Component: 4050  -  total flux cleaned = -0.0235116 Jy
! Total flux subtracted in 4096 components = -0.0235419 Jy
! Clean residual min=-0.001702 max=0.001560 Jy/beam
! Clean residual mean=0.000006 rms=0.000594 Jy/beam
! Combined flux in latest and established models = 4.11923 Jy
keep
! Adding 302 model components to the UV plane model.
! The established model now contains 446 components and 4.11923 Jy
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
wmap   /image_prod/ug002_uva//J0854+2006_S_a_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.326 x 10.89 at -18.51 degrees (North through East)
! Clean map  min=-0.010513  max=3.5303 Jy/beam
! Writing clean map to FITS file: /image_prod/ug002_uva//J0854+2006_S_a_map.fits
wmodel /image_prod/ug002_uva//J0854+2006_S_a_map.mod
! Writing 446 model components to file: /image_prod/ug002_uva//J0854+2006_S_a_map.mod
wobs   /image_prod/ug002_uva//J0854+2006_S_a_uvs.fits
! Writing UV FITS file: /image_prod/ug002_uva//J0854+2006_S_a_uvs.fits
wwins  /image_prod/ug002_uva//J0854+2006_S_a_map.win
! wwins: Wrote 4 windows to /image_prod/ug002_uva//J0854+2006_S_a_map.win
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 3.53034 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00058406 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.0496321 0.0496321 0.0992641 0.198528 0.397056 0.794113 1.58823 3.17645 6.3529 12.7058 25.4116 50.8232
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0017 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       3.53 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  6044.48 
quit
! Quitting program
! Log file /image_prod/ug002_uva//J0854+2006_S_a_dfm.log closed on Wed Mar  4 22:38:15 2020