! Started logfile: /vlbi/uf001/uva//J1209+433A_S_dfm.log on Mon Jan 14 01:41:02 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//J1209+433A_S_uva.fits
! Reading UV FITS file: /vlbi/uf001/uva//J1209+433A_S_uva.fits
! AN table 1: 36 integrations on 10 of 10 possible baselines.
! AN table 2: 229 integrations on 6 of 6 possible baselines.
! AN table 3: 37 integrations on 1 of 1 possible baselines.
! Apparent sampling: 0.541502 visibilities/baseline/integration-bin.
! Found source: J1209+433A
! 
! 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 2877 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 J1209+433A_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 556 telescope corrections were flagged in sub-array 2.
!  A total of 74 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 2.
!  A total of 556 telescope corrections were flagged in sub-array 2.
!  A total of 74 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 3.
!  A total of 556 telescope corrections were flagged in sub-array 2.
!  A total of 74 telescope corrections were flagged in sub-array 3.
! 
! Fit before self-cal, rms=1.026160Jy  sigma=21.285511
! Fit after  self-cal, rms=0.926511Jy  sigma=19.420008
! 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=10.04 mas, bmaj=27.91 mas, bpa=-75.4 degrees
! Estimated noise=0.884706 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.0567195 Jy
! Component: 100  -  total flux cleaned = 0.0872006 Jy
! Total flux subtracted in 100 components = 0.0872006 Jy
! Clean residual min=-0.017094 max=0.017193 Jy/beam
! Clean residual mean=-0.000402 rms=0.008505 Jy/beam
! Combined flux in latest and established models = 0.0872006 Jy
! Performing phase self-cal
! Adding 65 model components to the UV plane model.
! The established model now contains 65 components and 0.0872006 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.148248Jy  sigma=2.059603
! Fit after  self-cal, rms=0.146456Jy  sigma=2.043451
! 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=12.4 mas, bmaj=30.27 mas, bpa=-75.51 degrees
! Estimated noise=0.697844 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.0257394 Jy
! Component: 100  -  total flux cleaned = 0.0483809 Jy
! Component: 150  -  total flux cleaned = 0.0660143 Jy
! Component: 200  -  total flux cleaned = 0.0736453 Jy
! Total flux subtracted in 200 components = 0.0736453 Jy
! Clean residual min=-0.013625 max=0.013491 Jy/beam
! Clean residual mean=-0.000854 rms=0.005711 Jy/beam
! Combined flux in latest and established models = 0.160846 Jy
! Performing phase self-cal
! Adding 103 model components to the UV plane model.
! The established model now contains 166 components and 0.160846 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.130042Jy  sigma=1.532001
! Fit after  self-cal, rms=0.127558Jy  sigma=1.507526
! Inverting map 
print "============== Finished natural weighting clean =============="
! ============== Finished natural weighting clean ============== 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 12.4 x 30.27 at -75.51 degrees (North through East)
! Clean map  min=-0.015047  max=0.070236 Jy/beam
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.28     KP        1.07     LA        0.92     OV        0.81     
!   PT        0.66     
! 
!  A total of 216 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        1.42     KP        1.13     LA        0.90     PT        0.73     
! 
!  Telescope amplitude corrections in sub-array 3:
!   MK        0.00*    OV        0.00*    
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.34     KP        1.31     LA        0.81     OV        0.73     
!   PT        0.88     
! 
!  A total of 216 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        1.24     KP        1.13     LA        0.89     PT        1.03     
! 
!  Telescope amplitude corrections in sub-array 3:
!   MK        0.00*    OV        0.00*    
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.13     KP        1.24     LA        0.93     OV        0.75     
!   PT        0.90     
! 
!  A total of 216 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   FD        1.17     KP        1.06     LA        0.89     PT        0.84     
! 
!  Telescope amplitude corrections in sub-array 3:
!   MK        0.00*    OV        0.00*    
! 
! 
! Fit before self-cal, rms=0.127558Jy  sigma=1.507526
! Fit after  self-cal, rms=0.132458Jy  sigma=1.475922
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=11.47 mas, bmaj=27.56 mas, bpa=-81.88 degrees
! Estimated noise=0.793298 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.00233213 Jy
! Component: 100  -  total flux cleaned = 0.0020844 Jy
! Component: 150  -  total flux cleaned = 0.00194626 Jy
! Component: 200  -  total flux cleaned = 0.00314652 Jy
! Total flux subtracted in 200 components = 0.00314652 Jy
! Clean residual min=-0.001597 max=0.001598 Jy/beam
! Clean residual mean=-0.000035 rms=0.000904 Jy/beam
! Combined flux in latest and established models = 0.163993 Jy
! Performing phase self-cal
! Adding 86 model components to the UV plane model.
! The established model now contains 252 components and 0.163993 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.131471Jy  sigma=1.473190
! Fit after  self-cal, rms=0.131455Jy  sigma=1.470097
! 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.131455Jy  sigma=1.470097
! Fit after  self-cal, rms=0.140425Jy  sigma=1.468325
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=11.58 mas, bmaj=27.14 mas, bpa=-82.83 degrees
! Estimated noise=0.791912 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00015198 Jy
! Component: 100  -  total flux cleaned = 0.00105655 Jy
! Component: 150  -  total flux cleaned = 0.0017755 Jy
! Component: 200  -  total flux cleaned = 0.00225978 Jy
! Total flux subtracted in 200 components = 0.00225978 Jy
! Clean residual min=-0.001542 max=0.001561 Jy/beam
! Clean residual mean=-0.000052 rms=0.000869 Jy/beam
! Combined flux in latest and established models = 0.166252 Jy
selfcal
! Performing phase self-cal
! Adding 110 model components to the UV plane model.
! The established model now contains 359 components and 0.166252 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.139583Jy  sigma=1.467462
! Fit after  self-cal, rms=0.139688Jy  sigma=1.465841
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.000160335 Jy
! Component: 100  -  total flux cleaned = -0.000259174 Jy
! Component: 150  -  total flux cleaned = -0.000572309 Jy
! Component: 200  -  total flux cleaned = -0.00125231 Jy
! Total flux subtracted in 200 components = -0.00125231 Jy
! Clean residual min=-0.001223 max=0.001213 Jy/beam
! Clean residual mean=-0.000038 rms=0.000656 Jy/beam
! Combined flux in latest and established models = 0.165 Jy
! Performing phase self-cal
! Adding 107 model components to the UV plane model.
! The established model now contains 466 components and 0.165 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.139044Jy  sigma=1.465791
! Fit after  self-cal, rms=0.139019Jy  sigma=1.465195
! 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.139019Jy  sigma=1.465195
! Fit after  self-cal, rms=0.157558Jy  sigma=1.462801
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=11.59 mas, bmaj=27.09 mas, bpa=-83.2 degrees
! Estimated noise=0.79091 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000766801 Jy
! Component: 100  -  total flux cleaned = -0.00142762 Jy
! Component: 150  -  total flux cleaned = -0.00203936 Jy
! Component: 200  -  total flux cleaned = -0.00317507 Jy
! Total flux subtracted in 200 components = -0.00317507 Jy
! Clean residual min=-0.001545 max=0.001528 Jy/beam
! Clean residual mean=-0.000046 rms=0.000895 Jy/beam
! Combined flux in latest and established models = 0.161825 Jy
selfcal
! Performing phase self-cal
! Adding 119 model components to the UV plane model.
! The established model now contains 580 components and 0.161825 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.156690Jy  sigma=1.462609
! Fit after  self-cal, rms=0.156761Jy  sigma=1.461637
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.000225312 Jy
! Component: 100  -  total flux cleaned = 0.000230016 Jy
! Component: 150  -  total flux cleaned = -0.000627822 Jy
! Component: 200  -  total flux cleaned = -0.00127015 Jy
! Total flux subtracted in 200 components = -0.00127015 Jy
! Clean residual min=-0.001324 max=0.001305 Jy/beam
! Clean residual mean=-0.000035 rms=0.000697 Jy/beam
! Combined flux in latest and established models = 0.160555 Jy
! Performing phase self-cal
! Adding 91 model components to the UV plane model.
! The established model now contains 663 components and 0.160555 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.156067Jy  sigma=1.461623
! Fit after  self-cal, rms=0.156016Jy  sigma=1.460638
! 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.156016Jy  sigma=1.460638
! Fit after  self-cal, rms=0.174392Jy  sigma=1.458515
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=13.02 mas, bmaj=27.54 mas, bpa=-86.07 degrees
! Estimated noise=0.783493 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -2.90377e-05 Jy
! Component: 100  -  total flux cleaned = -0.000571798 Jy
! Component: 150  -  total flux cleaned = -0.00210787 Jy
! Component: 200  -  total flux cleaned = -0.0032675 Jy
! Total flux subtracted in 200 components = -0.0032675 Jy
! Clean residual min=-0.001564 max=0.001548 Jy/beam
! Clean residual mean=-0.000048 rms=0.000911 Jy/beam
! Combined flux in latest and established models = 0.157287 Jy
selfcal
! Performing phase self-cal
! Adding 109 model components to the UV plane model.
! The established model now contains 752 components and 0.157287 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.173510Jy  sigma=1.458074
! Fit after  self-cal, rms=0.173476Jy  sigma=1.457277
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.000392604 Jy
! Component: 100  -  total flux cleaned = 0.000407669 Jy
! Component: 150  -  total flux cleaned = 0.000410973 Jy
! Component: 200  -  total flux cleaned = -0.000139813 Jy
! Total flux subtracted in 200 components = -0.000139813 Jy
! Clean residual min=-0.001289 max=0.001277 Jy/beam
! Clean residual mean=-0.000037 rms=0.000729 Jy/beam
! Combined flux in latest and established models = 0.157147 Jy
! Performing phase self-cal
! Adding 102 model components to the UV plane model.
! The established model now contains 831 components and 0.157147 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.172785Jy  sigma=1.457088
! Fit after  self-cal, rms=0.172761Jy  sigma=1.456380
! 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=11.57 mas, bmaj=27.22 mas, bpa=-83.39 degrees
! Estimated noise=0.786714 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 11.57 x 27.22 at -83.39 degrees (North through East)
! Clean map  min=-0.0062454  max=0.054507 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=9.897 mas, bmaj=24.29 mas, bpa=-82.58 degrees
! Estimated noise=0.928976 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00247457 Jy
! Component: 100  -  total flux cleaned = -0.0053333 Jy
! Component: 150  -  total flux cleaned = -0.000371002 Jy
! Component: 200  -  total flux cleaned = -0.00142349 Jy
! Total flux subtracted in 200 components = -0.00142349 Jy
! Clean residual min=-0.009542 max=0.009560 Jy/beam
! Clean residual mean=-0.000814 rms=0.004558 Jy/beam
! Combined flux in latest and established models = -0.00142349 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 95 model components to the UV plane model.
! The established model now contains 95 components and -0.00142349 Jy
! Inverting map and beam 
! Estimated beam: bmin=11.57 mas, bmaj=27.22 mas, bpa=-83.39 degrees
! Estimated noise=0.786714 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.0116715 Jy
! Component: 100  -  total flux cleaned = -0.0177132 Jy
! Component: 150  -  total flux cleaned = -0.0195122 Jy
! Component: 200  -  total flux cleaned = -0.0173521 Jy
! Total flux subtracted in 200 components = -0.0173521 Jy
! Clean residual min=-0.006418 max=0.006510 Jy/beam
! Clean residual mean=-0.000534 rms=0.003500 Jy/beam
! Combined flux in latest and established models = -0.0187756 Jy
! Adding 125 model components to the UV plane model.
! The established model now contains 216 components and -0.0187756 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=13.02 mas, bmaj=27.54 mas, bpa=-86.07 degrees
! Estimated noise=0.783493 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00194154 Jy
! Component: 100  -  total flux cleaned = -0.000839961 Jy
! Component: 150  -  total flux cleaned = -0.00174807 Jy
! Component: 200  -  total flux cleaned = -0.00309615 Jy
! Total flux subtracted in 200 components = -0.00309615 Jy
! Clean residual min=-0.003127 max=0.003117 Jy/beam
! Clean residual mean=-0.000172 rms=0.001816 Jy/beam
! Combined flux in latest and established models = -0.0218718 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.000557232 Jy
! Component: 100  -  total flux cleaned = -0.00109729 Jy
! Component: 150  -  total flux cleaned = -0.000546223 Jy
! Component: 200  -  total flux cleaned = -0.000357628 Jy
! Total flux subtracted in 200 components = -0.000357628 Jy
! Clean residual min=-0.003176 max=0.003195 Jy/beam
! Clean residual mean=-0.000147 rms=0.001491 Jy/beam
! Combined flux in latest and established models = -0.0222294 Jy
! Adding 229 model components to the UV plane model.
! The established model now contains 444 components and -0.0222294 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=11.57 mas, bmaj=27.22 mas, bpa=-83.39 degrees
! Estimated noise=0.786714 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000885812 Jy
! Component: 100  -  total flux cleaned = -0.00044085 Jy
! Component: 150  -  total flux cleaned = -0.000300805 Jy
! Component: 200  -  total flux cleaned = -0.000298042 Jy
! Total flux subtracted in 200 components = -0.000298042 Jy
! Clean residual min=-0.002305 max=0.002328 Jy/beam
! Clean residual mean=-0.000083 rms=0.001137 Jy/beam
! Combined flux in latest and established models = -0.0225274 Jy
selfcal
! Performing phase self-cal
! Adding 88 model components to the UV plane model.
! The established model now contains 529 components and -0.0225274 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Fit before self-cal, rms=0.174250Jy  sigma=1.463795
! Fit after  self-cal, rms=0.174190Jy  sigma=1.460695
wmodel J1209+433A_S_map.mod
! Writing 529 model components to file: J1209+433A_S_map.mod
wobs   J1209+433A_S_uvs.fits
! Writing UV FITS file: J1209+433A_S_uvs.fits
wwins  J1209+433A_S_map.win
! wwins: No CLEAN windows to write.
! Error occured in command: wwins
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
! 
wmap   /vlbi/uf001/uva//J1209+433A_S_map.fits
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 11.57 x 27.22 at -83.39 degrees (North through East)
! Clean map  min=-0.018231  max=0.01885 Jy/beam
! Writing clean map to FITS file: /vlbi/uf001/uva//J1209+433A_S_map.fits
wmodel /vlbi/uf001/uva//J1209+433A_S_map.mod
! Writing 529 model components to file: /vlbi/uf001/uva//J1209+433A_S_map.mod
wobs   /vlbi/uf001/uva//J1209+433A_S_uvs.fits
! Writing UV FITS file: /vlbi/uf001/uva//J1209+433A_S_uvs.fits
wwins  /vlbi/uf001/uva//J1209+433A_S_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.0188497 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00110825 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -17.6381 17.6381 35.2763 70.5526
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0033 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       0.018 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  17.008 
quit
! Quitting program
! Log file /vlbi/uf001/uva//J1209+433A_S_dfm.log closed on Mon Jan 14 01:41:06 2019