! Started logfile: J1640+1220_S_2015_01_23_pus_uvs.log on Sun Jan 8 13:00:23 2017 obs J1640+1220_S_2015_01_23_pus_uva.fits ! Reading UV FITS file: J1640+1220_S_2015_01_23_pus_uva.fits ! AN table 1: 57 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.948928 visibilities/baseline/integration-bin. ! Found source: J1640+1220 ! ! 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.23588e+09 3.2e+07 1 3.2e+07 ! 02 2 2.26788e+09 3.2e+07 1 3.2e+07 ! 03 3 2.29988e+09 3.2e+07 1 3.2e+07 ! 04 4 2.36388e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 338 lines of history. ! ! Reading 9736 visibilities. select 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 print "field_size =",field_size ! field_size = 2048 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J1640+1220_S_2015_01_23_pus] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 2048x2048 pixels with 0.500x0.500 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 17 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 17 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 17 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 17 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.196551Jy sigma=0.005778 ! Fit after self-cal, rms=1.164389Jy sigma=0.005208 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. 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. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.725 mas, bmaj=6.495 mas, bpa=-2.838 degrees ! Estimated noise=989.165 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.437686 Jy ! Component: 100 - total flux cleaned = 0.587424 Jy ! Total flux subtracted in 100 components = 0.587424 Jy ! Clean residual min=-0.034600 max=0.097268 Jy/beam ! Clean residual mean=-0.000008 rms=0.011272 Jy/beam ! Combined flux in latest and established models = 0.587424 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.587424 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.167017Jy sigma=0.002885 ! Fit after self-cal, rms=1.166988Jy sigma=0.002884 ! Inverting map ! Added new window around map position (-4.5, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0989248 Jy ! Component: 100 - total flux cleaned = 0.163996 Jy ! Total flux subtracted in 100 components = 0.163996 Jy ! Clean residual min=-0.028244 max=0.046661 Jy/beam ! Clean residual mean=-0.000004 rms=0.008218 Jy/beam ! Combined flux in latest and established models = 0.75142 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 38 components and 0.75142 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.138301Jy sigma=0.002230 ! Fit after self-cal, rms=1.137816Jy sigma=0.002209 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.077 mas, bmaj=7.14 mas, bpa=-3.925 degrees ! Estimated noise=746.097 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0832301 Jy ! Component: 100 - total flux cleaned = 0.137111 Jy ! Component: 150 - total flux cleaned = 0.175912 Jy ! Component: 200 - total flux cleaned = 0.204825 Jy ! Total flux subtracted in 200 components = 0.204825 Jy ! Clean residual min=-0.033478 max=0.049831 Jy/beam ! Clean residual mean=-0.000002 rms=0.010586 Jy/beam ! Combined flux in latest and established models = 0.956246 Jy ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 89 components and 0.956245 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.114892Jy sigma=0.001804 ! Fit after self-cal, rms=1.114370Jy sigma=0.001779 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 0.49 HN 0.99 KP 1.03 ! LA 0.92 MK 1.02 NL 1.08 OV 1.04 ! PT 0.96 SC 1.00 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 0.59 HN 1.07 KP 0.95 ! LA 0.90 MK 1.03 NL 1.27 OV 1.08 ! PT 0.96 SC 1.07 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 0.79 HN 1.07 KP 0.84 ! LA 0.74 MK 1.06 NL 1.36 OV 1.01 ! PT 0.94 SC 1.06 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.13 FD 1.36 HN 1.04 KP 0.79 ! LA 0.79 MK 1.02 NL 1.02 OV 1.08 ! PT 0.98 SC 1.00 ! ! ! Fit before self-cal, rms=1.114370Jy sigma=0.001779 ! Fit after self-cal, rms=0.749433Jy sigma=0.001376 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.131 mas, bmaj=7.274 mas, bpa=-3.759 degrees ! Estimated noise=749.882 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00360698 Jy ! Component: 100 - total flux cleaned = 0.00466693 Jy ! Component: 150 - total flux cleaned = 0.00527549 Jy ! Component: 200 - total flux cleaned = 0.00543583 Jy ! Total flux subtracted in 200 components = 0.00543583 Jy ! Clean residual min=-0.023034 max=0.033067 Jy/beam ! Clean residual mean=0.000023 rms=0.007916 Jy/beam ! Combined flux in latest and established models = 0.961681 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 103 components and 0.961681 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.749205Jy sigma=0.001368 ! Fit after self-cal, rms=0.749199Jy sigma=0.001364 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.749199Jy sigma=0.001364 ! Fit after self-cal, rms=0.510348Jy sigma=0.001220 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.133 mas, bmaj=7.248 mas, bpa=-3.684 degrees ! Estimated noise=752.002 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00103952 Jy ! Component: 100 - total flux cleaned = 0.00153081 Jy ! Component: 150 - total flux cleaned = 0.00178819 Jy ! Component: 200 - total flux cleaned = 0.00201906 Jy ! Total flux subtracted in 200 components = 0.00201906 Jy ! Clean residual min=-0.021148 max=0.029709 Jy/beam ! Clean residual mean=0.000007 rms=0.006339 Jy/beam ! Combined flux in latest and established models = 0.9637 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 117 components and 0.9637 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.510255Jy sigma=0.001217 ! Fit after self-cal, rms=0.510269Jy sigma=0.001216 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-109, 94). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0307642 Jy ! Component: 100 - total flux cleaned = 0.0487701 Jy ! Component: 150 - total flux cleaned = 0.062224 Jy ! Component: 200 - total flux cleaned = 0.072676 Jy ! Total flux subtracted in 200 components = 0.072676 Jy ! Clean residual min=-0.019332 max=0.021351 Jy/beam ! Clean residual mean=0.000006 rms=0.004934 Jy/beam ! Combined flux in latest and established models = 1.03638 Jy ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 190 components and 1.03638 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.509081Jy sigma=0.001098 ! Fit after self-cal, rms=0.508882Jy sigma=0.001091 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.508882Jy sigma=0.001091 ! Fit after self-cal, rms=0.508822Jy sigma=0.001082 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.14 mas, bmaj=7.257 mas, bpa=-3.642 degrees ! Estimated noise=750.36 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0128014 Jy ! Component: 100 - total flux cleaned = 0.0227221 Jy ! Component: 150 - total flux cleaned = 0.0303908 Jy ! Component: 200 - total flux cleaned = 0.0353587 Jy ! Total flux subtracted in 200 components = 0.0353587 Jy ! Clean residual min=-0.017383 max=0.019931 Jy/beam ! Clean residual mean=0.000006 rms=0.004468 Jy/beam ! Combined flux in latest and established models = 1.07173 Jy selfcal ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 229 components and 1.07174 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.508712Jy sigma=0.001051 ! Fit after self-cal, rms=0.508612Jy sigma=0.001049 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (116, -42.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0188797 Jy ! Component: 100 - total flux cleaned = 0.0294471 Jy ! Component: 150 - total flux cleaned = 0.0373449 Jy ! Component: 200 - total flux cleaned = 0.0439617 Jy ! Total flux subtracted in 200 components = 0.0439617 Jy ! Clean residual min=-0.015539 max=0.016267 Jy/beam ! Clean residual mean=0.000003 rms=0.003807 Jy/beam ! Combined flux in latest and established models = 1.1157 Jy ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 295 components and 1.1157 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.508246Jy sigma=0.001001 ! Fit after self-cal, rms=0.508180Jy sigma=0.000999 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.508180Jy sigma=0.000999 ! Fit after self-cal, rms=0.511885Jy sigma=0.000993 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.191 mas, bmaj=10.98 mas, bpa=-4.728 degrees ! Estimated noise=978.436 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0162457 Jy ! Component: 100 - total flux cleaned = 0.027883 Jy ! Component: 150 - total flux cleaned = 0.0364814 Jy ! Component: 200 - total flux cleaned = 0.041055 Jy ! Total flux subtracted in 200 components = 0.041055 Jy ! Clean residual min=-0.021564 max=0.027597 Jy/beam ! Clean residual mean=0.000003 rms=0.006627 Jy/beam ! Combined flux in latest and established models = 1.15675 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 334 components and 1.15675 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.512043Jy sigma=0.000978 ! Fit after self-cal, rms=0.511979Jy sigma=0.000976 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00779409 Jy ! Component: 100 - total flux cleaned = 0.0115384 Jy ! Component: 150 - total flux cleaned = 0.014841 Jy ! Component: 200 - total flux cleaned = 0.0175736 Jy ! Total flux subtracted in 200 components = 0.0175736 Jy ! Clean residual min=-0.021553 max=0.026371 Jy/beam ! Clean residual mean=0.000000 rms=0.006429 Jy/beam ! Combined flux in latest and established models = 1.17433 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 340 components and 1.17433 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.512113Jy sigma=0.000974 ! Fit after self-cal, rms=0.512083Jy sigma=0.000973 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=3.144 mas, bmaj=7.269 mas, bpa=-3.612 degrees ! Estimated noise=748.39 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.144 x 7.269 at -3.612 degrees (North through East) ! Clean map min=-0.0137 max=0.58587 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.512083Jy sigma=0.000973 ! Fit after self-cal, rms=0.648271Jy sigma=0.000922 ! Inverting map and beam ! Estimated beam: bmin=3.145 mas, bmaj=7.279 mas, bpa=-3.652 degrees ! Estimated noise=746.705 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000603012 Jy ! Component: 100 - total flux cleaned = -0.000304979 Jy ! Component: 150 - total flux cleaned = 0.000357297 Jy ! Component: 200 - total flux cleaned = 0.000961639 Jy ! Total flux subtracted in 200 components = 0.000961639 Jy ! Clean residual min=-0.013638 max=0.012921 Jy/beam ! Clean residual mean=0.000001 rms=0.003283 Jy/beam ! Combined flux in latest and established models = 1.17529 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 359 components and 1.17529 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.648206Jy sigma=0.000919 ! Fit after self-cal, rms=0.648210Jy sigma=0.000919 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000466659 Jy ! Component: 100 - total flux cleaned = 0.00101559 Jy ! Component: 150 - total flux cleaned = 0.00132349 Jy ! Component: 200 - total flux cleaned = 0.00161177 Jy ! Total flux subtracted in 200 components = 0.00161177 Jy ! Clean residual min=-0.013655 max=0.012888 Jy/beam ! Clean residual mean=0.000000 rms=0.003273 Jy/beam ! Combined flux in latest and established models = 1.1769 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 367 components and 1.1769 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.648161Jy sigma=0.000917 ! Fit after self-cal, rms=0.648159Jy sigma=0.000917 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.648159Jy sigma=0.000917 ! Fit after self-cal, rms=0.654210Jy sigma=0.000916 print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. 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 clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.768 mas, bmaj=6.584 mas, bpa=-2.372 degrees ! Estimated noise=1004.15 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.442165 Jy ! Component: 100 - total flux cleaned = 0.593847 Jy ! Component: 150 - total flux cleaned = 0.687828 Jy ! Component: 200 - total flux cleaned = 0.749772 Jy ! Total flux subtracted in 200 components = 0.749772 Jy ! Clean residual min=-0.021722 max=0.033901 Jy/beam ! Clean residual mean=-0.000003 rms=0.005996 Jy/beam ! Combined flux in latest and established models = 0.749772 Jy uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 25 model components to the UV plane model. ! The established model now contains 25 components and 0.749772 Jy ! Inverting map and beam ! Estimated beam: bmin=3.149 mas, bmaj=7.286 mas, bpa=-3.594 degrees ! Estimated noise=747.261 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.06242 Jy ! Component: 100 - total flux cleaned = 0.104824 Jy ! Component: 150 - total flux cleaned = 0.137958 Jy ! Component: 200 - total flux cleaned = 0.16563 Jy ! Total flux subtracted in 200 components = 0.16563 Jy ! Clean residual min=-0.019822 max=0.023802 Jy/beam ! Clean residual mean=-0.000001 rms=0.006592 Jy/beam ! Combined flux in latest and established models = 0.915401 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 97 components and 0.915401 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.2 mas, bmaj=11.01 mas, bpa=-4.701 degrees ! Estimated noise=974.352 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0496843 Jy ! Component: 100 - total flux cleaned = 0.087349 Jy ! Component: 150 - total flux cleaned = 0.117112 Jy ! Component: 200 - total flux cleaned = 0.141214 Jy ! Total flux subtracted in 200 components = 0.141214 Jy ! Clean residual min=-0.029388 max=0.030554 Jy/beam ! Clean residual mean=-0.000005 rms=0.008900 Jy/beam ! Combined flux in latest and established models = 1.05662 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0200558 Jy ! Component: 100 - total flux cleaned = 0.0369142 Jy ! Component: 150 - total flux cleaned = 0.0511493 Jy ! Component: 200 - total flux cleaned = 0.0632588 Jy ! Total flux subtracted in 200 components = 0.0632588 Jy ! Clean residual min=-0.023896 max=0.026771 Jy/beam ! Clean residual mean=-0.000005 rms=0.007185 Jy/beam ! Combined flux in latest and established models = 1.11987 Jy ! Adding 133 model components to the UV plane model. ! The established model now contains 219 components and 1.11987 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.149 mas, bmaj=7.286 mas, bpa=-3.594 degrees ! Estimated noise=747.261 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00327996 Jy ! Component: 100 - total flux cleaned = 0.00726205 Jy ! Component: 150 - total flux cleaned = 0.0113526 Jy ! Component: 200 - total flux cleaned = 0.0152767 Jy ! Total flux subtracted in 200 components = 0.0152767 Jy ! Clean residual min=-0.015099 max=0.014130 Jy/beam ! Clean residual mean=-0.000002 rms=0.003513 Jy/beam ! Combined flux in latest and established models = 1.13515 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 237 components and 1.13515 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.654118Jy sigma=0.000941 ! Fit after self-cal, rms=0.653954Jy sigma=0.000936 wmodel J1640+1220_S_2015_01_23_pus_map.mod ! Writing 237 model components to file: J1640+1220_S_2015_01_23_pus_map.mod wobs J1640+1220_S_2015_01_23_pus_uvs.fits ! Writing UV FITS file: J1640+1220_S_2015_01_23_pus_uvs.fits wwins J1640+1220_S_2015_01_23_pus_map.win ! wwins: Wrote 4 windows to J1640+1220_S_2015_01_23_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=8192 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0004691 Jy ! Component: 100 - total flux cleaned = 0.00168775 Jy ! Component: 150 - total flux cleaned = 0.00259813 Jy ! Component: 200 - total flux cleaned = 0.00367989 Jy ! Component: 250 - total flux cleaned = 0.00409408 Jy ! Component: 300 - total flux cleaned = 0.00469015 Jy ! Component: 350 - total flux cleaned = 0.00544719 Jy ! Component: 400 - total flux cleaned = 0.00617835 Jy ! Component: 450 - total flux cleaned = 0.00688699 Jy ! Component: 500 - total flux cleaned = 0.00740205 Jy ! Component: 550 - total flux cleaned = 0.0075758 Jy ! Component: 600 - total flux cleaned = 0.00839332 Jy ! Component: 650 - total flux cleaned = 0.00933745 Jy ! Component: 700 - total flux cleaned = 0.0102529 Jy ! Component: 750 - total flux cleaned = 0.0112957 Jy ! Component: 800 - total flux cleaned = 0.012312 Jy ! Component: 850 - total flux cleaned = 0.0131591 Jy ! Component: 900 - total flux cleaned = 0.0142629 Jy ! Component: 950 - total flux cleaned = 0.0150739 Jy ! Component: 1000 - total flux cleaned = 0.0167859 Jy ! Component: 1050 - total flux cleaned = 0.0178158 Jy ! Component: 1100 - total flux cleaned = 0.0185707 Jy ! Component: 1150 - total flux cleaned = 0.0195546 Jy ! Component: 1200 - total flux cleaned = 0.0203986 Jy ! Component: 1250 - total flux cleaned = 0.0211084 Jy ! Component: 1300 - total flux cleaned = 0.021803 Jy ! Component: 1350 - total flux cleaned = 0.0224847 Jy ! Component: 1400 - total flux cleaned = 0.0230393 Jy ! Component: 1450 - total flux cleaned = 0.0238014 Jy ! Component: 1500 - total flux cleaned = 0.0245458 Jy ! Component: 1550 - total flux cleaned = 0.024859 Jy ! Component: 1600 - total flux cleaned = 0.0255772 Jy ! Component: 1650 - total flux cleaned = 0.0258844 Jy ! Component: 1700 - total flux cleaned = 0.0268752 Jy ! Component: 1750 - total flux cleaned = 0.0272631 Jy ! Component: 1800 - total flux cleaned = 0.0274554 Jy ! Component: 1850 - total flux cleaned = 0.0282058 Jy ! Component: 1900 - total flux cleaned = 0.0286669 Jy ! Component: 1950 - total flux cleaned = 0.0289386 Jy ! Component: 2000 - total flux cleaned = 0.0292078 Jy ! Component: 2050 - total flux cleaned = 0.0299969 Jy ! Component: 2100 - total flux cleaned = 0.030254 Jy ! Component: 2150 - total flux cleaned = 0.0306784 Jy ! Component: 2200 - total flux cleaned = 0.031096 Jy ! Component: 2250 - total flux cleaned = 0.0312601 Jy ! Component: 2300 - total flux cleaned = 0.0315022 Jy ! Component: 2350 - total flux cleaned = 0.0314219 Jy ! Component: 2400 - total flux cleaned = 0.0315005 Jy ! Component: 2450 - total flux cleaned = 0.0316559 Jy ! Component: 2500 - total flux cleaned = 0.0316572 Jy ! Component: 2550 - total flux cleaned = 0.0316578 Jy ! Component: 2600 - total flux cleaned = 0.0316572 Jy ! Component: 2650 - total flux cleaned = 0.0315846 Jy ! Component: 2700 - total flux cleaned = 0.0313712 Jy ! Component: 2750 - total flux cleaned = 0.0315823 Jy ! Component: 2800 - total flux cleaned = 0.0312361 Jy ! Component: 2850 - total flux cleaned = 0.0313049 Jy ! Component: 2900 - total flux cleaned = 0.0313059 Jy ! Component: 2950 - total flux cleaned = 0.0312405 Jy ! Component: 3000 - total flux cleaned = 0.0312418 Jy ! Component: 3050 - total flux cleaned = 0.031113 Jy ! Component: 3100 - total flux cleaned = 0.0310497 Jy ! Component: 3150 - total flux cleaned = 0.0313017 Jy ! Component: 3200 - total flux cleaned = 0.0310528 Jy ! Component: 3250 - total flux cleaned = 0.0309927 Jy ! Component: 3300 - total flux cleaned = 0.031053 Jy ! Component: 3350 - total flux cleaned = 0.0309345 Jy ! Component: 3400 - total flux cleaned = 0.0306393 Jy ! Component: 3450 - total flux cleaned = 0.0306981 Jy ! Component: 3500 - total flux cleaned = 0.0305838 Jy ! Component: 3550 - total flux cleaned = 0.0304708 Jy ! Component: 3600 - total flux cleaned = 0.0303589 Jy ! Component: 3650 - total flux cleaned = 0.0304152 Jy ! Component: 3700 - total flux cleaned = 0.03014 Jy ! Component: 3750 - total flux cleaned = 0.0300324 Jy ! Component: 3800 - total flux cleaned = 0.0298716 Jy ! Component: 3850 - total flux cleaned = 0.0298184 Jy ! Component: 3900 - total flux cleaned = 0.0295033 Jy ! Component: 3950 - total flux cleaned = 0.0294522 Jy ! Component: 4000 - total flux cleaned = 0.0295547 Jy ! Component: 4050 - total flux cleaned = 0.029403 Jy ! Component: 4100 - total flux cleaned = 0.0294033 Jy ! Component: 4150 - total flux cleaned = 0.029205 Jy ! Component: 4200 - total flux cleaned = 0.0292542 Jy ! Component: 4250 - total flux cleaned = 0.0290598 Jy ! Component: 4300 - total flux cleaned = 0.0289644 Jy ! Component: 4350 - total flux cleaned = 0.0287738 Jy ! Component: 4400 - total flux cleaned = 0.0285381 Jy ! Component: 4450 - total flux cleaned = 0.028632 Jy ! Component: 4500 - total flux cleaned = 0.0283529 Jy ! Component: 4550 - total flux cleaned = 0.0281687 Jy ! Component: 4600 - total flux cleaned = 0.0279863 Jy ! Component: 4650 - total flux cleaned = 0.0279411 Jy ! Component: 4700 - total flux cleaned = 0.027852 Jy ! Component: 4750 - total flux cleaned = 0.0278082 Jy ! Component: 4800 - total flux cleaned = 0.0276321 Jy ! Component: 4850 - total flux cleaned = 0.0274575 Jy ! Component: 4900 - total flux cleaned = 0.0272842 Jy ! Component: 4950 - total flux cleaned = 0.0272412 Jy ! Component: 5000 - total flux cleaned = 0.0270277 Jy ! Component: 5050 - total flux cleaned = 0.0269437 Jy ! Component: 5100 - total flux cleaned = 0.0268187 Jy ! Component: 5150 - total flux cleaned = 0.0266522 Jy ! Component: 5200 - total flux cleaned = 0.0266111 Jy ! Component: 5250 - total flux cleaned = 0.0264883 Jy ! Component: 5300 - total flux cleaned = 0.0262045 Jy ! Component: 5350 - total flux cleaned = 0.0260839 Jy ! Component: 5400 - total flux cleaned = 0.0260438 Jy ! Component: 5450 - total flux cleaned = 0.0259644 Jy ! Component: 5500 - total flux cleaned = 0.0257677 Jy ! Component: 5550 - total flux cleaned = 0.0258073 Jy ! Component: 5600 - total flux cleaned = 0.0257684 Jy ! Component: 5650 - total flux cleaned = 0.0254994 Jy ! Component: 5700 - total flux cleaned = 0.0254989 Jy ! Component: 5750 - total flux cleaned = 0.025272 Jy ! Component: 5800 - total flux cleaned = 0.0253094 Jy ! Component: 5850 - total flux cleaned = 0.0250859 Jy ! Component: 5900 - total flux cleaned = 0.0250488 Jy ! Component: 5950 - total flux cleaned = 0.0251958 Jy ! Component: 6000 - total flux cleaned = 0.0249766 Jy ! Component: 6050 - total flux cleaned = 0.0248679 Jy ! Component: 6100 - total flux cleaned = 0.024724 Jy ! Component: 6150 - total flux cleaned = 0.0248316 Jy ! Component: 6200 - total flux cleaned = 0.0247247 Jy ! Component: 6250 - total flux cleaned = 0.0246893 Jy ! Component: 6300 - total flux cleaned = 0.0247243 Jy ! Component: 6350 - total flux cleaned = 0.0245494 Jy ! Component: 6400 - total flux cleaned = 0.0245151 Jy ! Component: 6450 - total flux cleaned = 0.0245498 Jy ! Component: 6500 - total flux cleaned = 0.0245838 Jy ! Component: 6550 - total flux cleaned = 0.0245495 Jy ! Component: 6600 - total flux cleaned = 0.0243465 Jy ! Component: 6650 - total flux cleaned = 0.0243469 Jy ! Component: 6700 - total flux cleaned = 0.0243134 Jy ! Component: 6750 - total flux cleaned = 0.0243464 Jy ! Component: 6800 - total flux cleaned = 0.0243462 Jy ! Component: 6850 - total flux cleaned = 0.0242802 Jy ! Component: 6900 - total flux cleaned = 0.0241169 Jy ! Component: 6950 - total flux cleaned = 0.0241494 Jy ! Component: 7000 - total flux cleaned = 0.0240523 Jy ! Component: 7050 - total flux cleaned = 0.0240201 Jy ! Component: 7100 - total flux cleaned = 0.0239238 Jy ! Component: 7150 - total flux cleaned = 0.0238602 Jy ! Component: 7200 - total flux cleaned = 0.0238284 Jy ! Component: 7250 - total flux cleaned = 0.0235763 Jy ! Component: 7300 - total flux cleaned = 0.0236706 Jy ! Component: 7350 - total flux cleaned = 0.0234835 Jy ! Component: 7400 - total flux cleaned = 0.0235142 Jy ! Component: 7450 - total flux cleaned = 0.0234527 Jy ! Component: 7500 - total flux cleaned = 0.0234524 Jy ! Component: 7550 - total flux cleaned = 0.0234216 Jy ! Component: 7600 - total flux cleaned = 0.0233302 Jy ! Component: 7650 - total flux cleaned = 0.0232696 Jy ! Component: 7700 - total flux cleaned = 0.0232996 Jy ! Component: 7750 - total flux cleaned = 0.0231795 Jy ! Component: 7800 - total flux cleaned = 0.0231197 Jy ! Component: 7850 - total flux cleaned = 0.0231793 Jy ! Component: 7900 - total flux cleaned = 0.0231201 Jy ! Component: 7950 - total flux cleaned = 0.0230908 Jy ! Component: 8000 - total flux cleaned = 0.0230613 Jy ! Component: 8050 - total flux cleaned = 0.0230318 Jy ! Component: 8100 - total flux cleaned = 0.0230609 Jy ! Component: 8150 - total flux cleaned = 0.0229158 Jy ! Total flux subtracted in 8192 components = 0.0229446 Jy ! Clean residual min=-0.001480 max=0.001740 Jy/beam ! Clean residual mean=-0.000001 rms=0.000500 Jy/beam ! Combined flux in latest and established models = 1.1581 Jy keep ! Adding 815 model components to the UV plane model. ! The established model now contains 1052 components and 1.1581 Jy ![Exited script file: muppet] wmap "J1640+1220_S_2015_01_23_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.149 x 7.286 at -3.594 degrees (North through East) ! Clean map min=-0.0070159 max=0.59453 Jy/beam ! Writing clean map to FITS file: J1640+1220_S_2015_01_23_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.594533 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000480696 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.242558 0.242558 0.485116 0.970232 1.94046 3.88093 7.76186 15.5237 31.0474 62.0949 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0014 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.594 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.149 x 7.286 at -3.594 degrees (North through East) ! Clean map min=-0.0070159 max=0.59453 Jy/beam ! Error occured in command: mapplot ! ! ! ! ! ! 4 consecutive EOF's on stdin - 10 will exit program ! ! 5 consecutive EOF's on stdin - 10 will exit program ! ! 6 consecutive EOF's on stdin - 10 will exit program ! ! 7 consecutive EOF's on stdin - 10 will exit program ! ! 8 consecutive EOF's on stdin - 10 will exit program ! ! 9 consecutive EOF's on stdin - 10 will exit program ! ! 10 consecutive EOF's on stdin - 10 will exit program ! Log file J1640+1220_S_2015_01_23_pus_uvs.log closed on Sun Jan 8 13:00:52 2017