! Started logfile: J1753+2848_S_2011_06_28_pus_uvs.log on Tue Apr 24 22:44:08 2012 obs J1753+2848_S_2011_06_28_pus_uva.fits ! Reading UV FITS file: J1753+2848_S_2011_06_28_pus_uva.fits ! AN table 1: 283 integrations on 153 of 153 possible baselines. ! AN table 2: 99 integrations on 153 of 153 possible baselines. ! AN table 3: 11 integrations on 153 of 153 possible baselines. ! Apparent sampling: 0.102779 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1753+2848 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23697e+09 8e+06 1 8e+06 ! 02 2 2.26697e+09 8e+06 1 8e+06 ! 03 3 2.35697e+09 8e+06 1 8e+06 ! 04 4 2.37697e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 4589 lines of history. ! ! Reading 24720 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.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J1753+2848_S_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.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 3672 telescope corrections were flagged in sub-array 1. ! A total of 1595 telescope corrections were flagged in sub-array 2. ! A total of 198 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 3672 telescope corrections were flagged in sub-array 1. ! A total of 1595 telescope corrections were flagged in sub-array 2. ! A total of 198 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 3639 telescope corrections were flagged in sub-array 1. ! A total of 1540 telescope corrections were flagged in sub-array 2. ! A total of 198 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 3639 telescope corrections were flagged in sub-array 1. ! A total of 1540 telescope corrections were flagged in sub-array 2. ! A total of 198 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.643294Jy sigma=0.006322 ! Fit after self-cal, rms=0.406454Jy sigma=0.003096 ! 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=1.819 mas, bmaj=2.686 mas, bpa=-0.05389 degrees ! Estimated noise=959.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 = 0.51624 Jy ! Component: 100 - total flux cleaned = 0.634786 Jy ! Total flux subtracted in 100 components = 0.634786 Jy ! Clean residual min=-0.013351 max=0.044174 Jy/beam ! Clean residual mean=0.000235 rms=0.004778 Jy/beam ! Combined flux in latest and established models = 0.634786 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.634786 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.218151Jy sigma=0.001612 ! Fit after self-cal, rms=0.218196Jy sigma=0.001612 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0496939 Jy ! Component: 100 - total flux cleaned = 0.0784271 Jy ! Total flux subtracted in 100 components = 0.0784271 Jy ! Clean residual min=-0.011158 max=0.015091 Jy/beam ! Clean residual mean=0.000014 rms=0.003054 Jy/beam ! Combined flux in latest and established models = 0.713213 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 21 components and 0.713213 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.208745Jy sigma=0.001401 ! Fit after self-cal, rms=0.208696Jy sigma=0.001401 ! 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=2.135 mas, bmaj=3.098 mas, bpa=1.489 degrees ! Estimated noise=550.951 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.0206718 Jy ! Component: 100 - total flux cleaned = 0.0321518 Jy ! Component: 150 - total flux cleaned = 0.0389869 Jy ! Component: 200 - total flux cleaned = 0.0421457 Jy ! Total flux subtracted in 200 components = 0.0421457 Jy ! Clean residual min=-0.010090 max=0.010752 Jy/beam ! Clean residual mean=-0.000366 rms=0.002344 Jy/beam ! Combined flux in latest and established models = 0.755359 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 37 components and 0.755359 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207608Jy sigma=0.001359 ! Fit after self-cal, rms=0.207597Jy sigma=0.001359 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 39 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.03 ! HN 1.04 KK 0.99 KP 0.97 MC 1.00 ! MK 0.92 NL 0.98 NY 1.01 OV 0.90 ! PT 1.00 SC 1.04 SH 1.00* TS 1.00* ! WF 1.06 WZ 1.03 ! ! A total of 63 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* FT 1.00* HH 1.09 ! HN 1.00* KK 1.00* KP 1.00* MC 1.08 ! MK 1.00* NL 1.00* NY 0.97 OV 1.00* ! PT 1.00* SC 1.00* SH 0.97 TS 1.00* ! WF 1.00* WZ 1.00* ! ! 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 39 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.98 FT 1.00* HH 0.97 ! HN 1.02 KK 1.01 KP 1.00 MC 0.99 ! MK 0.91 NL 0.94 NY 1.03 OV 0.95 ! PT 0.97 SC 1.02 SH 1.00* TS 1.00* ! WF 1.01 WZ 0.99 ! ! A total of 63 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* FT 1.00* HH 1.01 ! HN 1.00* KK 1.00* KP 1.00* MC 0.94 ! MK 1.00* NL 1.00* NY 1.06 OV 1.00* ! PT 1.00* SC 1.00* SH 0.99 TS 1.00* ! WF 1.00* WZ 1.00* ! ! 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 72 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.00 FT 1.00* HH 0.98 ! HN 1.04 KK 1.10 KP 1.01 MC 1.05 ! MK 0.93 NL 1.01 NY 1.05 OV 0.97 ! PT 1.01 SC 1.04 SH 1.00* TS 1.00* ! WF 1.02 WZ 1.03 ! ! A total of 84 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* FT 1.00* HH 1.07 ! HN 1.00* KK 1.00* KP 1.00* MC 1.05 ! MK 1.00* NL 1.00* NY 1.09 OV 1.00* ! PT 1.00* SC 1.00* SH 1.09 TS 1.01 ! WF 1.00* WZ 1.00* ! ! 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 72 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.99 FT 1.00* HH 0.98 ! HN 1.02 KK 1.06 KP 1.00 MC 0.93 ! MK 0.96 NL 0.99 NY 1.09 OV 0.96 ! PT 0.99 SC 1.04 SH 1.00* TS 1.00* ! WF 1.02 WZ 1.04 ! ! A total of 84 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* FT 1.00* HH 1.07 ! HN 1.00* KK 1.00* KP 1.00* MC 1.05 ! MK 1.00* NL 1.00* NY 1.07 OV 1.00* ! PT 1.00* SC 1.00* SH 0.94 TS 1.00 ! WF 1.00* WZ 1.00* ! ! 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.207597Jy sigma=0.001359 ! Fit after self-cal, rms=0.209343Jy sigma=0.001314 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.154 mas, bmaj=3.171 mas, bpa=2.519 degrees ! Estimated noise=553.272 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.00377956 Jy ! Component: 100 - total flux cleaned = -0.00529465 Jy ! Component: 150 - total flux cleaned = -0.00571909 Jy ! Component: 200 - total flux cleaned = -0.00587464 Jy ! Total flux subtracted in 200 components = -0.00587464 Jy ! Clean residual min=-0.007310 max=0.007467 Jy/beam ! Clean residual mean=-0.000016 rms=0.001776 Jy/beam ! Combined flux in latest and established models = 0.749484 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 42 components and 0.749484 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.209223Jy sigma=0.001314 ! Fit after self-cal, rms=0.209325Jy sigma=0.001314 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.209325Jy sigma=0.001314 ! Fit after self-cal, rms=0.207470Jy sigma=0.001284 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.157 mas, bmaj=3.167 mas, bpa=2.187 degrees ! Estimated noise=554.166 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000250573 Jy ! Component: 100 - total flux cleaned = -0.000321792 Jy ! Component: 150 - total flux cleaned = -0.000322136 Jy ! Component: 200 - total flux cleaned = -0.000458889 Jy ! Total flux subtracted in 200 components = -0.000458889 Jy ! Clean residual min=-0.003500 max=0.006464 Jy/beam ! Clean residual mean=0.000257 rms=0.000955 Jy/beam ! Combined flux in latest and established models = 0.749025 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 44 components and 0.749025 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207400Jy sigma=0.001284 ! Fit after self-cal, rms=0.207436Jy sigma=0.001284 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 (4.5, -2.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00536481 Jy ! Component: 100 - total flux cleaned = 0.00695983 Jy ! Component: 150 - total flux cleaned = 0.00740558 Jy ! Component: 200 - total flux cleaned = 0.00756026 Jy ! Total flux subtracted in 200 components = 0.00756026 Jy ! Clean residual min=-0.003209 max=0.004368 Jy/beam ! Clean residual mean=0.000161 rms=0.000870 Jy/beam ! Combined flux in latest and established models = 0.756585 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 58 components and 0.756585 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207259Jy sigma=0.001282 ! Fit after self-cal, rms=0.207248Jy sigma=0.001282 ! Inverting map ! Added new window around map position (17, 14). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00420774 Jy ! Component: 100 - total flux cleaned = 0.00676393 Jy ! Component: 150 - total flux cleaned = 0.00778239 Jy ! Component: 200 - total flux cleaned = 0.00818244 Jy ! Total flux subtracted in 200 components = 0.00818244 Jy ! Clean residual min=-0.003038 max=0.003068 Jy/beam ! Clean residual mean=0.000070 rms=0.000783 Jy/beam ! Combined flux in latest and established models = 0.764768 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 73 components and 0.764768 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207140Jy sigma=0.001280 ! Fit after self-cal, rms=0.207133Jy sigma=0.001280 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207133Jy sigma=0.001280 ! Fit after self-cal, rms=0.206410Jy sigma=0.001279 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.156 mas, bmaj=3.162 mas, bpa=2.157 degrees ! Estimated noise=554.495 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000927253 Jy ! Component: 100 - total flux cleaned = 0.00113532 Jy ! Component: 150 - total flux cleaned = 0.00131897 Jy ! Component: 200 - total flux cleaned = 0.00137396 Jy ! Total flux subtracted in 200 components = 0.00137396 Jy ! Clean residual min=-0.002853 max=0.003042 Jy/beam ! Clean residual mean=0.000079 rms=0.000719 Jy/beam ! Combined flux in latest and established models = 0.766142 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 82 components and 0.766142 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206348Jy sigma=0.001279 ! Fit after self-cal, rms=0.206347Jy sigma=0.001279 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 (27.5, 17). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00181285 Jy ! Component: 100 - total flux cleaned = 0.00237845 Jy ! Component: 150 - total flux cleaned = 0.00248881 Jy ! Component: 200 - total flux cleaned = 0.00259003 Jy ! Total flux subtracted in 200 components = 0.00259003 Jy ! Clean residual min=-0.002628 max=0.002705 Jy/beam ! Clean residual mean=0.000054 rms=0.000697 Jy/beam ! Combined flux in latest and established models = 0.768732 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 95 components and 0.768732 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206285Jy sigma=0.001278 ! Fit after self-cal, rms=0.206280Jy sigma=0.001278 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206280Jy sigma=0.001278 ! Fit after self-cal, rms=0.206254Jy sigma=0.001278 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.528 mas, bmaj=8.743 mas, bpa=-5.606 degrees ! Estimated noise=785.389 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00204333 Jy ! Component: 100 - total flux cleaned = 0.00310773 Jy ! Component: 150 - total flux cleaned = 0.00347084 Jy ! Component: 200 - total flux cleaned = 0.00353143 Jy ! Total flux subtracted in 200 components = 0.00353143 Jy ! Clean residual min=-0.002542 max=0.002825 Jy/beam ! Clean residual mean=0.000086 rms=0.000925 Jy/beam ! Combined flux in latest and established models = 0.772263 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 110 components and 0.772263 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206253Jy sigma=0.001278 ! Fit after self-cal, rms=0.206253Jy sigma=0.001278 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.0004124 Jy ! Component: 100 - total flux cleaned = 0.000466199 Jy ! Component: 150 - total flux cleaned = 0.000491309 Jy ! Component: 200 - total flux cleaned = 0.000557506 Jy ! Total flux subtracted in 200 components = 0.000557506 Jy ! Clean residual min=-0.002503 max=0.002929 Jy/beam ! Clean residual mean=0.000068 rms=0.000925 Jy/beam ! Combined flux in latest and established models = 0.77282 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 112 components and 0.77282 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206257Jy sigma=0.001278 ! Fit after self-cal, rms=0.206259Jy sigma=0.001278 ! 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=2.156 mas, bmaj=3.161 mas, bpa=2.143 degrees ! Estimated noise=554.46 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.156 x 3.161 at 2.143 degrees (North through East) ! Clean map min=-0.0025773 max=0.6719 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.206259Jy sigma=0.001278 ! Fit after self-cal, rms=0.199796Jy sigma=0.001174 ! Inverting map and beam ! Estimated beam: bmin=2.158 mas, bmaj=3.16 mas, bpa=2.15 degrees ! Estimated noise=554.516 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000302942 Jy ! Component: 100 - total flux cleaned = -0.000361793 Jy ! Component: 150 - total flux cleaned = -0.000464353 Jy ! Component: 200 - total flux cleaned = -0.000321816 Jy ! Total flux subtracted in 200 components = -0.000321816 Jy ! Clean residual min=-0.002539 max=0.002425 Jy/beam ! Clean residual mean=0.000043 rms=0.000659 Jy/beam ! Combined flux in latest and established models = 0.772499 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 127 components and 0.772498 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199752Jy sigma=0.001174 ! Fit after self-cal, rms=0.199756Jy sigma=0.001174 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 9.17355e-05 Jy ! Component: 100 - total flux cleaned = 0.000221931 Jy ! Component: 150 - total flux cleaned = 0.000262484 Jy ! Component: 200 - total flux cleaned = 0.000303469 Jy ! Total flux subtracted in 200 components = 0.000303469 Jy ! Clean residual min=-0.002465 max=0.002420 Jy/beam ! Clean residual mean=0.000039 rms=0.000655 Jy/beam ! Combined flux in latest and established models = 0.772802 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 136 components and 0.772802 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199730Jy sigma=0.001173 ! Fit after self-cal, rms=0.199731Jy sigma=0.001173 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199731Jy sigma=0.001173 ! Fit after self-cal, rms=0.199728Jy sigma=0.001173 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=1.83 mas, bmaj=2.725 mas, bpa=-0.106 degrees ! Estimated noise=948 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.515125 Jy ! Component: 100 - total flux cleaned = 0.632297 Jy ! Component: 150 - total flux cleaned = 0.68038 Jy ! Component: 200 - total flux cleaned = 0.708059 Jy ! Total flux subtracted in 200 components = 0.708059 Jy ! Clean residual min=-0.006262 max=0.014168 Jy/beam ! Clean residual mean=0.000182 rms=0.001804 Jy/beam ! Combined flux in latest and established models = 0.708059 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 17 model components to the UV plane model. ! The established model now contains 17 components and 0.708059 Jy ! Inverting map and beam ! Estimated beam: bmin=2.157 mas, bmaj=3.158 mas, bpa=2.118 degrees ! Estimated noise=554.738 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.0209053 Jy ! Component: 100 - total flux cleaned = 0.0326775 Jy ! Component: 150 - total flux cleaned = 0.0401963 Jy ! Component: 200 - total flux cleaned = 0.0453539 Jy ! Total flux subtracted in 200 components = 0.0453539 Jy ! Clean residual min=-0.003200 max=0.003831 Jy/beam ! Clean residual mean=0.000260 rms=0.000811 Jy/beam ! Combined flux in latest and established models = 0.753413 Jy ! Adding 28 model components to the UV plane model. ! The established model now contains 37 components and 0.753413 Jy ! Inverting map ! Added new window around map position (6.5, 6.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00431074 Jy ! Component: 100 - total flux cleaned = 0.00772057 Jy ! Component: 150 - total flux cleaned = 0.0103222 Jy ! Component: 200 - total flux cleaned = 0.0120745 Jy ! Total flux subtracted in 200 components = 0.0120745 Jy ! Clean residual min=-0.002888 max=0.003164 Jy/beam ! Clean residual mean=0.000121 rms=0.000699 Jy/beam ! Combined flux in latest and established models = 0.765487 Jy ! Adding 37 model components to the UV plane model. ! The established model now contains 56 components and 0.765487 Jy ! Inverting map ! Added new window around map position (25.5, -45). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00281286 Jy ! Component: 100 - total flux cleaned = 0.00389182 Jy ! Component: 150 - total flux cleaned = 0.00457352 Jy ! Component: 200 - total flux cleaned = 0.0052064 Jy ! Total flux subtracted in 200 components = 0.0052064 Jy ! Clean residual min=-0.002770 max=0.002750 Jy/beam ! Clean residual mean=0.000063 rms=0.000661 Jy/beam ! Combined flux in latest and established models = 0.770694 Jy ! Adding 37 model components to the UV plane model. ! The established model now contains 74 components and 0.770694 Jy ! Inverting map 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 ! Estimated beam: bmin=1.989 mas, bmaj=2.686 mas, bpa=7.39 degrees ! Estimated noise=779.128 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000611689 Jy ! Component: 100 - total flux cleaned = 0.0014602 Jy ! Component: 150 - total flux cleaned = 0.00194347 Jy ! Component: 200 - total flux cleaned = 0.0022135 Jy ! Total flux subtracted in 200 components = 0.0022135 Jy ! Clean residual min=-0.003657 max=0.003384 Jy/beam ! Clean residual mean=0.000029 rms=0.000869 Jy/beam ! Combined flux in latest and established models = 0.772907 Jy clean ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000337711 Jy ! Component: 100 - total flux cleaned = 0.000499263 Jy ! Component: 150 - total flux cleaned = 0.000653671 Jy ! Component: 200 - total flux cleaned = 0.000879594 Jy ! Total flux subtracted in 200 components = 0.000879594 Jy ! Clean residual min=-0.003493 max=0.003263 Jy/beam ! Clean residual mean=0.000018 rms=0.000847 Jy/beam ! Combined flux in latest and established models = 0.773787 Jy selfcal ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 107 components and 0.773787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199650Jy sigma=0.001173 ! Fit after self-cal, rms=0.199655Jy sigma=0.001173 wmodel J1753+2848_S_2011_06_28_pus_map.mod ! Writing 107 model components to file: J1753+2848_S_2011_06_28_pus_map.mod wobs J1753+2848_S_2011_06_28_pus_uvs.fits ! Writing UV FITS file: J1753+2848_S_2011_06_28_pus_uvs.fits wwins J1753+2848_S_2011_06_28_pus_map.win ! wwins: Wrote 6 windows to J1753+2848_S_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 ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000512937 Jy ! Component: 100 - total flux cleaned = -0.000861601 Jy ! Component: 150 - total flux cleaned = -0.00091596 Jy ! Component: 200 - total flux cleaned = -0.000968611 Jy ! Component: 250 - total flux cleaned = -0.000915795 Jy ! Component: 300 - total flux cleaned = -0.000814242 Jy ! Component: 350 - total flux cleaned = -0.000763708 Jy ! Component: 400 - total flux cleaned = -0.000616761 Jy ! Component: 450 - total flux cleaned = -0.000519936 Jy ! Component: 500 - total flux cleaned = -0.000377856 Jy ! Component: 550 - total flux cleaned = -0.000331069 Jy ! Component: 600 - total flux cleaned = -0.000193644 Jy ! Component: 650 - total flux cleaned = -5.85728e-05 Jy ! Component: 700 - total flux cleaned = 3.03144e-05 Jy ! Component: 750 - total flux cleaned = 0.000205528 Jy ! Component: 800 - total flux cleaned = 0.000248029 Jy ! Component: 850 - total flux cleaned = 0.000376585 Jy ! Component: 900 - total flux cleaned = 0.000419283 Jy ! Component: 950 - total flux cleaned = 0.000586786 Jy ! Component: 1000 - total flux cleaned = 0.000711014 Jy ! Component: 1050 - total flux cleaned = 0.000874119 Jy ! Component: 1100 - total flux cleaned = 0.000873759 Jy ! Component: 1150 - total flux cleaned = 0.000994353 Jy ! Component: 1200 - total flux cleaned = 0.00119245 Jy ! Component: 1250 - total flux cleaned = 0.00115295 Jy ! Component: 1300 - total flux cleaned = 0.00126968 Jy ! Component: 1350 - total flux cleaned = 0.00127002 Jy ! Component: 1400 - total flux cleaned = 0.00146117 Jy ! Component: 1450 - total flux cleaned = 0.0014233 Jy ! Component: 1500 - total flux cleaned = 0.00161087 Jy ! Component: 1550 - total flux cleaned = 0.00161064 Jy ! Component: 1600 - total flux cleaned = 0.00161073 Jy ! Component: 1650 - total flux cleaned = 0.00168408 Jy ! Component: 1700 - total flux cleaned = 0.0017202 Jy ! Component: 1750 - total flux cleaned = 0.00172007 Jy ! Component: 1800 - total flux cleaned = 0.00175591 Jy ! Component: 1850 - total flux cleaned = 0.00186226 Jy ! Component: 1900 - total flux cleaned = 0.00175614 Jy ! Component: 1950 - total flux cleaned = 0.00189634 Jy ! Component: 2000 - total flux cleaned = 0.00189625 Jy ! Total flux subtracted in 2048 components = 0.00196541 Jy ! Clean residual min=-0.002007 max=0.002241 Jy/beam ! Clean residual mean=0.000006 rms=0.000625 Jy/beam ! Combined flux in latest and established models = 0.775752 Jy keep ! Adding 179 model components to the UV plane model. ! The established model now contains 286 components and 0.775752 Jy ![Exited script file: muppet] wmap "J1753+2848_S_2011_06_28_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.989 x 2.686 at 7.39 degrees (North through East) ! Clean map min=-0.0025318 max=0.66266 Jy/beam ! Writing clean map to FITS file: J1753+2848_S_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 0.662665 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000624192 device "J1753+2848_S_2011_06_28_pus\_map.eps/VPS" ! Attempting to open device: 'J1753+2848_S_2011_06_28_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.282583 0.282583 0.565165 1.13033 2.26066 4.52132 9.04265 18.0853 36.1706 72.3412 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0018 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.662 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.989 x 2.686 at 7.39 degrees (North through East) ! Clean map min=-0.0025318 max=0.66266 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1753+2848_S_2011_06_28_pus\_map.gif/vgif" ! Attempting to open device: 'J1753+2848_S_2011_06_28_pus_map.gif/vgif' mapplot cln