! Started logfile: difmap.log_242 on Mon Nov 11 09:13:26 2019 ![@/tmp/difmap__15021.1] ! Obeying initialization file: /tmp/difmap__15021.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__15021.2] ![Exited script file: /tmp/difmap__15021.2] obs J2219+1806_S_uva.fits ! Reading UV FITS file: J2219+1806_S_uva.fits ! AN table 1: 29 integrations on 21 of 21 possible baselines. ! AN table 2: 15 integrations on 3 of 3 possible baselines. ! AN table 3: 19 integrations on 28 of 28 possible baselines. ! AN table 4: 18 integrations on 15 of 15 possible baselines. ! Apparent sampling: 0.940247 visibilities/baseline/integration-bin. ! Found source: J2219+1806 ! ! 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.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 79 lines of history. ! ! Reading 5476 visibilities. ![Exited script file: /tmp/difmap__15021.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 11 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.011628Jy sigma=15.656699 ! Fit after self-cal, rms=0.883927Jy sigma=13.768539 ! 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 11 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.207 mas, bmaj=6.296 mas, bpa=-9.457 degrees ! Estimated noise=1.17335 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.074781 Jy ! Component: 100 - total flux cleaned = 0.102023 Jy ! Component: 150 - total flux cleaned = 0.110714 Jy ! Component: 200 - total flux cleaned = 0.11216 Jy ! Total flux subtracted in 200 components = 0.11216 Jy ! Clean residual min=-0.004615 max=0.004667 Jy/beam ! Clean residual mean=-0.000004 rms=0.001260 Jy/beam ! Combined flux in latest and established models = 0.11216 Jy mapl map,true ! ! Move the cursor into the plot window and press 'H' for help ! Note that 7 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.141594Jy sigma=2.111874 ! Fit after self-cal, rms=0.141594Jy sigma=2.111874 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.207 mas, bmaj=6.296 mas, bpa=-9.457 degrees ! Estimated noise=1.17335 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0464559 Jy ! Component: 100 - total flux cleaned = 0.0745619 Jy ! Component: 150 - total flux cleaned = 0.0915663 Jy ! Component: 200 - total flux cleaned = 0.101863 Jy ! Component: 250 - total flux cleaned = 0.10811 Jy ! Component: 300 - total flux cleaned = 0.1119 Jy ! Component: 350 - total flux cleaned = 0.114356 Jy ! Component: 400 - total flux cleaned = 0.116206 Jy ! Component: 450 - total flux cleaned = 0.117624 Jy ! Component: 500 - total flux cleaned = 0.118714 Jy ! Component: 550 - total flux cleaned = 0.119552 Jy ! Component: 600 - total flux cleaned = 0.120199 Jy ! Component: 650 - total flux cleaned = 0.120698 Jy ! Component: 700 - total flux cleaned = 0.121084 Jy ! Component: 750 - total flux cleaned = 0.121249 Jy ! Component: 800 - total flux cleaned = 0.121261 Jy ! Total flux subtracted in 800 components = 0.121261 Jy ! Clean residual min=-0.006763 max=0.005311 Jy/beam ! Clean residual mean=0.000004 rms=0.001273 Jy/beam ! Combined flux in latest and established models = 0.121261 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 11 model components to the UV plane model. ! The established model now contains 11 components and 0.121261 Jy ! Inverting map and beam ! Estimated beam: bmin=3.064 mas, bmaj=7.922 mas, bpa=-10.33 degrees ! Estimated noise=0.626496 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000344762 Jy ! Component: 100 - total flux cleaned = 0.000374999 Jy ! Component: 150 - total flux cleaned = 0.000374414 Jy ! Component: 200 - total flux cleaned = 0.000373962 Jy ! Component: 250 - total flux cleaned = 0.000373607 Jy ! Component: 300 - total flux cleaned = 0.000367705 Jy ! Component: 350 - total flux cleaned = 0.000367696 Jy ! Component: 400 - total flux cleaned = 0.000372975 Jy ! Component: 450 - total flux cleaned = 0.000383145 Jy ! Component: 500 - total flux cleaned = 0.000388027 Jy ! Component: 550 - total flux cleaned = 0.000397522 Jy ! Component: 600 - total flux cleaned = 0.000402033 Jy ! Component: 650 - total flux cleaned = 0.00040648 Jy ! Component: 700 - total flux cleaned = 0.000410806 Jy ! Component: 750 - total flux cleaned = 0.000410858 Jy ! Component: 800 - total flux cleaned = 0.000415098 Jy ! Component: 850 - total flux cleaned = 0.000419314 Jy ! Component: 900 - total flux cleaned = 0.0004191 Jy ! Component: 950 - total flux cleaned = 0.000419094 Jy ! Component: 1000 - total flux cleaned = 0.000415037 Jy ! Component: 1050 - total flux cleaned = 0.000410904 Jy ! Component: 1100 - total flux cleaned = 0.000406864 Jy ! Component: 1150 - total flux cleaned = 0.00040282 Jy ! Component: 1200 - total flux cleaned = 0.000402725 Jy ! Total flux subtracted in 1200 components = 0.000402725 Jy ! Clean residual min=-0.002959 max=0.003473 Jy/beam ! Clean residual mean=-0.000018 rms=0.000783 Jy/beam ! Combined flux in latest and established models = 0.121663 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 8 model components to the UV plane model. ! The established model now contains 16 components and 0.121663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.076318Jy sigma=1.006623 ! Fit after self-cal, rms=0.076178Jy sigma=1.005148 selfcal true,true,1.e6 ! Performing amp+phase self-cal over 1e+06 minute time intervals ! ! Correcting IF 1. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 36 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.076178Jy sigma=1.005148 ! Fit after self-cal, rms=0.074873Jy sigma=0.996717 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.139763Jy sigma=2.109135 ! Fit after self-cal, rms=0.139763Jy sigma=2.109135 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.186 mas, bmaj=6.236 mas, bpa=-9.586 degrees ! Estimated noise=1.19098 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0459808 Jy ! Component: 100 - total flux cleaned = 0.0737995 Jy ! Component: 150 - total flux cleaned = 0.0906325 Jy ! Component: 200 - total flux cleaned = 0.10084 Jy ! Component: 250 - total flux cleaned = 0.107034 Jy ! Component: 300 - total flux cleaned = 0.11083 Jy ! Component: 350 - total flux cleaned = 0.113518 Jy ! Component: 400 - total flux cleaned = 0.11555 Jy ! Component: 450 - total flux cleaned = 0.117103 Jy ! Component: 500 - total flux cleaned = 0.1183 Jy ! Component: 550 - total flux cleaned = 0.119224 Jy ! Component: 600 - total flux cleaned = 0.119939 Jy ! Component: 650 - total flux cleaned = 0.120493 Jy ! Component: 700 - total flux cleaned = 0.120924 Jy ! Component: 750 - total flux cleaned = 0.121054 Jy ! Component: 800 - total flux cleaned = 0.121108 Jy ! Total flux subtracted in 800 components = 0.121108 Jy ! Clean residual min=-0.005875 max=0.005700 Jy/beam ! Clean residual mean=0.000002 rms=0.001170 Jy/beam ! Combined flux in latest and established models = 0.121108 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.121108 Jy ! Inverting map and beam ! Estimated beam: bmin=3.029 mas, bmaj=7.77 mas, bpa=-10.07 degrees ! Estimated noise=0.634302 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000383092 Jy ! Component: 100 - total flux cleaned = 0.000614864 Jy ! Component: 150 - total flux cleaned = 0.000755088 Jy ! Component: 200 - total flux cleaned = 0.000839925 Jy ! Component: 250 - total flux cleaned = 0.000879759 Jy ! Component: 300 - total flux cleaned = 0.000879719 Jy ! Component: 350 - total flux cleaned = 0.000879773 Jy ! Component: 400 - total flux cleaned = 0.000881536 Jy ! Component: 450 - total flux cleaned = 0.000881614 Jy ! Component: 500 - total flux cleaned = 0.000881657 Jy ! Component: 550 - total flux cleaned = 0.000883284 Jy ! Component: 600 - total flux cleaned = 0.000883213 Jy ! Component: 650 - total flux cleaned = 0.000883248 Jy ! Component: 700 - total flux cleaned = 0.000884791 Jy ! Component: 750 - total flux cleaned = 0.000883293 Jy ! Component: 800 - total flux cleaned = 0.000881806 Jy ! Component: 850 - total flux cleaned = 0.000881811 Jy ! Component: 900 - total flux cleaned = 0.000880339 Jy ! Component: 950 - total flux cleaned = 0.00088034 Jy ! Component: 1000 - total flux cleaned = 0.000878878 Jy ! Component: 1050 - total flux cleaned = 0.000877434 Jy ! Component: 1100 - total flux cleaned = 0.000877424 Jy ! Component: 1150 - total flux cleaned = 0.000875981 Jy ! Component: 1200 - total flux cleaned = 0.000875981 Jy ! Total flux subtracted in 1200 components = 0.000875981 Jy ! Clean residual min=-0.002271 max=0.003122 Jy/beam ! Clean residual mean=-0.000007 rms=0.000648 Jy/beam ! Combined flux in latest and established models = 0.121984 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 14 components and 0.121984 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074862Jy sigma=0.996732 ! Fit after self-cal, rms=0.074888Jy sigma=0.996603 selfcal true,true,10. ! Performing amp+phase self-cal over 10 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074888Jy sigma=0.996603 ! Fit after self-cal, rms=0.074754Jy sigma=0.996342 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=3.021 mas, bmaj=7.712 mas, bpa=-10.08 degrees ! Estimated noise=0.635272 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.021 x 7.712 at -10.08 degrees (North through East) ! Clean map min=-0.002258 max=0.11853 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 (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00129979 Jy ! Component: 100 - total flux cleaned = 0.00191257 Jy ! Component: 150 - total flux cleaned = 0.0021505 Jy ! Component: 200 - total flux cleaned = 0.00218739 Jy ! Component: 250 - total flux cleaned = 0.00218803 Jy ! Component: 300 - total flux cleaned = 0.00218849 Jy ! Component: 350 - total flux cleaned = 0.00215613 Jy ! Component: 400 - total flux cleaned = 0.00209428 Jy ! Component: 450 - total flux cleaned = 0.00209435 Jy ! Component: 500 - total flux cleaned = 0.00212324 Jy ! Component: 550 - total flux cleaned = 0.00215103 Jy ! Component: 600 - total flux cleaned = 0.00217848 Jy ! Component: 650 - total flux cleaned = 0.00215212 Jy ! Component: 700 - total flux cleaned = 0.00207424 Jy ! Component: 750 - total flux cleaned = 0.00207417 Jy ! Component: 800 - total flux cleaned = 0.00209941 Jy ! Component: 850 - total flux cleaned = 0.00207556 Jy ! Component: 900 - total flux cleaned = 0.00209974 Jy ! Component: 950 - total flux cleaned = 0.00202953 Jy ! Component: 1000 - total flux cleaned = 0.00207607 Jy ! Component: 1050 - total flux cleaned = 0.00200772 Jy ! Component: 1100 - total flux cleaned = 0.00203023 Jy ! Component: 1150 - total flux cleaned = 0.00203054 Jy ! Component: 1200 - total flux cleaned = 0.00205224 Jy ! Component: 1250 - total flux cleaned = 0.00213771 Jy ! Component: 1300 - total flux cleaned = 0.00205323 Jy ! Component: 1350 - total flux cleaned = 0.00209516 Jy ! Component: 1400 - total flux cleaned = 0.00211549 Jy ! Component: 1450 - total flux cleaned = 0.00211514 Jy ! Component: 1500 - total flux cleaned = 0.00213504 Jy ! Component: 1550 - total flux cleaned = 0.00207599 Jy ! Component: 1600 - total flux cleaned = 0.00215434 Jy ! Component: 1650 - total flux cleaned = 0.00217331 Jy ! Component: 1700 - total flux cleaned = 0.00221151 Jy ! Component: 1750 - total flux cleaned = 0.0022114 Jy ! Component: 1800 - total flux cleaned = 0.00226709 Jy ! Component: 1850 - total flux cleaned = 0.00224878 Jy ! Component: 1900 - total flux cleaned = 0.00224907 Jy ! Component: 1950 - total flux cleaned = 0.00233923 Jy ! Component: 2000 - total flux cleaned = 0.00233929 Jy ! Component: 2050 - total flux cleaned = 0.00237442 Jy ! Component: 2100 - total flux cleaned = 0.00232191 Jy ! Component: 2150 - total flux cleaned = 0.00235658 Jy ! Component: 2200 - total flux cleaned = 0.00240802 Jy ! Component: 2250 - total flux cleaned = 0.00244182 Jy ! Component: 2300 - total flux cleaned = 0.00239135 Jy ! Component: 2350 - total flux cleaned = 0.00244154 Jy ! Component: 2400 - total flux cleaned = 0.00242525 Jy ! Component: 2450 - total flux cleaned = 0.00245809 Jy ! Component: 2500 - total flux cleaned = 0.0025231 Jy ! Component: 2550 - total flux cleaned = 0.00255514 Jy ! Component: 2600 - total flux cleaned = 0.00253925 Jy ! Component: 2650 - total flux cleaned = 0.00260271 Jy ! Component: 2700 - total flux cleaned = 0.00255582 Jy ! Component: 2750 - total flux cleaned = 0.00258707 Jy ! Component: 2800 - total flux cleaned = 0.00261794 Jy ! Component: 2850 - total flux cleaned = 0.00266402 Jy ! Component: 2900 - total flux cleaned = 0.00266391 Jy ! Component: 2950 - total flux cleaned = 0.00269411 Jy ! Component: 3000 - total flux cleaned = 0.00269401 Jy ! Component: 3050 - total flux cleaned = 0.0026792 Jy ! Component: 3100 - total flux cleaned = 0.00273835 Jy ! Component: 3150 - total flux cleaned = 0.002753 Jy ! Component: 3200 - total flux cleaned = 0.00273843 Jy ! Component: 3250 - total flux cleaned = 0.00276753 Jy ! Component: 3300 - total flux cleaned = 0.00272444 Jy ! Component: 3350 - total flux cleaned = 0.00273874 Jy ! Component: 3400 - total flux cleaned = 0.00276706 Jy ! Component: 3450 - total flux cleaned = 0.00276702 Jy ! Component: 3500 - total flux cleaned = 0.00276697 Jy ! Component: 3550 - total flux cleaned = 0.00278101 Jy ! Component: 3600 - total flux cleaned = 0.00280852 Jy ! Component: 3650 - total flux cleaned = 0.00276734 Jy ! Component: 3700 - total flux cleaned = 0.00278112 Jy ! Component: 3750 - total flux cleaned = 0.00279472 Jy ! Component: 3800 - total flux cleaned = 0.00279472 Jy ! Component: 3850 - total flux cleaned = 0.00279472 Jy ! Component: 3900 - total flux cleaned = 0.00282141 Jy ! Component: 3950 - total flux cleaned = 0.00278163 Jy ! Component: 4000 - total flux cleaned = 0.00276859 Jy ! Component: 4050 - total flux cleaned = 0.00276861 Jy ! Total flux subtracted in 4096 components = 0.00271643 Jy ! Clean residual min=-0.000954 max=0.000904 Jy/beam ! Clean residual mean=0.000000 rms=0.000231 Jy/beam ! Combined flux in latest and established models = 0.1247 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: 3.021 x 7.712 at -10.08 degrees (North through East) ! Clean map min=-0.00095409 max=0.11847 Jy/beam ! Clean map min=-0.0017787 max=0.11847 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 J2219+1806_S ! Writing UV FITS file: J2219+1806_S.uvf ! Writing 277 model components to file: J2219+1806_S.mod ! wwins: Wrote 3 windows to J2219+1806_S.win ! Adding 263 model components to the UV plane model. ! The established model now contains 277 components and 0.1247 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.021 x 7.712 at -10.08 degrees (North through East) ! Clean map min=-0.0016664 max=0.11845 Jy/beam ! Writing clean map to FITS file: J2219+1806_S.fits ! Writing difmap environment to: J2219+1806_S.par quit ! Quitting program ! Log file difmap.log_242 closed on Mon Nov 11 09:14:11 2019