! Started logfile: J1517-2422_S_2000_03_13_pus_uvs.log on Tue Apr 2 14:48:30 2013 obs J1517-2422_S_2000_03_13_pus_uva.fits ! Reading UV FITS file: J1517-2422_S_2000_03_13_pus_uva.fits ! AN table 1: 78 integrations on 190 of 190 possible baselines. ! AN table 2: 27 integrations on 190 of 190 possible baselines. ! AN table 3: 91 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.0950054 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1517-2422 ! ! 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.22497e+09 8e+06 1 8e+06 ! 02 2 2.24497e+09 8e+06 1 8e+06 ! 03 3 2.33497e+09 8e+06 1 8e+06 ! 04 4 2.36497e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 2247 lines of history. ! ! Reading 14152 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 J1517-2422_S_2000_03_13_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 1168 telescope corrections were flagged in sub-array 1. ! A total of 424 telescope corrections were flagged in sub-array 2. ! A total of 1315 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 1168 telescope corrections were flagged in sub-array 1. ! A total of 424 telescope corrections were flagged in sub-array 2. ! A total of 1315 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 1168 telescope corrections were flagged in sub-array 1. ! A total of 424 telescope corrections were flagged in sub-array 2. ! A total of 1316 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 1168 telescope corrections were flagged in sub-array 1. ! A total of 424 telescope corrections were flagged in sub-array 2. ! A total of 1315 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.449501Jy sigma=0.011080 ! Fit after self-cal, rms=1.047306Jy sigma=0.007975 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.265 mas, bmaj=7.452 mas, bpa=-2.01 degrees ! Estimated noise=1263.48 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.22396 Jy ! Component: 100 - total flux cleaned = 1.63884 Jy ! Total flux subtracted in 100 components = 1.63884 Jy ! Clean residual min=-0.053720 max=0.174471 Jy/beam ! Clean residual mean=0.001138 rms=0.018691 Jy/beam ! Combined flux in latest and established models = 1.63884 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.63884 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.613714Jy sigma=0.004188 ! Fit after self-cal, rms=0.612135Jy sigma=0.004178 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.194915 Jy ! Component: 100 - total flux cleaned = 0.307297 Jy ! Total flux subtracted in 100 components = 0.307297 Jy ! Clean residual min=-0.037538 max=0.085879 Jy/beam ! Clean residual mean=0.000707 rms=0.012029 Jy/beam ! Combined flux in latest and established models = 1.94613 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 37 components and 1.94613 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.483951Jy sigma=0.002946 ! Fit after self-cal, rms=0.483204Jy sigma=0.002940 ! Inverting map ! Added new window around map position (3.5, -8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0888309 Jy ! Component: 100 - total flux cleaned = 0.14958 Jy ! Total flux subtracted in 100 components = 0.14958 Jy ! Clean residual min=-0.032681 max=0.076531 Jy/beam ! Clean residual mean=0.000475 rms=0.009950 Jy/beam ! Combined flux in latest and established models = 2.09571 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 63 components and 2.09571 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.443112Jy sigma=0.002546 ! Fit after self-cal, rms=0.441750Jy sigma=0.002535 ! Inverting map ! Added new window around map position (8, -17.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0714138 Jy ! Component: 100 - total flux cleaned = 0.118003 Jy ! Total flux subtracted in 100 components = 0.118003 Jy ! Clean residual min=-0.029948 max=0.049245 Jy/beam ! Clean residual mean=0.000264 rms=0.008194 Jy/beam ! Combined flux in latest and established models = 2.21372 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 83 components and 2.21372 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.416349Jy sigma=0.002273 ! Fit after self-cal, rms=0.414965Jy sigma=0.002260 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.293 mas, bmaj=10.25 mas, bpa=-3.738 degrees ! Estimated noise=776.689 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.0875368 Jy ! Component: 100 - total flux cleaned = 0.142018 Jy ! Component: 150 - total flux cleaned = 0.178576 Jy ! Component: 200 - total flux cleaned = 0.203928 Jy ! Total flux subtracted in 200 components = 0.203928 Jy ! Clean residual min=-0.055262 max=0.044316 Jy/beam ! Clean residual mean=-0.000910 rms=0.014318 Jy/beam ! Combined flux in latest and established models = 2.41764 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 131 components and 2.41764 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.390613Jy sigma=0.002050 ! Fit after self-cal, rms=0.388778Jy sigma=0.002028 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.05 GC 1.00* KP 1.04 ! LA 1.02 NL 0.75 OV 1.05 PT 1.01 ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 0.95 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 1.06 GC 1.00* KP 0.98 ! LA 1.04 NL 1.00* OV 1.05 PT 1.00* ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.05 NY 1.00* ! ON 1.00* TS 1.04 WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 1.02 GC 1.00* KP 0.97 ! LA 0.99 NL 0.99 OV 1.03 PT 0.99 ! SC 1.01 WF 1.08 HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 0.98 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.05 ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 1.03 GC 1.00* KP 0.99 ! LA 1.01 NL 0.77 OV 1.01 PT 1.01 ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 0.96 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.09 FD 1.03 GC 1.00* KP 0.98 ! LA 1.06 NL 1.00* OV 1.05 PT 1.00* ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.06 NY 1.00* ! ON 1.00* TS 1.04 WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.01 FD 0.99 GC 1.00* KP 0.95 ! LA 0.99 NL 1.01 OV 1.01 PT 1.00 ! SC 1.01 WF 1.00 HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.00 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.05 ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.00 GC 1.00* KP 1.03 ! LA 1.02 NL 0.73 OV 1.04 PT 1.00 ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 0.96 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 0.99 GC 1.00* KP 1.17 ! LA 1.07 NL 1.00* OV 1.07 PT 1.00* ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.07 NY 1.00* ! ON 1.00* TS 1.06 WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 0.97 GC 1.00* KP 1.05 ! LA 1.00 NL 0.98 OV 1.04 PT 0.97 ! SC 0.99 WF 1.79 HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.01 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.07 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.04 GC 1.00* KP 1.02 ! LA 1.01 NL 0.73 OV 1.01 PT 1.04 ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 0.96 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 1.03 GC 1.00* KP 1.05 ! LA 1.04 NL 1.00* OV 1.05 PT 1.00* ! SC 1.00* WF 1.00* HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.08 NY 1.00* ! ON 1.00* TS 1.05 WZ 1.00* HN 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.01 GC 1.00* KP 1.00 ! LA 1.00 NL 0.98 OV 1.02 PT 0.98 ! SC 1.01 WF 1.03 HH 1.00* KK 1.00* ! MA 1.00* MC 1.00* MK 1.02 NY 1.00* ! ON 1.00* TS 1.00* WZ 1.00* HN 1.04 ! ! ! Fit before self-cal, rms=0.388778Jy sigma=0.002028 ! Fit after self-cal, rms=0.336310Jy sigma=0.001518 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.346 mas, bmaj=10.28 mas, bpa=-3.687 degrees ! Estimated noise=782.292 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 (12.5, -20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0287309 Jy ! Component: 100 - total flux cleaned = 0.0306351 Jy ! Component: 150 - total flux cleaned = 0.0314896 Jy ! Component: 200 - total flux cleaned = 0.032257 Jy ! Total flux subtracted in 200 components = 0.032257 Jy ! Clean residual min=-0.014031 max=0.017225 Jy/beam ! Clean residual mean=0.000269 rms=0.004095 Jy/beam ! Combined flux in latest and established models = 2.4499 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 168 components and 2.4499 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.333823Jy sigma=0.001481 ! Fit after self-cal, rms=0.331963Jy sigma=0.001478 ! 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.331963Jy sigma=0.001478 ! Fit after self-cal, rms=0.353523Jy sigma=0.001418 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.349 mas, bmaj=10.28 mas, bpa=-3.685 degrees ! Estimated noise=781.64 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00736726 Jy ! Component: 100 - total flux cleaned = 0.00878322 Jy ! Component: 150 - total flux cleaned = 0.0097063 Jy ! Component: 200 - total flux cleaned = 0.0105339 Jy ! Total flux subtracted in 200 components = 0.0105339 Jy ! Clean residual min=-0.012378 max=0.014896 Jy/beam ! Clean residual mean=0.000464 rms=0.003566 Jy/beam ! Combined flux in latest and established models = 2.46043 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 185 components and 2.46043 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.352780Jy sigma=0.001408 ! Fit after self-cal, rms=0.352822Jy sigma=0.001407 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.00143161 Jy ! Component: 100 - total flux cleaned = 0.00246764 Jy ! Component: 150 - total flux cleaned = 0.00339228 Jy ! Component: 200 - total flux cleaned = 0.00402682 Jy ! Total flux subtracted in 200 components = 0.00402682 Jy ! Clean residual min=-0.011763 max=0.013934 Jy/beam ! Clean residual mean=0.000370 rms=0.003433 Jy/beam ! Combined flux in latest and established models = 2.46446 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 197 components and 2.46446 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.352465Jy sigma=0.001403 ! Fit after self-cal, rms=0.352399Jy sigma=0.001403 ! 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.352399Jy sigma=0.001403 ! Fit after self-cal, rms=0.323999Jy sigma=0.001293 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.356 mas, bmaj=10.31 mas, bpa=-3.718 degrees ! Estimated noise=781.988 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00356415 Jy ! Component: 100 - total flux cleaned = 0.00490205 Jy ! Component: 150 - total flux cleaned = 0.00546962 Jy ! Component: 200 - total flux cleaned = 0.00581127 Jy ! Total flux subtracted in 200 components = 0.00581127 Jy ! Clean residual min=-0.008203 max=0.011423 Jy/beam ! Clean residual mean=0.000287 rms=0.002814 Jy/beam ! Combined flux in latest and established models = 2.47027 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 217 components and 2.47027 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323727Jy sigma=0.001289 ! Fit after self-cal, rms=0.323656Jy sigma=0.001288 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.000680081 Jy ! Component: 100 - total flux cleaned = 0.00118698 Jy ! Component: 150 - total flux cleaned = 0.00149627 Jy ! Component: 200 - total flux cleaned = 0.00163971 Jy ! Total flux subtracted in 200 components = 0.00163971 Jy ! Clean residual min=-0.008082 max=0.011156 Jy/beam ! Clean residual mean=0.000248 rms=0.002758 Jy/beam ! Combined flux in latest and established models = 2.47191 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 223 components and 2.47191 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323477Jy sigma=0.001286 ! Fit after self-cal, rms=0.323417Jy sigma=0.001285 ! 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.323417Jy sigma=0.001285 ! Fit after self-cal, rms=0.323798Jy sigma=0.001285 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.495 mas, bmaj=18.54 mas, bpa=8.197 degrees ! Estimated noise=911.946 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00206556 Jy ! Component: 100 - total flux cleaned = 0.00289905 Jy ! Component: 150 - total flux cleaned = 0.00341419 Jy ! Component: 200 - total flux cleaned = 0.00391718 Jy ! Total flux subtracted in 200 components = 0.00391718 Jy ! Clean residual min=-0.008931 max=0.015625 Jy/beam ! Clean residual mean=0.000458 rms=0.004030 Jy/beam ! Combined flux in latest and established models = 2.47583 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 240 components and 2.47583 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323694Jy sigma=0.001284 ! Fit after self-cal, rms=0.323611Jy 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) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 2.70423e-05 Jy ! Component: 100 - total flux cleaned = 0.000555857 Jy ! Component: 150 - total flux cleaned = 0.000886567 Jy ! Component: 200 - total flux cleaned = 0.00115281 Jy ! Total flux subtracted in 200 components = 0.00115281 Jy ! Clean residual min=-0.008908 max=0.015353 Jy/beam ! Clean residual mean=0.000410 rms=0.003987 Jy/beam ! Combined flux in latest and established models = 2.47698 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 241 components and 2.47698 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323561Jy sigma=0.001283 ! Fit after self-cal, rms=0.323510Jy sigma=0.001283 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=3.361 mas, bmaj=10.32 mas, bpa=-3.716 degrees ! Estimated noise=782.151 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.361 x 10.32 at -3.716 degrees (North through East) ! Clean map min=-0.008469 max=1.7095 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.323510Jy sigma=0.001283 ! Fit after self-cal, rms=0.418792Jy sigma=0.001099 ! Inverting map and beam ! Estimated beam: bmin=3.364 mas, bmaj=10.32 mas, bpa=-3.711 degrees ! Estimated noise=781.486 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00158516 Jy ! Component: 100 - total flux cleaned = -0.00189564 Jy ! Component: 150 - total flux cleaned = -0.00189096 Jy ! Component: 200 - total flux cleaned = -0.0018841 Jy ! Total flux subtracted in 200 components = -0.0018841 Jy ! Clean residual min=-0.007897 max=0.010795 Jy/beam ! Clean residual mean=0.000241 rms=0.002631 Jy/beam ! Combined flux in latest and established models = 2.4751 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 248 components and 2.4751 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.418660Jy sigma=0.001097 ! Fit after self-cal, rms=0.418636Jy sigma=0.001096 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -3.88622e-05 Jy ! Component: 100 - total flux cleaned = -3.26965e-05 Jy ! Component: 150 - total flux cleaned = -2.93966e-05 Jy ! Component: 200 - total flux cleaned = 8.6788e-05 Jy ! Total flux subtracted in 200 components = 8.6788e-05 Jy ! Clean residual min=-0.007613 max=0.010698 Jy/beam ! Clean residual mean=0.000239 rms=0.002604 Jy/beam ! Combined flux in latest and established models = 2.47519 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 253 components and 2.47519 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.418546Jy sigma=0.001095 ! Fit after self-cal, rms=0.418516Jy sigma=0.001095 ! 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.418516Jy sigma=0.001095 ! Fit after self-cal, rms=0.419379Jy sigma=0.001094 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=2.293 mas, bmaj=7.475 mas, bpa=-2.071 degrees ! Estimated noise=1281.72 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.24043 Jy ! Component: 100 - total flux cleaned = 1.65935 Jy ! Component: 150 - total flux cleaned = 1.8548 Jy ! Component: 200 - total flux cleaned = 1.97096 Jy ! Total flux subtracted in 200 components = 1.97096 Jy ! Clean residual min=-0.019792 max=0.064378 Jy/beam ! Clean residual mean=0.000796 rms=0.009135 Jy/beam ! Combined flux in latest and established models = 1.97096 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 34 model components to the UV plane model. ! The established model now contains 34 components and 1.97096 Jy ! Inverting map and beam ! Estimated beam: bmin=3.369 mas, bmaj=10.33 mas, bpa=-3.714 degrees ! Estimated noise=781.485 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.165391 Jy ! Component: 100 - total flux cleaned = 0.266133 Jy ! Component: 150 - total flux cleaned = 0.332407 Jy ! Component: 200 - total flux cleaned = 0.37756 Jy ! Total flux subtracted in 200 components = 0.37756 Jy ! Clean residual min=-0.013254 max=0.024879 Jy/beam ! Clean residual mean=0.003131 rms=0.007316 Jy/beam ! Combined flux in latest and established models = 2.34852 Jy ! Adding 115 model components to the UV plane model. ! The established model now contains 138 components and 2.34852 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.5 mas, bmaj=18.55 mas, bpa=8.258 degrees ! Estimated noise=909.631 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0467176 Jy ! Component: 100 - total flux cleaned = 0.0742187 Jy ! Component: 150 - total flux cleaned = 0.091503 Jy ! Component: 200 - total flux cleaned = 0.102806 Jy ! Total flux subtracted in 200 components = 0.102806 Jy ! Clean residual min=-0.010828 max=0.016612 Jy/beam ! Clean residual mean=0.001339 rms=0.004703 Jy/beam ! Combined flux in latest and established models = 2.45133 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.00751141 Jy ! Component: 100 - total flux cleaned = 0.0125529 Jy ! Component: 150 - total flux cleaned = 0.0159512 Jy ! Component: 200 - total flux cleaned = 0.0182516 Jy ! Total flux subtracted in 200 components = 0.0182516 Jy ! Clean residual min=-0.009453 max=0.015744 Jy/beam ! Clean residual mean=0.000611 rms=0.004012 Jy/beam ! Combined flux in latest and established models = 2.46958 Jy ! Adding 142 model components to the UV plane model. ! The established model now contains 264 components and 2.46958 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.369 mas, bmaj=10.33 mas, bpa=-3.714 degrees ! Estimated noise=781.485 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000990359 Jy ! Component: 100 - total flux cleaned = 0.000189592 Jy ! Component: 150 - total flux cleaned = 0.000591501 Jy ! Component: 200 - total flux cleaned = 0.000941514 Jy ! Total flux subtracted in 200 components = 0.000941514 Jy ! Clean residual min=-0.007646 max=0.011057 Jy/beam ! Clean residual mean=0.000316 rms=0.002663 Jy/beam ! Combined flux in latest and established models = 2.47052 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 279 components and 2.47052 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.419304Jy sigma=0.001097 ! Fit after self-cal, rms=0.419031Jy sigma=0.001093 wmodel J1517-2422_S_2000_03_13_pus_map.mod ! Writing 279 model components to file: J1517-2422_S_2000_03_13_pus_map.mod wobs J1517-2422_S_2000_03_13_pus_uvs.fits ! Writing UV FITS file: J1517-2422_S_2000_03_13_pus_uvs.fits wwins J1517-2422_S_2000_03_13_pus_map.win ! wwins: Wrote 4 windows to J1517-2422_S_2000_03_13_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.0043842 Jy ! Component: 100 - total flux cleaned = 0.00582085 Jy ! Component: 150 - total flux cleaned = 0.00686022 Jy ! Component: 200 - total flux cleaned = 0.00746829 Jy ! Component: 250 - total flux cleaned = 0.00792492 Jy ! Component: 300 - total flux cleaned = 0.00856836 Jy ! Component: 350 - total flux cleaned = 0.00897343 Jy ! Component: 400 - total flux cleaned = 0.00936109 Jy ! Component: 450 - total flux cleaned = 0.00982475 Jy ! Component: 500 - total flux cleaned = 0.0103626 Jy ! Component: 550 - total flux cleaned = 0.0105354 Jy ! Component: 600 - total flux cleaned = 0.0106194 Jy ! Component: 650 - total flux cleaned = 0.0107017 Jy ! Component: 700 - total flux cleaned = 0.0108629 Jy ! Component: 750 - total flux cleaned = 0.0110197 Jy ! Component: 800 - total flux cleaned = 0.0113255 Jy ! Component: 850 - total flux cleaned = 0.0115487 Jy ! Component: 900 - total flux cleaned = 0.0116228 Jy ! Component: 950 - total flux cleaned = 0.0116239 Jy ! Component: 1000 - total flux cleaned = 0.0116937 Jy ! Component: 1050 - total flux cleaned = 0.0114861 Jy ! Component: 1100 - total flux cleaned = 0.0113519 Jy ! Component: 1150 - total flux cleaned = 0.0113517 Jy ! Component: 1200 - total flux cleaned = 0.0112861 Jy ! Component: 1250 - total flux cleaned = 0.0114147 Jy ! Component: 1300 - total flux cleaned = 0.0113508 Jy ! Component: 1350 - total flux cleaned = 0.0115388 Jy ! Component: 1400 - total flux cleaned = 0.0115377 Jy ! Component: 1450 - total flux cleaned = 0.0115372 Jy ! Component: 1500 - total flux cleaned = 0.011597 Jy ! Component: 1550 - total flux cleaned = 0.0116561 Jy ! Component: 1600 - total flux cleaned = 0.0116564 Jy ! Component: 1650 - total flux cleaned = 0.0116561 Jy ! Component: 1700 - total flux cleaned = 0.0116001 Jy ! Component: 1750 - total flux cleaned = 0.0116573 Jy ! Component: 1800 - total flux cleaned = 0.0117133 Jy ! Component: 1850 - total flux cleaned = 0.011768 Jy ! Component: 1900 - total flux cleaned = 0.0117682 Jy ! Component: 1950 - total flux cleaned = 0.0117686 Jy ! Component: 2000 - total flux cleaned = 0.011769 Jy ! Total flux subtracted in 2048 components = 0.0117689 Jy ! Clean residual min=-0.003457 max=0.003193 Jy/beam ! Clean residual mean=0.000062 rms=0.001137 Jy/beam ! Combined flux in latest and established models = 2.48229 Jy keep ! Adding 158 model components to the UV plane model. ! The established model now contains 437 components and 2.48229 Jy ![Exited script file: muppet] wmap "J1517-2422_S_2000_03_13_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.369 x 10.33 at -3.714 degrees (North through East) ! Clean map min=-0.005732 max=1.7113 Jy/beam ! Writing clean map to FITS file: J1517-2422_S_2000_03_13_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.71125 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.001116 device "J1517-2422_S_2000_03_13_pus\_map.eps/VPS" ! Attempting to open device: 'J1517-2422_S_2000_03_13_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.195646 0.195646 0.391293 0.782585 1.56517 3.13034 6.26068 12.5214 25.0427 50.0854 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0033 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.711 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.369 x 10.33 at -3.714 degrees (North through East) ! Clean map min=-0.005732 max=1.7113 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J1517-2422_S_2000_03_13_pus_uvs.log closed on Tue Apr 2 14:48:39 2013