! Started logfile: difmap.log_282 on Sat Feb  8 18:32:20 2020
![@/tmp/difmap__21047.1]
! Obeying initialization file: /tmp/difmap__21047.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__21047.2]
![Exited script file: /tmp/difmap__21047.2]
obs J1604+4705_X_uva.fits
! Reading UV FITS file: J1604+4705_X_uva.fits
! AN table 1: 16 integrations on 1 of 1 possible baselines.
! AN table 2: 16 integrations on 15 of 15 possible baselines.
! Apparent sampling: 0.980469 visibilities/baseline/integration-bin.
! Found source: J1604+4705
! 
! 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    7.392e+09      3.2e+07          1      3.2e+07
!  02        2    7.424e+09      3.2e+07          1      3.2e+07
!  03        3    7.456e+09      3.2e+07          1      3.2e+07
!  04        4    7.552e+09      3.2e+07          1      3.2e+07
!  05        5    7.744e+09      3.2e+07          1      3.2e+07
!  06        6    7.776e+09      3.2e+07          1      3.2e+07
!  07        7    7.808e+09      3.2e+07          1      3.2e+07
!  08        8     7.84e+09      3.2e+07          1      3.2e+07
! 
! Polarization(s): RR
! 
! Read 94 lines of history.
! 
! Reading 2008 visibilities.
![Exited script file: /tmp/difmap__21047.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.1
taper_size = 70.
mapsize field_size,field_cell
! Map grid = 1024x1024 pixels with 0.100x0.100 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 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 5.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 6.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 7.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 8.
!  A total of 32 telescope corrections were flagged in sub-array 1.
!  A total of 1 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=1.001425Jy  sigma=21.441787
! Fit after  self-cal, rms=0.977431Jy  sigma=21.171561
! 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 16 buffered edits.
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.1
taper_size = 70.
mapsize field_size,field_cell
! Map grid = 1024x1024 pixels with 0.100x0.100 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=0.978119Jy  sigma=21.219092
! Fit after  self-cal, rms=0.975224Jy  sigma=21.161998
! 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
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=2.371 mas, bmaj=2.633 mas, bpa=-39.95 degrees
! Estimated noise=1.05825 mJy/beam.
! clean: niter=200  gain=0.02  cutoff=0
! Component: 050  -  total flux cleaned = 0.0162241 Jy
! Component: 100  -  total flux cleaned = 0.0215454 Jy
! Component: 150  -  total flux cleaned = 0.0213247 Jy
! Component: 200  -  total flux cleaned = 0.0214951 Jy
! Total flux subtracted in 200 components = 0.0214951 Jy
! Clean residual min=-0.001955 max=0.001951 Jy/beam
! Clean residual mean=-0.000003 rms=0.000747 Jy/beam
! Combined flux in latest and established models = 0.0214951 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.039973Jy  sigma=0.851217
! Fit after  self-cal, rms=0.039973Jy  sigma=0.851217
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=2.371 mas, bmaj=2.633 mas, bpa=-39.95 degrees
! Estimated noise=1.05825 mJy/beam.
! clean: niter=800  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.0100775 Jy
! Component: 100  -  total flux cleaned = 0.0161765 Jy
! Component: 150  -  total flux cleaned = 0.01987 Jy
! Component: 200  -  total flux cleaned = 0.022108 Jy
! Component: 250  -  total flux cleaned = 0.0234644 Jy
! Component: 300  -  total flux cleaned = 0.0242869 Jy
! Component: 350  -  total flux cleaned = 0.0247256 Jy
! Component: 400  -  total flux cleaned = 0.0247745 Jy
! Component: 450  -  total flux cleaned = 0.0248041 Jy
! Component: 500  -  total flux cleaned = 0.0247909 Jy
! Component: 550  -  total flux cleaned = 0.0247513 Jy
! Component: 600  -  total flux cleaned = 0.0246999 Jy
! Component: 650  -  total flux cleaned = 0.02465 Jy
! Component: 700  -  total flux cleaned = 0.0246013 Jy
! Component: 750  -  total flux cleaned = 0.0245538 Jy
! Component: 800  -  total flux cleaned = 0.024496 Jy
! Total flux subtracted in 800 components = 0.024496 Jy
! Clean residual min=-0.003560 max=0.003025 Jy/beam
! Clean residual mean=-0.000006 rms=0.000977 Jy/beam
! Combined flux in latest and established models = 0.024496 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 18 model components to the UV plane model.
! The established model now contains 18 components and 0.024496 Jy
! Inverting map and beam 
! Estimated beam: bmin=2.547 mas, bmaj=2.72 mas, bpa=3.516 degrees
! Estimated noise=0.754909 mJy/beam.
! clean: niter=1200  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.000159652 Jy
! Component: 100  -  total flux cleaned = -0.000159305 Jy
! Component: 150  -  total flux cleaned = -0.0001591 Jy
! Component: 200  -  total flux cleaned = -0.000159018 Jy
! Component: 250  -  total flux cleaned = -0.000159016 Jy
! Component: 300  -  total flux cleaned = -0.000174558 Jy
! Component: 350  -  total flux cleaned = -0.000204627 Jy
! Component: 400  -  total flux cleaned = -0.000226218 Jy
! Component: 450  -  total flux cleaned = -0.000254067 Jy
! Component: 500  -  total flux cleaned = -0.000274369 Jy
! Component: 550  -  total flux cleaned = -0.000294167 Jy
! Component: 600  -  total flux cleaned = -0.000319533 Jy
! Component: 650  -  total flux cleaned = -0.000338105 Jy
! Component: 700  -  total flux cleaned = -0.000356274 Jy
! Component: 750  -  total flux cleaned = -0.000379823 Jy
! Component: 800  -  total flux cleaned = -0.000397088 Jy
! Component: 850  -  total flux cleaned = -0.000413968 Jy
! Component: 900  -  total flux cleaned = -0.000435851 Jy
! Component: 950  -  total flux cleaned = -0.000451895 Jy
! Component: 1000  -  total flux cleaned = -0.000467578 Jy
! Component: 1050  -  total flux cleaned = -0.000487908 Jy
! Component: 1100  -  total flux cleaned = -0.000502782 Jy
! Component: 1150  -  total flux cleaned = -0.000522084 Jy
! Component: 1200  -  total flux cleaned = -0.000536183 Jy
! Total flux subtracted in 1200 components = -0.000536183 Jy
! Clean residual min=-0.002500 max=0.002079 Jy/beam
! Clean residual mean=-0.000004 rms=0.000738 Jy/beam
! Combined flux in latest and established models = 0.0239598 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 12 model components to the UV plane model.
! The established model now contains 27 components and 0.0239598 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.030886Jy  sigma=0.653206
! Fit after  self-cal, rms=0.030318Jy  sigma=0.641259
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.030318Jy  sigma=0.641259
! Fit after  self-cal, rms=0.029847Jy  sigma=0.631082
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
! Inverting map and beam 
! Estimated beam: bmin=2.395 mas, bmaj=2.593 mas, bpa=-25.07 degrees
! Estimated noise=0.719817 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.395 x 2.593 at -25.07 degrees (North through East)
! Clean map  min=-0.0021312  max=0.025154 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mapcol grey
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
x = (field_size-8) * field_cell / 4 
addwin -x,x,-x,x 
clean (512),0.01
! Inverting map 
! clean: niter=512  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = -0.00070427 Jy
! Component: 100  -  total flux cleaned = -0.000765704 Jy
! Component: 150  -  total flux cleaned = -0.00067527 Jy
! Component: 200  -  total flux cleaned = -0.000575814 Jy
! Component: 250  -  total flux cleaned = -0.000522624 Jy
! Component: 300  -  total flux cleaned = -0.000440943 Jy
! Component: 350  -  total flux cleaned = -0.000324026 Jy
! Component: 400  -  total flux cleaned = -0.000257206 Jy
! Component: 450  -  total flux cleaned = -0.000159253 Jy
! Component: 500  -  total flux cleaned = -5.80142e-05 Jy
! Total flux subtracted in 512 components = -4.73054e-05 Jy
! Clean residual min=-0.000531 max=0.000532 Jy/beam
! Clean residual mean=0.000000 rms=0.000191 Jy/beam
! Combined flux in latest and established models = 0.0239125 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: 2.395 x 2.593 at -25.07 degrees (North through East)
! Clean map  min=-0.0005283  max=0.025107 Jy/beam
! Clean map  min=-0.0011368  max=0.025107 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 J1604+4705_X
! Writing UV FITS file: J1604+4705_X.uvf
! Writing 94 model components to file: J1604+4705_X.mod
! wwins: Wrote 2 windows to J1604+4705_X.win
! Adding 67 model components to the UV plane model.
! The established model now contains 94 components and 0.0239125 Jy
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.395 x 2.593 at -25.07 degrees (North through East)
! Clean map  min=-0.0011252  max=0.025109 Jy/beam
! Writing clean map to FITS file: J1604+4705_X.fits
! Writing difmap environment to: J1604+4705_X.par
quit
! Quitting program
! Log file difmap.log_282 closed on Sat Feb  8 18:32:48 2020