! Started logfile: J0922-3959_S_2014_12_20_pus_uvs.log on Wed Dec 28 13:30:21 2016 obs J0922-3959_S_2014_12_20_pus_uva.fits ! Reading UV FITS file: J0922-3959_S_2014_12_20_pus_uva.fits ! AN table 1: 15 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.582222 visibilities/baseline/integration-bin. ! 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 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 344 lines of history. ! ! Reading 1572 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 = 1024 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J0922-3959_S_2014_12_20_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 = 1024x1024 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 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=0.409096Jy sigma=0.002645 ! Fit after self-cal, rms=0.376391Jy sigma=0.002369 ! 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.988 mas, bmaj=18.73 mas, bpa=7.212 degrees ! Estimated noise=3167.86 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.595892 Jy ! Component: 100 - total flux cleaned = 0.802159 Jy ! Total flux subtracted in 100 components = 0.802159 Jy ! Clean residual min=-0.048818 max=0.111023 Jy/beam ! Clean residual mean=0.000004 rms=0.021442 Jy/beam ! Combined flux in latest and established models = 0.802159 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 23 components and 0.802159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.301923Jy sigma=0.002072 ! Fit after self-cal, rms=0.301915Jy sigma=0.002072 ! 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.753 mas, bmaj=20.27 mas, bpa=6.018 degrees ! Estimated noise=2598.89 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.176675 Jy ! Component: 100 - total flux cleaned = 0.26658 Jy ! Component: 150 - total flux cleaned = 0.317303 Jy ! Component: 200 - total flux cleaned = 0.347278 Jy ! Total flux subtracted in 200 components = 0.347278 Jy ! Clean residual min=-0.042487 max=0.034638 Jy/beam ! Clean residual mean=-0.000103 rms=0.012071 Jy/beam ! Combined flux in latest and established models = 1.14944 Jy ! Performing phase self-cal ! Adding 97 model components to the UV plane model. ! The established model now contains 112 components and 1.14944 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187177Jy sigma=0.001154 ! Fit after self-cal, rms=0.187009Jy sigma=0.001153 ! 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.00* FD 0.84 HN 1.07 KP 1.07 ! LA 0.95 MK 1.00* NL 1.09 OV 0.99 ! PT 1.05 SC 0.95 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.02 HN 1.18 KP 1.05 ! LA 0.94 MK 1.00* NL 1.09 OV 1.01 ! PT 1.01 SC 1.04 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.13 HN 1.10 KP 0.94 ! LA 0.75 MK 1.00* NL 1.11 OV 0.93 ! PT 1.02 SC 1.07 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00 HN 1.05 KP 0.90 ! LA 0.83 MK 1.00* NL 1.08 OV 0.98 ! PT 1.03 SC 1.07 ! ! ! Fit before self-cal, rms=0.187009Jy sigma=0.001153 ! Fit after self-cal, rms=0.142359Jy sigma=0.000851 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.796 mas, bmaj=20.83 mas, bpa=6.132 degrees ! Estimated noise=2611.46 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.0123319 Jy ! Component: 100 - total flux cleaned = -0.0138523 Jy ! Component: 150 - total flux cleaned = -0.0147398 Jy ! Component: 200 - total flux cleaned = -0.0157147 Jy ! Total flux subtracted in 200 components = -0.0157147 Jy ! Clean residual min=-0.022050 max=0.021384 Jy/beam ! Clean residual mean=-0.000106 rms=0.005569 Jy/beam ! Combined flux in latest and established models = 1.13372 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 132 components and 1.13372 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139513Jy sigma=0.000841 ! Fit after self-cal, rms=0.139558Jy sigma=0.000841 ! 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.139558Jy sigma=0.000841 ! Fit after self-cal, rms=0.140103Jy sigma=0.000837 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.803 mas, bmaj=20.96 mas, bpa=6.168 degrees ! Estimated noise=2620.09 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00755399 Jy ! Component: 100 - total flux cleaned = -0.00921152 Jy ! Component: 150 - total flux cleaned = -0.0101274 Jy ! Component: 200 - total flux cleaned = -0.010925 Jy ! Total flux subtracted in 200 components = -0.010925 Jy ! Clean residual min=-0.021724 max=0.020086 Jy/beam ! Clean residual mean=-0.000116 rms=0.005286 Jy/beam ! Combined flux in latest and established models = 1.1228 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 141 components and 1.1228 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138580Jy sigma=0.000833 ! Fit after self-cal, rms=0.138620Jy sigma=0.000833 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.00051364 Jy ! Component: 100 - total flux cleaned = -0.000998345 Jy ! Component: 150 - total flux cleaned = -0.00129282 Jy ! Component: 200 - total flux cleaned = -0.00142979 Jy ! Total flux subtracted in 200 components = -0.00142979 Jy ! Clean residual min=-0.021510 max=0.019851 Jy/beam ! Clean residual mean=-0.000117 rms=0.005249 Jy/beam ! Combined flux in latest and established models = 1.12137 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 152 components and 1.12137 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138234Jy sigma=0.000831 ! Fit after self-cal, rms=0.138237Jy sigma=0.000831 ! 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.138237Jy sigma=0.000831 ! Fit after self-cal, rms=0.139367Jy sigma=0.000828 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.815 mas, bmaj=21.1 mas, bpa=6.2 degrees ! Estimated noise=2632.08 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0051357 Jy ! Component: 100 - total flux cleaned = -0.00639864 Jy ! Component: 150 - total flux cleaned = -0.00710999 Jy ! Component: 200 - total flux cleaned = -0.00774802 Jy ! Total flux subtracted in 200 components = -0.00774802 Jy ! Clean residual min=-0.021253 max=0.018810 Jy/beam ! Clean residual mean=-0.000126 rms=0.005061 Jy/beam ! Combined flux in latest and established models = 1.11362 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 158 components and 1.11362 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138547Jy sigma=0.000825 ! Fit after self-cal, rms=0.138583Jy sigma=0.000825 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.000437938 Jy ! Component: 100 - total flux cleaned = -0.000575945 Jy ! Component: 150 - total flux cleaned = -0.000836558 Jy ! Component: 200 - total flux cleaned = -0.000954791 Jy ! Total flux subtracted in 200 components = -0.000954791 Jy ! Clean residual min=-0.021055 max=0.018701 Jy/beam ! Clean residual mean=-0.000126 rms=0.005030 Jy/beam ! Combined flux in latest and established models = 1.11267 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 162 components and 1.11266 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138334Jy sigma=0.000824 ! Fit after self-cal, rms=0.138335Jy sigma=0.000824 ! 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.138335Jy sigma=0.000824 ! Fit after self-cal, rms=0.139480Jy sigma=0.000822 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.194 mas, bmaj=31.05 mas, bpa=-3.325 degrees ! Estimated noise=2971.59 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00340186 Jy ! Component: 100 - total flux cleaned = -0.00362191 Jy ! Component: 150 - total flux cleaned = -0.00372399 Jy ! Component: 200 - total flux cleaned = -0.00391296 Jy ! Total flux subtracted in 200 components = -0.00391296 Jy ! Clean residual min=-0.015147 max=0.018106 Jy/beam ! Clean residual mean=-0.000219 rms=0.005951 Jy/beam ! Combined flux in latest and established models = 1.10875 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 170 components and 1.10875 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139093Jy sigma=0.000822 ! Fit after self-cal, rms=0.139114Jy sigma=0.000822 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.000104932 Jy ! Component: 100 - total flux cleaned = -0.000109679 Jy ! Component: 150 - total flux cleaned = -0.000107717 Jy ! Component: 200 - total flux cleaned = -3.17372e-05 Jy ! Total flux subtracted in 200 components = -3.17372e-05 Jy ! Clean residual min=-0.015056 max=0.018030 Jy/beam ! Clean residual mean=-0.000219 rms=0.005940 Jy/beam ! Combined flux in latest and established models = 1.10872 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 176 components and 1.10872 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139031Jy sigma=0.000822 ! Fit after self-cal, rms=0.139029Jy sigma=0.000822 ! 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.826 mas, bmaj=21.2 mas, bpa=6.221 degrees ! Estimated noise=2641.2 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.826 x 21.2 at 6.221 degrees (North through East) ! Clean map min=-0.019769 max=0.64539 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.139029Jy sigma=0.000822 ! Fit after self-cal, rms=0.132650Jy sigma=0.000754 ! Inverting map and beam ! Estimated beam: bmin=3.839 mas, bmaj=21.25 mas, bpa=6.167 degrees ! Estimated noise=2650.8 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00343122 Jy ! Component: 100 - total flux cleaned = -0.00407158 Jy ! Component: 150 - total flux cleaned = -0.00449925 Jy ! Component: 200 - total flux cleaned = -0.00476977 Jy ! Total flux subtracted in 200 components = -0.00476977 Jy ! Clean residual min=-0.020788 max=0.017741 Jy/beam ! Clean residual mean=-0.000136 rms=0.004897 Jy/beam ! Combined flux in latest and established models = 1.10395 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 187 components and 1.10395 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132239Jy sigma=0.000752 ! Fit after self-cal, rms=0.132242Jy sigma=0.000752 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000253981 Jy ! Component: 100 - total flux cleaned = -0.00049012 Jy ! Component: 150 - total flux cleaned = -0.000604735 Jy ! Component: 200 - total flux cleaned = -0.000719005 Jy ! Total flux subtracted in 200 components = -0.000719005 Jy ! Clean residual min=-0.020645 max=0.017646 Jy/beam ! Clean residual mean=-0.000136 rms=0.004869 Jy/beam ! Combined flux in latest and established models = 1.10323 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 188 components and 1.10323 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132052Jy sigma=0.000751 ! Fit after self-cal, rms=0.132054Jy sigma=0.000751 ! 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.132054Jy sigma=0.000751 ! Fit after self-cal, rms=0.133107Jy sigma=0.000750 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=3.011 mas, bmaj=19.47 mas, bpa=7.426 degrees ! Estimated noise=3320.58 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.622309 Jy ! Component: 100 - total flux cleaned = 0.822024 Jy ! Component: 150 - total flux cleaned = 0.928889 Jy ! Component: 200 - total flux cleaned = 0.990743 Jy ! Total flux subtracted in 200 components = 0.990743 Jy ! Clean residual min=-0.025680 max=0.046263 Jy/beam ! Clean residual mean=-0.000051 rms=0.008285 Jy/beam ! Combined flux in latest and established models = 0.990743 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 69 model components to the UV plane model. ! The established model now contains 69 components and 0.990743 Jy ! Inverting map and beam ! Estimated beam: bmin=3.849 mas, bmaj=21.32 mas, bpa=6.181 degrees ! Estimated noise=2657.87 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) ! Added new window around map position (4.5, -14.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.056241 Jy ! Component: 100 - total flux cleaned = 0.0825428 Jy ! Component: 150 - total flux cleaned = 0.0971848 Jy ! Component: 200 - total flux cleaned = 0.106578 Jy ! Total flux subtracted in 200 components = 0.106578 Jy ! Clean residual min=-0.019182 max=0.017262 Jy/beam ! Clean residual mean=-0.000129 rms=0.004668 Jy/beam ! Combined flux in latest and established models = 1.09732 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 133 components and 1.09732 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=8.248 mas, bmaj=31.47 mas, bpa=-3.468 degrees ! Estimated noise=2967.96 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00843356 Jy ! Component: 100 - total flux cleaned = 0.0105861 Jy ! Component: 150 - total flux cleaned = 0.0121141 Jy ! Component: 200 - total flux cleaned = 0.0133389 Jy ! Total flux subtracted in 200 components = 0.0133389 Jy ! Clean residual min=-0.015686 max=0.019433 Jy/beam ! Clean residual mean=-0.000190 rms=0.005591 Jy/beam ! Combined flux in latest and established models = 1.11066 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.000677335 Jy ! Component: 100 - total flux cleaned = 0.00110545 Jy ! Component: 150 - total flux cleaned = 0.00150117 Jy ! Component: 200 - total flux cleaned = 0.00178712 Jy ! Total flux subtracted in 200 components = 0.00178712 Jy ! Clean residual min=-0.015506 max=0.019623 Jy/beam ! Clean residual mean=-0.000183 rms=0.005530 Jy/beam ! Combined flux in latest and established models = 1.11245 Jy ! Adding 28 model components to the UV plane model. ! The established model now contains 155 components and 1.11245 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.849 mas, bmaj=21.32 mas, bpa=6.181 degrees ! Estimated noise=2657.87 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00134917 Jy ! Component: 100 - total flux cleaned = -0.00156214 Jy ! Component: 150 - total flux cleaned = -0.00155654 Jy ! Component: 200 - total flux cleaned = -0.0014049 Jy ! Total flux subtracted in 200 components = -0.0014049 Jy ! Clean residual min=-0.016403 max=0.016529 Jy/beam ! Clean residual mean=-0.000113 rms=0.004290 Jy/beam ! Combined flux in latest and established models = 1.11104 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 178 components and 1.11104 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130406Jy sigma=0.000731 ! Fit after self-cal, rms=0.130256Jy sigma=0.000729 wmodel J0922-3959_S_2014_12_20_pus_map.mod ! Writing 178 model components to file: J0922-3959_S_2014_12_20_pus_map.mod wobs J0922-3959_S_2014_12_20_pus_uvs.fits ! Writing UV FITS file: J0922-3959_S_2014_12_20_pus_uvs.fits wwins J0922-3959_S_2014_12_20_pus_map.win ! wwins: Wrote 2 windows to J0922-3959_S_2014_12_20_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=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00151043 Jy ! Component: 100 - total flux cleaned = 0.00178296 Jy ! Component: 150 - total flux cleaned = 0.00200904 Jy ! Component: 200 - total flux cleaned = 0.00221905 Jy ! Component: 250 - total flux cleaned = 0.00222295 Jy ! Component: 300 - total flux cleaned = 0.00273476 Jy ! Component: 350 - total flux cleaned = 0.00256971 Jy ! Component: 400 - total flux cleaned = 0.00226521 Jy ! Component: 450 - total flux cleaned = 0.00197403 Jy ! Component: 500 - total flux cleaned = 0.00155509 Jy ! Component: 550 - total flux cleaned = 0.00155693 Jy ! Component: 600 - total flux cleaned = 0.00142897 Jy ! Component: 650 - total flux cleaned = 0.00130392 Jy ! Component: 700 - total flux cleaned = 0.000944056 Jy ! Component: 750 - total flux cleaned = 0.000826492 Jy ! Component: 800 - total flux cleaned = 0.000825973 Jy ! Component: 850 - total flux cleaned = 0.000615732 Jy ! Component: 900 - total flux cleaned = 0.000826525 Jy ! Component: 950 - total flux cleaned = 0.000924092 Jy ! Component: 1000 - total flux cleaned = 0.000731123 Jy ! Component: 1050 - total flux cleaned = 0.000825661 Jy ! Component: 1100 - total flux cleaned = 0.000641664 Jy ! Component: 1150 - total flux cleaned = 0.000551731 Jy ! Component: 1200 - total flux cleaned = 0.000550197 Jy ! Component: 1250 - total flux cleaned = 0.000383981 Jy ! Component: 1300 - total flux cleaned = 0.000304971 Jy ! Component: 1350 - total flux cleaned = -1.13636e-05 Jy ! Component: 1400 - total flux cleaned = -0.000552682 Jy ! Component: 1450 - total flux cleaned = -0.000854562 Jy ! Component: 1500 - total flux cleaned = -0.00129913 Jy ! Component: 1550 - total flux cleaned = -0.00166296 Jy ! Component: 1600 - total flux cleaned = -0.00201965 Jy ! Component: 1650 - total flux cleaned = -0.0024406 Jy ! Component: 1700 - total flux cleaned = -0.00299054 Jy ! Component: 1750 - total flux cleaned = -0.00346253 Jy ! Component: 1800 - total flux cleaned = -0.00372705 Jy ! Component: 1850 - total flux cleaned = -0.00418527 Jy ! Component: 1900 - total flux cleaned = -0.00476602 Jy ! Component: 1950 - total flux cleaned = -0.00533917 Jy ! Component: 2000 - total flux cleaned = -0.00590709 Jy ! Component: 2050 - total flux cleaned = -0.00647011 Jy ! Component: 2100 - total flux cleaned = -0.00690494 Jy ! Component: 2150 - total flux cleaned = -0.00752089 Jy ! Component: 2200 - total flux cleaned = -0.00794903 Jy ! Component: 2250 - total flux cleaned = -0.00843577 Jy ! Component: 2300 - total flux cleaned = -0.00867748 Jy ! Component: 2350 - total flux cleaned = -0.00903898 Jy ! Component: 2400 - total flux cleaned = -0.00951783 Jy ! Component: 2450 - total flux cleaned = -0.00969686 Jy ! Component: 2500 - total flux cleaned = -0.0100539 Jy ! Component: 2550 - total flux cleaned = -0.0104106 Jy ! Component: 2600 - total flux cleaned = -0.0107071 Jy ! Component: 2650 - total flux cleaned = -0.0108257 Jy ! Component: 2700 - total flux cleaned = -0.0112411 Jy ! Component: 2750 - total flux cleaned = -0.0114793 Jy ! Component: 2800 - total flux cleaned = -0.0117173 Jy ! Component: 2850 - total flux cleaned = -0.0118965 Jy ! Component: 2900 - total flux cleaned = -0.0121961 Jy ! Component: 2950 - total flux cleaned = -0.0124971 Jy ! Component: 3000 - total flux cleaned = -0.0126178 Jy ! Component: 3050 - total flux cleaned = -0.012921 Jy ! Component: 3100 - total flux cleaned = -0.0130432 Jy ! Component: 3150 - total flux cleaned = -0.0131046 Jy ! Component: 3200 - total flux cleaned = -0.0131046 Jy ! Component: 3250 - total flux cleaned = -0.0133518 Jy ! Component: 3300 - total flux cleaned = -0.0134142 Jy ! Component: 3350 - total flux cleaned = -0.0136018 Jy ! Component: 3400 - total flux cleaned = -0.0137912 Jy ! Component: 3450 - total flux cleaned = -0.0138543 Jy ! Component: 3500 - total flux cleaned = -0.0141115 Jy ! Component: 3550 - total flux cleaned = -0.0142417 Jy ! Component: 3600 - total flux cleaned = -0.0143724 Jy ! Component: 3650 - total flux cleaned = -0.0145703 Jy ! Component: 3700 - total flux cleaned = -0.0147708 Jy ! Component: 3750 - total flux cleaned = -0.0148373 Jy ! Component: 3800 - total flux cleaned = -0.0151086 Jy ! Component: 3850 - total flux cleaned = -0.0152454 Jy ! Component: 3900 - total flux cleaned = -0.0153834 Jy ! Component: 3950 - total flux cleaned = -0.015664 Jy ! Component: 4000 - total flux cleaned = -0.0157355 Jy ! Component: 4050 - total flux cleaned = -0.0158068 Jy ! Total flux subtracted in 4096 components = -0.0159517 Jy ! Clean residual min=-0.004187 max=0.004674 Jy/beam ! Clean residual mean=-0.000042 rms=0.001249 Jy/beam ! Combined flux in latest and established models = 1.09509 Jy keep ! Adding 282 model components to the UV plane model. ! The established model now contains 460 components and 1.09509 Jy ![Exited script file: muppet] wmap "J0922-3959_S_2014_12_20_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.849 x 21.32 at 6.181 degrees (North through East) ! Clean map min=-0.0099353 max=0.67306 Jy/beam ! Writing clean map to FITS file: J0922-3959_S_2014_12_20_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.673064 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00116258 device "J0922-3959_S_2014_12_20_pus\_map.eps/VPS" ! Attempting to open device: 'J0922-3959_S_2014_12_20_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.518189 0.518189 1.03638 2.07276 4.14551 8.29103 16.5821 33.1641 66.3282 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0034 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.673 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.849 x 21.32 at 6.181 degrees (North through East) ! Clean map min=-0.0099353 max=0.67306 Jy/beam quit ! Quitting program ! Log file J0922-3959_S_2014_12_20_pus_uvs.log closed on Wed Dec 28 13:30:27 2016