! Started logfile: /scr/pima/rdv76_uvs/J1753+2848_X_dfm.log on Fri Jun 15 03:56:18 2018
float field_size
field_size =  1024
float field_cell
field_cell =   0.2000
float taper_size
taper_size =   70.0
float freq
freq       =    8.40599
obs     /scr/pima/rdv76_uvs/J1753+2848_X_uva.fits
! Reading UV FITS file: /scr/pima/rdv76_uvs/J1753+2848_X_uva.fits
! AN table 1: 16 integrations on 10 of 10 possible baselines.
! AN table 2: 32 integrations on 55 of 55 possible baselines.
! AN table 3: 31 integrations on 55 of 55 possible baselines.
! AN table 4: 160 integrations on 15 of 15 possible baselines.
! Apparent sampling: 0.443154 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J1753+2848
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  8.40599e+09        8e+06          1        8e+06
!  02        2  8.47599e+09        8e+06          1        8e+06
!  03        3  8.79099e+09        8e+06          1        8e+06
!  04        4  8.89599e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 79 lines of history.
! 
! Reading 10680 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 4 channels.
! Selecting polarization: RR,  channels: 1..4
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
print "field_size =", field_size
! field_size = 1024 
print "field_cell =", field_cell
! field_cell = 0.2 
print "taper_size =", taper_size
! taper_size = 70 
![@/opt64/bin/pima_mupet_01.dfm J1753+2848_X]
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.200x0.200 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 43 telescope corrections were flagged in sub-array 1.
!  A total of 169 telescope corrections were flagged in sub-array 2.
!  A total of 172 telescope corrections were flagged in sub-array 3.
!  A total of 729 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 2.
!  A total of 43 telescope corrections were flagged in sub-array 1.
!  A total of 169 telescope corrections were flagged in sub-array 2.
!  A total of 172 telescope corrections were flagged in sub-array 3.
!  A total of 729 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
!  A total of 43 telescope corrections were flagged in sub-array 1.
!  A total of 169 telescope corrections were flagged in sub-array 2.
!  A total of 172 telescope corrections were flagged in sub-array 3.
!  A total of 729 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 4.
!  A total of 43 telescope corrections were flagged in sub-array 1.
!  A total of 169 telescope corrections were flagged in sub-array 2.
!  A total of 172 telescope corrections were flagged in sub-array 3.
!  A total of 729 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=1.664256Jy  sigma=16.142484
! Fit after  self-cal, rms=0.424896Jy  sigma=4.164153
! 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 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.5356 mas, bmaj=0.9394 mas, bpa=0.01221 degrees
! Estimated noise=1.26127 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 = 1.01183 Jy
! Component: 100  -  total flux cleaned = 1.24563 Jy
! Total flux subtracted in 100 components = 1.24563 Jy
! Clean residual min=-0.028617 max=0.084773 Jy/beam
! Clean residual mean=0.000030 rms=0.006850 Jy/beam
! Combined flux in latest and established models = 1.24563 Jy
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 5 components and 1.24563 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.288800Jy  sigma=2.685122
! Fit after  self-cal, rms=0.288866Jy  sigma=2.685000
! Inverting map 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0942628 Jy
! Component: 100  -  total flux cleaned = 0.152999 Jy
! Total flux subtracted in 100 components = 0.152999 Jy
! Clean residual min=-0.022980 max=0.031576 Jy/beam
! Clean residual mean=0.000012 rms=0.005133 Jy/beam
! Combined flux in latest and established models = 1.39863 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 18 components and 1.39863 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.250830Jy  sigma=2.236794
! Fit after  self-cal, rms=0.250839Jy  sigma=2.236256
! Inverting map 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0394813 Jy
! Component: 100  -  total flux cleaned = 0.0658768 Jy
! Total flux subtracted in 100 components = 0.0658768 Jy
! Clean residual min=-0.021212 max=0.024858 Jy/beam
! Clean residual mean=0.000005 rms=0.004745 Jy/beam
! Combined flux in latest and established models = 1.46451 Jy
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 20 components and 1.46451 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.243742Jy  sigma=2.149495
! Fit after  self-cal, rms=0.243737Jy  sigma=2.149267
! Inverting map 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
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 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6407 mas, bmaj=1.093 mas, bpa=-1.073 degrees
! Estimated noise=0.966848 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.016851 Jy
! Component: 100  -  total flux cleaned = 0.0255807 Jy
! Component: 150  -  total flux cleaned = 0.0289147 Jy
! Component: 200  -  total flux cleaned = 0.0292357 Jy
! Total flux subtracted in 200 components = 0.0292357 Jy
! Clean residual min=-0.020162 max=0.026716 Jy/beam
! Clean residual mean=0.000011 rms=0.005637 Jy/beam
! Combined flux in latest and established models = 1.49374 Jy
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 23 components and 1.49374 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.242447Jy  sigma=2.132610
! Fit after  self-cal, rms=0.242409Jy  sigma=2.132337
! Inverting map 
! Your choice of large map pixels excluded 7.23% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
print "============== Finished natural weighting clean =============="
! ============== Finished natural weighting clean ============== 
! 
restore
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6407 x 1.093 at -1.073 degrees (North through East)
! Clean map  min=-0.017731  max=1.3282 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:
!   HN        1.70     NL        1.05     SC        0.95     WZ        0.97     
!   ZC        0.92     
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.99     FD        1.06     KP        0.95     LA        0.94     
!   NL        0.97     NY        1.00     OV        1.10     PT        0.99     
!   SC        0.98     WZ        1.00     ZC        0.94     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.99     FD        0.98     KK        1.00     KP        0.96     
!   LA        1.22     MK        0.98     NL        0.95     NY        1.01     
!   OV        0.97     PT        0.96     SC        0.96     
! 
!  A total of 96 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        1.01     KK        1.02     MK        1.00     NY        1.02     
!   WZ        1.08     ZC        0.96     
! 
! 
! Correcting IF 2.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   HN        1.79     NL        1.07     SC        0.94     WZ        0.94     
!   ZC        0.96     
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00     FD        1.07     KP        0.96     LA        0.96     
!   NL        0.98     NY        0.98     OV        1.11     PT        1.01     
!   SC        0.96     WZ        1.01     ZC        0.95     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.98     FD        0.99     KK        1.01     KP        0.95     
!   LA        1.23     MK        1.00     NL        0.97     NY        1.03     
!   OV        0.97     PT        0.95     SC        0.94     
! 
!  A total of 96 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        1.03     KK        1.06     MK        0.98     NY        1.01     
!   WZ        1.11     ZC        0.94     
! 
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   HN        1.64     NL        1.02     SC        0.99     WZ        0.98     
!   ZC        0.95     
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.99     FD        1.08     KP        0.96     LA        0.97     
!   NL        0.97     NY        1.00     OV        1.12     PT        1.03     
!   SC        0.99     WZ        1.02     ZC        0.95     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.98     FD        0.96     KK        0.96     KP        0.95     
!   LA        1.24     MK        1.01     NL        0.97     NY        1.03     
!   OV        0.98     PT        0.97     SC        0.94     
! 
!  A total of 96 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        1.08     KK        1.01     MK        0.97     NY        1.02     
!   WZ        1.14     ZC        0.91     
! 
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   HN        1.99     NL        1.03     SC        0.99     WZ        0.97     
!   ZC        0.98     
! 
!  A total of 12 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.98     FD        1.05     KP        0.96     LA        0.95     
!   NL        0.96     NY        1.03     OV        1.12     PT        1.01     
!   SC        0.99     WZ        1.02     ZC        0.96     
! 
!  A total of 3 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.99     FD        0.97     KK        0.99     KP        0.96     
!   LA        1.26     MK        0.99     NL        0.98     NY        1.12     
!   OV        0.97     PT        0.96     SC        0.98     
! 
!  A total of 96 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        1.07     KK        1.05     MK        0.99     NY        1.01     
!   WZ        1.10     ZC        0.92     
! 
! 
! Fit before self-cal, rms=0.242409Jy  sigma=2.132337
! Fit after  self-cal, rms=0.218148Jy  sigma=1.785433
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6209 mas, bmaj=1.102 mas, bpa=-1.246 degrees
! Estimated noise=1.00892 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.00779464 Jy
! Component: 100  -  total flux cleaned = 0.00922102 Jy
! Component: 150  -  total flux cleaned = 0.0092223 Jy
! Component: 200  -  total flux cleaned = 0.00935841 Jy
! Total flux subtracted in 200 components = 0.00935841 Jy
! Clean residual min=-0.013737 max=0.010601 Jy/beam
! Clean residual mean=0.000019 rms=0.002873 Jy/beam
! Combined flux in latest and established models = 1.5031 Jy
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 27 components and 1.5031 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.217775Jy  sigma=1.782406
! Fit after  self-cal, rms=0.218114Jy  sigma=1.781547
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
selfcal true, true, soltime1
! Performing amp+phase self-cal over 120 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.218114Jy  sigma=1.781547
! Fit after  self-cal, rms=0.213859Jy  sigma=1.724122
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6207 mas, bmaj=1.1 mas, bpa=-1.418 degrees
! Estimated noise=1.01109 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00636612 Jy
! Component: 100  -  total flux cleaned = 0.00666254 Jy
! Component: 150  -  total flux cleaned = 0.00678955 Jy
! Component: 200  -  total flux cleaned = 0.00679037 Jy
! Total flux subtracted in 200 components = 0.00679037 Jy
! Clean residual min=-0.009743 max=0.008242 Jy/beam
! Clean residual mean=0.000038 rms=0.002044 Jy/beam
! Combined flux in latest and established models = 1.50989 Jy
selfcal
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 28 components and 1.50989 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.213598Jy  sigma=1.722092
! Fit after  self-cal, rms=0.213618Jy  sigma=1.722001
flux_cutoff = imstat(rms) * dynam
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
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.000305248 Jy
! Component: 100  -  total flux cleaned = 9.34791e-05 Jy
! Component: 150  -  total flux cleaned = -9.29856e-07 Jy
! Component: 200  -  total flux cleaned = -1.31654e-06 Jy
! Total flux subtracted in 200 components = -1.31654e-06 Jy
! Clean residual min=-0.009669 max=0.008054 Jy/beam
! Clean residual mean=0.000038 rms=0.002039 Jy/beam
! Combined flux in latest and established models = 1.50989 Jy
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 29 components and 1.50989 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.213486Jy  sigma=1.721015
! Fit after  self-cal, rms=0.213482Jy  sigma=1.720970
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.213482Jy  sigma=1.720970
! Fit after  self-cal, rms=0.211443Jy  sigma=1.694999
clean clean_niter,clean_gain
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.621 mas, bmaj=1.098 mas, bpa=-1.535 degrees
! Estimated noise=1.01344 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00530536 Jy
! Component: 100  -  total flux cleaned = 0.00569353 Jy
! Component: 150  -  total flux cleaned = 0.00559397 Jy
! Component: 200  -  total flux cleaned = 0.00541309 Jy
! Total flux subtracted in 200 components = 0.00541309 Jy
! Clean residual min=-0.007678 max=0.007157 Jy/beam
! Clean residual mean=0.000047 rms=0.001714 Jy/beam
! Combined flux in latest and established models = 1.5153 Jy
selfcal
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 30 components and 1.5153 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211311Jy  sigma=1.694243
! Fit after  self-cal, rms=0.211330Jy  sigma=1.694114
flux_cutoff = imstat(rms) * dynam
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
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.000189873 Jy
! Component: 100  -  total flux cleaned = -0.000355071 Jy
! Component: 150  -  total flux cleaned = -0.000504929 Jy
! Component: 200  -  total flux cleaned = -0.000647547 Jy
! Total flux subtracted in 200 components = -0.000647547 Jy
! Clean residual min=-0.007597 max=0.007196 Jy/beam
! Clean residual mean=0.000047 rms=0.001710 Jy/beam
! Combined flux in latest and established models = 1.51466 Jy
! Performing phase self-cal
! Adding 11 model components to the UV plane model.
! The established model now contains 31 components and 1.51466 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211286Jy  sigma=1.693838
! Fit after  self-cal, rms=0.211289Jy  sigma=1.693815
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 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.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211289Jy  sigma=1.693815
! Fit after  self-cal, rms=0.211643Jy  sigma=1.692650
clean clean_niter,clean_gain
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=2.483 mas, bmaj=2.65 mas, bpa=-53.25 degrees
! Estimated noise=1.54522 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00617773 Jy
! Component: 100  -  total flux cleaned = 0.00813405 Jy
! Component: 150  -  total flux cleaned = 0.00865466 Jy
! Component: 200  -  total flux cleaned = 0.00870042 Jy
! Total flux subtracted in 200 components = 0.00870042 Jy
! Clean residual min=-0.010503 max=0.009929 Jy/beam
! Clean residual mean=0.000144 rms=0.002980 Jy/beam
! Combined flux in latest and established models = 1.52336 Jy
selfcal
! Performing phase self-cal
! Adding 11 model components to the UV plane model.
! The established model now contains 32 components and 1.52336 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211648Jy  sigma=1.692842
! Fit after  self-cal, rms=0.211613Jy  sigma=1.692695
flux_cutoff = imstat(rms) * dynam
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
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.000219971 Jy
! Component: 100  -  total flux cleaned = 0.000221496 Jy
! Component: 150  -  total flux cleaned = 0.000272096 Jy
! Component: 200  -  total flux cleaned = 0.000318031 Jy
! Total flux subtracted in 200 components = 0.000318031 Jy
! Clean residual min=-0.010520 max=0.009972 Jy/beam
! Clean residual mean=0.000143 rms=0.002981 Jy/beam
! Combined flux in latest and established models = 1.52367 Jy
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 33 components and 1.52367 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211651Jy  sigma=1.693005
! Fit after  self-cal, rms=0.211641Jy  sigma=1.692973
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
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 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6212 mas, bmaj=1.096 mas, bpa=-1.591 degrees
! Estimated noise=1.01557 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6212 x 1.096 at -1.591 degrees (North through East)
! Clean map  min=-0.0064874  max=1.3343 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.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.211641Jy  sigma=1.692973
! Fit after  self-cal, rms=0.189440Jy  sigma=1.488322
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6208 mas, bmaj=1.096 mas, bpa=-1.637 degrees
! Estimated noise=1.0157 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00129113 Jy
! Component: 100  -  total flux cleaned = 0.00130122 Jy
! Component: 150  -  total flux cleaned = 0.00120544 Jy
! Component: 200  -  total flux cleaned = 0.00103364 Jy
! Total flux subtracted in 200 components = 0.00103364 Jy
! Clean residual min=-0.007456 max=0.006361 Jy/beam
! Clean residual mean=0.000038 rms=0.001589 Jy/beam
! Combined flux in latest and established models = 1.52471 Jy
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 33 components and 1.52471 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.189344Jy  sigma=1.487700
! Fit after  self-cal, rms=0.189350Jy  sigma=1.487652
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 8.42501e-05 Jy
! Component: 100  -  total flux cleaned = -5.54221e-05 Jy
! Component: 150  -  total flux cleaned = -0.000182623 Jy
! Component: 200  -  total flux cleaned = -0.000303374 Jy
! Total flux subtracted in 200 components = -0.000303374 Jy
! Clean residual min=-0.007424 max=0.006389 Jy/beam
! Clean residual mean=0.000039 rms=0.001586 Jy/beam
! Combined flux in latest and established models = 1.5244 Jy
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 33 components and 1.5244 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.189330Jy  sigma=1.487637
! Fit after  self-cal, rms=0.189334Jy  sigma=1.487623
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.189334Jy  sigma=1.487623
! Fit after  self-cal, rms=0.189678Jy  sigma=1.487046
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 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.5225 mas, bmaj=0.9376 mas, bpa=-0.2762 degrees
! Estimated noise=1.31323 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.01665 Jy
! Component: 100  -  total flux cleaned = 1.25458 Jy
! Component: 150  -  total flux cleaned = 1.35379 Jy
! Component: 200  -  total flux cleaned = 1.41203 Jy
! Total flux subtracted in 200 components = 1.41203 Jy
! Clean residual min=-0.012272 max=0.030745 Jy/beam
! Clean residual mean=0.000030 rms=0.002939 Jy/beam
! Combined flux in latest and established models = 1.41203 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 1.41203 Jy
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6213 mas, bmaj=1.095 mas, bpa=-1.641 degrees
! Estimated noise=1.01753 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.0532732 Jy
! Component: 100  -  total flux cleaned = 0.0816298 Jy
! Component: 150  -  total flux cleaned = 0.0971146 Jy
! Component: 200  -  total flux cleaned = 0.105789 Jy
! Total flux subtracted in 200 components = 0.105789 Jy
! Clean residual min=-0.007548 max=0.007337 Jy/beam
! Clean residual mean=0.000043 rms=0.001633 Jy/beam
! Combined flux in latest and established models = 1.51782 Jy
! Adding 23 model components to the UV plane model.
! The established model now contains 28 components and 1.51782 Jy
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=2.481 mas, bmaj=2.649 mas, bpa=-53.69 degrees
! Estimated noise=1.5513 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00905398 Jy
! Component: 100  -  total flux cleaned = 0.0117225 Jy
! Component: 150  -  total flux cleaned = 0.0125085 Jy
! Component: 200  -  total flux cleaned = 0.0125622 Jy
! Total flux subtracted in 200 components = 0.0125622 Jy
! Clean residual min=-0.009690 max=0.009395 Jy/beam
! Clean residual mean=0.000122 rms=0.002627 Jy/beam
! Combined flux in latest and established models = 1.53038 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 = 2.57712e-05 Jy
! Component: 100  -  total flux cleaned = 7.61913e-05 Jy
! Component: 150  -  total flux cleaned = 0.000124435 Jy
! Component: 200  -  total flux cleaned = 0.000169898 Jy
! Total flux subtracted in 200 components = 0.000169898 Jy
! Clean residual min=-0.009685 max=0.009381 Jy/beam
! Clean residual mean=0.000121 rms=0.002626 Jy/beam
! Combined flux in latest and established models = 1.53055 Jy
! Adding 10 model components to the UV plane model.
! The established model now contains 33 components and 1.53055 Jy
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! Estimated beam: bmin=0.6213 mas, bmaj=1.095 mas, bpa=-1.641 degrees
! Estimated noise=1.01753 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00280182 Jy
! Component: 100  -  total flux cleaned = -0.00307055 Jy
! Component: 150  -  total flux cleaned = -0.00318359 Jy
! Component: 200  -  total flux cleaned = -0.00338231 Jy
! Total flux subtracted in 200 components = -0.00338231 Jy
! Clean residual min=-0.007685 max=0.006664 Jy/beam
! Clean residual mean=0.000039 rms=0.001594 Jy/beam
! Combined flux in latest and established models = 1.52717 Jy
selfcal
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 34 components and 1.52717 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.190448Jy  sigma=1.494388
! Fit after  self-cal, rms=0.190412Jy  sigma=1.494248
wmodel J1753+2848_X_map.mod
! Writing 34 model components to file: J1753+2848_X_map.mod
wobs   J1753+2848_X_uvs.fits
! Writing UV FITS file: J1753+2848_X_uvs.fits
wwins  J1753+2848_X_map.win
! wwins: Wrote 1 windows to J1753+2848_X_map.win
x = (field_size-8) * field_cell / 4
addwin -x,x,-x,x
clean (field_size*4),0.01
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! clean: niter=4096  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.00250324 Jy
! Component: 100  -  total flux cleaned = -0.00374986 Jy
! Component: 150  -  total flux cleaned = -0.00458633 Jy
! Component: 200  -  total flux cleaned = -0.0052831 Jy
! Component: 250  -  total flux cleaned = -0.00606788 Jy
! Component: 300  -  total flux cleaned = -0.00563069 Jy
! Component: 350  -  total flux cleaned = -0.00542003 Jy
! Component: 400  -  total flux cleaned = -0.00542175 Jy
! Component: 450  -  total flux cleaned = -0.00510928 Jy
! Component: 500  -  total flux cleaned = -0.00490153 Jy
! Component: 550  -  total flux cleaned = -0.00419256 Jy
! Component: 600  -  total flux cleaned = -0.00359284 Jy
! Component: 650  -  total flux cleaned = -0.00349394 Jy
! Component: 700  -  total flux cleaned = -0.0030035 Jy
! Component: 750  -  total flux cleaned = -0.00261669 Jy
! Component: 800  -  total flux cleaned = -0.00175574 Jy
! Component: 850  -  total flux cleaned = -0.00137762 Jy
! Component: 900  -  total flux cleaned = -0.000812597 Jy
! Component: 950  -  total flux cleaned = -0.000439133 Jy
! Component: 1000  -  total flux cleaned = -7.02106e-05 Jy
! Component: 1050  -  total flux cleaned = 0.000297025 Jy
! Component: 1100  -  total flux cleaned = 0.000930888 Jy
! Component: 1150  -  total flux cleaned = 0.00075251 Jy
! Component: 1200  -  total flux cleaned = 0.00128807 Jy
! Component: 1250  -  total flux cleaned = 0.00173061 Jy
! Component: 1300  -  total flux cleaned = 0.0015548 Jy
! Component: 1350  -  total flux cleaned = 0.00225203 Jy
! Component: 1400  -  total flux cleaned = 0.00233808 Jy
! Component: 1450  -  total flux cleaned = 0.00216591 Jy
! Component: 1500  -  total flux cleaned = 0.00250777 Jy
! Component: 1550  -  total flux cleaned = 0.00242419 Jy
! Component: 1600  -  total flux cleaned = 0.00267669 Jy
! Component: 1650  -  total flux cleaned = 0.00284474 Jy
! Component: 1700  -  total flux cleaned = 0.00301091 Jy
! Component: 1750  -  total flux cleaned = 0.00359009 Jy
! Component: 1800  -  total flux cleaned = 0.00375429 Jy
! Component: 1850  -  total flux cleaned = 0.00399779 Jy
! Component: 1900  -  total flux cleaned = 0.00464562 Jy
! Component: 1950  -  total flux cleaned = 0.00456423 Jy
! Component: 2000  -  total flux cleaned = 0.00496414 Jy
! Component: 2050  -  total flux cleaned = 0.00544221 Jy
! Component: 2100  -  total flux cleaned = 0.00536262 Jy
! Component: 2150  -  total flux cleaned = 0.0056773 Jy
! Component: 2200  -  total flux cleaned = 0.00630414 Jy
! Component: 2250  -  total flux cleaned = 0.00653719 Jy
! Component: 2300  -  total flux cleaned = 0.0069244 Jy
! Component: 2350  -  total flux cleaned = 0.00684715 Jy
! Component: 2400  -  total flux cleaned = 0.00723002 Jy
! Component: 2450  -  total flux cleaned = 0.00753498 Jy
! Component: 2500  -  total flux cleaned = 0.00806616 Jy
! Component: 2550  -  total flux cleaned = 0.008367 Jy
! Component: 2600  -  total flux cleaned = 0.00874262 Jy
! Component: 2650  -  total flux cleaned = 0.00874211 Jy
! Component: 2700  -  total flux cleaned = 0.00926284 Jy
! Component: 2750  -  total flux cleaned = 0.00926303 Jy
! Component: 2800  -  total flux cleaned = 0.00977833 Jy
! Component: 2850  -  total flux cleaned = 0.0100715 Jy
! Component: 2900  -  total flux cleaned = 0.0102175 Jy
! Component: 2950  -  total flux cleaned = 0.0104359 Jy
! Component: 3000  -  total flux cleaned = 0.0109416 Jy
! Component: 3050  -  total flux cleaned = 0.0112293 Jy
! Component: 3100  -  total flux cleaned = 0.0113722 Jy
! Component: 3150  -  total flux cleaned = 0.0117287 Jy
! Component: 3200  -  total flux cleaned = 0.0117997 Jy
! Component: 3250  -  total flux cleaned = 0.0122237 Jy
! Component: 3300  -  total flux cleaned = 0.0120828 Jy
! Component: 3350  -  total flux cleaned = 0.0126431 Jy
! Component: 3400  -  total flux cleaned = 0.0127821 Jy
! Component: 3450  -  total flux cleaned = 0.0131295 Jy
! Component: 3500  -  total flux cleaned = 0.0133373 Jy
! Component: 3550  -  total flux cleaned = 0.0136815 Jy
! Component: 3600  -  total flux cleaned = 0.0138876 Jy
! Component: 3650  -  total flux cleaned = 0.0146386 Jy
! Component: 3700  -  total flux cleaned = 0.015046 Jy
! Component: 3750  -  total flux cleaned = 0.0155203 Jy
! Component: 3800  -  total flux cleaned = 0.0153853 Jy
! Component: 3850  -  total flux cleaned = 0.0156536 Jy
! Component: 3900  -  total flux cleaned = 0.0160555 Jy
! Component: 3950  -  total flux cleaned = 0.0165888 Jy
! Component: 4000  -  total flux cleaned = 0.0169206 Jy
! Component: 4050  -  total flux cleaned = 0.0173832 Jy
! Total flux subtracted in 4096 components = 0.0177782 Jy
! Clean residual min=-0.004164 max=0.003730 Jy/beam
! Clean residual mean=0.000025 rms=0.001092 Jy/beam
! Combined flux in latest and established models = 1.54495 Jy
keep
! Adding 366 model components to the UV plane model.
! The established model now contains 400 components and 1.54495 Jy
![Exited script file: /opt64/bin/pima_mupet_01.dfm]
wmap   /scr/pima/rdv76_uvs/J1753+2848_X_map.fits
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6213 x 1.095 at -1.641 degrees (North through East)
! Clean map  min=-0.0057618  max=1.3374 Jy/beam
! Writing clean map to FITS file: /scr/pima/rdv76_uvs/J1753+2848_X_map.fits
wmodel /scr/pima/rdv76_uvs/J1753+2848_X_map.mod
! Writing 400 model components to file: /scr/pima/rdv76_uvs/J1753+2848_X_map.mod
wobs   /scr/pima/rdv76_uvs/J1753+2848_X_uvs.fits
! Writing UV FITS file: /scr/pima/rdv76_uvs/J1753+2848_X_uvs.fits
wwins  /scr/pima/rdv76_uvs/J1753+2848_X_map.win
! wwins: Wrote 2 windows to /scr/pima/rdv76_uvs/J1753+2848_X_map.win
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is", peak_flux
! The peak flux is 1.33742 
invert
! Inverting map 
! Your choice of large map pixels excluded 7.38% of the data.
!  The x-axis pixel size should ideally be below 0.1641 milli-arcsec
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00109016 
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.244536 0.244536 0.489072 0.978143 1.95629 3.91257 7.82514 15.6503 31.3006 62.6012
print "clev=      ", int(3*image_rms*10000)/10000
! clev=       0.0032 
print "peak=      ", int(1000*peak_flux)/1000
! peak=       1.337 
print "dyn_range= ", int(1000*peak_flux/image_rms)/1000
! dyn_range=  1226.81 
quit
! Quitting program
! Log file /scr/pima/rdv76_uvs/J1753+2848_X_dfm.log closed on Fri Jun 15 03:56:22 2018