! Started logfile: J1751+0939_X_2011_06_28_pus_uvs.log on Sat Apr 21 13:06:56 2012 obs J1751+0939_X_2011_06_28_pus_uva.fits ! Reading UV FITS file: J1751+0939_X_2011_06_28_pus_uva.fits ! AN table 1: 139 integrations on 153 of 153 possible baselines. ! AN table 2: 112 integrations on 153 of 153 possible baselines. ! AN table 3: 42 integrations on 153 of 153 possible baselines. ! Apparent sampling: 0.0718285 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1751+0939 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40997e+09 8e+06 1 8e+06 ! 02 2 8.47997e+09 8e+06 1 8e+06 ! 03 3 8.79497e+09 8e+06 1 8e+06 ! 04 4 8.89997e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 5249 lines of history. ! ! Reading 12880 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 = 512 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@muppet J1751+0939_X_2011_06_28_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 1717 telescope corrections were flagged in sub-array 1. ! A total of 1627 telescope corrections were flagged in sub-array 2. ! A total of 672 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 1717 telescope corrections were flagged in sub-array 1. ! A total of 1627 telescope corrections were flagged in sub-array 2. ! A total of 672 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 1717 telescope corrections were flagged in sub-array 1. ! A total of 1627 telescope corrections were flagged in sub-array 2. ! A total of 672 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 1717 telescope corrections were flagged in sub-array 1. ! A total of 1627 telescope corrections were flagged in sub-array 2. ! A total of 672 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=3.754114Jy sigma=0.035659 ! Fit after self-cal, rms=3.731781Jy sigma=0.035633 ! 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 ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.5871 mas, bmaj=0.9209 mas, bpa=2.274 degrees ! Estimated noise=1045.37 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 = 3.64296 Jy ! Component: 100 - total flux cleaned = 4.48887 Jy ! Total flux subtracted in 100 components = 4.48887 Jy ! Clean residual min=-0.085430 max=0.307725 Jy/beam ! Clean residual mean=0.000403 rms=0.022882 Jy/beam ! Combined flux in latest and established models = 4.48887 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 4.48887 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.688113Jy sigma=0.006563 ! Fit after self-cal, rms=0.688121Jy sigma=0.006563 ! Inverting map ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.318025 Jy ! Component: 100 - total flux cleaned = 0.495979 Jy ! Total flux subtracted in 100 components = 0.495979 Jy ! Clean residual min=-0.070888 max=0.093351 Jy/beam ! Clean residual mean=0.000099 rms=0.014336 Jy/beam ! Combined flux in latest and established models = 4.98485 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 17 components and 4.98485 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.508889Jy sigma=0.004102 ! Fit after self-cal, rms=0.508267Jy sigma=0.004099 ! Inverting map ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.109941 Jy ! Component: 100 - total flux cleaned = 0.180658 Jy ! Total flux subtracted in 100 components = 0.180658 Jy ! Clean residual min=-0.066911 max=0.068310 Jy/beam ! Clean residual mean=-0.000012 rms=0.012674 Jy/beam ! Combined flux in latest and established models = 5.16551 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 23 components and 5.16551 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.482645Jy sigma=0.003586 ! Fit after self-cal, rms=0.482297Jy sigma=0.003584 ! Inverting map ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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 ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7103 mas, bmaj=1.076 mas, bpa=0.9192 degrees ! Estimated noise=756.928 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.0783989 Jy ! Component: 100 - total flux cleaned = 0.118702 Jy ! Component: 150 - total flux cleaned = 0.141623 Jy ! Component: 200 - total flux cleaned = 0.151101 Jy ! Total flux subtracted in 200 components = 0.151101 Jy ! Clean residual min=-0.044007 max=0.053011 Jy/beam ! Clean residual mean=-0.000402 rms=0.010240 Jy/beam ! Combined flux in latest and established models = 5.31661 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 27 components and 5.31661 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.479393Jy sigma=0.003372 ! Fit after self-cal, rms=0.478995Jy sigma=0.003370 ! Inverting map ! Your choice of large map pixels excluded 5.45% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 90 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.00 FT 1.00* HH 1.06 ! HN 1.00 KK 0.99 KP 1.01 MC 1.01 ! MK 1.00 NL 0.99 NY 1.05 OV 0.96 ! PT 0.98 SC 0.99 SH 1.18 TS 1.04 ! WF 1.00 WZ 0.96 ! ! A total of 57 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 1.07 FT 1.00* HH 1.02 ! HN 1.27 KK 0.99 KP 1.05 MC 1.01 ! MK 0.99 NL 1.00* NY 1.06 OV 1.00 ! PT 1.01 SC 1.01 SH 0.96 TS 0.98 ! WF 1.06 WZ 0.98 ! ! A total of 84 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* FT 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* MC 0.00* ! MK 0.00* NL 0.00* NY 0.00* OV 0.00* ! PT 0.00* SC 0.00* SH 0.00* TS 0.00* ! WF 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! A total of 90 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.97 FT 1.00* HH 1.06 ! HN 1.01 KK 0.96 KP 0.99 MC 1.01 ! MK 0.97 NL 1.03 NY 1.02 OV 0.98 ! PT 1.00 SC 1.03 SH 1.20 TS 1.06 ! WF 0.98 WZ 0.95 ! ! A total of 57 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 1.03 FT 1.00* HH 1.01 ! HN 1.31 KK 0.98 KP 1.03 MC 1.01 ! MK 0.96 NL 1.00* NY 1.03 OV 1.02 ! PT 1.03 SC 1.04 SH 1.00 TS 0.99 ! WF 1.01 WZ 0.98 ! ! A total of 84 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* FT 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* MC 0.00* ! MK 0.00* NL 0.00* NY 0.00* OV 0.00* ! PT 0.00* SC 0.00* SH 0.00* TS 0.00* ! WF 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! A total of 90 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.00 FT 1.00* HH 1.11 ! HN 1.00 KK 1.02 KP 0.99 MC 1.01 ! MK 0.95 NL 0.97 NY 1.05 OV 0.97 ! PT 0.97 SC 1.03 SH 1.28 TS 1.02 ! WF 1.02 WZ 0.94 ! ! A total of 57 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 1.07 FT 1.00* HH 1.08 ! HN 1.30 KK 1.03 KP 1.03 MC 1.01 ! MK 0.94 NL 1.00* NY 1.05 OV 1.01 ! PT 1.01 SC 1.07 SH 1.05 TS 0.98 ! WF 1.05 WZ 0.98 ! ! A total of 84 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* FT 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* MC 0.00* ! MK 0.00* NL 0.00* NY 0.00* OV 0.00* ! PT 0.00* SC 0.00* SH 0.00* TS 0.00* ! WF 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! A total of 90 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00 FT 1.00* HH 1.07 ! HN 1.01 KK 1.02 KP 0.96 MC 0.98 ! MK 1.01 NL 1.00 NY 1.03 OV 1.01 ! PT 0.99 SC 1.04 SH 1.22 TS 1.02 ! WF 0.99 WZ 1.00 ! ! A total of 57 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 1.08 FT 1.00* HH 1.02 ! HN 1.29 KK 1.02 KP 1.00 MC 0.98 ! MK 1.00 NL 1.00* NY 1.05 OV 1.05 ! PT 1.02 SC 1.05 SH 0.98 TS 0.98 ! WF 1.03 WZ 1.03 ! ! A total of 84 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* FT 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* MC 0.00* ! MK 0.00* NL 0.00* NY 0.00* OV 0.00* ! PT 0.00* SC 0.00* SH 0.00* TS 0.00* ! WF 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.478995Jy sigma=0.003370 ! Fit after self-cal, rms=0.405270Jy sigma=0.002696 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.719 mas, bmaj=1.128 mas, bpa=-2.094 degrees ! Estimated noise=772.315 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 (2, 5.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0131196 Jy ! Component: 100 - total flux cleaned = 0.00757706 Jy ! Component: 150 - total flux cleaned = 0.0058797 Jy ! Component: 200 - total flux cleaned = 0.00517797 Jy ! Total flux subtracted in 200 components = 0.00517797 Jy ! Clean residual min=-0.031882 max=0.045043 Jy/beam ! Clean residual mean=0.001279 rms=0.008065 Jy/beam ! Combined flux in latest and established models = 5.32179 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 37 components and 5.32179 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.397500Jy sigma=0.002651 ! Fit after self-cal, rms=0.397310Jy sigma=0.002645 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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.397310Jy sigma=0.002645 ! Fit after self-cal, rms=0.271597Jy sigma=0.001745 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7193 mas, bmaj=1.125 mas, bpa=-2.171 degrees ! Estimated noise=770.982 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0112059 Jy ! Component: 100 - total flux cleaned = -0.0112056 Jy ! Component: 150 - total flux cleaned = -0.0112103 Jy ! Component: 200 - total flux cleaned = -0.0119179 Jy ! Total flux subtracted in 200 components = -0.0119179 Jy ! Clean residual min=-0.015652 max=0.014771 Jy/beam ! Clean residual mean=-0.000476 rms=0.003263 Jy/beam ! Combined flux in latest and established models = 5.30987 Jy selfcal ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 42 components and 5.30987 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.270068Jy sigma=0.001731 ! Fit after self-cal, rms=0.269880Jy sigma=0.001730 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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.00255004 Jy ! Component: 100 - total flux cleaned = -0.00222392 Jy ! Component: 150 - total flux cleaned = -0.00222462 Jy ! Component: 200 - total flux cleaned = -0.00222302 Jy ! Total flux subtracted in 200 components = -0.00222302 Jy ! Clean residual min=-0.015576 max=0.014363 Jy/beam ! Clean residual mean=-0.000437 rms=0.003197 Jy/beam ! Combined flux in latest and established models = 5.30765 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 45 components and 5.30764 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.269151Jy sigma=0.001723 ! Fit after self-cal, rms=0.268998Jy sigma=0.001723 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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.268998Jy sigma=0.001723 ! Fit after self-cal, rms=0.250055Jy sigma=0.001588 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7202 mas, bmaj=1.126 mas, bpa=-2.214 degrees ! Estimated noise=771.448 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00421407 Jy ! Component: 100 - total flux cleaned = -0.00529219 Jy ! Component: 150 - total flux cleaned = -0.00594915 Jy ! Component: 200 - total flux cleaned = -0.00655637 Jy ! Total flux subtracted in 200 components = -0.00655637 Jy ! Clean residual min=-0.011777 max=0.009550 Jy/beam ! Clean residual mean=-0.000127 rms=0.002116 Jy/beam ! Combined flux in latest and established models = 5.30109 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 48 components and 5.30109 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.249461Jy sigma=0.001581 ! Fit after self-cal, rms=0.249345Jy sigma=0.001581 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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.00145562 Jy ! Component: 100 - total flux cleaned = -0.0020043 Jy ! Component: 150 - total flux cleaned = -0.0025157 Jy ! Component: 200 - total flux cleaned = -0.00299874 Jy ! Total flux subtracted in 200 components = -0.00299874 Jy ! Clean residual min=-0.010312 max=0.009335 Jy/beam ! Clean residual mean=-0.000076 rms=0.002052 Jy/beam ! Combined flux in latest and established models = 5.29809 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 51 components and 5.29809 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.248907Jy sigma=0.001576 ! Fit after self-cal, rms=0.248795Jy sigma=0.001576 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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. ! ! Fit before self-cal, rms=0.248795Jy sigma=0.001576 ! Fit after self-cal, rms=0.248949Jy sigma=0.001573 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=2.102 mas, bmaj=2.649 mas, bpa=-17.78 degrees ! Estimated noise=1037.19 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00464675 Jy ! Component: 100 - total flux cleaned = -0.00514374 Jy ! Component: 150 - total flux cleaned = -0.00555927 Jy ! Component: 200 - total flux cleaned = -0.00574225 Jy ! Total flux subtracted in 200 components = -0.00574225 Jy ! Clean residual min=-0.009047 max=0.011689 Jy/beam ! Clean residual mean=0.000019 rms=0.002848 Jy/beam ! Combined flux in latest and established models = 5.29235 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 57 components and 5.29235 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.248772Jy sigma=0.001570 ! Fit after self-cal, rms=0.248530Jy sigma=0.001568 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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 = -2.28865e-05 Jy ! Component: 100 - total flux cleaned = -0.000242666 Jy ! Component: 150 - total flux cleaned = -0.000592869 Jy ! Component: 200 - total flux cleaned = -0.000890863 Jy ! Total flux subtracted in 200 components = -0.000890863 Jy ! Clean residual min=-0.008960 max=0.011256 Jy/beam ! Clean residual mean=0.000064 rms=0.002798 Jy/beam ! Combined flux in latest and established models = 5.29146 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 60 components and 5.29146 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.248359Jy sigma=0.001566 ! Fit after self-cal, rms=0.248196Jy sigma=0.001565 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec 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 ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7207 mas, bmaj=1.128 mas, bpa=-2.286 degrees ! Estimated noise=772.01 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7207 x 1.128 at -2.286 degrees (North through East) ! Clean map min=-0.0081193 max=4.7753 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.248196Jy sigma=0.001565 ! Fit after self-cal, rms=0.210719Jy sigma=0.001283 ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7207 mas, bmaj=1.128 mas, bpa=-2.29 degrees ! Estimated noise=771.502 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00210016 Jy ! Component: 100 - total flux cleaned = -0.0033115 Jy ! Component: 150 - total flux cleaned = -0.00469373 Jy ! Component: 200 - total flux cleaned = -0.00571687 Jy ! Total flux subtracted in 200 components = -0.00571687 Jy ! Clean residual min=-0.010178 max=0.008347 Jy/beam ! Clean residual mean=0.000044 rms=0.001831 Jy/beam ! Combined flux in latest and established models = 5.28574 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 61 components and 5.28574 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210371Jy sigma=0.001279 ! Fit after self-cal, rms=0.210345Jy sigma=0.001279 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00106609 Jy ! Component: 100 - total flux cleaned = -0.00151942 Jy ! Component: 150 - total flux cleaned = -0.00194694 Jy ! Component: 200 - total flux cleaned = -0.00234368 Jy ! Total flux subtracted in 200 components = -0.00234368 Jy ! Clean residual min=-0.008847 max=0.008367 Jy/beam ! Clean residual mean=0.000085 rms=0.001796 Jy/beam ! Combined flux in latest and established models = 5.2834 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 62 components and 5.2834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210174Jy sigma=0.001276 ! Fit after self-cal, rms=0.210136Jy sigma=0.001276 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210136Jy sigma=0.001276 ! Fit after self-cal, rms=0.210085Jy sigma=0.001274 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 ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.5999 mas, bmaj=0.9698 mas, bpa=-1.935 degrees ! Estimated noise=1076.55 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.71122 Jy ! Component: 100 - total flux cleaned = 4.54695 Jy ! Component: 150 - total flux cleaned = 4.83572 Jy ! Component: 200 - total flux cleaned = 4.99134 Jy ! Total flux subtracted in 200 components = 4.99134 Jy ! Clean residual min=-0.024531 max=0.078076 Jy/beam ! Clean residual mean=0.000177 rms=0.006243 Jy/beam ! Combined flux in latest and established models = 4.99134 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 16 model components to the UV plane model. ! The established model now contains 16 components and 4.99134 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.7211 mas, bmaj=1.13 mas, bpa=-2.367 degrees ! Estimated noise=772.016 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.121942 Jy ! Component: 100 - total flux cleaned = 0.18391 Jy ! Component: 150 - total flux cleaned = 0.217819 Jy ! Component: 200 - total flux cleaned = 0.238671 Jy ! Total flux subtracted in 200 components = 0.238671 Jy ! Clean residual min=-0.006909 max=0.020618 Jy/beam ! Clean residual mean=0.000953 rms=0.002283 Jy/beam ! Combined flux in latest and established models = 5.23001 Jy ! Adding 19 model components to the UV plane model. ! The established model now contains 27 components and 5.23002 Jy ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Added new window around map position (-0.4, 1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0177959 Jy ! Component: 100 - total flux cleaned = 0.0235006 Jy ! Component: 150 - total flux cleaned = 0.0279251 Jy ! Component: 200 - total flux cleaned = 0.0317665 Jy ! Total flux subtracted in 200 components = 0.0317665 Jy ! Clean residual min=-0.006926 max=0.011378 Jy/beam ! Clean residual mean=0.000399 rms=0.001864 Jy/beam ! Combined flux in latest and established models = 5.26178 Jy ! Adding 20 model components to the UV plane model. ! The established model now contains 40 components and 5.26178 Jy ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Added new window around map position (1.4, 2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0092336 Jy ! Component: 100 - total flux cleaned = 0.0114947 Jy ! Component: 150 - total flux cleaned = 0.0135184 Jy ! Component: 200 - total flux cleaned = 0.0150988 Jy ! Total flux subtracted in 200 components = 0.0150988 Jy ! Clean residual min=-0.006891 max=0.006971 Jy/beam ! Clean residual mean=0.000135 rms=0.001693 Jy/beam ! Combined flux in latest and established models = 5.27688 Jy ! Adding 21 model components to the UV plane model. ! The established model now contains 51 components and 5.27688 Jy ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec uvw 0,2 ! Uniform weighting is not currently selected. ! Amplitude error weighting is not currently selected. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! Estimated beam: bmin=0.6713 mas, bmaj=1.008 mas, bpa=1.087 degrees ! Estimated noise=996.005 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00200379 Jy ! Component: 100 - total flux cleaned = -0.00169399 Jy ! Component: 150 - total flux cleaned = -0.000525125 Jy ! Component: 200 - total flux cleaned = 0.000314762 Jy ! Total flux subtracted in 200 components = 0.000314762 Jy ! Clean residual min=-0.008114 max=0.008135 Jy/beam ! Clean residual mean=0.000056 rms=0.002078 Jy/beam ! Combined flux in latest and established models = 5.27719 Jy clean ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000520664 Jy ! Component: 100 - total flux cleaned = 0.000777089 Jy ! Component: 150 - total flux cleaned = 0.0010251 Jy ! Component: 200 - total flux cleaned = 0.00126006 Jy ! Total flux subtracted in 200 components = 0.00126006 Jy ! Clean residual min=-0.007909 max=0.007966 Jy/beam ! Clean residual mean=0.000022 rms=0.002038 Jy/beam ! Combined flux in latest and established models = 5.27845 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 67 components and 5.27845 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.213762Jy sigma=0.001287 ! Fit after self-cal, rms=0.213090Jy sigma=0.001283 wmodel J1751+0939_X_2011_06_28_pus_map.mod ! Writing 67 model components to file: J1751+0939_X_2011_06_28_pus_map.mod wobs J1751+0939_X_2011_06_28_pus_uvs.fits ! Writing UV FITS file: J1751+0939_X_2011_06_28_pus_uvs.fits wwins J1751+0939_X_2011_06_28_pus_map.win ! wwins: Wrote 4 windows to J1751+0939_X_2011_06_28_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000620182 Jy ! Component: 100 - total flux cleaned = 0.000898695 Jy ! Component: 150 - total flux cleaned = 0.000899922 Jy ! Component: 200 - total flux cleaned = 0.000900962 Jy ! Component: 250 - total flux cleaned = 0.00102876 Jy ! Component: 300 - total flux cleaned = 0.00115228 Jy ! Component: 350 - total flux cleaned = 0.00103044 Jy ! Component: 400 - total flux cleaned = 0.00115195 Jy ! Component: 450 - total flux cleaned = 0.00115105 Jy ! Component: 500 - total flux cleaned = 0.000917667 Jy ! Component: 550 - total flux cleaned = 0.000918983 Jy ! Component: 600 - total flux cleaned = 0.000919327 Jy ! Component: 650 - total flux cleaned = 0.000807722 Jy ! Component: 700 - total flux cleaned = 0.000917599 Jy ! Component: 750 - total flux cleaned = 0.000808228 Jy ! Component: 800 - total flux cleaned = 0.000807968 Jy ! Component: 850 - total flux cleaned = 0.000807469 Jy ! Component: 900 - total flux cleaned = 0.000596386 Jy ! Component: 950 - total flux cleaned = 0.00059697 Jy ! Component: 1000 - total flux cleaned = 0.000597103 Jy ! Component: 1050 - total flux cleaned = 0.00049425 Jy ! Component: 1100 - total flux cleaned = 0.000495094 Jy ! Component: 1150 - total flux cleaned = 0.000293701 Jy ! Component: 1200 - total flux cleaned = 0.000394259 Jy ! Component: 1250 - total flux cleaned = 0.000394499 Jy ! Component: 1300 - total flux cleaned = 0.000394375 Jy ! Component: 1350 - total flux cleaned = 0.000393809 Jy ! Component: 1400 - total flux cleaned = 0.000199556 Jy ! Component: 1450 - total flux cleaned = 0.000199524 Jy ! Component: 1500 - total flux cleaned = 0.000199158 Jy ! Component: 1550 - total flux cleaned = 0.000104113 Jy ! Component: 1600 - total flux cleaned = 1.04255e-05 Jy ! Component: 1650 - total flux cleaned = 0.000103392 Jy ! Component: 1700 - total flux cleaned = 1.07304e-05 Jy ! Component: 1750 - total flux cleaned = 1.08223e-05 Jy ! Component: 1800 - total flux cleaned = 0.000101748 Jy ! Component: 1850 - total flux cleaned = 1.11126e-05 Jy ! Component: 1900 - total flux cleaned = 0.000100766 Jy ! Component: 1950 - total flux cleaned = 0.000190142 Jy ! Component: 2000 - total flux cleaned = 0.000100457 Jy ! Total flux subtracted in 2048 components = 1.19826e-05 Jy ! Clean residual min=-0.004952 max=0.005527 Jy/beam ! Clean residual mean=-0.000019 rms=0.001579 Jy/beam ! Combined flux in latest and established models = 5.27847 Jy keep ! Adding 179 model components to the UV plane model. ! The established model now contains 243 components and 5.27847 Jy ![Exited script file: muppet] wmap "J1751+0939_X_2011_06_28_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6713 x 1.008 at 1.087 degrees (North through East) ! Clean map min=-0.0070371 max=4.7415 Jy/beam ! Writing clean map to FITS file: J1751+0939_X_2011_06_28_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 4.74147 invert ! Inverting map ! Your choice of large map pixels excluded 5.47% of the data. ! The x-axis pixel size should ideally be below 0.1484 milli-arcsec ! The y-axis pixel size should ideally be below 0.1779 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00157752 device "J1751+0939_X_2011_06_28_pus\_map.eps/VPS" ! Attempting to open device: 'J1751+0939_X_2011_06_28_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.0998117 0.0998117 0.199623 0.399247 0.798494 1.59699 3.19398 6.38795 12.7759 25.5518 51.1036 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0047 print "peak=",int(1000*peak_flux)/1000 ! peak= 4.741 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.6713 x 1.008 at 1.087 degrees (North through East) ! Clean map min=-0.0070371 max=4.7415 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1751+0939_X_2011_06_28_pus\_map.gif/vgif" ! Attempting to open device: 'J1751+0939_X_2011_06_28_pus_map.gif/vgif' mapplot cln