! Started logfile: difmap.log_200 on Sun Dec 31 00:01:35 2017 ![@/tmp/difmap__28972.1] ! Obeying initialization file: /tmp/difmap__28972.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__28972.2] ![Exited script file: /tmp/difmap__28972.2] obs J0452+1403_C_a_uvm.fits ! Reading UV FITS file: J0452+1403_C_a_uvm.fits ! AN table 1: 17 integrations on 28 of 28 possible baselines. ! AN table 2: 17 integrations on 1 of 1 possible baselines. ! AN table 3: 18 integrations on 3 of 3 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J0452+1403 ! ! 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 4376 visibilities. ![Exited script file: /tmp/difmap__28972.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 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 5. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 6. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 7. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 8. ! A total of 34 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.000864Jy sigma=30.452219 ! Fit after self-cal, rms=0.986873Jy sigma=30.227190 ! 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 4 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=1.587 mas, bmaj=5.637 mas, bpa=-16.53 degrees ! Estimated noise=0.473354 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.00668088 Jy ! Component: 100 - total flux cleaned = 0.00908489 Jy ! Component: 150 - total flux cleaned = 0.0095055 Jy ! Component: 200 - total flux cleaned = 0.00955558 Jy ! Total flux subtracted in 200 components = 0.00955558 Jy ! Clean residual min=-0.001218 max=0.001220 Jy/beam ! Clean residual mean=-0.000004 rms=0.000419 Jy/beam ! Combined flux in latest and established models = 0.00955559 Jy mapl map,true ! ! Move the cursor into the plot window and press 'H' for help ! Note that 31 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.030472Jy sigma=0.798725 ! Fit after self-cal, rms=0.030472Jy sigma=0.798725 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=1.587 mas, bmaj=5.637 mas, bpa=-16.53 degrees ! Estimated noise=0.473354 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00414835 Jy ! Component: 100 - total flux cleaned = 0.00666127 Jy ! Component: 150 - total flux cleaned = 0.0081839 Jy ! Component: 200 - total flux cleaned = 0.00910871 Jy ! Component: 250 - total flux cleaned = 0.00967266 Jy ! Component: 300 - total flux cleaned = 0.0100257 Jy ! Component: 350 - total flux cleaned = 0.010249 Jy ! Component: 400 - total flux cleaned = 0.0103945 Jy ! Component: 450 - total flux cleaned = 0.0104951 Jy ! Component: 500 - total flux cleaned = 0.0105652 Jy ! Component: 550 - total flux cleaned = 0.010614 Jy ! Component: 600 - total flux cleaned = 0.010648 Jy ! Component: 650 - total flux cleaned = 0.0106717 Jy ! Component: 700 - total flux cleaned = 0.0106797 Jy ! Component: 750 - total flux cleaned = 0.0106811 Jy ! Component: 800 - total flux cleaned = 0.0106824 Jy ! Total flux subtracted in 800 components = 0.0106824 Jy ! Clean residual min=-0.001849 max=0.001573 Jy/beam ! Clean residual mean=-0.000004 rms=0.000480 Jy/beam ! Combined flux in latest and established models = 0.0106824 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 8 model components to the UV plane model. ! The established model now contains 8 components and 0.0106824 Jy ! Inverting map and beam ! Estimated beam: bmin=2.09 mas, bmaj=5.633 mas, bpa=-15.22 degrees ! Estimated noise=0.355431 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000210001 Jy ! Component: 100 - total flux cleaned = 0.000343515 Jy ! Component: 150 - total flux cleaned = 0.000431422 Jy ! Component: 200 - total flux cleaned = 0.000489301 Jy ! Component: 250 - total flux cleaned = 0.000527409 Jy ! Component: 300 - total flux cleaned = 0.000536615 Jy ! Component: 350 - total flux cleaned = 0.000537704 Jy ! Component: 400 - total flux cleaned = 0.000538792 Jy ! Component: 450 - total flux cleaned = 0.000539852 Jy ! Component: 500 - total flux cleaned = 0.000540925 Jy ! Component: 550 - total flux cleaned = 0.000541998 Jy ! Component: 600 - total flux cleaned = 0.000544123 Jy ! Component: 650 - total flux cleaned = 0.00054519 Jy ! Component: 700 - total flux cleaned = 0.000546251 Jy ! Component: 750 - total flux cleaned = 0.000547318 Jy ! Component: 800 - total flux cleaned = 0.000548354 Jy ! Component: 850 - total flux cleaned = 0.000549398 Jy ! Component: 900 - total flux cleaned = 0.000550443 Jy ! Component: 950 - total flux cleaned = 0.000552499 Jy ! Component: 1000 - total flux cleaned = 0.000553515 Jy ! Component: 1050 - total flux cleaned = 0.000554542 Jy ! Component: 1100 - total flux cleaned = 0.00055557 Jy ! Component: 1150 - total flux cleaned = 0.000556591 Jy ! Component: 1200 - total flux cleaned = 0.000557608 Jy ! Total flux subtracted in 1200 components = 0.000557608 Jy ! Clean residual min=-0.001810 max=0.002008 Jy/beam ! Clean residual mean=0.000001 rms=0.000531 Jy/beam ! Combined flux in latest and established models = 0.01124 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 3 model components to the UV plane model. ! The established model now contains 8 components and 0.01124 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.027411Jy sigma=0.724132 ! Fit after self-cal, rms=0.026005Jy sigma=0.667657 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 3. ! ! Correcting IF 2. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 5. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 6. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 7. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 8. ! A total of 54 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.026005Jy sigma=0.667657 ! Fit after self-cal, rms=0.018892Jy sigma=0.651956 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.022105Jy sigma=0.794088 ! Fit after self-cal, rms=0.022105Jy sigma=0.794088 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=1.623 mas, bmaj=5.65 mas, bpa=-16.4 degrees ! Estimated noise=0.377149 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00430018 Jy ! Component: 100 - total flux cleaned = 0.00690337 Jy ! Component: 150 - total flux cleaned = 0.00848063 Jy ! Component: 200 - total flux cleaned = 0.00943861 Jy ! Component: 250 - total flux cleaned = 0.0100215 Jy ! Component: 300 - total flux cleaned = 0.0103761 Jy ! Component: 350 - total flux cleaned = 0.0105918 Jy ! Component: 400 - total flux cleaned = 0.0107231 Jy ! Component: 450 - total flux cleaned = 0.0108031 Jy ! Component: 500 - total flux cleaned = 0.0108519 Jy ! Component: 550 - total flux cleaned = 0.0108818 Jy ! Component: 600 - total flux cleaned = 0.0109007 Jy ! Component: 650 - total flux cleaned = 0.0109133 Jy ! Component: 700 - total flux cleaned = 0.0109217 Jy ! Component: 750 - total flux cleaned = 0.0109273 Jy ! Component: 800 - total flux cleaned = 0.0109311 Jy ! Total flux subtracted in 800 components = 0.0109311 Jy ! Clean residual min=-0.000964 max=0.001075 Jy/beam ! Clean residual mean=-0.000000 rms=0.000275 Jy/beam ! Combined flux in latest and established models = 0.0109311 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 8 model components to the UV plane model. ! The established model now contains 8 components and 0.0109311 Jy ! Inverting map and beam ! Estimated beam: bmin=2.163 mas, bmaj=5.568 mas, bpa=-14.46 degrees ! Estimated noise=0.279899 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 9.88988e-05 Jy ! Component: 100 - total flux cleaned = 0.000158733 Jy ! Component: 150 - total flux cleaned = 0.000194933 Jy ! Component: 200 - total flux cleaned = 0.000202165 Jy ! Component: 250 - total flux cleaned = 0.000202123 Jy ! Component: 300 - total flux cleaned = 0.000202097 Jy ! Component: 350 - total flux cleaned = 0.000202084 Jy ! Component: 400 - total flux cleaned = 0.000202079 Jy ! Component: 450 - total flux cleaned = 0.000202079 Jy ! Component: 500 - total flux cleaned = 0.000202082 Jy ! Component: 550 - total flux cleaned = 0.000202085 Jy ! Component: 600 - total flux cleaned = 0.000202089 Jy ! Component: 650 - total flux cleaned = 0.000202093 Jy ! Component: 700 - total flux cleaned = 0.000202097 Jy ! Component: 750 - total flux cleaned = 0.000202101 Jy ! Component: 800 - total flux cleaned = 0.000202104 Jy ! Component: 850 - total flux cleaned = 0.000202108 Jy ! Component: 900 - total flux cleaned = 0.000202111 Jy ! Component: 950 - total flux cleaned = 0.000202114 Jy ! Component: 1000 - total flux cleaned = 0.000202117 Jy ! Component: 1050 - total flux cleaned = 0.00020212 Jy ! Component: 1100 - total flux cleaned = 0.000202123 Jy ! Component: 1150 - total flux cleaned = 0.000202126 Jy ! Component: 1200 - total flux cleaned = 0.000202129 Jy ! Total flux subtracted in 1200 components = 0.000202129 Jy ! Clean residual min=-0.000794 max=0.000877 Jy/beam ! Clean residual mean=0.000001 rms=0.000212 Jy/beam ! Combined flux in latest and established models = 0.0111333 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 proj ! Adding 5 model components to the UV plane model. ! The established model now contains 12 components and 0.0111333 Jy ! Using default options string "m3" ! 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 ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.163 x 5.568 at -14.46 degrees (North through East) ! Clean map min=-0.00078326 max=0.011093 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.000183325 Jy ! Component: 100 - total flux cleaned = 0.00028643 Jy ! Component: 150 - total flux cleaned = 0.000240549 Jy ! Component: 200 - total flux cleaned = 0.000218673 Jy ! Component: 250 - total flux cleaned = 0.000218139 Jy ! Component: 300 - total flux cleaned = 0.00021835 Jy ! Component: 350 - total flux cleaned = 0.000282353 Jy ! Component: 400 - total flux cleaned = 0.000308683 Jy ! Component: 450 - total flux cleaned = 0.000376172 Jy ! Component: 500 - total flux cleaned = 0.000440697 Jy ! Total flux subtracted in 512 components = 0.000448581 Jy ! Clean residual min=-0.000392 max=0.000458 Jy/beam ! Clean residual mean=0.000001 rms=0.000137 Jy/beam ! Combined flux in latest and established models = 0.0115819 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.163 x 5.568 at -14.46 degrees (North through East) ! Clean map min=-0.00038843 max=0.011142 Jy/beam ! Clean map min=-0.00068173 max=0.011142 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 J0452+1403_C_a ! Writing UV FITS file: J0452+1403_C_a.uvf ! Writing 72 model components to file: J0452+1403_C_a.mod ! wwins: Wrote 2 windows to J0452+1403_C_a.win ! Adding 60 model components to the UV plane model. ! The established model now contains 72 components and 0.0115819 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.163 x 5.568 at -14.46 degrees (North through East) ! Clean map min=-0.00068914 max=0.011145 Jy/beam ! Writing clean map to FITS file: J0452+1403_C_a.fits ! Writing difmap environment to: J0452+1403_C_a.par quit ! Quitting program ! Log file difmap.log_200 closed on Sun Dec 31 00:02:01 2017