! Started logfile: difmap.log_59 on Mon Dec 11 00:57:07 2017 ![@/tmp/difmap__19781.1] ! Obeying initialization file: /tmp/difmap__19781.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__19781.2] ![Exited script file: /tmp/difmap__19781.2] obs J2331-3811_S_uva.fits ! Reading UV FITS file: J2331-3811_S_uva.fits ! AN table 1: 4 integrations on 3 of 3 possible baselines. ! AN table 2: 2 integrations on 3 of 3 possible baselines. ! AN table 3: 10 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.842105 visibilities/baseline/integration-bin. ! Found source: J2331-3811 ! ! 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.23275e+09 8e+06 1 8e+06 ! 02 2 2.26275e+09 8e+06 1 8e+06 ! 03 3 2.35275e+09 8e+06 1 8e+06 ! 04 4 2.37275e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 74 lines of history. ! ! Reading 768 visibilities. ![Exited script file: /tmp/difmap__19781.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 6 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.168291Jy sigma=14.472680 ! Fit after self-cal, rms=0.450987Jy sigma=5.442328 ! 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 9 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=2.57 mas, bmaj=6.7 mas, bpa=20.13 degrees ! Estimated noise=2.79834 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.345289 Jy ! Component: 100 - total flux cleaned = 0.472462 Jy ! Component: 150 - total flux cleaned = 0.523588 Jy ! Component: 200 - total flux cleaned = 0.551852 Jy ! Total flux subtracted in 200 components = 0.551852 Jy ! Clean residual min=-0.013455 max=0.022740 Jy/beam ! Clean residual mean=-0.000206 rms=0.003775 Jy/beam ! Combined flux in latest and established models = 0.551852 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. ! A total of 3 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=0.564219Jy sigma=7.011868 ! Fit after self-cal, rms=0.564017Jy sigma=7.014475 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.57 mas, bmaj=6.698 mas, bpa=20.12 degrees ! Estimated noise=2.81967 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.214488 Jy ! Component: 100 - total flux cleaned = 0.344254 Jy ! Component: 150 - total flux cleaned = 0.422794 Jy ! Component: 200 - total flux cleaned = 0.471648 Jy ! Component: 250 - total flux cleaned = 0.502787 Jy ! Component: 300 - total flux cleaned = 0.522639 Jy ! Component: 350 - total flux cleaned = 0.535302 Jy ! Component: 400 - total flux cleaned = 0.54342 Jy ! Component: 450 - total flux cleaned = 0.548723 Jy ! Component: 500 - total flux cleaned = 0.552213 Jy ! Component: 550 - total flux cleaned = 0.553705 Jy ! Component: 600 - total flux cleaned = 0.553787 Jy ! Component: 650 - total flux cleaned = 0.55386 Jy ! Component: 700 - total flux cleaned = 0.553921 Jy ! Component: 750 - total flux cleaned = 0.553977 Jy ! Component: 800 - total flux cleaned = 0.553977 Jy ! Total flux subtracted in 800 components = 0.553977 Jy ! Clean residual min=-0.017842 max=0.029289 Jy/beam ! Clean residual mean=-0.000190 rms=0.003972 Jy/beam ! Combined flux in latest and established models = 0.553978 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 16 model components to the UV plane model. ! The established model now contains 16 components and 0.553978 Jy ! Inverting map and beam ! Estimated beam: bmin=3.798 mas, bmaj=10.18 mas, bpa=15.07 degrees ! Estimated noise=2.08593 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00620289 Jy ! Component: 100 - total flux cleaned = 0.0102276 Jy ! Component: 150 - total flux cleaned = 0.0128885 Jy ! Component: 200 - total flux cleaned = 0.0146573 Jy ! Component: 250 - total flux cleaned = 0.0158331 Jy ! Component: 300 - total flux cleaned = 0.0166147 Jy ! Component: 350 - total flux cleaned = 0.0169314 Jy ! Component: 400 - total flux cleaned = 0.0169515 Jy ! Component: 450 - total flux cleaned = 0.0169713 Jy ! Component: 500 - total flux cleaned = 0.0169908 Jy ! Component: 550 - total flux cleaned = 0.01701 Jy ! Component: 600 - total flux cleaned = 0.0170288 Jy ! Component: 650 - total flux cleaned = 0.0170475 Jy ! Component: 700 - total flux cleaned = 0.0170659 Jy ! Component: 750 - total flux cleaned = 0.0170842 Jy ! Component: 800 - total flux cleaned = 0.0171022 Jy ! Component: 850 - total flux cleaned = 0.017137 Jy ! Component: 900 - total flux cleaned = 0.0171539 Jy ! Component: 950 - total flux cleaned = 0.0171707 Jy ! Component: 1000 - total flux cleaned = 0.0171872 Jy ! Component: 1050 - total flux cleaned = 0.0172036 Jy ! Component: 1100 - total flux cleaned = 0.0172197 Jy ! Component: 1150 - total flux cleaned = 0.0172357 Jy ! Component: 1200 - total flux cleaned = 0.0172513 Jy ! Total flux subtracted in 1200 components = 0.0172513 Jy ! Clean residual min=-0.008931 max=0.022335 Jy/beam ! Clean residual mean=-0.000296 rms=0.003259 Jy/beam ! Combined flux in latest and established models = 0.57123 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 20 components and 0.57123 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.076372Jy sigma=0.862321 ! Fit after self-cal, rms=0.076318Jy sigma=0.861270 selfcal true,true,1.e6 ! Performing amp+phase self-cal over 1e+06 minute time intervals ! ! Correcting IF 1. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.076318Jy sigma=0.861270 ! Fit after self-cal, rms=0.063046Jy sigma=0.687306 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.571267Jy sigma=7.022190 ! Fit after self-cal, rms=0.571267Jy sigma=7.022190 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.406 mas, bmaj=16.39 mas, bpa=8.195 degrees ! Estimated noise=3.04486 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.224138 Jy ! Component: 100 - total flux cleaned = 0.359744 Jy ! Component: 150 - total flux cleaned = 0.441785 Jy ! Component: 200 - total flux cleaned = 0.491421 Jy ! Component: 250 - total flux cleaned = 0.521452 Jy ! Component: 300 - total flux cleaned = 0.53968 Jy ! Component: 350 - total flux cleaned = 0.550776 Jy ! Component: 400 - total flux cleaned = 0.557538 Jy ! Component: 450 - total flux cleaned = 0.561713 Jy ! Component: 500 - total flux cleaned = 0.564337 Jy ! Component: 550 - total flux cleaned = 0.566006 Jy ! Component: 600 - total flux cleaned = 0.567069 Jy ! Component: 650 - total flux cleaned = 0.567536 Jy ! Component: 700 - total flux cleaned = 0.567536 Jy ! Component: 750 - total flux cleaned = 0.567535 Jy ! Component: 800 - total flux cleaned = 0.567535 Jy ! Total flux subtracted in 800 components = 0.567535 Jy ! Clean residual min=-0.011492 max=0.011542 Jy/beam ! Clean residual mean=-0.000077 rms=0.002731 Jy/beam ! Combined flux in latest and established models = 0.567535 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 9 model components to the UV plane model. ! The established model now contains 9 components and 0.567535 Jy ! Inverting map and beam ! Estimated beam: bmin=4.733 mas, bmaj=22.3 mas, bpa=7.474 degrees ! Estimated noise=2.19698 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00107945 Jy ! Component: 100 - total flux cleaned = 0.00173252 Jy ! Component: 150 - total flux cleaned = 0.00212763 Jy ! Component: 200 - total flux cleaned = 0.00236668 Jy ! Component: 250 - total flux cleaned = 0.00238151 Jy ! Component: 300 - total flux cleaned = 0.00238177 Jy ! Component: 350 - total flux cleaned = 0.0023819 Jy ! Component: 400 - total flux cleaned = 0.00238196 Jy ! Component: 450 - total flux cleaned = 0.00238197 Jy ! Component: 500 - total flux cleaned = 0.00238195 Jy ! Component: 550 - total flux cleaned = 0.00238191 Jy ! Component: 600 - total flux cleaned = 0.00238173 Jy ! Component: 650 - total flux cleaned = 0.0023779 Jy ! Component: 700 - total flux cleaned = 0.00237796 Jy ! Component: 750 - total flux cleaned = 0.00237797 Jy ! Component: 800 - total flux cleaned = 0.0023778 Jy ! Component: 850 - total flux cleaned = 0.00237766 Jy ! Component: 900 - total flux cleaned = 0.00237737 Jy ! Component: 950 - total flux cleaned = 0.00237715 Jy ! Component: 1000 - total flux cleaned = 0.00237684 Jy ! Component: 1050 - total flux cleaned = 0.00237656 Jy ! Component: 1100 - total flux cleaned = 0.00237628 Jy ! Component: 1150 - total flux cleaned = 0.00237647 Jy ! Component: 1200 - total flux cleaned = 0.00237661 Jy ! Total flux subtracted in 1200 components = 0.00237661 Jy ! Clean residual min=-0.005867 max=0.006697 Jy/beam ! Clean residual mean=-0.000084 rms=0.001614 Jy/beam ! Combined flux in latest and established models = 0.569912 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 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.733 x 22.3 at 7.474 degrees (North through East) ! Clean map min=-0.0057994 max=0.56722 Jy/beam ! Clean map min=-0.0057994 max=0.56973 Jy/beam ! ! Move the cursor into the plot window and press 'H' for help mapcol grey ! Mapplot colormap: grey, contrast: 1 brightness: 0.5. uncalib selfcal ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 11 components and 0.569912 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063309Jy sigma=0.689113 ! Fit after self-cal, rms=0.063398Jy sigma=0.688857 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.063398Jy sigma=0.688857 ! Fit after self-cal, rms=0.063485Jy sigma=0.688576 x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map and beam ! Estimated beam: bmin=4.745 mas, bmaj=22.37 mas, bpa=7.47 degrees ! Estimated noise=2.19745 mJy/beam. ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000114185 Jy ! Component: 100 - total flux cleaned = -0.0011937 Jy ! Component: 150 - total flux cleaned = -0.00200287 Jy ! Component: 200 - total flux cleaned = -0.00260283 Jy ! Component: 250 - total flux cleaned = -0.00276978 Jy ! Component: 300 - total flux cleaned = -0.00277053 Jy ! Component: 350 - total flux cleaned = -0.00238656 Jy ! Component: 400 - total flux cleaned = -0.00208935 Jy ! Component: 450 - total flux cleaned = -0.00201728 Jy ! Component: 500 - total flux cleaned = -0.00180543 Jy ! Component: 550 - total flux cleaned = -0.00159889 Jy ! Component: 600 - total flux cleaned = -0.00133124 Jy ! Component: 650 - total flux cleaned = -0.00107 Jy ! Component: 700 - total flux cleaned = -0.000687077 Jy ! Component: 750 - total flux cleaned = -0.000437518 Jy ! Component: 800 - total flux cleaned = -0.000375957 Jy ! Component: 850 - total flux cleaned = -7.50541e-05 Jy ! Component: 900 - total flux cleaned = 0.000101502 Jy ! Component: 950 - total flux cleaned = 0.000331689 Jy ! Component: 1000 - total flux cleaned = 0.000388165 Jy ! Component: 1050 - total flux cleaned = 0.000553945 Jy ! Component: 1100 - total flux cleaned = 0.000716549 Jy ! Component: 1150 - total flux cleaned = 0.000715978 Jy ! Component: 1200 - total flux cleaned = 0.000716606 Jy ! Component: 1250 - total flux cleaned = 0.000767536 Jy ! Component: 1300 - total flux cleaned = 0.000666414 Jy ! Component: 1350 - total flux cleaned = 0.000716768 Jy ! Component: 1400 - total flux cleaned = 0.000667927 Jy ! Component: 1450 - total flux cleaned = 0.000715915 Jy ! Component: 1500 - total flux cleaned = 0.000715481 Jy ! Component: 1550 - total flux cleaned = 0.000668074 Jy ! Component: 1600 - total flux cleaned = 0.00052776 Jy ! Component: 1650 - total flux cleaned = 0.000296923 Jy ! Component: 1700 - total flux cleaned = 6.87886e-05 Jy ! Component: 1750 - total flux cleaned = -0.000381417 Jy ! Component: 1800 - total flux cleaned = -0.000514374 Jy ! Component: 1850 - total flux cleaned = -0.000911232 Jy ! Component: 1900 - total flux cleaned = -0.000998776 Jy ! Component: 1950 - total flux cleaned = -0.00125793 Jy ! Component: 2000 - total flux cleaned = -0.00155708 Jy ! Component: 2050 - total flux cleaned = -0.00185284 Jy ! Component: 2100 - total flux cleaned = -0.00197835 Jy ! Component: 2150 - total flux cleaned = -0.00226813 Jy ! Component: 2200 - total flux cleaned = -0.00230906 Jy ! Component: 2250 - total flux cleaned = -0.00259367 Jy ! Component: 2300 - total flux cleaned = -0.00279402 Jy ! Component: 2350 - total flux cleaned = -0.00299239 Jy ! Component: 2400 - total flux cleaned = -0.00307078 Jy ! Component: 2450 - total flux cleaned = -0.0032258 Jy ! Component: 2500 - total flux cleaned = -0.00337963 Jy ! Component: 2550 - total flux cleaned = -0.00341857 Jy ! Component: 2600 - total flux cleaned = -0.00372227 Jy ! Component: 2650 - total flux cleaned = -0.0039102 Jy ! Component: 2700 - total flux cleaned = -0.00405976 Jy ! Component: 2750 - total flux cleaned = -0.00417065 Jy ! Component: 2800 - total flux cleaned = -0.00442732 Jy ! Component: 2850 - total flux cleaned = -0.0045001 Jy ! Component: 2900 - total flux cleaned = -0.00468064 Jy ! Component: 2950 - total flux cleaned = -0.00482358 Jy ! Component: 3000 - total flux cleaned = -0.00492996 Jy ! Component: 3050 - total flux cleaned = -0.00503583 Jy ! Component: 3100 - total flux cleaned = -0.00524501 Jy ! Component: 3150 - total flux cleaned = -0.00534895 Jy ! Component: 3200 - total flux cleaned = -0.00548631 Jy ! Component: 3250 - total flux cleaned = -0.00548652 Jy ! Component: 3300 - total flux cleaned = -0.00579044 Jy ! Component: 3350 - total flux cleaned = -0.00585721 Jy ! Component: 3400 - total flux cleaned = -0.00599056 Jy ! Component: 3450 - total flux cleaned = -0.00602368 Jy ! Component: 3500 - total flux cleaned = -0.00622007 Jy ! Component: 3550 - total flux cleaned = -0.00622017 Jy ! Component: 3600 - total flux cleaned = -0.00631671 Jy ! Component: 3650 - total flux cleaned = -0.00638078 Jy ! Component: 3700 - total flux cleaned = -0.00641251 Jy ! Component: 3750 - total flux cleaned = -0.00650713 Jy ! Component: 3800 - total flux cleaned = -0.00660105 Jy ! Component: 3850 - total flux cleaned = -0.00672555 Jy ! Component: 3900 - total flux cleaned = -0.00675644 Jy ! Component: 3950 - total flux cleaned = -0.00681783 Jy ! Component: 4000 - total flux cleaned = -0.00690943 Jy ! Component: 4050 - total flux cleaned = -0.00684875 Jy ! Total flux subtracted in 4096 components = -0.00693927 Jy ! Clean residual min=-0.001503 max=0.001595 Jy/beam ! Clean residual mean=-0.000029 rms=0.000655 Jy/beam ! Combined flux in latest and established models = 0.562973 Jy levs=0.4,0.8,1.6,3.2,6.4,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.745 x 22.37 at 7.47 degrees (North through East) ! Clean map min=-0.0014929 max=0.56949 Jy/beam ! Clean map min=-0.0074639 max=0.56949 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. save J2331-3811_S ! Writing UV FITS file: J2331-3811_S.uvf ! Writing 265 model components to file: J2331-3811_S.mod ! wwins: Wrote 2 windows to J2331-3811_S.win ! Adding 254 model components to the UV plane model. ! The established model now contains 265 components and 0.562973 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.745 x 22.37 at 7.47 degrees (North through East) ! Clean map min=-0.007479 max=0.56946 Jy/beam ! Writing clean map to FITS file: J2331-3811_S.fits ! Writing difmap environment to: J2331-3811_S.par quit ! Quitting program ! Log file difmap.log_59 closed on Mon Dec 11 00:57:36 2017