! Started logfile: difmap.log_82 on Sun Jul 31 17:47:08 2016 ![@/tmp/difmap__5720.1] ! Obeying initialization file: /tmp/difmap__5720.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__5720.2] ![Exited script file: /tmp/difmap__5720.2] obs J0628-2805_S_uva.fits ! Reading UV FITS file: J0628-2805_S_uva.fits ! AN table 1: 28 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.762897 visibilities/baseline/integration-bin. ! Found source: J0628-2805 ! ! 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 2.22e+09 3.2e+07 1 3.2e+07 ! 02 2 2.252e+09 3.2e+07 1 3.2e+07 ! 03 3 2.284e+09 3.2e+07 1 3.2e+07 ! 04 4 2.348e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 54 lines of history. ! ! Reading 3076 visibilities. ![Exited script file: /tmp/difmap__5720.1] sele I ! Polarization I is unavailable. ! 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 field_size=1024 field_cell=0.4 taper_size = 20. mapsize field_size,field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 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 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.008284Jy sigma=16.916478 ! Fit after self-cal, rms=0.905892Jy sigma=15.218242 ! 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 2 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=3.089 mas, bmaj=13.02 mas, bpa=18.05 degrees ! Estimated noise=1.18362 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.0613527 Jy ! Component: 100 - total flux cleaned = 0.0836955 Jy ! Component: 150 - total flux cleaned = 0.0889263 Jy ! Component: 200 - total flux cleaned = 0.0913754 Jy ! Total flux subtracted in 200 components = 0.0913754 Jy ! Clean residual min=-0.005561 max=0.005585 Jy/beam ! Clean residual mean=0.000021 rms=0.001758 Jy/beam ! Combined flux in latest and established models = 0.0913754 Jy mapl map,true ! ! Move the cursor into the plot window and press 'H' for help ! Note that 17 components lie outside the plot. 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. ! ! Fit before self-cal, rms=0.127717Jy sigma=2.092460 ! Fit after self-cal, rms=0.127717Jy sigma=2.092460 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=3.089 mas, bmaj=13.02 mas, bpa=18.05 degrees ! Estimated noise=1.18362 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0381139 Jy ! Component: 100 - total flux cleaned = 0.061173 Jy ! Component: 150 - total flux cleaned = 0.0751239 Jy ! Component: 200 - total flux cleaned = 0.0835643 Jy ! Component: 250 - total flux cleaned = 0.0886739 Jy ! Component: 300 - total flux cleaned = 0.0917704 Jy ! Component: 350 - total flux cleaned = 0.0936829 Jy ! Component: 400 - total flux cleaned = 0.0949579 Jy ! Component: 450 - total flux cleaned = 0.0958477 Jy ! Component: 500 - total flux cleaned = 0.0964826 Jy ! Component: 550 - total flux cleaned = 0.0969416 Jy ! Component: 600 - total flux cleaned = 0.0972767 Jy ! Component: 650 - total flux cleaned = 0.0975237 Jy ! Component: 700 - total flux cleaned = 0.097706 Jy ! Component: 750 - total flux cleaned = 0.0977762 Jy ! Component: 800 - total flux cleaned = 0.0977927 Jy ! Total flux subtracted in 800 components = 0.0977927 Jy ! Clean residual min=-0.007788 max=0.006949 Jy/beam ! Clean residual mean=0.000035 rms=0.001867 Jy/beam ! Combined flux in latest and established models = 0.0977927 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 14 model components to the UV plane model. ! The established model now contains 14 components and 0.0977927 Jy ! Inverting map and beam ! Estimated beam: bmin=4.049 mas, bmaj=16.18 mas, bpa=16.75 degrees ! Estimated noise=0.844177 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000286711 Jy ! Component: 100 - total flux cleaned = 0.000460252 Jy ! Component: 150 - total flux cleaned = 0.000476126 Jy ! Component: 200 - total flux cleaned = 0.000486473 Jy ! Component: 250 - total flux cleaned = 0.000491575 Jy ! Component: 300 - total flux cleaned = 0.000501626 Jy ! Component: 350 - total flux cleaned = 0.000506603 Jy ! Component: 400 - total flux cleaned = 0.000516377 Jy ! Component: 450 - total flux cleaned = 0.00052119 Jy ! Component: 500 - total flux cleaned = 0.000530737 Jy ! Component: 550 - total flux cleaned = 0.00054013 Jy ! Component: 600 - total flux cleaned = 0.000544769 Jy ! Component: 650 - total flux cleaned = 0.000553985 Jy ! Component: 700 - total flux cleaned = 0.000558627 Jy ! Component: 750 - total flux cleaned = 0.000567626 Jy ! Component: 800 - total flux cleaned = 0.000572092 Jy ! Component: 850 - total flux cleaned = 0.00058093 Jy ! Component: 900 - total flux cleaned = 0.00058968 Jy ! Component: 950 - total flux cleaned = 0.000594 Jy ! Component: 1000 - total flux cleaned = 0.000602553 Jy ! Component: 1050 - total flux cleaned = 0.000606871 Jy ! Component: 1100 - total flux cleaned = 0.000615387 Jy ! Component: 1150 - total flux cleaned = 0.000623712 Jy ! Component: 1200 - total flux cleaned = 0.000627865 Jy ! Total flux subtracted in 1200 components = 0.000627865 Jy ! Clean residual min=-0.006366 max=0.007088 Jy/beam ! Clean residual mean=0.000060 rms=0.001811 Jy/beam ! Combined flux in latest and established models = 0.0984205 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 6 model components to the UV plane model. ! The established model now contains 17 components and 0.0984205 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.082208Jy sigma=1.306649 ! Fit after self-cal, rms=0.081750Jy sigma=1.295539 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. ! ! Fit before self-cal, rms=0.081750Jy sigma=1.295539 ! Fit after self-cal, rms=0.082138Jy sigma=1.270849 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. ! ! Fit before self-cal, rms=0.127418Jy sigma=2.092460 ! Fit after self-cal, rms=0.127418Jy sigma=2.092460 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=3.08 mas, bmaj=13.24 mas, bpa=18.25 degrees ! Estimated noise=1.18103 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0381045 Jy ! Component: 100 - total flux cleaned = 0.061158 Jy ! Component: 150 - total flux cleaned = 0.0751055 Jy ! Component: 200 - total flux cleaned = 0.0835476 Jy ! Component: 250 - total flux cleaned = 0.0886627 Jy ! Component: 300 - total flux cleaned = 0.091762 Jy ! Component: 350 - total flux cleaned = 0.0936399 Jy ! Component: 400 - total flux cleaned = 0.0948093 Jy ! Component: 450 - total flux cleaned = 0.0956448 Jy ! Component: 500 - total flux cleaned = 0.0962783 Jy ! Component: 550 - total flux cleaned = 0.0967646 Jy ! Component: 600 - total flux cleaned = 0.0971391 Jy ! Component: 650 - total flux cleaned = 0.0974283 Jy ! Component: 700 - total flux cleaned = 0.0976517 Jy ! Component: 750 - total flux cleaned = 0.0978243 Jy ! Component: 800 - total flux cleaned = 0.0978872 Jy ! Total flux subtracted in 800 components = 0.0978872 Jy ! Clean residual min=-0.005838 max=0.005966 Jy/beam ! Clean residual mean=0.000011 rms=0.001474 Jy/beam ! Combined flux in latest and established models = 0.0978872 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 13 model components to the UV plane model. ! The established model now contains 13 components and 0.0978872 Jy ! Inverting map and beam ! Estimated beam: bmin=4.092 mas, bmaj=16.82 mas, bpa=17.05 degrees ! Estimated noise=0.831801 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000262226 Jy ! Component: 100 - total flux cleaned = 0.000420983 Jy ! Component: 150 - total flux cleaned = 0.000507691 Jy ! Component: 200 - total flux cleaned = 0.000514637 Jy ! Component: 250 - total flux cleaned = 0.00051812 Jy ! Component: 300 - total flux cleaned = 0.000524791 Jy ! Component: 350 - total flux cleaned = 0.000531357 Jy ! Component: 400 - total flux cleaned = 0.000537768 Jy ! Component: 450 - total flux cleaned = 0.000540907 Jy ! Component: 500 - total flux cleaned = 0.00054712 Jy ! Component: 550 - total flux cleaned = 0.000553157 Jy ! Component: 600 - total flux cleaned = 0.000556131 Jy ! Component: 650 - total flux cleaned = 0.000561996 Jy ! Component: 700 - total flux cleaned = 0.000567674 Jy ! Component: 750 - total flux cleaned = 0.000570492 Jy ! Component: 800 - total flux cleaned = 0.000576005 Jy ! Component: 850 - total flux cleaned = 0.000581362 Jy ! Component: 900 - total flux cleaned = 0.000584027 Jy ! Component: 950 - total flux cleaned = 0.00058921 Jy ! Component: 1000 - total flux cleaned = 0.000594266 Jy ! Component: 1050 - total flux cleaned = 0.000599217 Jy ! Component: 1100 - total flux cleaned = 0.000601657 Jy ! Component: 1150 - total flux cleaned = 0.000606426 Jy ! Component: 1200 - total flux cleaned = 0.000611079 Jy ! Total flux subtracted in 1200 components = 0.000611079 Jy ! Clean residual min=-0.004082 max=0.004022 Jy/beam ! Clean residual mean=0.000016 rms=0.001063 Jy/beam ! Combined flux in latest and established models = 0.0984983 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 (field_size*4),0.01 ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000366148 Jy ! Component: 100 - total flux cleaned = -0.00129411 Jy ! Component: 150 - total flux cleaned = -0.00196882 Jy ! Component: 200 - total flux cleaned = -0.00196942 Jy ! Component: 250 - total flux cleaned = -0.00213533 Jy ! Component: 300 - total flux cleaned = -0.00234701 Jy ! Component: 350 - total flux cleaned = -0.00255036 Jy ! Component: 400 - total flux cleaned = -0.00264818 Jy ! Component: 450 - total flux cleaned = -0.00288518 Jy ! Component: 500 - total flux cleaned = -0.00283915 Jy ! Component: 550 - total flux cleaned = -0.00274908 Jy ! Component: 600 - total flux cleaned = -0.00261685 Jy ! Component: 650 - total flux cleaned = -0.00248814 Jy ! Component: 700 - total flux cleaned = -0.00240343 Jy ! Component: 750 - total flux cleaned = -0.00228027 Jy ! Component: 800 - total flux cleaned = -0.00219971 Jy ! Component: 850 - total flux cleaned = -0.00212086 Jy ! Component: 900 - total flux cleaned = -0.00200396 Jy ! Component: 950 - total flux cleaned = -0.00192747 Jy ! Component: 1000 - total flux cleaned = -0.00185239 Jy ! Component: 1050 - total flux cleaned = -0.00185286 Jy ! Component: 1100 - total flux cleaned = -0.00181666 Jy ! Component: 1150 - total flux cleaned = -0.00185192 Jy ! Component: 1200 - total flux cleaned = -0.00178127 Jy ! Component: 1250 - total flux cleaned = -0.00174635 Jy ! Component: 1300 - total flux cleaned = -0.00164393 Jy ! Component: 1350 - total flux cleaned = -0.00164403 Jy ! Component: 1400 - total flux cleaned = -0.00164353 Jy ! Component: 1450 - total flux cleaned = -0.00174135 Jy ! Component: 1500 - total flux cleaned = -0.00177313 Jy ! Component: 1550 - total flux cleaned = -0.00180438 Jy ! Component: 1600 - total flux cleaned = -0.00186712 Jy ! Component: 1650 - total flux cleaned = -0.00199132 Jy ! Component: 1700 - total flux cleaned = -0.00202213 Jy ! Component: 1750 - total flux cleaned = -0.00202214 Jy ! Component: 1800 - total flux cleaned = -0.00214193 Jy ! Component: 1850 - total flux cleaned = -0.00217132 Jy ! Component: 1900 - total flux cleaned = -0.00220058 Jy ! Component: 1950 - total flux cleaned = -0.00231598 Jy ! Component: 2000 - total flux cleaned = -0.00231596 Jy ! Component: 2050 - total flux cleaned = -0.00231604 Jy ! Component: 2100 - total flux cleaned = -0.00239951 Jy ! Component: 2150 - total flux cleaned = -0.00234425 Jy ! Component: 2200 - total flux cleaned = -0.00242572 Jy ! Component: 2250 - total flux cleaned = -0.00239848 Jy ! Component: 2300 - total flux cleaned = -0.00245195 Jy ! Component: 2350 - total flux cleaned = -0.00237332 Jy ! Component: 2400 - total flux cleaned = -0.0024254 Jy ! Component: 2450 - total flux cleaned = -0.00245124 Jy ! Component: 2500 - total flux cleaned = -0.00234903 Jy ! Component: 2550 - total flux cleaned = -0.00239991 Jy ! Component: 2600 - total flux cleaned = -0.00237468 Jy ! Component: 2650 - total flux cleaned = -0.00244924 Jy ! Component: 2700 - total flux cleaned = -0.00244907 Jy ! Component: 2750 - total flux cleaned = -0.00254664 Jy ! Component: 2800 - total flux cleaned = -0.00261904 Jy ! Component: 2850 - total flux cleaned = -0.00266693 Jy ! Component: 2900 - total flux cleaned = -0.0027142 Jy ! Component: 2950 - total flux cleaned = -0.00271427 Jy ! Component: 3000 - total flux cleaned = -0.00280775 Jy ! Component: 3050 - total flux cleaned = -0.00276135 Jy ! Component: 3100 - total flux cleaned = -0.00283045 Jy ! Component: 3150 - total flux cleaned = -0.00283049 Jy ! Component: 3200 - total flux cleaned = -0.00287592 Jy ! Component: 3250 - total flux cleaned = -0.00289841 Jy ! Component: 3300 - total flux cleaned = -0.00294273 Jy ! Component: 3350 - total flux cleaned = -0.00292034 Jy ! Component: 3400 - total flux cleaned = -0.00289823 Jy ! Component: 3450 - total flux cleaned = -0.00296362 Jy ! Component: 3500 - total flux cleaned = -0.00289866 Jy ! Component: 3550 - total flux cleaned = -0.00296313 Jy ! Component: 3600 - total flux cleaned = -0.00294194 Jy ! Component: 3650 - total flux cleaned = -0.00300528 Jy ! Component: 3700 - total flux cleaned = -0.00300515 Jy ! Component: 3750 - total flux cleaned = -0.00300515 Jy ! Component: 3800 - total flux cleaned = -0.00300517 Jy ! Component: 3850 - total flux cleaned = -0.00304631 Jy ! Component: 3900 - total flux cleaned = -0.00306673 Jy ! Component: 3950 - total flux cleaned = -0.00312781 Jy ! Component: 4000 - total flux cleaned = -0.00314808 Jy ! Component: 4050 - total flux cleaned = -0.00310803 Jy ! Total flux subtracted in 4096 components = -0.00314776 Jy ! Clean residual min=-0.001101 max=0.001141 Jy/beam ! Clean residual mean=0.000001 rms=0.000402 Jy/beam ! Combined flux in latest and established models = 0.0953505 Jy levs=0.75,1.5,3.0,6.0,12,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: 4.092 x 16.82 at 17.05 degrees (North through East) ! Clean map min=-0.0011011 max=0.097027 Jy/beam ! Clean map min=-0.0029856 max=0.097424 Jy/beam ! ! Move the cursor into the plot window and press 'H' for help mapcol grey ! Mapplot colormap: grey, contrast: 1 brightness: 0.5. 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 ! ! 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 ! Adding 266 model components to the UV plane model. ! The established model now contains 276 components and 0.0953505 Jy ! Inverting map ! clean: niter=512 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000375475 Jy ! Component: 100 - total flux cleaned = -0.000538268 Jy ! Component: 150 - total flux cleaned = -0.000676444 Jy ! Component: 200 - total flux cleaned = -0.000772988 Jy ! Component: 250 - total flux cleaned = -0.000791944 Jy ! Component: 300 - total flux cleaned = -0.000829601 Jy ! Component: 350 - total flux cleaned = -0.000829461 Jy ! Component: 400 - total flux cleaned = -0.00081105 Jy ! Component: 450 - total flux cleaned = -0.000720737 Jy ! Component: 500 - total flux cleaned = -0.00075681 Jy ! Total flux subtracted in 512 components = -0.000721253 Jy ! Clean residual min=-0.000889 max=0.001033 Jy/beam ! Clean residual mean=-0.000002 rms=0.000348 Jy/beam ! Combined flux in latest and established models = 0.0946293 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: 4.092 x 16.82 at 17.05 degrees (North through East) ! Clean map min=-0.0025757 max=0.097385 Jy/beam ! Clean map min=-0.0025757 max=0.097385 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 J0628-2805_S ! Writing UV FITS file: J0628-2805_S.uvf ! Writing 346 model components to file: J0628-2805_S.mod ! wwins: Wrote 3 windows to J0628-2805_S.win ! Adding 70 model components to the UV plane model. ! The established model now contains 324 components and 0.0946293 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.092 x 16.82 at 17.05 degrees (North through East) ! Clean map min=-0.0025895 max=0.097386 Jy/beam ! Writing clean map to FITS file: J0628-2805_S.fits ! Writing difmap environment to: J0628-2805_S.par quit ! Quitting program ! Log file difmap.log_82 closed on Sun Jul 31 17:47:46 2016