! Started logfile: J2258-2758_X_2014_08_05_pus_uvs.log on Sat Nov 12 21:17:27 2016 obs J2258-2758_X_2014_08_05_pus_uva.fits ! Reading UV FITS file: J2258-2758_X_2014_08_05_pus_uva.fits ! AN table 1: 16 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.604167 visibilities/baseline/integration-bin. ! Found source: J2258-2758 ! ! 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 414 lines of history. ! ! Reading 5220 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 J2258-2758_X_2014_08_05_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 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. ! ! Correcting IF 5. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 9. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 10. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 11. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 12. ! A total of 34 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=0.921679Jy sigma=0.009369 ! Fit after self-cal, rms=0.884710Jy sigma=0.008959 ! 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.7925 mas, bmaj=6.891 mas, bpa=30.98 degrees ! Estimated noise=1347.88 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 = 1.32875 Jy ! Component: 100 - total flux cleaned = 1.64525 Jy ! Total flux subtracted in 100 components = 1.64525 Jy ! Clean residual min=-0.114925 max=0.135671 Jy/beam ! Clean residual mean=0.000180 rms=0.031516 Jy/beam ! Combined flux in latest and established models = 1.64525 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 1.64525 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.403379Jy sigma=0.003882 ! Fit after self-cal, rms=0.403397Jy sigma=0.003882 ! 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=1.041 mas, bmaj=6.645 mas, bpa=30.06 degrees ! Estimated noise=964.054 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.191432 Jy ! Component: 100 - total flux cleaned = 0.298697 Jy ! Component: 150 - total flux cleaned = 0.352968 Jy ! Component: 200 - total flux cleaned = 0.380014 Jy ! Total flux subtracted in 200 components = 0.380014 Jy ! Clean residual min=-0.123944 max=0.076563 Jy/beam ! Clean residual mean=-0.000343 rms=0.030795 Jy/beam ! Combined flux in latest and established models = 2.02527 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 65 components and 2.02527 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.282345Jy sigma=0.002750 ! Fit after self-cal, rms=0.282140Jy sigma=0.002748 ! 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 0.91 FD 1.09 HN 1.00* KP 0.92 ! LA 1.01 MK 1.00* NL 1.03 OV 0.98 ! PT 1.02 SC 1.02 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.09 HN 1.00* KP 0.87 ! LA 1.02 MK 1.00* NL 1.03 OV 0.98 ! PT 0.97 SC 1.03 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.91 FD 1.10 HN 1.00* KP 0.83 ! LA 1.05 MK 1.00* NL 1.02 OV 0.98 ! PT 0.95 SC 1.03 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.91 FD 1.08 HN 1.00* KP 0.79 ! LA 1.06 MK 1.00* NL 1.03 OV 0.97 ! PT 0.96 SC 1.03 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.08 HN 1.00* KP 0.75 ! LA 1.05 MK 1.00* NL 1.04 OV 0.98 ! PT 0.97 SC 1.03 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.14 HN 1.00* KP 0.75 ! LA 1.06 MK 1.00* NL 1.05 OV 1.00 ! PT 0.96 SC 1.00 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.22 HN 1.00* KP 0.80 ! LA 1.15 MK 1.00* NL 1.10 OV 1.06 ! PT 1.04 SC 1.11 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.15 HN 1.00* KP 0.75 ! LA 1.05 MK 1.00* NL 1.05 OV 1.02 ! PT 0.97 SC 1.05 ! ! ! Correcting IF 9. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.12 HN 1.00* KP 0.77 ! LA 1.08 MK 1.00* NL 1.06 OV 1.02 ! PT 0.98 SC 1.00 ! ! ! Correcting IF 10. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.11 HN 1.00* KP 0.80 ! LA 1.10 MK 1.00* NL 1.07 OV 1.02 ! PT 0.98 SC 1.00 ! ! ! Correcting IF 11. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.12 HN 1.00* KP 0.84 ! LA 1.11 MK 1.00* NL 1.07 OV 1.00 ! PT 0.99 SC 1.04 ! ! ! Correcting IF 12. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.12 HN 1.00* KP 0.87 ! LA 1.09 MK 1.00* NL 1.06 OV 1.01 ! PT 1.00 SC 1.02 ! ! ! Fit before self-cal, rms=0.282140Jy sigma=0.002748 ! Fit after self-cal, rms=0.096904Jy sigma=0.000951 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.049 mas, bmaj=6.726 mas, bpa=29.96 degrees ! Estimated noise=952.627 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.0159429 Jy ! Component: 100 - total flux cleaned = -0.0159711 Jy ! Component: 150 - total flux cleaned = -0.0145833 Jy ! Component: 200 - total flux cleaned = -0.0125556 Jy ! Total flux subtracted in 200 components = -0.0125556 Jy ! Clean residual min=-0.023233 max=0.012417 Jy/beam ! Clean residual mean=-0.000066 rms=0.004531 Jy/beam ! Combined flux in latest and established models = 2.01271 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 84 components and 2.01271 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.090556Jy sigma=0.000887 ! Fit after self-cal, rms=0.090453Jy sigma=0.000885 ! 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.090453Jy sigma=0.000885 ! Fit after self-cal, rms=0.090243Jy sigma=0.000882 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.048 mas, bmaj=6.74 mas, bpa=29.96 degrees ! Estimated noise=952.253 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00349965 Jy ! Component: 100 - total flux cleaned = -0.00247637 Jy ! Component: 150 - total flux cleaned = -0.00186213 Jy ! Component: 200 - total flux cleaned = -0.001319 Jy ! Total flux subtracted in 200 components = -0.001319 Jy ! Clean residual min=-0.019582 max=0.010697 Jy/beam ! Clean residual mean=-0.000054 rms=0.003755 Jy/beam ! Combined flux in latest and established models = 2.01139 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 96 components and 2.01139 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.087794Jy sigma=0.000859 ! Fit after self-cal, rms=0.087759Jy sigma=0.000858 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.000999087 Jy ! Component: 100 - total flux cleaned = 0.000780531 Jy ! Component: 150 - total flux cleaned = 0.000578423 Jy ! Component: 200 - total flux cleaned = 0.000394678 Jy ! Total flux subtracted in 200 components = 0.000394678 Jy ! Clean residual min=-0.017723 max=0.010045 Jy/beam ! Clean residual mean=-0.000053 rms=0.003514 Jy/beam ! Combined flux in latest and established models = 2.01179 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 104 components and 2.01179 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.086731Jy sigma=0.000849 ! Fit after self-cal, rms=0.086714Jy sigma=0.000849 ! 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.086714Jy sigma=0.000849 ! Fit after self-cal, rms=0.086414Jy sigma=0.000846 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.047 mas, bmaj=6.751 mas, bpa=29.97 degrees ! Estimated noise=951.475 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00399652 Jy ! Component: 100 - total flux cleaned = -0.00380362 Jy ! Component: 150 - total flux cleaned = -0.00397501 Jy ! Component: 200 - total flux cleaned = -0.00380892 Jy ! Total flux subtracted in 200 components = -0.00380892 Jy ! Clean residual min=-0.014923 max=0.009487 Jy/beam ! Clean residual mean=-0.000040 rms=0.003076 Jy/beam ! Combined flux in latest and established models = 2.00798 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 115 components and 2.00798 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.085567Jy sigma=0.000837 ! Fit after self-cal, rms=0.085557Jy sigma=0.000837 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.00014724 Jy ! Component: 100 - total flux cleaned = 0.000146267 Jy ! Component: 150 - total flux cleaned = 0.000145626 Jy ! Component: 200 - total flux cleaned = 0.000144638 Jy ! Total flux subtracted in 200 components = 0.000144638 Jy ! Clean residual min=-0.014054 max=0.009143 Jy/beam ! Clean residual mean=-0.000040 rms=0.002972 Jy/beam ! Combined flux in latest and established models = 2.00812 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 121 components and 2.00812 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.085173Jy sigma=0.000834 ! Fit after self-cal, rms=0.085168Jy sigma=0.000834 ! 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.085168Jy sigma=0.000834 ! Fit after self-cal, rms=0.085028Jy sigma=0.000832 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.512 mas, bmaj=8.499 mas, bpa=18.32 degrees ! Estimated noise=1059.28 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00151174 Jy ! Component: 100 - total flux cleaned = -0.0012272 Jy ! Component: 150 - total flux cleaned = -0.000839041 Jy ! Component: 200 - total flux cleaned = -0.000562533 Jy ! Total flux subtracted in 200 components = -0.000562533 Jy ! Clean residual min=-0.013104 max=0.009639 Jy/beam ! Clean residual mean=-0.000053 rms=0.003379 Jy/beam ! Combined flux in latest and established models = 2.00756 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 137 components and 2.00756 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.084863Jy sigma=0.000830 ! Fit after self-cal, rms=0.084856Jy sigma=0.000830 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.000361257 Jy ! Component: 100 - total flux cleaned = -0.000235747 Jy ! Component: 150 - total flux cleaned = -6.62857e-05 Jy ! Component: 200 - total flux cleaned = 9.4733e-05 Jy ! Total flux subtracted in 200 components = 9.4733e-05 Jy ! Clean residual min=-0.012787 max=0.009685 Jy/beam ! Clean residual mean=-0.000053 rms=0.003324 Jy/beam ! Combined flux in latest and established models = 2.00766 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 146 components and 2.00766 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.084789Jy sigma=0.000829 ! Fit after self-cal, rms=0.084789Jy sigma=0.000829 ! 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=1.047 mas, bmaj=6.761 mas, bpa=29.97 degrees ! Estimated noise=950.923 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.047 x 6.761 at 29.97 degrees (North through East) ! Clean map min=-0.011547 max=1.7513 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. ! ! 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.084789Jy sigma=0.000829 ! Fit after self-cal, rms=0.075417Jy sigma=0.000730 ! Inverting map and beam ! Estimated beam: bmin=1.046 mas, bmaj=6.771 mas, bpa=29.97 degrees ! Estimated noise=950.269 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00436315 Jy ! Component: 100 - total flux cleaned = -0.00562661 Jy ! Component: 150 - total flux cleaned = -0.00562607 Jy ! Component: 200 - total flux cleaned = -0.00562609 Jy ! Total flux subtracted in 200 components = -0.00562609 Jy ! Clean residual min=-0.010834 max=0.008835 Jy/beam ! Clean residual mean=-0.000030 rms=0.002503 Jy/beam ! Combined flux in latest and established models = 2.00203 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 150 components and 2.00203 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.074992Jy sigma=0.000726 ! Fit after self-cal, rms=0.074989Jy sigma=0.000726 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 9.58977e-05 Jy ! Component: 100 - total flux cleaned = 0.000182299 Jy ! Component: 150 - total flux cleaned = 0.000183979 Jy ! Component: 200 - total flux cleaned = 0.000342855 Jy ! Total flux subtracted in 200 components = 0.000342855 Jy ! Clean residual min=-0.010415 max=0.008833 Jy/beam ! Clean residual mean=-0.000030 rms=0.002454 Jy/beam ! Combined flux in latest and established models = 2.00237 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 155 components and 2.00237 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.074803Jy sigma=0.000724 ! Fit after self-cal, rms=0.074803Jy sigma=0.000724 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Correcting IF 9. ! ! Correcting IF 10. ! ! Correcting IF 11. ! ! Correcting IF 12. ! ! Fit before self-cal, rms=0.074803Jy sigma=0.000724 ! Fit after self-cal, rms=0.074737Jy sigma=0.000723 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.79 mas, bmaj=7.064 mas, bpa=30.88 degrees ! Estimated noise=1341.22 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.32203 Jy ! Component: 100 - total flux cleaned = 1.63259 Jy ! Component: 150 - total flux cleaned = 1.75752 Jy ! Component: 200 - total flux cleaned = 1.82505 Jy ! Total flux subtracted in 200 components = 1.82505 Jy ! Clean residual min=-0.018889 max=0.039230 Jy/beam ! Clean residual mean=0.000112 rms=0.008420 Jy/beam ! Combined flux in latest and established models = 1.82505 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 35 model components to the UV plane model. ! The established model now contains 35 components and 1.82505 Jy ! Inverting map and beam ! Estimated beam: bmin=1.046 mas, bmaj=6.779 mas, bpa=29.98 degrees ! Estimated noise=949.838 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 (3.8, -15.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0630049 Jy ! Component: 100 - total flux cleaned = 0.100578 Jy ! Component: 150 - total flux cleaned = 0.125212 Jy ! Component: 200 - total flux cleaned = 0.142314 Jy ! Total flux subtracted in 200 components = 0.142314 Jy ! Clean residual min=-0.009479 max=0.013982 Jy/beam ! Clean residual mean=0.000020 rms=0.003727 Jy/beam ! Combined flux in latest and established models = 1.96736 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 100 components and 1.96736 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.514 mas, bmaj=8.513 mas, bpa=18.37 degrees ! Estimated noise=1058.05 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0165446 Jy ! Component: 100 - total flux cleaned = 0.0255526 Jy ! Component: 150 - total flux cleaned = 0.0284377 Jy ! Component: 200 - total flux cleaned = 0.030569 Jy ! Total flux subtracted in 200 components = 0.030569 Jy ! Clean residual min=-0.008206 max=0.009655 Jy/beam ! Clean residual mean=-0.000025 rms=0.002711 Jy/beam ! Combined flux in latest and established models = 1.99793 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.00157308 Jy ! Component: 100 - total flux cleaned = 0.00234479 Jy ! Component: 150 - total flux cleaned = 0.00218582 Jy ! Component: 200 - total flux cleaned = 0.00172761 Jy ! Total flux subtracted in 200 components = 0.00172761 Jy ! Clean residual min=-0.007353 max=0.008019 Jy/beam ! Clean residual mean=-0.000027 rms=0.002205 Jy/beam ! Combined flux in latest and established models = 1.99966 Jy ! Adding 61 model components to the UV plane model. ! The established model now contains 154 components and 1.99966 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.046 mas, bmaj=6.779 mas, bpa=29.98 degrees ! Estimated noise=949.838 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.001085 Jy ! Component: 100 - total flux cleaned = 0.00108065 Jy ! Component: 150 - total flux cleaned = 0.000930128 Jy ! Component: 200 - total flux cleaned = 0.000793673 Jy ! Total flux subtracted in 200 components = 0.000793673 Jy ! Clean residual min=-0.005347 max=0.006021 Jy/beam ! Clean residual mean=-0.000023 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 2.00045 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 184 components and 2.00045 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.073408Jy sigma=0.000709 ! Fit after self-cal, rms=0.073245Jy sigma=0.000707 wmodel J2258-2758_X_2014_08_05_pus_map.mod ! Writing 184 model components to file: J2258-2758_X_2014_08_05_pus_map.mod wobs J2258-2758_X_2014_08_05_pus_uvs.fits ! Writing UV FITS file: J2258-2758_X_2014_08_05_pus_uvs.fits wwins J2258-2758_X_2014_08_05_pus_map.win ! wwins: Wrote 2 windows to J2258-2758_X_2014_08_05_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.00101688 Jy ! Component: 100 - total flux cleaned = 0.00102198 Jy ! Component: 150 - total flux cleaned = 0.00149486 Jy ! Component: 200 - total flux cleaned = 0.00199901 Jy ! Component: 250 - total flux cleaned = 0.00240524 Jy ! Component: 300 - total flux cleaned = 0.00273237 Jy ! Component: 350 - total flux cleaned = 0.00298939 Jy ! Component: 400 - total flux cleaned = 0.00323234 Jy ! Component: 450 - total flux cleaned = 0.00341676 Jy ! Component: 500 - total flux cleaned = 0.00368186 Jy ! Component: 550 - total flux cleaned = 0.00389303 Jy ! Component: 600 - total flux cleaned = 0.00401347 Jy ! Component: 650 - total flux cleaned = 0.00405229 Jy ! Component: 700 - total flux cleaned = 0.00427651 Jy ! Component: 750 - total flux cleaned = 0.0044924 Jy ! Component: 800 - total flux cleaned = 0.00463101 Jy ! Component: 850 - total flux cleaned = 0.00466436 Jy ! Component: 900 - total flux cleaned = 0.00476248 Jy ! Component: 950 - total flux cleaned = 0.00479355 Jy ! Component: 1000 - total flux cleaned = 0.00476296 Jy ! Component: 1050 - total flux cleaned = 0.00491217 Jy ! Component: 1100 - total flux cleaned = 0.00494092 Jy ! Component: 1150 - total flux cleaned = 0.00505407 Jy ! Component: 1200 - total flux cleaned = 0.00516395 Jy ! Component: 1250 - total flux cleaned = 0.00527086 Jy ! Component: 1300 - total flux cleaned = 0.00534915 Jy ! Component: 1350 - total flux cleaned = 0.00539987 Jy ! Component: 1400 - total flux cleaned = 0.00539955 Jy ! Component: 1450 - total flux cleaned = 0.00542396 Jy ! Component: 1500 - total flux cleaned = 0.00540011 Jy ! Component: 1550 - total flux cleaned = 0.00547076 Jy ! Component: 1600 - total flux cleaned = 0.00549402 Jy ! Component: 1650 - total flux cleaned = 0.00549397 Jy ! Component: 1700 - total flux cleaned = 0.00549378 Jy ! Component: 1750 - total flux cleaned = 0.00547112 Jy ! Component: 1800 - total flux cleaned = 0.00549326 Jy ! Component: 1850 - total flux cleaned = 0.00549336 Jy ! Component: 1900 - total flux cleaned = 0.00553723 Jy ! Component: 1950 - total flux cleaned = 0.00555892 Jy ! Component: 2000 - total flux cleaned = 0.00555886 Jy ! Total flux subtracted in 2048 components = 0.00555905 Jy ! Clean residual min=-0.001182 max=0.001280 Jy/beam ! Clean residual mean=-0.000018 rms=0.000414 Jy/beam ! Combined flux in latest and established models = 2.00601 Jy keep ! Adding 128 model components to the UV plane model. ! The established model now contains 312 components and 2.00601 Jy ![Exited script file: muppet] wmap "J2258-2758_X_2014_08_05_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.046 x 6.779 at 29.98 degrees (North through East) ! Clean map min=-0.0050973 max=1.7457 Jy/beam ! Writing clean map to FITS file: J2258-2758_X_2014_08_05_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.74568 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000393335 device "J2258-2758_X_2014_08_05_pus\_map.eps/VPS" ! Attempting to open device: 'J2258-2758_X_2014_08_05_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.0675957 0.0675957 0.135191 0.270383 0.540765 1.08153 2.16306 4.32612 8.65225 17.3045 34.609 69.218 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0011 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.745 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.046 x 6.779 at 29.98 degrees (North through East) ! Clean map min=-0.0050973 max=1.7457 Jy/beam quit ! Quitting program ! Log file J2258-2758_X_2014_08_05_pus_uvs.log closed on Sat Nov 12 21:17:29 2016