! Started logfile: J0016-0015_X_2017_03_27_pus_uvs.log on Sat Apr 27 21:56:08 2019 obs J0016-0015_X_2017_03_27_pus_uva.fits ! Reading UV FITS file: J0016-0015_X_2017_03_27_pus_uva.fits ! AN table 1: 8 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.866667 visibilities/baseline/integration-bin. ! Found source: J0016-0015 ! ! There are 12 IFs, and a total of 12 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.44388e+09 3.2e+07 1 3.2e+07 ! 02 2 8.47588e+09 3.2e+07 1 3.2e+07 ! 03 3 8.50788e+09 3.2e+07 1 3.2e+07 ! 04 4 8.53988e+09 3.2e+07 1 3.2e+07 ! 05 5 8.60388e+09 3.2e+07 1 3.2e+07 ! 06 6 8.63588e+09 3.2e+07 1 3.2e+07 ! 07 7 8.69988e+09 3.2e+07 1 3.2e+07 ! 08 8 8.73188e+09 3.2e+07 1 3.2e+07 ! 09 9 8.79588e+09 3.2e+07 1 3.2e+07 ! 10 10 8.82788e+09 3.2e+07 1 3.2e+07 ! 11 11 8.85988e+09 3.2e+07 1 3.2e+07 ! 12 12 8.89188e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 1844 lines of history. ! ! Reading 3744 visibilities. select I ! Polarization I is unavailable. ! Selecting polarization: RR, channels: 1..12 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: 5..5 ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 ! Reading IF 9 channels: 9..9 ! Reading IF 10 channels: 10..10 ! Reading IF 11 channels: 11..11 ! Reading IF 12 channels: 12..12 print "field_size =",field_size ! field_size = 512 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@muppet J0016-0015_X_2017_03_27_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 = 512x512 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 9. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 10. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 11. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 12. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=0.641131Jy sigma=0.005250 ! Fit after self-cal, rms=0.629222Jy sigma=0.005156 ! 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=0.7396 mas, bmaj=1.803 mas, bpa=-2.407 degrees ! Estimated noise=2101.08 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.282859 Jy ! Component: 100 - total flux cleaned = 0.353669 Jy ! Total flux subtracted in 100 components = 0.353669 Jy ! Clean residual min=-0.012933 max=0.027813 Jy/beam ! Clean residual mean=0.000005 rms=0.004855 Jy/beam ! Combined flux in latest and established models = 0.353669 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.353669 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.103781Jy sigma=0.000819 ! Fit after self-cal, rms=0.103762Jy sigma=0.000818 ! 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=0.9056 mas, bmaj=1.997 mas, bpa=-0.2023 degrees ! Estimated noise=1423.47 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.0384341 Jy ! Component: 100 - total flux cleaned = 0.0572868 Jy ! Component: 150 - total flux cleaned = 0.0673969 Jy ! Component: 200 - total flux cleaned = 0.0733223 Jy ! Total flux subtracted in 200 components = 0.0733223 Jy ! Clean residual min=-0.009065 max=0.011380 Jy/beam ! Clean residual mean=0.000027 rms=0.002571 Jy/beam ! Combined flux in latest and established models = 0.426992 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 43 components and 0.426992 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.089372Jy sigma=0.000702 ! Fit after self-cal, rms=0.089344Jy sigma=0.000701 ! 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.97 HN 1.01 KP 0.96 ! LA 0.94 MK 1.00 NL 0.96 OV 1.05 ! PT 1.02 SC 1.02 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.98 HN 0.99 KP 0.99 ! LA 0.99 MK 1.01 NL 0.99 OV 1.04 ! PT 0.98 SC 0.99 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 0.94 HN 0.96 KP 0.91 ! LA 1.01 MK 1.04 NL 0.99 OV 1.07 ! PT 0.98 SC 1.02 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.93 HN 1.07 KP 0.95 ! LA 1.00 MK 0.99 NL 0.95 OV 1.02 ! PT 1.03 SC 1.04 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.94 HN 0.93 KP 0.95 ! LA 0.97 MK 1.08 NL 1.00 OV 1.02 ! PT 1.01 SC 1.00 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.95 HN 0.96 KP 0.95 ! LA 0.96 MK 1.00 NL 1.01 OV 1.05 ! PT 1.01 SC 0.99 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.93 HN 0.94 KP 0.93 ! LA 1.02 MK 1.02 NL 0.98 OV 1.09 ! PT 1.02 SC 1.00 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.03 HN 0.93 KP 0.97 ! LA 0.95 MK 1.02 NL 0.99 OV 1.05 ! PT 1.03 SC 0.96 ! ! ! Correcting IF 9. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 0.96 HN 0.87 KP 0.94 ! LA 1.00 MK 0.97 NL 0.97 OV 1.07 ! PT 1.05 SC 0.97 ! ! ! Correcting IF 10. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.97 HN 0.99 KP 0.85 ! LA 0.98 MK 1.01 NL 0.86 OV 1.17 ! PT 1.10 SC 0.92 ! ! ! Correcting IF 11. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 FD 1.02 HN 1.04 KP 0.97 ! LA 1.05 MK 1.07 NL 0.91 OV 1.20 ! PT 1.16 SC 1.08 ! ! ! Correcting IF 12. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 0.95 HN 0.95 KP 0.93 ! LA 1.03 MK 1.03 NL 0.92 OV 1.20 ! PT 1.09 SC 1.05 ! ! ! Fit before self-cal, rms=0.089344Jy sigma=0.000701 ! Fit after self-cal, rms=0.082988Jy sigma=0.000654 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9108 mas, bmaj=2.003 mas, bpa=-0.1916 degrees ! Estimated noise=1415.66 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.00223064 Jy ! Component: 100 - total flux cleaned = 0.00254289 Jy ! Component: 150 - total flux cleaned = 0.00274809 Jy ! Component: 200 - total flux cleaned = 0.00293534 Jy ! Total flux subtracted in 200 components = 0.00293534 Jy ! Clean residual min=-0.004184 max=0.005515 Jy/beam ! Clean residual mean=0.000006 rms=0.001271 Jy/beam ! Combined flux in latest and established models = 0.429927 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 48 components and 0.429927 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082842Jy sigma=0.000653 ! Fit after self-cal, rms=0.082816Jy sigma=0.000653 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082816Jy sigma=0.000653 ! Fit after self-cal, rms=0.082753Jy sigma=0.000653 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9098 mas, bmaj=2 mas, bpa=-0.1808 degrees ! Estimated noise=1415.66 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00168801 Jy ! Component: 100 - total flux cleaned = 0.00186817 Jy ! Component: 150 - total flux cleaned = 0.00208053 Jy ! Component: 200 - total flux cleaned = 0.00227168 Jy ! Total flux subtracted in 200 components = 0.00227168 Jy ! Clean residual min=-0.004118 max=0.005272 Jy/beam ! Clean residual mean=0.000006 rms=0.001218 Jy/beam ! Combined flux in latest and established models = 0.432199 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 50 components and 0.432199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082668Jy sigma=0.000652 ! Fit after self-cal, rms=0.082657Jy sigma=0.000652 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 (9.6, 2.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00421394 Jy ! Component: 100 - total flux cleaned = 0.00463301 Jy ! Component: 150 - total flux cleaned = 0.00476587 Jy ! Component: 200 - total flux cleaned = 0.00477362 Jy ! Total flux subtracted in 200 components = 0.00477362 Jy ! Clean residual min=-0.003682 max=0.005005 Jy/beam ! Clean residual mean=0.000006 rms=0.001116 Jy/beam ! Combined flux in latest and established models = 0.436972 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 63 components and 0.436972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082412Jy sigma=0.000650 ! Fit after self-cal, rms=0.082410Jy sigma=0.000650 ! Inverting map ! Added new window around map position (-10, 9). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00465947 Jy ! Component: 100 - total flux cleaned = 0.00619378 Jy ! Component: 150 - total flux cleaned = 0.00695247 Jy ! Component: 200 - total flux cleaned = 0.00742303 Jy ! Total flux subtracted in 200 components = 0.00742303 Jy ! Clean residual min=-0.003496 max=0.003451 Jy/beam ! Clean residual mean=-0.000003 rms=0.000942 Jy/beam ! Combined flux in latest and established models = 0.444395 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 80 components and 0.444395 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082103Jy sigma=0.000648 ! Fit after self-cal, rms=0.082097Jy sigma=0.000648 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.082097Jy sigma=0.000648 ! Fit after self-cal, rms=0.081935Jy sigma=0.000647 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9085 mas, bmaj=1.997 mas, bpa=-0.1467 degrees ! Estimated noise=1414.22 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019202 Jy ! Component: 100 - total flux cleaned = 0.0025231 Jy ! Component: 150 - total flux cleaned = 0.00285176 Jy ! Component: 200 - total flux cleaned = 0.00309966 Jy ! Total flux subtracted in 200 components = 0.00309966 Jy ! Clean residual min=-0.003466 max=0.003399 Jy/beam ! Clean residual mean=-0.000002 rms=0.000886 Jy/beam ! Combined flux in latest and established models = 0.447495 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 92 components and 0.447495 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081845Jy sigma=0.000647 ! Fit after self-cal, rms=0.081843Jy sigma=0.000647 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.000651358 Jy ! Component: 100 - total flux cleaned = 0.000702017 Jy ! Component: 150 - total flux cleaned = 0.000787284 Jy ! Component: 200 - total flux cleaned = 0.00078699 Jy ! Total flux subtracted in 200 components = 0.00078699 Jy ! Clean residual min=-0.003376 max=0.003332 Jy/beam ! Clean residual mean=-0.000003 rms=0.000865 Jy/beam ! Combined flux in latest and established models = 0.448282 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 102 components and 0.448282 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081796Jy sigma=0.000647 ! Fit after self-cal, rms=0.081796Jy sigma=0.000647 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081796Jy sigma=0.000647 ! Fit after self-cal, rms=0.081770Jy sigma=0.000646 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.173 mas, bmaj=3.529 mas, bpa=12.89 degrees ! Estimated noise=1769.66 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000805376 Jy ! Component: 100 - total flux cleaned = 0.00089256 Jy ! Component: 150 - total flux cleaned = 0.000894168 Jy ! Component: 200 - total flux cleaned = 0.000860463 Jy ! Total flux subtracted in 200 components = 0.000860463 Jy ! Clean residual min=-0.004420 max=0.003484 Jy/beam ! Clean residual mean=-0.000007 rms=0.001117 Jy/beam ! Combined flux in latest and established models = 0.449142 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 113 components and 0.449142 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081748Jy sigma=0.000646 ! Fit after self-cal, rms=0.081748Jy sigma=0.000646 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 = -7.75706e-05 Jy ! Component: 100 - total flux cleaned = -0.000108552 Jy ! Component: 150 - total flux cleaned = -0.000107799 Jy ! Component: 200 - total flux cleaned = -0.000106729 Jy ! Total flux subtracted in 200 components = -0.000106729 Jy ! Clean residual min=-0.004328 max=0.003399 Jy/beam ! Clean residual mean=-0.000007 rms=0.001097 Jy/beam ! Combined flux in latest and established models = 0.449036 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 117 components and 0.449036 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081739Jy sigma=0.000646 ! Fit after self-cal, rms=0.081740Jy sigma=0.000646 ! 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=0.9082 mas, bmaj=1.996 mas, bpa=-0.1287 degrees ! Estimated noise=1414.36 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9082 x 1.996 at -0.1287 degrees (North through East) ! Clean map min=-0.0036413 max=0.36844 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 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=0.7398 mas, bmaj=1.801 mas, bpa=-2.429 degrees ! Estimated noise=2092.5 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.281633 Jy ! Component: 100 - total flux cleaned = 0.352218 Jy ! Component: 150 - total flux cleaned = 0.382633 Jy ! Component: 200 - total flux cleaned = 0.400201 Jy ! Total flux subtracted in 200 components = 0.400201 Jy ! Clean residual min=-0.007507 max=0.009966 Jy/beam ! Clean residual mean=-0.000002 rms=0.002318 Jy/beam ! Combined flux in latest and established models = 0.400201 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.400201 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9082 mas, bmaj=1.996 mas, bpa=-0.1287 degrees ! Estimated noise=1414.36 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 (1, -0.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0145218 Jy ! Component: 100 - total flux cleaned = 0.0231935 Jy ! Component: 150 - total flux cleaned = 0.0291347 Jy ! Component: 200 - total flux cleaned = 0.0336037 Jy ! Total flux subtracted in 200 components = 0.0336037 Jy ! Clean residual min=-0.004050 max=0.003603 Jy/beam ! Clean residual mean=0.000010 rms=0.001082 Jy/beam ! Combined flux in latest and established models = 0.433805 Jy ! Adding 46 model components to the UV plane model. ! The established model now contains 65 components and 0.433805 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.173 mas, bmaj=3.529 mas, bpa=12.89 degrees ! Estimated noise=1769.66 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00559747 Jy ! Component: 100 - total flux cleaned = 0.00899794 Jy ! Component: 150 - total flux cleaned = 0.0109161 Jy ! Component: 200 - total flux cleaned = 0.0120448 Jy ! Total flux subtracted in 200 components = 0.0120448 Jy ! Clean residual min=-0.004677 max=0.003394 Jy/beam ! Clean residual mean=-0.000004 rms=0.001119 Jy/beam ! Combined flux in latest and established models = 0.44585 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.000811065 Jy ! Component: 100 - total flux cleaned = 0.00133601 Jy ! Component: 150 - total flux cleaned = 0.00185063 Jy ! Component: 200 - total flux cleaned = 0.00231236 Jy ! Total flux subtracted in 200 components = 0.00231236 Jy ! Clean residual min=-0.004351 max=0.003371 Jy/beam ! Clean residual mean=-0.000008 rms=0.001059 Jy/beam ! Combined flux in latest and established models = 0.448162 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 122 components and 0.448162 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.9082 mas, bmaj=1.996 mas, bpa=-0.1287 degrees ! Estimated noise=1414.36 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000525359 Jy ! Component: 100 - total flux cleaned = 0.00129541 Jy ! Component: 150 - total flux cleaned = 0.00165594 Jy ! Component: 200 - total flux cleaned = 0.00198119 Jy ! Total flux subtracted in 200 components = 0.00198119 Jy ! Clean residual min=-0.003306 max=0.002745 Jy/beam ! Clean residual mean=-0.000006 rms=0.000814 Jy/beam ! Combined flux in latest and established models = 0.450144 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 139 components and 0.450144 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.081743Jy sigma=0.000646 ! Fit after self-cal, rms=0.081730Jy sigma=0.000646 wmodel J0016-0015_X_2017_03_27_pus_map.mod ! Writing 139 model components to file: J0016-0015_X_2017_03_27_pus_map.mod wobs J0016-0015_X_2017_03_27_pus_uvs.fits ! Writing UV FITS file: J0016-0015_X_2017_03_27_pus_uvs.fits wwins J0016-0015_X_2017_03_27_pus_map.win ! wwins: Wrote 4 windows to J0016-0015_X_2017_03_27_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=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00136303 Jy ! Component: 100 - total flux cleaned = -0.00171944 Jy ! Component: 150 - total flux cleaned = -0.00218478 Jy ! Component: 200 - total flux cleaned = -0.00280043 Jy ! Component: 250 - total flux cleaned = -0.00321667 Jy ! Component: 300 - total flux cleaned = -0.00345474 Jy ! Component: 350 - total flux cleaned = -0.00360597 Jy ! Component: 400 - total flux cleaned = -0.00382363 Jy ! Component: 450 - total flux cleaned = -0.00382195 Jy ! Component: 500 - total flux cleaned = -0.00392517 Jy ! Component: 550 - total flux cleaned = -0.00405862 Jy ! Component: 600 - total flux cleaned = -0.00415633 Jy ! Component: 650 - total flux cleaned = -0.00437757 Jy ! Component: 700 - total flux cleaned = -0.00446929 Jy ! Component: 750 - total flux cleaned = -0.00449958 Jy ! Component: 800 - total flux cleaned = -0.00455834 Jy ! Component: 850 - total flux cleaned = -0.00464522 Jy ! Component: 900 - total flux cleaned = -0.00472993 Jy ! Component: 950 - total flux cleaned = -0.00467485 Jy ! Component: 1000 - total flux cleaned = -0.00467561 Jy ! Component: 1050 - total flux cleaned = -0.00467554 Jy ! Component: 1100 - total flux cleaned = -0.00464941 Jy ! Component: 1150 - total flux cleaned = -0.00459778 Jy ! Component: 1200 - total flux cleaned = -0.00452199 Jy ! Component: 1250 - total flux cleaned = -0.00447219 Jy ! Component: 1300 - total flux cleaned = -0.00439917 Jy ! Component: 1350 - total flux cleaned = -0.0042798 Jy ! Component: 1400 - total flux cleaned = -0.0042099 Jy ! Component: 1450 - total flux cleaned = -0.00409438 Jy ! Component: 1500 - total flux cleaned = -0.0041173 Jy ! Component: 1550 - total flux cleaned = -0.00404981 Jy ! Component: 1600 - total flux cleaned = -0.00402738 Jy ! Component: 1650 - total flux cleaned = -0.00391785 Jy ! Component: 1700 - total flux cleaned = -0.00385372 Jy ! Component: 1750 - total flux cleaned = -0.00379008 Jy ! Component: 1800 - total flux cleaned = -0.0037269 Jy ! Component: 1850 - total flux cleaned = -0.00368548 Jy ! Component: 1900 - total flux cleaned = -0.00358361 Jy ! Component: 1950 - total flux cleaned = -0.00350335 Jy ! Component: 2000 - total flux cleaned = -0.00356316 Jy ! Total flux subtracted in 2048 components = -0.00346443 Jy ! Clean residual min=-0.001008 max=0.001082 Jy/beam ! Clean residual mean=-0.000000 rms=0.000357 Jy/beam ! Combined flux in latest and established models = 0.446679 Jy keep ! Adding 120 model components to the UV plane model. ! The established model now contains 259 components and 0.446679 Jy ![Exited script file: muppet] wmap "J0016-0015_X_2017_03_27_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9082 x 1.996 at -0.1287 degrees (North through East) ! Clean map min=-0.0025895 max=0.36816 Jy/beam ! Writing clean map to FITS file: J0016-0015_X_2017_03_27_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.368156 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000354546 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.28891 0.28891 0.577819 1.15564 2.31128 4.62255 9.24511 18.4902 36.9804 73.9608 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.001 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.368 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9082 x 1.996 at -0.1287 degrees (North through East) ! Clean map min=-0.0025895 max=0.36816 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 J0016-0015_X_2017_03_27_pus_uvs.log closed on Sat Apr 27 21:56:09 2019