! Started logfile: /vlbi/uf001/uva//J2301+3726_S_dfm.log on Mon Jan 14 06:30:48 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//J2301+3726_S_uva.fits
! Reading UV FITS file: /vlbi/uf001/uva//J2301+3726_S_uva.fits
! AN table 1: 42 integrations on 28 of 28 possible baselines.
! Apparent sampling: 0.789116 visibilities/baseline/integration-bin.
! Found source: J2301+3726
! 
! 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 2784 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 J2301+3726_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 83 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 83 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 83 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=1.187052Jy  sigma=170.047500
! Fit after  self-cal, rms=0.611488Jy  sigma=165.402144
! 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.156 mas, bmaj=15.73 mas, bpa=-0.9008 degrees
! Estimated noise=0.2055 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.109698 Jy
! Component: 100  -  total flux cleaned = 0.189006 Jy
! Total flux subtracted in 100 components = 0.189006 Jy
! Clean residual min=-0.046592 max=0.061121 Jy/beam
! Clean residual mean=0.000000 rms=0.012374 Jy/beam
! Combined flux in latest and established models = 0.189006 Jy
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 12 components and 0.189006 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.550915Jy  sigma=14.129481
! Fit after  self-cal, rms=0.522891Jy  sigma=13.716373
! 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=2.337 mas, bmaj=15.3 mas, bpa=-0.1836 degrees
! Estimated noise=0.214058 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 (-3, 2).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.125534 Jy
! Component: 100  -  total flux cleaned = 0.219458 Jy
! Component: 150  -  total flux cleaned = 0.291112 Jy
! Component: 200  -  total flux cleaned = 0.346616 Jy
! Total flux subtracted in 200 components = 0.346616 Jy
! Clean residual min=-0.053779 max=0.071054 Jy/beam
! Clean residual mean=0.000015 rms=0.016875 Jy/beam
! Combined flux in latest and established models = 0.535622 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 26 components and 0.535622 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.292695Jy  sigma=16.557200
! Fit after  self-cal, rms=0.283729Jy  sigma=16.485501
! Inverting map 
print "============== Finished natural weighting clean =============="
! ============== Finished natural weighting clean ============== 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.337 x 15.3 at -0.1836 degrees (North through East)
! Clean map  min=-0.051301  max=0.27261 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.85     HN        0.79     KP        0.84     LA        0.76     
!   MK       127.71     NL        0.88     OV        0.77     PT        0.83     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.83     HN        0.83     KP        0.84     LA        0.75     
!   MK       128.75     NL        0.85     OV        0.73     PT        0.83     
! 
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   FD        0.84     HN        0.81     KP        0.84     LA        0.77     
!   MK       176.54     NL        0.84     OV        0.76     PT        0.82     
! 
! 
! Fit before self-cal, rms=0.283729Jy  sigma=16.485501
! Fit after  self-cal, rms=0.144891Jy  sigma=3.119635
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=5.348 mas, bmaj=11.13 mas, bpa=-16.78 degrees
! Estimated noise=0.598407 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.0179527 Jy
! Component: 100  -  total flux cleaned = 0.0239538 Jy
! Component: 150  -  total flux cleaned = 0.0228615 Jy
! Component: 200  -  total flux cleaned = 0.0212306 Jy
! Total flux subtracted in 200 components = 0.0212306 Jy
! Clean residual min=-0.036302 max=0.015658 Jy/beam
! Clean residual mean=0.000008 rms=0.004538 Jy/beam
! Combined flux in latest and established models = 0.556853 Jy
! Performing phase self-cal
! Adding 13 model components to the UV plane model.
! The established model now contains 38 components and 0.556853 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.121549Jy  sigma=1.878297
! Fit after  self-cal, rms=0.116373Jy  sigma=1.751131
! 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.116373Jy  sigma=1.751131
! Fit after  self-cal, rms=0.174645Jy  sigma=1.436818
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=5.282 mas, bmaj=11.37 mas, bpa=-15.43 degrees
! Estimated noise=0.595839 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 3.63635e-05 Jy
! Component: 100  -  total flux cleaned = -0.000672784 Jy
! Component: 150  -  total flux cleaned = 0.000224962 Jy
! Component: 200  -  total flux cleaned = 0.000768004 Jy
! Total flux subtracted in 200 components = 0.000768004 Jy
! Clean residual min=-0.014128 max=0.008988 Jy/beam
! Clean residual mean=0.000008 rms=0.002804 Jy/beam
! Combined flux in latest and established models = 0.557621 Jy
selfcal
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 59 components and 0.557621 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.173298Jy  sigma=1.305516
! Fit after  self-cal, rms=0.173162Jy  sigma=1.298728
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 = -1.7196e-05 Jy
! Component: 100  -  total flux cleaned = 0.000220366 Jy
! Component: 150  -  total flux cleaned = 0.000667045 Jy
! Component: 200  -  total flux cleaned = 0.00066672 Jy
! Total flux subtracted in 200 components = 0.00066672 Jy
! Clean residual min=-0.013160 max=0.007908 Jy/beam
! Clean residual mean=0.000008 rms=0.002541 Jy/beam
! Combined flux in latest and established models = 0.558287 Jy
! Performing phase self-cal
! Adding 25 model components to the UV plane model.
! The established model now contains 74 components and 0.558287 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.172668Jy  sigma=1.237080
! Fit after  self-cal, rms=0.172542Jy  sigma=1.234822
! 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.172542Jy  sigma=1.234822
! Fit after  self-cal, rms=0.187808Jy  sigma=1.223267
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=5.289 mas, bmaj=11.32 mas, bpa=-15.49 degrees
! Estimated noise=0.595333 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00121936 Jy
! Component: 100  -  total flux cleaned = 0.0022394 Jy
! Component: 150  -  total flux cleaned = 0.00297745 Jy
! Component: 200  -  total flux cleaned = 0.00331235 Jy
! Total flux subtracted in 200 components = 0.00331235 Jy
! Clean residual min=-0.012221 max=0.007453 Jy/beam
! Clean residual mean=0.000008 rms=0.002302 Jy/beam
! Combined flux in latest and established models = 0.5616 Jy
selfcal
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 82 components and 0.5616 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.187416Jy  sigma=1.174388
! Fit after  self-cal, rms=0.187300Jy  sigma=1.172619
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.78136e-05 Jy
! Component: 100  -  total flux cleaned = 0.000576269 Jy
! Component: 150  -  total flux cleaned = 0.000873516 Jy
! Component: 200  -  total flux cleaned = 0.00101674 Jy
! Total flux subtracted in 200 components = 0.00101674 Jy
! Clean residual min=-0.011365 max=0.007240 Jy/beam
! Clean residual mean=0.000008 rms=0.002161 Jy/beam
! Combined flux in latest and established models = 0.562616 Jy
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 89 components and 0.562616 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.187102Jy  sigma=1.141965
! Fit after  self-cal, rms=0.187020Jy  sigma=1.141020
! 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.187020Jy  sigma=1.141020
! Fit after  self-cal, rms=0.197479Jy  sigma=1.137878
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.021 mas, bmaj=13.44 mas, bpa=-18.4 degrees
! Estimated noise=0.655326 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000968175 Jy
! Component: 100  -  total flux cleaned = 0.00145939 Jy
! Component: 150  -  total flux cleaned = 0.00172405 Jy
! Component: 200  -  total flux cleaned = 0.00183818 Jy
! Total flux subtracted in 200 components = 0.00183818 Jy
! Clean residual min=-0.008673 max=0.006276 Jy/beam
! Clean residual mean=0.000005 rms=0.001888 Jy/beam
! Combined flux in latest and established models = 0.564454 Jy
selfcal
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 100 components and 0.564455 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.197426Jy  sigma=1.132250
! Fit after  self-cal, rms=0.197394Jy  sigma=1.131132
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.000136994 Jy
! Component: 100  -  total flux cleaned = 0.000380168 Jy
! Component: 150  -  total flux cleaned = 0.000328378 Jy
! Component: 200  -  total flux cleaned = 0.000325823 Jy
! Total flux subtracted in 200 components = 0.000325823 Jy
! Clean residual min=-0.008597 max=0.006275 Jy/beam
! Clean residual mean=0.000006 rms=0.001873 Jy/beam
! Combined flux in latest and established models = 0.56478 Jy
! Performing phase self-cal
! Adding 24 model components to the UV plane model.
! The established model now contains 110 components and 0.56478 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.197349Jy  sigma=1.125925
! Fit after  self-cal, rms=0.197318Jy  sigma=1.124965
! 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=5.303 mas, bmaj=11.29 mas, bpa=-15.64 degrees
! Estimated noise=0.595657 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 5.303 x 11.29 at -15.64 degrees (North through East)
! Clean map  min=-0.01114  max=0.43531 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=4.65 mas, bmaj=10.27 mas, bpa=-15.94 degrees
! Estimated noise=0.807525 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.336016 Jy
! Component: 100  -  total flux cleaned = 0.427568 Jy
! Component: 150  -  total flux cleaned = 0.476612 Jy
! Component: 200  -  total flux cleaned = 0.505586 Jy
! Total flux subtracted in 200 components = 0.505586 Jy
! Clean residual min=-0.014065 max=0.018647 Jy/beam
! Clean residual mean=-0.000009 rms=0.003294 Jy/beam
! Combined flux in latest and established models = 0.505587 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 27 model components to the UV plane model.
! The established model now contains 27 components and 0.505587 Jy
! Inverting map and beam 
! Estimated beam: bmin=5.303 mas, bmaj=11.29 mas, bpa=-15.64 degrees
! Estimated noise=0.595657 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)
! Added new window around map position (3, -3.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0229834 Jy
! Component: 100  -  total flux cleaned = 0.0345408 Jy
! Component: 150  -  total flux cleaned = 0.0390529 Jy
! Component: 200  -  total flux cleaned = 0.0424204 Jy
! Total flux subtracted in 200 components = 0.0424204 Jy
! Clean residual min=-0.008030 max=0.008364 Jy/beam
! Clean residual mean=0.000001 rms=0.002002 Jy/beam
! Combined flux in latest and established models = 0.548007 Jy
! Adding 41 model components to the UV plane model.
! The established model now contains 66 components and 0.548007 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.021 mas, bmaj=13.44 mas, bpa=-18.4 degrees
! Estimated noise=0.655326 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00667228 Jy
! Component: 100  -  total flux cleaned = 0.00929227 Jy
! Component: 150  -  total flux cleaned = 0.0107297 Jy
! Component: 200  -  total flux cleaned = 0.0111887 Jy
! Total flux subtracted in 200 components = 0.0111887 Jy
! Clean residual min=-0.007835 max=0.006424 Jy/beam
! Clean residual mean=0.000003 rms=0.001896 Jy/beam
! Combined flux in latest and established models = 0.559196 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.000404949 Jy
! Component: 100  -  total flux cleaned = 0.000521647 Jy
! Component: 150  -  total flux cleaned = 0.000739105 Jy
! Component: 200  -  total flux cleaned = 0.000842245 Jy
! Total flux subtracted in 200 components = 0.000842245 Jy
! Clean residual min=-0.007604 max=0.005986 Jy/beam
! Clean residual mean=0.000003 rms=0.001787 Jy/beam
! Combined flux in latest and established models = 0.560038 Jy
! Adding 87 model components to the UV plane model.
! The established model now contains 145 components and 0.560038 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=5.303 mas, bmaj=11.29 mas, bpa=-15.64 degrees
! Estimated noise=0.595657 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000392085 Jy
! Component: 100  -  total flux cleaned = 0.00119478 Jy
! Component: 150  -  total flux cleaned = 0.00166095 Jy
! Component: 200  -  total flux cleaned = 0.00176554 Jy
! Total flux subtracted in 200 components = 0.00176554 Jy
! Clean residual min=-0.007185 max=0.005849 Jy/beam
! Clean residual mean=0.000004 rms=0.001571 Jy/beam
! Combined flux in latest and established models = 0.561803 Jy
selfcal
! Performing phase self-cal
! Adding 22 model components to the UV plane model.
! The established model now contains 163 components and 0.561803 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.202655Jy  sigma=1.025169
! Fit after  self-cal, rms=0.197389Jy  sigma=1.015818
wmodel J2301+3726_S_map.mod
! Writing 163 model components to file: J2301+3726_S_map.mod
wobs   J2301+3726_S_uvs.fits
! Writing UV FITS file: J2301+3726_S_uvs.fits
wwins  J2301+3726_S_map.win
! wwins: Wrote 3 windows to J2301+3726_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 = -0.0018699 Jy
! Component: 100  -  total flux cleaned = -0.00206507 Jy
! Component: 150  -  total flux cleaned = -0.00261863 Jy
! Component: 200  -  total flux cleaned = -0.00304857 Jy
! Component: 250  -  total flux cleaned = -0.00369726 Jy
! Component: 300  -  total flux cleaned = -0.00431099 Jy
! Component: 350  -  total flux cleaned = -0.005044 Jy
! Component: 400  -  total flux cleaned = -0.00560055 Jy
! Component: 450  -  total flux cleaned = -0.00586991 Jy
! Component: 500  -  total flux cleaned = -0.00625679 Jy
! Component: 550  -  total flux cleaned = -0.00662944 Jy
! Component: 600  -  total flux cleaned = -0.00680847 Jy
! Component: 650  -  total flux cleaned = -0.00715554 Jy
! Component: 700  -  total flux cleaned = -0.00765723 Jy
! Component: 750  -  total flux cleaned = -0.00798008 Jy
! Component: 800  -  total flux cleaned = -0.00829371 Jy
! Component: 850  -  total flux cleaned = -0.0085977 Jy
! Component: 900  -  total flux cleaned = -0.00879402 Jy
! Component: 950  -  total flux cleaned = -0.00903301 Jy
! Component: 1000  -  total flux cleaned = -0.00922097 Jy
! Component: 1050  -  total flux cleaned = -0.00958592 Jy
! Component: 1100  -  total flux cleaned = -0.00971876 Jy
! Component: 1150  -  total flux cleaned = -0.0098915 Jy
! Component: 1200  -  total flux cleaned = -0.010103 Jy
! Component: 1250  -  total flux cleaned = -0.0102681 Jy
! Component: 1300  -  total flux cleaned = -0.0104702 Jy
! Component: 1350  -  total flux cleaned = -0.0105888 Jy
! Component: 1400  -  total flux cleaned = -0.0107448 Jy
! Component: 1450  -  total flux cleaned = -0.0108595 Jy
! Component: 1500  -  total flux cleaned = -0.0110461 Jy
! Component: 1550  -  total flux cleaned = -0.0110823 Jy
! Component: 1600  -  total flux cleaned = -0.0111916 Jy
! Component: 1650  -  total flux cleaned = -0.0112624 Jy
! Component: 1700  -  total flux cleaned = -0.0112968 Jy
! Component: 1750  -  total flux cleaned = -0.0113649 Jy
! Component: 1800  -  total flux cleaned = -0.011432 Jy
! Component: 1850  -  total flux cleaned = -0.0113667 Jy
! Component: 1900  -  total flux cleaned = -0.0114317 Jy
! Component: 1950  -  total flux cleaned = -0.011463 Jy
! Component: 2000  -  total flux cleaned = -0.0114315 Jy
! Component: 2050  -  total flux cleaned = -0.0114621 Jy
! Component: 2100  -  total flux cleaned = -0.0114921 Jy
! Component: 2150  -  total flux cleaned = -0.0114624 Jy
! Component: 2200  -  total flux cleaned = -0.0114042 Jy
! Component: 2250  -  total flux cleaned = -0.0114327 Jy
! Component: 2300  -  total flux cleaned = -0.011348 Jy
! Component: 2350  -  total flux cleaned = -0.0114316 Jy
! Component: 2400  -  total flux cleaned = -0.011404 Jy
! Component: 2450  -  total flux cleaned = -0.0114312 Jy
! Component: 2500  -  total flux cleaned = -0.0115907 Jy
! Component: 2550  -  total flux cleaned = -0.0116426 Jy
! Component: 2600  -  total flux cleaned = -0.0116677 Jy
! Component: 2650  -  total flux cleaned = -0.0117448 Jy
! Component: 2700  -  total flux cleaned = -0.011896 Jy
! Component: 2750  -  total flux cleaned = -0.011946 Jy
! Component: 2800  -  total flux cleaned = -0.0119951 Jy
! Component: 2850  -  total flux cleaned = -0.0120922 Jy
! Component: 2900  -  total flux cleaned = -0.0122123 Jy
! Component: 2950  -  total flux cleaned = -0.0122123 Jy
! Component: 3000  -  total flux cleaned = -0.0122823 Jy
! Component: 3050  -  total flux cleaned = -0.0123287 Jy
! Component: 3100  -  total flux cleaned = -0.0123512 Jy
! Component: 3150  -  total flux cleaned = -0.0124648 Jy
! Component: 3200  -  total flux cleaned = -0.0125321 Jy
! Component: 3250  -  total flux cleaned = -0.0125988 Jy
! Component: 3300  -  total flux cleaned = -0.0126869 Jy
! Component: 3350  -  total flux cleaned = -0.0128611 Jy
! Component: 3400  -  total flux cleaned = -0.0128828 Jy
! Component: 3450  -  total flux cleaned = -0.0130536 Jy
! Component: 3500  -  total flux cleaned = -0.0131386 Jy
! Component: 3550  -  total flux cleaned = -0.0133067 Jy
! Component: 3600  -  total flux cleaned = -0.0134107 Jy
! Component: 3650  -  total flux cleaned = -0.0135142 Jy
! Component: 3700  -  total flux cleaned = -0.0136988 Jy
! Component: 3750  -  total flux cleaned = -0.0137191 Jy
! Component: 3800  -  total flux cleaned = -0.0138203 Jy
! Component: 3850  -  total flux cleaned = -0.0139407 Jy
! Component: 3900  -  total flux cleaned = -0.0140604 Jy
! Component: 3950  -  total flux cleaned = -0.0141793 Jy
! Component: 4000  -  total flux cleaned = -0.014258 Jy
! Component: 4050  -  total flux cleaned = -0.0143557 Jy
! Total flux subtracted in 4096 components = -0.0144334 Jy
! Clean residual min=-0.001105 max=0.001097 Jy/beam
! Clean residual mean=0.000007 rms=0.000378 Jy/beam
! Combined flux in latest and established models = 0.54737 Jy
keep
! Adding 370 model components to the UV plane model.
! The established model now contains 533 components and 0.54737 Jy
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
wmap   /vlbi/uf001/uva//J2301+3726_S_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 5.303 x 11.29 at -15.64 degrees (North through East)
! Clean map  min=-0.0074059  max=0.44766 Jy/beam
! Writing clean map to FITS file: /vlbi/uf001/uva//J2301+3726_S_map.fits
wmodel /vlbi/uf001/uva//J2301+3726_S_map.mod
! Writing 533 model components to file: /vlbi/uf001/uva//J2301+3726_S_map.mod
wobs   /vlbi/uf001/uva//J2301+3726_S_uvs.fits
! Writing UV FITS file: /vlbi/uf001/uva//J2301+3726_S_uvs.fits
wwins  /vlbi/uf001/uva//J2301+3726_S_map.win
! wwins: Wrote 4 windows to /vlbi/uf001/uva//J2301+3726_S_map.win
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 0.447657 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000375946 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.251943 0.251943 0.503886 1.00777 2.01554 4.03109 8.06217 16.1243 32.2487 64.4974
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0011 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       0.447 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  1190.75 
quit
! Quitting program
! Log file /vlbi/uf001/uva//J2301+3726_S_dfm.log closed on Mon Jan 14 06:30:52 2019