! Started logfile: difmap.log_319 on Tue Jan 2 19:54:25 2018 ![@/tmp/difmap__4708.1] ! Obeying initialization file: /tmp/difmap__4708.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__4708.2] ![Exited script file: /tmp/difmap__4708.2] obs J0940+5812_C_a_uvm.fits ! Reading UV FITS file: J0940+5812_C_a_uvm.fits ! AN table 1: 39 integrations on 6 of 6 possible baselines. ! AN table 2: 18 integrations on 3 of 3 possible baselines. ! Apparent sampling: 0.868056 visibilities/baseline/integration-bin. ! Found source: J0940+5812 ! ! 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 2000 visibilities. ![Exited script file: /tmp/difmap__4708.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. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 44 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.001072Jy sigma=34.423182 ! Fit after self-cal, rms=0.983380Jy sigma=33.919784 ! 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 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 clean 200, 0.02 ! Inverting map and beam ! Estimated beam: bmin=13.08 mas, bmaj=19.67 mas, bpa=17.34 degrees ! Estimated noise=0.580024 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.0104514 Jy ! Component: 100 - total flux cleaned = 0.0144217 Jy ! Component: 150 - total flux cleaned = 0.0168219 Jy ! Component: 200 - total flux cleaned = 0.0182025 Jy ! Total flux subtracted in 200 components = 0.0182025 Jy ! Clean residual min=-0.001425 max=0.001431 Jy/beam ! Clean residual mean=-0.000171 rms=0.000844 Jy/beam ! Combined flux in latest and established models = 0.0182025 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.026814Jy sigma=0.819808 ! Fit after self-cal, rms=0.026814Jy sigma=0.819808 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=13.08 mas, bmaj=19.67 mas, bpa=17.34 degrees ! Estimated noise=0.580024 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00649257 Jy ! Component: 100 - total flux cleaned = 0.0104208 Jy ! Component: 150 - total flux cleaned = 0.0127978 Jy ! Component: 200 - total flux cleaned = 0.0143547 Jy ! Component: 250 - total flux cleaned = 0.015575 Jy ! Component: 300 - total flux cleaned = 0.0165719 Jy ! Component: 350 - total flux cleaned = 0.0174102 Jy ! Component: 400 - total flux cleaned = 0.0181229 Jy ! Component: 450 - total flux cleaned = 0.0187309 Jy ! Component: 500 - total flux cleaned = 0.0192536 Jy ! Component: 550 - total flux cleaned = 0.0197059 Jy ! Component: 600 - total flux cleaned = 0.0200999 Jy ! Component: 650 - total flux cleaned = 0.0204484 Jy ! Component: 700 - total flux cleaned = 0.0207597 Jy ! Component: 750 - total flux cleaned = 0.0210387 Jy ! Component: 800 - total flux cleaned = 0.0212888 Jy ! Total flux subtracted in 800 components = 0.0212888 Jy ! Clean residual min=-0.001531 max=0.001910 Jy/beam ! Clean residual mean=-0.000120 rms=0.000700 Jy/beam ! Combined flux in latest and established models = 0.0212888 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 199 model components to the UV plane model. ! The established model now contains 199 components and 0.0212888 Jy ! Inverting map and beam ! Estimated beam: bmin=12.4 mas, bmaj=17.33 mas, bpa=29.15 degrees ! Estimated noise=0.492813 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000251165 Jy ! Component: 100 - total flux cleaned = 0.000424502 Jy ! Component: 150 - total flux cleaned = 0.00057349 Jy ! Component: 200 - total flux cleaned = 0.000706984 Jy ! Component: 250 - total flux cleaned = 0.000827438 Jy ! Component: 300 - total flux cleaned = 0.000889801 Jy ! Component: 350 - total flux cleaned = 0.000932445 Jy ! Component: 400 - total flux cleaned = 0.000969484 Jy ! Component: 450 - total flux cleaned = 0.00100927 Jy ! Component: 500 - total flux cleaned = 0.00104765 Jy ! Component: 550 - total flux cleaned = 0.00108107 Jy ! Component: 600 - total flux cleaned = 0.00111349 Jy ! Component: 650 - total flux cleaned = 0.00114497 Jy ! Component: 700 - total flux cleaned = 0.00117555 Jy ! Component: 750 - total flux cleaned = 0.00120527 Jy ! Component: 800 - total flux cleaned = 0.00123411 Jy ! Component: 850 - total flux cleaned = 0.00126216 Jy ! Component: 900 - total flux cleaned = 0.0012864 Jy ! Component: 950 - total flux cleaned = 0.00130428 Jy ! Component: 1000 - total flux cleaned = 0.00131903 Jy ! Component: 1050 - total flux cleaned = 0.00133367 Jy ! Component: 1100 - total flux cleaned = 0.0013482 Jy ! Component: 1150 - total flux cleaned = 0.0013626 Jy ! Component: 1200 - total flux cleaned = 0.00137695 Jy ! Total flux subtracted in 1200 components = 0.00137695 Jy ! Clean residual min=-0.001260 max=0.001887 Jy/beam ! Clean residual mean=-0.000067 rms=0.000574 Jy/beam ! Combined flux in latest and established models = 0.0226658 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.026814Jy sigma=0.819808 ! Fit after self-cal, rms=0.026814Jy sigma=0.819808 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=13.08 mas, bmaj=19.67 mas, bpa=17.34 degrees ! Estimated noise=0.580024 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00649257 Jy ! Component: 100 - total flux cleaned = 0.0104208 Jy ! Component: 150 - total flux cleaned = 0.0127978 Jy ! Component: 200 - total flux cleaned = 0.0143557 Jy ! Component: 250 - total flux cleaned = 0.0155803 Jy ! Component: 300 - total flux cleaned = 0.0165796 Jy ! Component: 350 - total flux cleaned = 0.0174169 Jy ! Component: 400 - total flux cleaned = 0.0181338 Jy ! Component: 450 - total flux cleaned = 0.0187689 Jy ! Component: 500 - total flux cleaned = 0.01934 Jy ! Component: 550 - total flux cleaned = 0.0198568 Jy ! Component: 600 - total flux cleaned = 0.0203295 Jy ! Component: 650 - total flux cleaned = 0.0207641 Jy ! Component: 700 - total flux cleaned = 0.0211642 Jy ! Component: 750 - total flux cleaned = 0.0215326 Jy ! Component: 800 - total flux cleaned = 0.0218722 Jy ! Total flux subtracted in 800 components = 0.0218722 Jy ! Clean residual min=-0.001349 max=0.001599 Jy/beam ! Clean residual mean=-0.000125 rms=0.000632 Jy/beam ! Combined flux in latest and established models = 0.0218722 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 202 model components to the UV plane model. ! The established model now contains 202 components and 0.0218722 Jy ! Inverting map and beam ! Estimated beam: bmin=12.4 mas, bmaj=17.33 mas, bpa=29.15 degrees ! Estimated noise=0.492813 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000296502 Jy ! Component: 100 - total flux cleaned = 0.000557594 Jy ! Component: 150 - total flux cleaned = 0.00079988 Jy ! Component: 200 - total flux cleaned = 0.0010275 Jy ! Component: 250 - total flux cleaned = 0.0012419 Jy ! Component: 300 - total flux cleaned = 0.00144389 Jy ! Component: 350 - total flux cleaned = 0.00163422 Jy ! Component: 400 - total flux cleaned = 0.00181361 Jy ! Component: 450 - total flux cleaned = 0.00198269 Jy ! Component: 500 - total flux cleaned = 0.00214213 Jy ! Component: 550 - total flux cleaned = 0.00229245 Jy ! Component: 600 - total flux cleaned = 0.00243422 Jy ! Component: 650 - total flux cleaned = 0.00256793 Jy ! Component: 700 - total flux cleaned = 0.00269405 Jy ! Component: 750 - total flux cleaned = 0.00276173 Jy ! Component: 800 - total flux cleaned = 0.00282293 Jy ! Component: 850 - total flux cleaned = 0.00288706 Jy ! Component: 900 - total flux cleaned = 0.0029451 Jy ! Component: 950 - total flux cleaned = 0.00300603 Jy ! Component: 1000 - total flux cleaned = 0.00306116 Jy ! Component: 1050 - total flux cleaned = 0.00311497 Jy ! Component: 1100 - total flux cleaned = 0.00315948 Jy ! Component: 1150 - total flux cleaned = 0.00319919 Jy ! Component: 1200 - total flux cleaned = 0.003242 Jy ! Total flux subtracted in 1200 components = 0.003242 Jy ! Clean residual min=-0.001166 max=0.001225 Jy/beam ! Clean residual mean=-0.000051 rms=0.000407 Jy/beam ! Combined flux in latest and established models = 0.0251142 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 65 model components to the UV plane model. ! The established model now contains 264 components and 0.0251142 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.020596Jy sigma=0.605133 ! Fit after self-cal, rms=0.019047Jy sigma=0.542435 selfcal true,true,1.e6 ! Performing amp+phase self-cal over 1e+06 minute time intervals ! ! Correcting IF 1. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 5. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 6. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 7. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 8. ! A total of 54 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.019047Jy sigma=0.542435 ! Fit after self-cal, rms=0.016835Jy sigma=0.526859 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=11.35 mas, bmaj=15.47 mas, bpa=36.86 degrees ! Estimated noise=0.488946 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.35 x 15.47 at 36.86 degrees (North through East) ! Clean map min=-0.00035036 max=0.015637 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 = 1.67508e-05 Jy ! Component: 100 - total flux cleaned = 4.50786e-05 Jy ! Component: 150 - total flux cleaned = 8.16474e-05 Jy ! Component: 200 - total flux cleaned = 0.000106288 Jy ! Component: 250 - total flux cleaned = 0.000121792 Jy ! Component: 300 - total flux cleaned = 0.000110313 Jy ! Component: 350 - total flux cleaned = 0.000101987 Jy ! Component: 400 - total flux cleaned = 9.37044e-05 Jy ! Component: 450 - total flux cleaned = 8.54107e-05 Jy ! Component: 500 - total flux cleaned = 7.7126e-05 Jy ! Total flux subtracted in 512 components = 7.43879e-05 Jy ! Clean residual min=-0.000142 max=0.000140 Jy/beam ! Clean residual mean=0.000001 rms=0.000047 Jy/beam ! Combined flux in latest and established models = 0.0251886 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: 11.35 x 15.47 at 36.86 degrees (North through East) ! Clean map min=-0.00014188 max=0.015604 Jy/beam ! Clean map min=-0.00024738 max=0.015604 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 J0940+5812_C_a ! Writing UV FITS file: J0940+5812_C_a.uvf ! Writing 335 model components to file: J0940+5812_C_a.mod ! wwins: Wrote 3 windows to J0940+5812_C_a.win ! Adding 71 model components to the UV plane model. ! The established model now contains 335 components and 0.0251886 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 11.35 x 15.47 at 36.86 degrees (North through East) ! Clean map min=-0.00019599 max=0.015604 Jy/beam ! Writing clean map to FITS file: J0940+5812_C_a.fits ! Writing difmap environment to: J0940+5812_C_a.par quit ! Quitting program ! Log file difmap.log_319 closed on Tue Jan 2 19:54:46 2018