! Started logfile: difmap.log_137 on Sat Feb  1 09:07:39 2020
![@/tmp/difmap__6860.1]
! Obeying initialization file: /tmp/difmap__6860.1
![@/home/lpetrov/config/difmap.cnf]
device /XW
! Attempting to open device: '/XW'
vflags="fb3"
integer clean_niter
clean_niter = 100
float   clean_gain
clean_gain = 0.03
float   dynam
dynam = 6.0
float   soltime1
soltime1 = 800000
float   soltime2
soltime2 = 15
float   soltime3
soltime3 = 1
float   thresh
thresh = 0.5
float   win_mult
win_mult = 1.8
float   taper_size
taper_size = 70.
integer field_size
field_size = 1024
float   field_cell
field_cell = 0.1
float old_peak
float new_peak
float flux_cutoff
![Exited script file: /home/lpetrov/config/difmap.cnf]
![@/tmp/difmap__6860.2]
![Exited script file: /tmp/difmap__6860.2]
obs J1007+1616_C_uva.fits
! Reading UV FITS file: J1007+1616_C_uva.fits
! AN table 1: 16 integrations on 10 of 10 possible baselines.
! Apparent sampling: 1 visibilities/baseline/integration-bin.
! Found source: J1007+1616
! 
! There are 8 IFs, and a total of 8 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1    4.128e+09      3.2e+07          1      3.2e+07
!  02        2     4.16e+09      3.2e+07          1      3.2e+07
!  03        3    4.192e+09      3.2e+07          1      3.2e+07
!  04        4    4.224e+09      3.2e+07          1      3.2e+07
!  05        5    4.416e+09      3.2e+07          1      3.2e+07
!  06        6    4.512e+09      3.2e+07          1      3.2e+07
!  07        7    4.544e+09      3.2e+07          1      3.2e+07
!  08        8    4.576e+09      3.2e+07          1      3.2e+07
! 
! Polarization(s): RR
! 
! Read 94 lines of history.
! 
! Reading 1280 visibilities.
![Exited script file: /tmp/difmap__6860.1]
sele I
! Polarization I is unavailable.
! Selecting polarization: RR,  channels: 1..8
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
! Reading IF 5 channels: 5..5
! Reading IF 6 channels: 6..6
! Reading IF 7 channels: 7..7
! Reading IF 8 channels: 8..8
field_size=1024
field_cell=0.2
taper_size = 50.
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.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=1.000562Jy  sigma=23.989338
! Fit after  self-cal, rms=0.980139Jy  sigma=23.499900
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
radplot
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 7 buffered edits.
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.
clrmod true,true
clean  200, 0.02
! Inverting map and beam 
! Estimated beam: bmin=7.787 mas, bmaj=16.37 mas, bpa=43.05 degrees
! Estimated noise=0.902601 mJy/beam.
! clean: niter=200  gain=0.02  cutoff=0
! Component: 050  -  total flux cleaned = 0.0125922 Jy
! Component: 100  -  total flux cleaned = 0.0171794 Jy
! Component: 150  -  total flux cleaned = 0.0173632 Jy
! Component: 200  -  total flux cleaned = 0.0171011 Jy
! Total flux subtracted in 200 components = 0.0171011 Jy
! Clean residual min=-0.001495 max=0.001499 Jy/beam
! Clean residual mean=-0.000014 rms=0.000735 Jy/beam
! Combined flux in latest and established models = 0.0171011 Jy
mapl  map,true
! 
! Move the cursor into the plot window and press 'H' for help
uvaver 555
! Averaging into 555 second bins.
! Selecting polarization: RR,  channels: 1..8
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
! Reading IF 5 channels: 5..5
! Reading IF 6 channels: 6..6
! Reading IF 7 channels: 7..7
! Reading IF 8 channels: 8..8
proj
! Adding 56 model components to the UV plane model.
! The established model now contains 56 components and 0.0171011 Jy
! Using default options string "m3"
! Move the cursor into the plot window and press 'H' for help
! Applying 5 buffered edits.
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.
clrmod true,true
! clrmod: Cleared the established and tentative models.
clean  200, 0.02
! Inverting map and beam 
! Estimated beam: bmin=7.777 mas, bmaj=16.34 mas, bpa=43.05 degrees
! Estimated noise=0.92468 mJy/beam.
! clean: niter=200  gain=0.02  cutoff=0
! Component: 050  -  total flux cleaned = 0.0119619 Jy
! Component: 100  -  total flux cleaned = 0.0163199 Jy
! Component: 150  -  total flux cleaned = 0.0181036 Jy
! Component: 200  -  total flux cleaned = 0.0192129 Jy
! Total flux subtracted in 200 components = 0.0192129 Jy
! Clean residual min=-0.002018 max=0.001267 Jy/beam
! Clean residual mean=-0.000014 rms=0.000629 Jy/beam
! Combined flux in latest and established models = 0.0192129 Jy
mapl  map,true
! 
! Move the cursor into the plot window and press 'H' for help
clrmod true,true
! clrmod: Cleared the established and tentative models.
selfcal
! Performing phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.019801Jy  sigma=1.914104
! Fit after  self-cal, rms=0.019801Jy  sigma=1.914104
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  800, 0.01
! Inverting map and beam 
! Estimated beam: bmin=7.777 mas, bmaj=16.34 mas, bpa=43.05 degrees
! Estimated noise=0.92468 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00743076 Jy
! Component: 100  -  total flux cleaned = 0.0119269 Jy
! Component: 150  -  total flux cleaned = 0.0146477 Jy
! Component: 200  -  total flux cleaned = 0.0162943 Jy
! Component: 250  -  total flux cleaned = 0.0173185 Jy
! Component: 300  -  total flux cleaned = 0.0180828 Jy
! Component: 350  -  total flux cleaned = 0.0186943 Jy
! Component: 400  -  total flux cleaned = 0.0191943 Jy
! Component: 450  -  total flux cleaned = 0.0196061 Jy
! Component: 500  -  total flux cleaned = 0.0199454 Jy
! Component: 550  -  total flux cleaned = 0.0202252 Jy
! Component: 600  -  total flux cleaned = 0.0204561 Jy
! Component: 650  -  total flux cleaned = 0.0206468 Jy
! Component: 700  -  total flux cleaned = 0.0208044 Jy
! Component: 750  -  total flux cleaned = 0.0208907 Jy
! Component: 800  -  total flux cleaned = 0.0209211 Jy
! Total flux subtracted in 800 components = 0.0209211 Jy
! Clean residual min=-0.001841 max=0.001308 Jy/beam
! Clean residual mean=-0.000009 rms=0.000569 Jy/beam
! Combined flux in latest and established models = 0.0209212 Jy
uvw 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean 1200, 0.01
! Adding 53 model components to the UV plane model.
! The established model now contains 53 components and 0.0209212 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.038 mas, bmaj=15.04 mas, bpa=37.68 degrees
! Estimated noise=0.851747 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 2.32743e-05 Jy
! Component: 100  -  total flux cleaned = 7.5684e-05 Jy
! Component: 150  -  total flux cleaned = 0.000125129 Jy
! Component: 200  -  total flux cleaned = 0.000165358 Jy
! Component: 250  -  total flux cleaned = 0.000210265 Jy
! Component: 300  -  total flux cleaned = 0.000253354 Jy
! Component: 350  -  total flux cleaned = 0.000294596 Jy
! Component: 400  -  total flux cleaned = 0.000334206 Jy
! Component: 450  -  total flux cleaned = 0.000372327 Jy
! Component: 500  -  total flux cleaned = 0.00040401 Jy
! Component: 550  -  total flux cleaned = 0.000439888 Jy
! Component: 600  -  total flux cleaned = 0.000474645 Jy
! Component: 650  -  total flux cleaned = 0.000508321 Jy
! Component: 700  -  total flux cleaned = 0.000536348 Jy
! Component: 750  -  total flux cleaned = 0.000568108 Jy
! Component: 800  -  total flux cleaned = 0.000598888 Jy
! Component: 850  -  total flux cleaned = 0.00062873 Jy
! Component: 900  -  total flux cleaned = 0.000657635 Jy
! Component: 950  -  total flux cleaned = 0.000685654 Jy
! Component: 1000  -  total flux cleaned = 0.000712828 Jy
! Component: 1050  -  total flux cleaned = 0.000735437 Jy
! Component: 1100  -  total flux cleaned = 0.000761079 Jy
! Component: 1150  -  total flux cleaned = 0.000785929 Jy
! Component: 1200  -  total flux cleaned = 0.000810018 Jy
! Total flux subtracted in 1200 components = 0.000810018 Jy
! Clean residual min=-0.001599 max=0.001134 Jy/beam
! Clean residual mean=-0.000015 rms=0.000457 Jy/beam
! Combined flux in latest and established models = 0.0217312 Jy
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapl map,true
! 
! Move the cursor into the plot window and press 'H' for help
clrmod true,true
! clrmod: Cleared the established and tentative models.
selfcal
! Performing phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.019801Jy  sigma=1.914104
! Fit after  self-cal, rms=0.019801Jy  sigma=1.914104
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  800, 0.01
! Inverting map and beam 
! Estimated beam: bmin=7.777 mas, bmaj=16.34 mas, bpa=43.05 degrees
! Estimated noise=0.92468 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00743076 Jy
! Component: 100  -  total flux cleaned = 0.0119269 Jy
! Component: 150  -  total flux cleaned = 0.0146477 Jy
! Component: 200  -  total flux cleaned = 0.0162943 Jy
! Component: 250  -  total flux cleaned = 0.0173185 Jy
! Component: 300  -  total flux cleaned = 0.0180846 Jy
! Component: 350  -  total flux cleaned = 0.0187027 Jy
! Component: 400  -  total flux cleaned = 0.0192117 Jy
! Component: 450  -  total flux cleaned = 0.0196328 Jy
! Component: 500  -  total flux cleaned = 0.0199815 Jy
! Component: 550  -  total flux cleaned = 0.0202704 Jy
! Component: 600  -  total flux cleaned = 0.02051 Jy
! Component: 650  -  total flux cleaned = 0.0207088 Jy
! Component: 700  -  total flux cleaned = 0.0208738 Jy
! Component: 750  -  total flux cleaned = 0.0210108 Jy
! Component: 800  -  total flux cleaned = 0.0211246 Jy
! Total flux subtracted in 800 components = 0.0211246 Jy
! Clean residual min=-0.001858 max=0.001336 Jy/beam
! Clean residual mean=-0.000011 rms=0.000572 Jy/beam
! Combined flux in latest and established models = 0.0211247 Jy
uvw 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean 1200, 0.01
! Adding 52 model components to the UV plane model.
! The established model now contains 52 components and 0.0211247 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.038 mas, bmaj=15.04 mas, bpa=37.68 degrees
! Estimated noise=0.851747 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.000198141 Jy
! Component: 100  -  total flux cleaned = 0.000248919 Jy
! Component: 150  -  total flux cleaned = 0.000291401 Jy
! Component: 200  -  total flux cleaned = 0.000337372 Jy
! Component: 250  -  total flux cleaned = 0.000375476 Jy
! Component: 300  -  total flux cleaned = 0.000411786 Jy
! Component: 350  -  total flux cleaned = 0.000451399 Jy
! Component: 400  -  total flux cleaned = 0.000484219 Jy
! Component: 450  -  total flux cleaned = 0.000519792 Jy
! Component: 500  -  total flux cleaned = 0.000549413 Jy
! Component: 550  -  total flux cleaned = 0.000577714 Jy
! Component: 600  -  total flux cleaned = 0.000608355 Jy
! Component: 650  -  total flux cleaned = 0.000633734 Jy
! Component: 700  -  total flux cleaned = 0.000658064 Jy
! Component: 750  -  total flux cleaned = 0.000684449 Jy
! Component: 800  -  total flux cleaned = 0.000706329 Jy
! Component: 850  -  total flux cleaned = 0.000730012 Jy
! Component: 900  -  total flux cleaned = 0.000749845 Jy
! Component: 950  -  total flux cleaned = 0.000768698 Jy
! Component: 1000  -  total flux cleaned = 0.000789112 Jy
! Component: 1050  -  total flux cleaned = 0.000806067 Jy
! Component: 1100  -  total flux cleaned = 0.000824572 Jy
! Component: 1150  -  total flux cleaned = 0.000839907 Jy
! Component: 1200  -  total flux cleaned = 0.000854485 Jy
! Total flux subtracted in 1200 components = 0.000854485 Jy
! Clean residual min=-0.001554 max=0.001191 Jy/beam
! Clean residual mean=-0.000016 rms=0.000456 Jy/beam
! Combined flux in latest and established models = 0.0219791 Jy
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapl map,true
! 
! Move the cursor into the plot window and press 'H' for help
uncalib
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 54 components and 0.0219791 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.006855Jy  sigma=0.653615
! Fit after  self-cal, rms=0.006523Jy  sigma=0.623170
selfcal true,true,1.e6
! Performing amp+phase self-cal over 1e+06 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.006523Jy  sigma=0.623170
! Fit after  self-cal, rms=0.006255Jy  sigma=0.508737
clrmod true,true
! clrmod: Cleared the established and tentative models.
selfcal
! Performing phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.019943Jy  sigma=1.914104
! Fit after  self-cal, rms=0.019943Jy  sigma=1.914104
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  800, 0.01
! Inverting map and beam 
! Estimated beam: bmin=7.966 mas, bmaj=16.59 mas, bpa=42.85 degrees
! Estimated noise=0.93716 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00757414 Jy
! Component: 100  -  total flux cleaned = 0.0121573 Jy
! Component: 150  -  total flux cleaned = 0.0149306 Jy
! Component: 200  -  total flux cleaned = 0.0166089 Jy
! Component: 250  -  total flux cleaned = 0.0176274 Jy
! Component: 300  -  total flux cleaned = 0.0183372 Jy
! Component: 350  -  total flux cleaned = 0.0189021 Jy
! Component: 400  -  total flux cleaned = 0.0193608 Jy
! Component: 450  -  total flux cleaned = 0.0197364 Jy
! Component: 500  -  total flux cleaned = 0.0200444 Jy
! Component: 550  -  total flux cleaned = 0.0202977 Jy
! Component: 600  -  total flux cleaned = 0.0205068 Jy
! Component: 650  -  total flux cleaned = 0.0206797 Jy
! Component: 700  -  total flux cleaned = 0.0208229 Jy
! Component: 750  -  total flux cleaned = 0.0209416 Jy
! Component: 800  -  total flux cleaned = 0.0210363 Jy
! Total flux subtracted in 800 components = 0.0210363 Jy
! Clean residual min=-0.001041 max=0.000791 Jy/beam
! Clean residual mean=0.000003 rms=0.000304 Jy/beam
! Combined flux in latest and established models = 0.0210363 Jy
uvw 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean 1200, 0.01
! Adding 51 model components to the UV plane model.
! The established model now contains 51 components and 0.0210363 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.438 mas, bmaj=15.72 mas, bpa=37.84 degrees
! Estimated noise=0.841818 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -9.30269e-05 Jy
! Component: 100  -  total flux cleaned = -9.30443e-05 Jy
! Component: 150  -  total flux cleaned = -9.30785e-05 Jy
! Component: 200  -  total flux cleaned = -8.10895e-05 Jy
! Component: 250  -  total flux cleaned = -6.36594e-05 Jy
! Component: 300  -  total flux cleaned = -4.41135e-05 Jy
! Component: 350  -  total flux cleaned = -2.54055e-05 Jy
! Component: 400  -  total flux cleaned = -7.47011e-06 Jy
! Component: 450  -  total flux cleaned = 9.7136e-06 Jy
! Component: 500  -  total flux cleaned = 2.6178e-05 Jy
! Component: 550  -  total flux cleaned = 3.97242e-05 Jy
! Component: 600  -  total flux cleaned = 5.4885e-05 Jy
! Component: 650  -  total flux cleaned = 6.9432e-05 Jy
! Component: 700  -  total flux cleaned = 8.33906e-05 Jy
! Component: 750  -  total flux cleaned = 9.67599e-05 Jy
! Component: 800  -  total flux cleaned = 0.00010955 Jy
! Component: 850  -  total flux cleaned = 0.000120085 Jy
! Component: 900  -  total flux cleaned = 0.000131892 Jy
! Component: 950  -  total flux cleaned = 0.000143225 Jy
! Component: 1000  -  total flux cleaned = 0.000154082 Jy
! Component: 1050  -  total flux cleaned = 0.000164482 Jy
! Component: 1100  -  total flux cleaned = 0.00017305 Jy
! Component: 1150  -  total flux cleaned = 0.000182641 Jy
! Component: 1200  -  total flux cleaned = 0.00019184 Jy
! Total flux subtracted in 1200 components = 0.00019184 Jy
! Clean residual min=-0.000827 max=0.000744 Jy/beam
! Clean residual mean=-0.000000 rms=0.000267 Jy/beam
! Combined flux in latest and established models = 0.0212281 Jy
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapl map,true
! 
! Move the cursor into the plot window and press 'H' for help
proj
! Adding 5 model components to the UV plane model.
! The established model now contains 52 components and 0.0212281 Jy
! Using default options string "m3"
! Move the cursor into the plot window and press 'H' for help
clrmod true,true
! clrmod: Cleared the established and tentative models.
selfcal
! Performing phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.019943Jy  sigma=1.914104
! Fit after  self-cal, rms=0.019943Jy  sigma=1.914104
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  800, 0.01
! Inverting map and beam 
! Estimated beam: bmin=7.966 mas, bmaj=16.59 mas, bpa=42.85 degrees
! Estimated noise=0.93716 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00757414 Jy
! Component: 100  -  total flux cleaned = 0.0121573 Jy
! Component: 150  -  total flux cleaned = 0.0149306 Jy
! Component: 200  -  total flux cleaned = 0.0166089 Jy
! Component: 250  -  total flux cleaned = 0.0176274 Jy
! Component: 300  -  total flux cleaned = 0.0183372 Jy
! Component: 350  -  total flux cleaned = 0.0189021 Jy
! Component: 400  -  total flux cleaned = 0.0193608 Jy
! Component: 450  -  total flux cleaned = 0.0197364 Jy
! Component: 500  -  total flux cleaned = 0.0200444 Jy
! Component: 550  -  total flux cleaned = 0.0202977 Jy
! Component: 600  -  total flux cleaned = 0.0205068 Jy
! Component: 650  -  total flux cleaned = 0.0206797 Jy
! Component: 700  -  total flux cleaned = 0.0208229 Jy
! Component: 750  -  total flux cleaned = 0.0209416 Jy
! Component: 800  -  total flux cleaned = 0.0210363 Jy
! Total flux subtracted in 800 components = 0.0210363 Jy
! Clean residual min=-0.001041 max=0.000791 Jy/beam
! Clean residual mean=0.000003 rms=0.000304 Jy/beam
! Combined flux in latest and established models = 0.0210363 Jy
uvw 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean 1200, 0.01
! Adding 51 model components to the UV plane model.
! The established model now contains 51 components and 0.0210363 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.438 mas, bmaj=15.72 mas, bpa=37.84 degrees
! Estimated noise=0.841818 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -9.30269e-05 Jy
! Component: 100  -  total flux cleaned = -9.30443e-05 Jy
! Component: 150  -  total flux cleaned = -9.30785e-05 Jy
! Component: 200  -  total flux cleaned = -8.10895e-05 Jy
! Component: 250  -  total flux cleaned = -6.36594e-05 Jy
! Component: 300  -  total flux cleaned = -4.41135e-05 Jy
! Component: 350  -  total flux cleaned = -2.54055e-05 Jy
! Component: 400  -  total flux cleaned = -7.47011e-06 Jy
! Component: 450  -  total flux cleaned = 9.7136e-06 Jy
! Component: 500  -  total flux cleaned = 2.6178e-05 Jy
! Component: 550  -  total flux cleaned = 3.97242e-05 Jy
! Component: 600  -  total flux cleaned = 5.4885e-05 Jy
! Component: 650  -  total flux cleaned = 6.9432e-05 Jy
! Component: 700  -  total flux cleaned = 8.33906e-05 Jy
! Component: 750  -  total flux cleaned = 9.67599e-05 Jy
! Component: 800  -  total flux cleaned = 0.00010955 Jy
! Component: 850  -  total flux cleaned = 0.000120085 Jy
! Component: 900  -  total flux cleaned = 0.000131892 Jy
! Component: 950  -  total flux cleaned = 0.000143225 Jy
! Component: 1000  -  total flux cleaned = 0.000154082 Jy
! Component: 1050  -  total flux cleaned = 0.000164482 Jy
! Component: 1100  -  total flux cleaned = 0.00017305 Jy
! Component: 1150  -  total flux cleaned = 0.000182641 Jy
! Component: 1200  -  total flux cleaned = 0.00019184 Jy
! Total flux subtracted in 1200 components = 0.00019184 Jy
! Clean residual min=-0.000827 max=0.000744 Jy/beam
! Clean residual mean=-0.000000 rms=0.000267 Jy/beam
! Combined flux in latest and established models = 0.0212281 Jy
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapl map,true
! 
! Move the cursor into the plot window and press 'H' for help
clrmod true,true
! clrmod: Cleared the established and tentative models.
selfcal
! Performing phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.019943Jy  sigma=1.914104
! Fit after  self-cal, rms=0.019943Jy  sigma=1.914104
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  800, 0.01
! Inverting map and beam 
! Estimated beam: bmin=7.966 mas, bmaj=16.59 mas, bpa=42.85 degrees
! Estimated noise=0.93716 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00757414 Jy
! Component: 100  -  total flux cleaned = 0.0121573 Jy
! Component: 150  -  total flux cleaned = 0.0149306 Jy
! Component: 200  -  total flux cleaned = 0.0166089 Jy
! Component: 250  -  total flux cleaned = 0.0176274 Jy
! Component: 300  -  total flux cleaned = 0.0183386 Jy
! Component: 350  -  total flux cleaned = 0.0189154 Jy
! Component: 400  -  total flux cleaned = 0.0193958 Jy
! Component: 450  -  total flux cleaned = 0.0197991 Jy
! Component: 500  -  total flux cleaned = 0.0201387 Jy
! Component: 550  -  total flux cleaned = 0.020425 Jy
! Component: 600  -  total flux cleaned = 0.0206667 Jy
! Component: 650  -  total flux cleaned = 0.0208708 Jy
! Component: 700  -  total flux cleaned = 0.021045 Jy
! Component: 750  -  total flux cleaned = 0.0211952 Jy
! Component: 800  -  total flux cleaned = 0.0213257 Jy
! Total flux subtracted in 800 components = 0.0213257 Jy
! Clean residual min=-0.001011 max=0.000713 Jy/beam
! Clean residual mean=-0.000001 rms=0.000273 Jy/beam
! Combined flux in latest and established models = 0.0213256 Jy
uvw 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean 1200, 0.01
! Adding 67 model components to the UV plane model.
! The established model now contains 67 components and 0.0213256 Jy
! Inverting map and beam 
! Estimated beam: bmin=9.438 mas, bmaj=15.72 mas, bpa=37.84 degrees
! Estimated noise=0.841818 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.000104021 Jy
! Component: 100  -  total flux cleaned = 0.000185905 Jy
! Component: 150  -  total flux cleaned = 0.00025312 Jy
! Component: 200  -  total flux cleaned = 0.000310906 Jy
! Component: 250  -  total flux cleaned = 0.000360726 Jy
! Component: 300  -  total flux cleaned = 0.000390615 Jy
! Component: 350  -  total flux cleaned = 0.00039718 Jy
! Component: 400  -  total flux cleaned = 0.000402006 Jy
! Component: 450  -  total flux cleaned = 0.000408276 Jy
! Component: 500  -  total flux cleaned = 0.000412893 Jy
! Component: 550  -  total flux cleaned = 0.000418867 Jy
! Component: 600  -  total flux cleaned = 0.000418863 Jy
! Component: 650  -  total flux cleaned = 0.000421732 Jy
! Component: 700  -  total flux cleaned = 0.000423143 Jy
! Component: 750  -  total flux cleaned = 0.000424533 Jy
! Component: 800  -  total flux cleaned = 0.000425929 Jy
! Component: 850  -  total flux cleaned = 0.000427307 Jy
! Component: 900  -  total flux cleaned = 0.000430006 Jy
! Component: 950  -  total flux cleaned = 0.00043132 Jy
! Component: 1000  -  total flux cleaned = 0.000432627 Jy
! Component: 1050  -  total flux cleaned = 0.000433928 Jy
! Component: 1100  -  total flux cleaned = 0.000435219 Jy
! Component: 1150  -  total flux cleaned = 0.000436487 Jy
! Component: 1200  -  total flux cleaned = 0.000438985 Jy
! Total flux subtracted in 1200 components = 0.000438985 Jy
! Clean residual min=-0.000806 max=0.000681 Jy/beam
! Clean residual mean=-0.000002 rms=0.000248 Jy/beam
! Combined flux in latest and established models = 0.0217646 Jy
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapl map,true
! 
! Move the cursor into the plot window and press 'H' for help
x = (field_size-8) * field_cell / 4 
addwin -x,x,-x,x 
clean (512),0.01
! clean: niter=512  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.000133988 Jy
! Component: 100  -  total flux cleaned = -9.08887e-05 Jy
! Component: 150  -  total flux cleaned = -0.000100104 Jy
! Component: 200  -  total flux cleaned = -0.000143636 Jy
! Component: 250  -  total flux cleaned = -0.000127719 Jy
! Component: 300  -  total flux cleaned = -0.000113232 Jy
! Component: 350  -  total flux cleaned = -0.000119874 Jy
! Component: 400  -  total flux cleaned = -0.000139962 Jy
! Component: 450  -  total flux cleaned = -0.00018495 Jy
! Component: 500  -  total flux cleaned = -0.000216205 Jy
! Total flux subtracted in 512 components = -0.000216232 Jy
! Clean residual min=-0.000307 max=0.000310 Jy/beam
! Clean residual mean=-0.000003 rms=0.000118 Jy/beam
! Combined flux in latest and established models = 0.0215484 Jy
levs=1.563,3.125,6.25,12.5,25.0,50.0
mapcol color
! Mapplot colormap: rainbow, contrast: 1 brightness: 0.5.
mapl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 9.438 x 15.72 at 37.84 degrees (North through East)
! Clean map  min=-0.00030674  max=0.019105 Jy/beam
! Clean map  min=-0.00061231  max=0.019253 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
save J1007+1616_C
! Writing UV FITS file: J1007+1616_C.uvf
! Writing 158 model components to file: J1007+1616_C.mod
! wwins: Wrote 18 windows to J1007+1616_C.win
! Adding 91 model components to the UV plane model.
! The established model now contains 157 components and 0.0215484 Jy
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 9.438 x 15.72 at 37.84 degrees (North through East)
! Clean map  min=-0.0006301  max=0.019256 Jy/beam
! Writing clean map to FITS file: J1007+1616_C.fits
! Writing difmap environment to: J1007+1616_C.par
quit
! Quitting program
! Log file difmap.log_137 closed on Sat Feb  1 09:09:14 2020