! Started logfile: difmap.log_23 on Sat Jul 28 10:35:33 2018 ![@/tmp/difmap__34289.1] ! Obeying initialization file: /tmp/difmap__34289.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__34289.2] ![Exited script file: /tmp/difmap__34289.2] obs J0922-3959_X_uva.fits ! Reading UV FITS file: J0922-3959_X_uva.fits ! AN table 1: 15 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.32381 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0922-3959 ! ! 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 8.36575e+09 8e+06 1 8e+06 ! 02 2 8.44575e+09 8e+06 1 8e+06 ! 03 3 8.80575e+09 8e+06 1 8e+06 ! 04 4 8.92575e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 74 lines of history. ! ! Reading 544 visibilities. ![Exited script file: /tmp/difmap__34289.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.1 taper_size = 70. mapsize field_size,field_cell ! Map grid = 1024x1024 pixels with 0.100x0.100 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 76 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 76 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 76 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 76 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.297132Jy sigma=18.066639 ! Fit after self-cal, rms=0.375119Jy sigma=5.488460 ! 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 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=0.6844 mas, bmaj=2.246 mas, bpa=-8.324 degrees ! Estimated noise=3.25905 mJy/beam. ! clean: niter=200 gain=0.02 cutoff=0 ! Component: 050 - total flux cleaned = 0.426496 Jy ! Component: 100 - total flux cleaned = 0.636817 Jy ! Component: 150 - total flux cleaned = 0.748376 Jy ! Component: 200 - total flux cleaned = 0.814779 Jy ! Total flux subtracted in 200 components = 0.814779 Jy ! Clean residual min=-0.027950 max=0.055447 Jy/beam ! Clean residual mean=0.000720 rms=0.010177 Jy/beam ! Combined flux in latest and established models = 0.814779 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. ! ! Fit before self-cal, rms=0.831818Jy sigma=11.343837 ! Fit after self-cal, rms=0.831818Jy sigma=11.343837 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=0.6844 mas, bmaj=2.246 mas, bpa=-8.324 degrees ! Estimated noise=3.25905 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.257311 Jy ! Component: 100 - total flux cleaned = 0.425087 Jy ! Component: 150 - total flux cleaned = 0.546712 Jy ! Component: 200 - total flux cleaned = 0.63527 Jy ! Component: 250 - total flux cleaned = 0.699773 Jy ! Component: 300 - total flux cleaned = 0.747042 Jy ! Component: 350 - total flux cleaned = 0.782687 Jy ! Component: 400 - total flux cleaned = 0.81106 Jy ! Component: 450 - total flux cleaned = 0.834624 Jy ! Component: 500 - total flux cleaned = 0.854395 Jy ! Component: 550 - total flux cleaned = 0.871196 Jy ! Component: 600 - total flux cleaned = 0.885565 Jy ! Component: 650 - total flux cleaned = 0.897975 Jy ! Component: 700 - total flux cleaned = 0.908798 Jy ! Component: 750 - total flux cleaned = 0.918241 Jy ! Component: 800 - total flux cleaned = 0.926484 Jy ! Total flux subtracted in 800 components = 0.926484 Jy ! Clean residual min=-0.022188 max=0.043386 Jy/beam ! Clean residual mean=0.000357 rms=0.007253 Jy/beam ! Combined flux in latest and established models = 0.926484 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 52 model components to the UV plane model. ! The established model now contains 52 components and 0.926484 Jy ! Inverting map and beam ! Estimated beam: bmin=0.858 mas, bmaj=2.909 mas, bpa=-11.6 degrees ! Estimated noise=2.44516 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0175051 Jy ! Component: 100 - total flux cleaned = 0.030466 Jy ! Component: 150 - total flux cleaned = 0.0402754 Jy ! Component: 200 - total flux cleaned = 0.0479657 Jy ! Component: 250 - total flux cleaned = 0.0540083 Jy ! Component: 300 - total flux cleaned = 0.0587561 Jy ! Component: 350 - total flux cleaned = 0.0624866 Jy ! Component: 400 - total flux cleaned = 0.0654179 Jy ! Component: 450 - total flux cleaned = 0.0677213 Jy ! Component: 500 - total flux cleaned = 0.0695312 Jy ! Component: 550 - total flux cleaned = 0.0709542 Jy ! Component: 600 - total flux cleaned = 0.0714801 Jy ! Component: 650 - total flux cleaned = 0.0716153 Jy ! Component: 700 - total flux cleaned = 0.0717038 Jy ! Component: 750 - total flux cleaned = 0.0717906 Jy ! Component: 800 - total flux cleaned = 0.0718761 Jy ! Component: 850 - total flux cleaned = 0.0720027 Jy ! Component: 900 - total flux cleaned = 0.0720858 Jy ! Component: 950 - total flux cleaned = 0.0721677 Jy ! Component: 1000 - total flux cleaned = 0.0722482 Jy ! Component: 1050 - total flux cleaned = 0.0723276 Jy ! Component: 1100 - total flux cleaned = 0.0724441 Jy ! Component: 1150 - total flux cleaned = 0.0725209 Jy ! Component: 1200 - total flux cleaned = 0.0725957 Jy ! Total flux subtracted in 1200 components = 0.0725957 Jy ! Clean residual min=-0.016032 max=0.029409 Jy/beam ! Clean residual mean=0.000383 rms=0.005626 Jy/beam ! Combined flux in latest and established models = 0.999079 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 15 model components to the UV plane model. ! The established model now contains 57 components and 0.999079 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093915Jy sigma=1.214822 ! Fit after self-cal, rms=0.089000Jy sigma=1.144851 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.089000Jy sigma=1.144851 ! Fit after self-cal, rms=0.074548Jy sigma=0.948165 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.838062Jy sigma=11.343837 ! Fit after self-cal, rms=0.838062Jy sigma=11.343837 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=0.6899 mas, bmaj=2.252 mas, bpa=-7.965 degrees ! Estimated noise=3.39637 mJy/beam. ! clean: niter=800 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.267743 Jy ! Component: 100 - total flux cleaned = 0.44155 Jy ! Component: 150 - total flux cleaned = 0.566795 Jy ! Component: 200 - total flux cleaned = 0.657489 Jy ! Component: 250 - total flux cleaned = 0.72325 Jy ! Component: 300 - total flux cleaned = 0.770977 Jy ! Component: 350 - total flux cleaned = 0.80639 Jy ! Component: 400 - total flux cleaned = 0.834496 Jy ! Component: 450 - total flux cleaned = 0.857586 Jy ! Component: 500 - total flux cleaned = 0.876699 Jy ! Component: 550 - total flux cleaned = 0.892668 Jy ! Component: 600 - total flux cleaned = 0.906046 Jy ! Component: 650 - total flux cleaned = 0.917284 Jy ! Component: 700 - total flux cleaned = 0.926747 Jy ! Component: 750 - total flux cleaned = 0.934739 Jy ! Component: 800 - total flux cleaned = 0.941497 Jy ! Total flux subtracted in 800 components = 0.941497 Jy ! Clean residual min=-0.015797 max=0.021398 Jy/beam ! Clean residual mean=0.000309 rms=0.004415 Jy/beam ! Combined flux in latest and established models = 0.941497 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 63 model components to the UV plane model. ! The established model now contains 63 components and 0.941497 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8927 mas, bmaj=2.965 mas, bpa=-10.97 degrees ! Estimated noise=2.48798 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0117398 Jy ! Component: 100 - total flux cleaned = 0.0204363 Jy ! Component: 150 - total flux cleaned = 0.0269736 Jy ! Component: 200 - total flux cleaned = 0.031919 Jy ! Component: 250 - total flux cleaned = 0.0357077 Jy ! Component: 300 - total flux cleaned = 0.0386185 Jy ! Component: 350 - total flux cleaned = 0.0408549 Jy ! Component: 400 - total flux cleaned = 0.0425749 Jy ! Component: 450 - total flux cleaned = 0.0439029 Jy ! Component: 500 - total flux cleaned = 0.0449301 Jy ! Component: 550 - total flux cleaned = 0.045725 Jy ! Component: 600 - total flux cleaned = 0.0463402 Jy ! Component: 650 - total flux cleaned = 0.046603 Jy ! Component: 700 - total flux cleaned = 0.0466595 Jy ! Component: 750 - total flux cleaned = 0.0467148 Jy ! Component: 800 - total flux cleaned = 0.0467868 Jy ! Component: 850 - total flux cleaned = 0.0468397 Jy ! Component: 900 - total flux cleaned = 0.0468906 Jy ! Component: 950 - total flux cleaned = 0.0469413 Jy ! Component: 1000 - total flux cleaned = 0.0469906 Jy ! Component: 1050 - total flux cleaned = 0.0470391 Jy ! Component: 1100 - total flux cleaned = 0.0470864 Jy ! Component: 1150 - total flux cleaned = 0.0471329 Jy ! Component: 1200 - total flux cleaned = 0.0471784 Jy ! Total flux subtracted in 1200 components = 0.0471784 Jy ! Clean residual min=-0.010808 max=0.013468 Jy/beam ! Clean residual mean=0.000525 rms=0.002937 Jy/beam ! Combined flux in latest and established models = 0.988675 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: 0.8927 x 2.965 at -10.97 degrees (North through East) ! Clean map min=-0.010808 max=0.69899 Jy/beam ! Clean map min=-0.010808 max=0.71379 Jy/beam ! ! Move the cursor into the plot window and press 'H' for help mapcol grey ! Mapplot colormap: grey, contrast: 1 brightness: 0.5. 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.838062Jy sigma=11.343837 ! Fit after self-cal, rms=0.838062Jy sigma=11.343837 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 1200, 0.01 ! Inverting map and beam ! Estimated beam: bmin=0.6899 mas, bmaj=2.252 mas, bpa=-7.965 degrees ! Estimated noise=3.39637 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.267743 Jy ! Component: 100 - total flux cleaned = 0.44155 Jy ! Component: 150 - total flux cleaned = 0.566795 Jy ! Component: 200 - total flux cleaned = 0.657489 Jy ! Component: 250 - total flux cleaned = 0.72325 Jy ! Component: 300 - total flux cleaned = 0.770977 Jy ! Component: 350 - total flux cleaned = 0.80639 Jy ! Component: 400 - total flux cleaned = 0.834496 Jy ! Component: 450 - total flux cleaned = 0.857586 Jy ! Component: 500 - total flux cleaned = 0.876699 Jy ! Component: 550 - total flux cleaned = 0.892668 Jy ! Component: 600 - total flux cleaned = 0.906046 Jy ! Component: 650 - total flux cleaned = 0.917284 Jy ! Component: 700 - total flux cleaned = 0.926747 Jy ! Component: 750 - total flux cleaned = 0.934739 Jy ! Component: 800 - total flux cleaned = 0.941497 Jy ! Component: 850 - total flux cleaned = 0.947238 Jy ! Component: 900 - total flux cleaned = 0.952144 Jy ! Component: 950 - total flux cleaned = 0.956344 Jy ! Component: 1000 - total flux cleaned = 0.959942 Jy ! Component: 1050 - total flux cleaned = 0.963031 Jy ! Component: 1100 - total flux cleaned = 0.965686 Jy ! Component: 1150 - total flux cleaned = 0.967967 Jy ! Component: 1200 - total flux cleaned = 0.969929 Jy ! Total flux subtracted in 1200 components = 0.969929 Jy ! Clean residual min=-0.015419 max=0.016789 Jy/beam ! Clean residual mean=0.000215 rms=0.003986 Jy/beam ! Combined flux in latest and established models = 0.969929 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 2000, 0.01 ! Adding 68 model components to the UV plane model. ! The established model now contains 68 components and 0.969929 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8927 mas, bmaj=2.965 mas, bpa=-10.97 degrees ! Estimated noise=2.48798 mJy/beam. ! clean: niter=2000 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00512958 Jy ! Component: 100 - total flux cleaned = 0.00883083 Jy ! Component: 150 - total flux cleaned = 0.0115568 Jy ! Component: 200 - total flux cleaned = 0.0136034 Jy ! Component: 250 - total flux cleaned = 0.0151593 Jy ! Component: 300 - total flux cleaned = 0.0163431 Jy ! Component: 350 - total flux cleaned = 0.0172478 Jy ! Component: 400 - total flux cleaned = 0.0179407 Jy ! Component: 450 - total flux cleaned = 0.0184339 Jy ! Component: 500 - total flux cleaned = 0.0185077 Jy ! Component: 550 - total flux cleaned = 0.0185435 Jy ! Component: 600 - total flux cleaned = 0.0185964 Jy ! Component: 650 - total flux cleaned = 0.0186482 Jy ! Component: 700 - total flux cleaned = 0.0186988 Jy ! Component: 750 - total flux cleaned = 0.0187645 Jy ! Component: 800 - total flux cleaned = 0.0188124 Jy ! Component: 850 - total flux cleaned = 0.0188594 Jy ! Component: 900 - total flux cleaned = 0.0189053 Jy ! Component: 950 - total flux cleaned = 0.0189503 Jy ! Component: 1000 - total flux cleaned = 0.0189943 Jy ! Component: 1050 - total flux cleaned = 0.0190374 Jy ! Component: 1100 - total flux cleaned = 0.01908 Jy ! Component: 1150 - total flux cleaned = 0.0191216 Jy ! Component: 1200 - total flux cleaned = 0.0191621 Jy ! Component: 1250 - total flux cleaned = 0.0192023 Jy ! Component: 1300 - total flux cleaned = 0.0192417 Jy ! Component: 1350 - total flux cleaned = 0.0192804 Jy ! Component: 1400 - total flux cleaned = 0.0193056 Jy ! Component: 1450 - total flux cleaned = 0.0193435 Jy ! Component: 1500 - total flux cleaned = 0.0193681 Jy ! Component: 1550 - total flux cleaned = 0.0193925 Jy ! Component: 1600 - total flux cleaned = 0.0194289 Jy ! Component: 1650 - total flux cleaned = 0.0194526 Jy ! Component: 1700 - total flux cleaned = 0.0194763 Jy ! Component: 1750 - total flux cleaned = 0.0195115 Jy ! Component: 1800 - total flux cleaned = 0.0195344 Jy ! Component: 1850 - total flux cleaned = 0.0195574 Jy ! Component: 1900 - total flux cleaned = 0.0195798 Jy ! Component: 1950 - total flux cleaned = 0.0195911 Jy ! Component: 2000 - total flux cleaned = 0.0196133 Jy ! Total flux subtracted in 2000 components = 0.0196133 Jy ! Clean residual min=-0.010735 max=0.013436 Jy/beam ! Clean residual mean=0.000518 rms=0.002927 Jy/beam ! Combined flux in latest and established models = 0.989542 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 29 model components to the UV plane model. ! The established model now contains 86 components and 0.989542 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.071954Jy sigma=0.921306 ! Fit after self-cal, rms=0.070740Jy sigma=0.911600 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.070740Jy sigma=0.911600 ! Fit after self-cal, rms=0.065733Jy sigma=0.837795 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.840158Jy sigma=11.343837 ! Fit after self-cal, rms=0.840158Jy sigma=11.343837 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 1200, 0.01 ! Inverting map and beam ! Estimated beam: bmin=0.6882 mas, bmaj=2.238 mas, bpa=-8.173 degrees ! Estimated noise=3.40885 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.268434 Jy ! Component: 100 - total flux cleaned = 0.443798 Jy ! Component: 150 - total flux cleaned = 0.570222 Jy ! Component: 200 - total flux cleaned = 0.661786 Jy ! Component: 250 - total flux cleaned = 0.728179 Jy ! Component: 300 - total flux cleaned = 0.77641 Jy ! Component: 350 - total flux cleaned = 0.812251 Jy ! Component: 400 - total flux cleaned = 0.840522 Jy ! Component: 450 - total flux cleaned = 0.863654 Jy ! Component: 500 - total flux cleaned = 0.882764 Jy ! Component: 550 - total flux cleaned = 0.89867 Jy ! Component: 600 - total flux cleaned = 0.911956 Jy ! Component: 650 - total flux cleaned = 0.923071 Jy ! Component: 700 - total flux cleaned = 0.932389 Jy ! Component: 750 - total flux cleaned = 0.94021 Jy ! Component: 800 - total flux cleaned = 0.946801 Jy ! Component: 850 - total flux cleaned = 0.952366 Jy ! Component: 900 - total flux cleaned = 0.957068 Jy ! Component: 950 - total flux cleaned = 0.961055 Jy ! Component: 1000 - total flux cleaned = 0.964443 Jy ! Component: 1050 - total flux cleaned = 0.967327 Jy ! Component: 1100 - total flux cleaned = 0.969791 Jy ! Component: 1150 - total flux cleaned = 0.971896 Jy ! Component: 1200 - total flux cleaned = 0.973696 Jy ! Total flux subtracted in 1200 components = 0.973696 Jy ! Clean residual min=-0.012770 max=0.013475 Jy/beam ! Clean residual mean=0.000136 rms=0.003262 Jy/beam ! Combined flux in latest and established models = 0.973696 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 2000, 0.01 ! Adding 71 model components to the UV plane model. ! The established model now contains 71 components and 0.973696 Jy ! Inverting map and beam ! Estimated beam: bmin=0.894 mas, bmaj=2.951 mas, bpa=-11.27 degrees ! Estimated noise=2.48846 mJy/beam. ! clean: niter=2000 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00361096 Jy ! Component: 100 - total flux cleaned = 0.00612144 Jy ! Component: 150 - total flux cleaned = 0.00794128 Jy ! Component: 200 - total flux cleaned = 0.00930729 Jy ! Component: 250 - total flux cleaned = 0.0103349 Jy ! Component: 300 - total flux cleaned = 0.0111091 Jy ! Component: 350 - total flux cleaned = 0.0116925 Jy ! Component: 400 - total flux cleaned = 0.0121322 Jy ! Component: 450 - total flux cleaned = 0.0122205 Jy ! Component: 500 - total flux cleaned = 0.0122624 Jy ! Component: 550 - total flux cleaned = 0.0122897 Jy ! Component: 600 - total flux cleaned = 0.0123423 Jy ! Component: 650 - total flux cleaned = 0.0123678 Jy ! Component: 700 - total flux cleaned = 0.0124049 Jy ! Component: 750 - total flux cleaned = 0.0124408 Jy ! Component: 800 - total flux cleaned = 0.0124756 Jy ! Component: 850 - total flux cleaned = 0.0124981 Jy ! Component: 900 - total flux cleaned = 0.0125308 Jy ! Component: 950 - total flux cleaned = 0.0125622 Jy ! Component: 1000 - total flux cleaned = 0.0125826 Jy ! Component: 1050 - total flux cleaned = 0.0126026 Jy ! Component: 1100 - total flux cleaned = 0.0126317 Jy ! Component: 1150 - total flux cleaned = 0.0126505 Jy ! Component: 1200 - total flux cleaned = 0.012669 Jy ! Component: 1250 - total flux cleaned = 0.0126959 Jy ! Component: 1300 - total flux cleaned = 0.0127133 Jy ! Component: 1350 - total flux cleaned = 0.0127304 Jy ! Component: 1400 - total flux cleaned = 0.0127389 Jy ! Component: 1450 - total flux cleaned = 0.0127389 Jy ! Component: 1500 - total flux cleaned = 0.0127392 Jy ! Component: 1550 - total flux cleaned = 0.0127471 Jy ! Component: 1600 - total flux cleaned = 0.0127551 Jy ! Component: 1650 - total flux cleaned = 0.0127552 Jy ! Component: 1700 - total flux cleaned = 0.0127632 Jy ! Component: 1750 - total flux cleaned = 0.0127709 Jy ! Component: 1800 - total flux cleaned = 0.0127709 Jy ! Component: 1850 - total flux cleaned = 0.0127788 Jy ! Component: 1900 - total flux cleaned = 0.0127864 Jy ! Component: 1950 - total flux cleaned = 0.012794 Jy ! Component: 2000 - total flux cleaned = 0.012794 Jy ! Total flux subtracted in 2000 components = 0.012794 Jy ! Clean residual min=-0.008518 max=0.009887 Jy/beam ! Clean residual mean=0.000318 rms=0.002202 Jy/beam ! Combined flux in latest and established models = 0.98649 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 29 model components to the UV plane model. ! The established model now contains 93 components and 0.98649 Jy ! Using default options string "m3" ! Move the cursor into the plot window and press 'H' for help ! Applying 2 buffered edits. uncalib selfcal ! Performing phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063481Jy sigma=0.801076 ! Fit after self-cal, rms=0.063048Jy sigma=0.793837 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.063048Jy sigma=0.793837 ! Fit after self-cal, rms=0.063246Jy sigma=0.788998 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.842117Jy sigma=11.350198 ! Fit after self-cal, rms=0.842117Jy sigma=11.350198 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 1200, 0.01 ! Inverting map and beam ! Estimated beam: bmin=0.6907 mas, bmaj=2.243 mas, bpa=-8.069 degrees ! Estimated noise=3.4442 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.270299 Jy ! Component: 100 - total flux cleaned = 0.446619 Jy ! Component: 150 - total flux cleaned = 0.57352 Jy ! Component: 200 - total flux cleaned = 0.665292 Jy ! Component: 250 - total flux cleaned = 0.731741 Jy ! Component: 300 - total flux cleaned = 0.779907 Jy ! Component: 350 - total flux cleaned = 0.815589 Jy ! Component: 400 - total flux cleaned = 0.843691 Jy ! Component: 450 - total flux cleaned = 0.866631 Jy ! Component: 500 - total flux cleaned = 0.885533 Jy ! Component: 550 - total flux cleaned = 0.901222 Jy ! Component: 600 - total flux cleaned = 0.914282 Jy ! Component: 650 - total flux cleaned = 0.925177 Jy ! Component: 700 - total flux cleaned = 0.934278 Jy ! Component: 750 - total flux cleaned = 0.941894 Jy ! Component: 800 - total flux cleaned = 0.948281 Jy ! Component: 850 - total flux cleaned = 0.953649 Jy ! Component: 900 - total flux cleaned = 0.958171 Jy ! Component: 950 - total flux cleaned = 0.961983 Jy ! Component: 1000 - total flux cleaned = 0.965199 Jy ! Component: 1050 - total flux cleaned = 0.967924 Jy ! Component: 1100 - total flux cleaned = 0.970234 Jy ! Component: 1150 - total flux cleaned = 0.972198 Jy ! Component: 1200 - total flux cleaned = 0.97387 Jy ! Total flux subtracted in 1200 components = 0.97387 Jy ! Clean residual min=-0.012045 max=0.012169 Jy/beam ! Clean residual mean=0.000128 rms=0.003138 Jy/beam ! Combined flux in latest and established models = 0.97387 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 2000, 0.01 ! Adding 67 model components to the UV plane model. ! The established model now contains 67 components and 0.97387 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9091 mas, bmaj=2.963 mas, bpa=-11.14 degrees ! Estimated noise=2.50637 mJy/beam. ! clean: niter=2000 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00318663 Jy ! Component: 100 - total flux cleaned = 0.00540803 Jy ! Component: 150 - total flux cleaned = 0.00703361 Jy ! Component: 200 - total flux cleaned = 0.00823474 Jy ! Component: 250 - total flux cleaned = 0.00912467 Jy ! Component: 300 - total flux cleaned = 0.00978917 Jy ! Component: 350 - total flux cleaned = 0.0102869 Jy ! Component: 400 - total flux cleaned = 0.0106598 Jy ! Component: 450 - total flux cleaned = 0.0107103 Jy ! Component: 500 - total flux cleaned = 0.0107465 Jy ! Component: 550 - total flux cleaned = 0.0107815 Jy ! Component: 600 - total flux cleaned = 0.0108153 Jy ! Component: 650 - total flux cleaned = 0.0108479 Jy ! Component: 700 - total flux cleaned = 0.0108792 Jy ! Component: 750 - total flux cleaned = 0.0109096 Jy ! Component: 800 - total flux cleaned = 0.0109389 Jy ! Component: 850 - total flux cleaned = 0.0109672 Jy ! Component: 900 - total flux cleaned = 0.0109946 Jy ! Component: 950 - total flux cleaned = 0.0110035 Jy ! Component: 1000 - total flux cleaned = 0.0110035 Jy ! Component: 1050 - total flux cleaned = 0.0110121 Jy ! Component: 1100 - total flux cleaned = 0.011012 Jy ! Component: 1150 - total flux cleaned = 0.0110204 Jy ! Component: 1200 - total flux cleaned = 0.0110204 Jy ! Component: 1250 - total flux cleaned = 0.0110287 Jy ! Component: 1300 - total flux cleaned = 0.0110285 Jy ! Component: 1350 - total flux cleaned = 0.0110287 Jy ! Component: 1400 - total flux cleaned = 0.0110286 Jy ! Component: 1450 - total flux cleaned = 0.0110362 Jy ! Component: 1500 - total flux cleaned = 0.0110363 Jy ! Component: 1550 - total flux cleaned = 0.0110363 Jy ! Component: 1600 - total flux cleaned = 0.0110438 Jy ! Component: 1650 - total flux cleaned = 0.0110438 Jy ! Component: 1700 - total flux cleaned = 0.0110438 Jy ! Component: 1750 - total flux cleaned = 0.0110439 Jy ! Component: 1800 - total flux cleaned = 0.0110509 Jy ! Component: 1850 - total flux cleaned = 0.0110509 Jy ! Component: 1900 - total flux cleaned = 0.011051 Jy ! Component: 1950 - total flux cleaned = 0.0110577 Jy ! Component: 2000 - total flux cleaned = 0.0110578 Jy ! Total flux subtracted in 2000 components = 0.0110578 Jy ! Clean residual min=-0.007538 max=0.007647 Jy/beam ! Clean residual mean=0.000304 rms=0.002044 Jy/beam ! Combined flux in latest and established models = 0.984928 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. ! ! Fit before self-cal, rms=0.842117Jy sigma=11.350198 ! Fit after self-cal, rms=0.842117Jy sigma=11.350198 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 1200, 0.01 ! Inverting map and beam ! Estimated beam: bmin=0.6907 mas, bmaj=2.243 mas, bpa=-8.069 degrees ! Estimated noise=3.4442 mJy/beam. ! clean: niter=1200 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.270299 Jy ! Component: 100 - total flux cleaned = 0.446619 Jy ! Component: 150 - total flux cleaned = 0.57352 Jy ! Component: 200 - total flux cleaned = 0.665292 Jy ! Component: 250 - total flux cleaned = 0.731741 Jy ! Component: 300 - total flux cleaned = 0.779907 Jy ! Component: 350 - total flux cleaned = 0.815589 Jy ! Component: 400 - total flux cleaned = 0.843688 Jy ! Component: 450 - total flux cleaned = 0.866452 Jy ! Component: 500 - total flux cleaned = 0.885184 Jy ! Component: 550 - total flux cleaned = 0.900759 Jy ! Component: 600 - total flux cleaned = 0.913766 Jy ! Component: 650 - total flux cleaned = 0.924642 Jy ! Component: 700 - total flux cleaned = 0.933752 Jy ! Component: 750 - total flux cleaned = 0.941404 Jy ! Component: 800 - total flux cleaned = 0.94787 Jy ! Component: 850 - total flux cleaned = 0.953378 Jy ! Component: 900 - total flux cleaned = 0.95809 Jy ! Component: 950 - total flux cleaned = 0.962127 Jy ! Component: 1000 - total flux cleaned = 0.965586 Jy ! Component: 1050 - total flux cleaned = 0.968553 Jy ! Component: 1100 - total flux cleaned = 0.971109 Jy ! Component: 1150 - total flux cleaned = 0.973313 Jy ! Component: 1200 - total flux cleaned = 0.975213 Jy ! Total flux subtracted in 1200 components = 0.975213 Jy ! Clean residual min=-0.011818 max=0.012259 Jy/beam ! Clean residual mean=0.000126 rms=0.003162 Jy/beam ! Combined flux in latest and established models = 0.975214 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 2000, 0.01 ! Adding 86 model components to the UV plane model. ! The established model now contains 86 components and 0.975214 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9091 mas, bmaj=2.963 mas, bpa=-11.14 degrees ! Estimated noise=2.50637 mJy/beam. ! clean: niter=2000 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0032881 Jy ! Component: 100 - total flux cleaned = 0.00563932 Jy ! Component: 150 - total flux cleaned = 0.00733228 Jy ! Component: 200 - total flux cleaned = 0.00855549 Jy ! Component: 250 - total flux cleaned = 0.00945182 Jy ! Component: 300 - total flux cleaned = 0.0101098 Jy ! Component: 350 - total flux cleaned = 0.0103209 Jy ! Component: 400 - total flux cleaned = 0.0103591 Jy ! Component: 450 - total flux cleaned = 0.0103589 Jy ! Component: 500 - total flux cleaned = 0.010359 Jy ! Component: 550 - total flux cleaned = 0.0103592 Jy ! Component: 600 - total flux cleaned = 0.0103589 Jy ! Component: 650 - total flux cleaned = 0.0103425 Jy ! Component: 700 - total flux cleaned = 0.0103426 Jy ! Component: 750 - total flux cleaned = 0.0103422 Jy ! Component: 800 - total flux cleaned = 0.0103427 Jy ! Component: 850 - total flux cleaned = 0.0103426 Jy ! Component: 900 - total flux cleaned = 0.0103276 Jy ! Component: 950 - total flux cleaned = 0.0103278 Jy ! Component: 1000 - total flux cleaned = 0.0103277 Jy ! Component: 1050 - total flux cleaned = 0.0103274 Jy ! Component: 1100 - total flux cleaned = 0.0103275 Jy ! Component: 1150 - total flux cleaned = 0.0103138 Jy ! Component: 1200 - total flux cleaned = 0.0103139 Jy ! Component: 1250 - total flux cleaned = 0.0103139 Jy ! Component: 1300 - total flux cleaned = 0.010314 Jy ! Component: 1350 - total flux cleaned = 0.0103141 Jy ! Component: 1400 - total flux cleaned = 0.0103138 Jy ! Component: 1450 - total flux cleaned = 0.0103019 Jy ! Component: 1500 - total flux cleaned = 0.0103018 Jy ! Component: 1550 - total flux cleaned = 0.0103021 Jy ! Component: 1600 - total flux cleaned = 0.0103023 Jy ! Component: 1650 - total flux cleaned = 0.0103026 Jy ! Component: 1700 - total flux cleaned = 0.0103023 Jy ! Component: 1750 - total flux cleaned = 0.0103026 Jy ! Component: 1800 - total flux cleaned = 0.0103028 Jy ! Component: 1850 - total flux cleaned = 0.0103132 Jy ! Component: 1900 - total flux cleaned = 0.0103135 Jy ! Component: 1950 - total flux cleaned = 0.0103137 Jy ! Component: 2000 - total flux cleaned = 0.0103137 Jy ! Total flux subtracted in 2000 components = 0.0103137 Jy ! Clean residual min=-0.007431 max=0.007545 Jy/beam ! Clean residual mean=0.000301 rms=0.002063 Jy/beam ! Combined flux in latest and established models = 0.985528 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 19 model components to the UV plane model. ! The established model now contains 98 components and 0.985528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063337Jy sigma=0.789822 ! Fit after self-cal, rms=0.063273Jy sigma=0.788823 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.063273Jy sigma=0.788823 ! Fit after self-cal, rms=0.063387Jy sigma=0.788008 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=0.9134 mas, bmaj=2.977 mas, bpa=-11.09 degrees ! Estimated noise=2.51061 mJy/beam. ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00183539 Jy ! Component: 100 - total flux cleaned = 0.00350161 Jy ! Component: 150 - total flux cleaned = 0.00482622 Jy ! Component: 200 - total flux cleaned = 0.00619812 Jy ! Component: 250 - total flux cleaned = 0.00761701 Jy ! Component: 300 - total flux cleaned = 0.0092912 Jy ! Component: 350 - total flux cleaned = 0.0106988 Jy ! Component: 400 - total flux cleaned = 0.0120556 Jy ! Component: 450 - total flux cleaned = 0.0130823 Jy ! Component: 500 - total flux cleaned = 0.0141638 Jy ! Component: 550 - total flux cleaned = 0.0151221 Jy ! Component: 600 - total flux cleaned = 0.0161337 Jy ! Component: 650 - total flux cleaned = 0.0170297 Jy ! Component: 700 - total flux cleaned = 0.0177411 Jy ! Component: 750 - total flux cleaned = 0.0185868 Jy ! Component: 800 - total flux cleaned = 0.0194088 Jy ! Component: 850 - total flux cleaned = 0.0202093 Jy ! Component: 900 - total flux cleaned = 0.0208496 Jy ! Component: 950 - total flux cleaned = 0.0214746 Jy ! Component: 1000 - total flux cleaned = 0.0220846 Jy ! Component: 1050 - total flux cleaned = 0.022482 Jy ! Component: 1100 - total flux cleaned = 0.0228716 Jy ! Component: 1150 - total flux cleaned = 0.0231897 Jy ! Component: 1200 - total flux cleaned = 0.023564 Jy ! Component: 1250 - total flux cleaned = 0.0237477 Jy ! Component: 1300 - total flux cleaned = 0.0240495 Jy ! Component: 1350 - total flux cleaned = 0.0242871 Jy ! Component: 1400 - total flux cleaned = 0.024579 Jy ! Component: 1450 - total flux cleaned = 0.0249234 Jy ! Component: 1500 - total flux cleaned = 0.025092 Jy ! Component: 1550 - total flux cleaned = 0.0251474 Jy ! Component: 1600 - total flux cleaned = 0.0253121 Jy ! Component: 1650 - total flux cleaned = 0.0255823 Jy ! Component: 1700 - total flux cleaned = 0.025635 Jy ! Component: 1750 - total flux cleaned = 0.0257405 Jy ! Component: 1800 - total flux cleaned = 0.0259489 Jy ! Component: 1850 - total flux cleaned = 0.0261544 Jy ! Component: 1900 - total flux cleaned = 0.0263061 Jy ! Component: 1950 - total flux cleaned = 0.0263569 Jy ! Component: 2000 - total flux cleaned = 0.0265054 Jy ! Component: 2050 - total flux cleaned = 0.0266035 Jy ! Component: 2100 - total flux cleaned = 0.0266995 Jy ! Component: 2150 - total flux cleaned = 0.0266523 Jy ! Component: 2200 - total flux cleaned = 0.0268412 Jy ! Component: 2250 - total flux cleaned = 0.0269802 Jy ! Component: 2300 - total flux cleaned = 0.0269351 Jy ! Component: 2350 - total flux cleaned = 0.0270265 Jy ! Component: 2400 - total flux cleaned = 0.0271616 Jy ! Component: 2450 - total flux cleaned = 0.0272503 Jy ! Component: 2500 - total flux cleaned = 0.0274696 Jy ! Component: 2550 - total flux cleaned = 0.0275126 Jy ! Component: 2600 - total flux cleaned = 0.0275559 Jy ! Component: 2650 - total flux cleaned = 0.0276408 Jy ! Component: 2700 - total flux cleaned = 0.0275994 Jy ! Component: 2750 - total flux cleaned = 0.0277661 Jy ! Component: 2800 - total flux cleaned = 0.027766 Jy ! Component: 2850 - total flux cleaned = 0.0278069 Jy ! Component: 2900 - total flux cleaned = 0.027927 Jy ! Component: 2950 - total flux cleaned = 0.027887 Jy ! Component: 3000 - total flux cleaned = 0.0279663 Jy ! Component: 3050 - total flux cleaned = 0.028083 Jy ! Component: 3100 - total flux cleaned = 0.0280439 Jy ! Component: 3150 - total flux cleaned = 0.028006 Jy ! Component: 3200 - total flux cleaned = 0.0282342 Jy ! Component: 3250 - total flux cleaned = 0.0281964 Jy ! Component: 3300 - total flux cleaned = 0.0281965 Jy ! Component: 3350 - total flux cleaned = 0.0283078 Jy ! Component: 3400 - total flux cleaned = 0.0283811 Jy ! Component: 3450 - total flux cleaned = 0.0283814 Jy ! Component: 3500 - total flux cleaned = 0.0284896 Jy ! Component: 3550 - total flux cleaned = 0.0284896 Jy ! Component: 3600 - total flux cleaned = 0.0284897 Jy ! Component: 3650 - total flux cleaned = 0.0285957 Jy ! Component: 3700 - total flux cleaned = 0.0286303 Jy ! Component: 3750 - total flux cleaned = 0.0285955 Jy ! Component: 3800 - total flux cleaned = 0.0288373 Jy ! Component: 3850 - total flux cleaned = 0.0287685 Jy ! Component: 3900 - total flux cleaned = 0.0288707 Jy ! Component: 3950 - total flux cleaned = 0.0289383 Jy ! Component: 4000 - total flux cleaned = 0.0290055 Jy ! Component: 4050 - total flux cleaned = 0.0290723 Jy ! Total flux subtracted in 4096 components = 0.0291388 Jy ! Clean residual min=-0.001944 max=0.001741 Jy/beam ! Clean residual mean=0.000037 rms=0.000646 Jy/beam ! Combined flux in latest and established models = 1.01467 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: 0.9134 x 2.977 at -11.09 degrees (North through East) ! Clean map min=-0.0021424 max=0.72016 Jy/beam ! Clean map min=-0.0057455 max=0.72016 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=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 ! ! Move the cursor into the plot window and press 'H' for help mapcol grey ! Mapplot colormap: grey, contrast: 1 brightness: 0.5. save J0922-3959_X ! Writing UV FITS file: J0922-3959_X.uvf ! Writing 365 model components to file: J0922-3959_X.mod ! wwins: Wrote 9 windows to J0922-3959_X.win ! Adding 267 model components to the UV plane model. ! The established model now contains 365 components and 1.01467 Jy ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9134 x 2.977 at -11.09 degrees (North through East) ! Clean map min=-0.0056714 max=0.72016 Jy/beam ! Writing clean map to FITS file: J0922-3959_X.fits ! Writing difmap environment to: J0922-3959_X.par quit ! Quitting program ! Log file difmap.log_23 closed on Sat Jul 28 10:36:39 2018