! Started logfile: J2146-1525_S_2011_06_28_pus_uvs.log on Tue Apr 24 22:46:30 2012 obs J2146-1525_S_2011_06_28_pus_uva.fits ! Reading UV FITS file: J2146-1525_S_2011_06_28_pus_uva.fits ! AN table 1: 203 integrations on 153 of 153 possible baselines. ! AN table 2: 131 integrations on 153 of 153 possible baselines. ! AN table 3: 225 integrations on 153 of 153 possible baselines. ! AN table 4: 391 integrations on 153 of 153 possible baselines. ! AN table 5: 106 integrations on 153 of 153 possible baselines. ! Apparent sampling: 0.0744269 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2146-1525 ! ! 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 48100 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 J2146-1525_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 3006 telescope corrections were flagged in sub-array 1. ! A total of 2241 telescope corrections were flagged in sub-array 2. ! A total of 3855 telescope corrections were flagged in sub-array 3. ! A total of 4766 telescope corrections were flagged in sub-array 4. ! A total of 1908 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 3006 telescope corrections were flagged in sub-array 1. ! A total of 2241 telescope corrections were flagged in sub-array 2. ! A total of 3855 telescope corrections were flagged in sub-array 3. ! A total of 4766 telescope corrections were flagged in sub-array 4. ! A total of 1908 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 2994 telescope corrections were flagged in sub-array 1. ! A total of 2241 telescope corrections were flagged in sub-array 2. ! A total of 3821 telescope corrections were flagged in sub-array 3. ! A total of 4766 telescope corrections were flagged in sub-array 4. ! A total of 1863 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 2994 telescope corrections were flagged in sub-array 1. ! A total of 2241 telescope corrections were flagged in sub-array 2. ! A total of 3821 telescope corrections were flagged in sub-array 3. ! A total of 4766 telescope corrections were flagged in sub-array 4. ! A total of 1863 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.690325Jy sigma=0.006523 ! Fit after self-cal, rms=0.652457Jy sigma=0.006367 ! 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.777 mas, bmaj=3.801 mas, bpa=9.415 degrees ! Estimated noise=828.124 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.281923 Jy ! Component: 100 - total flux cleaned = 0.349701 Jy ! Total flux subtracted in 100 components = 0.349701 Jy ! Clean residual min=-0.009922 max=0.027684 Jy/beam ! Clean residual mean=0.000058 rms=0.002928 Jy/beam ! Combined flux in latest and established models = 0.349701 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.349701 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199988Jy sigma=0.001687 ! Fit after self-cal, rms=0.199949Jy sigma=0.001687 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0305322 Jy ! Component: 100 - total flux cleaned = 0.0465071 Jy ! Total flux subtracted in 100 components = 0.0465071 Jy ! Clean residual min=-0.007093 max=0.014624 Jy/beam ! Clean residual mean=0.000081 rms=0.001983 Jy/beam ! Combined flux in latest and established models = 0.396208 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 12 components and 0.396209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.193448Jy sigma=0.001586 ! Fit after self-cal, rms=0.193409Jy sigma=0.001585 ! Inverting map ! Added new window around map position (8.5, 5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0149177 Jy ! Component: 100 - total flux cleaned = 0.024139 Jy ! Total flux subtracted in 100 components = 0.024139 Jy ! Clean residual min=-0.006691 max=0.010036 Jy/beam ! Clean residual mean=0.000094 rms=0.001694 Jy/beam ! Combined flux in latest and established models = 0.420347 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 20 components and 0.420347 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.191219Jy sigma=0.001553 ! Fit after self-cal, rms=0.191076Jy sigma=0.001551 ! 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.257 mas, bmaj=4.639 mas, bpa=9.228 degrees ! Estimated noise=348.252 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 (5.5, 4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0215973 Jy ! Component: 100 - total flux cleaned = 0.0341725 Jy ! Component: 150 - total flux cleaned = 0.042487 Jy ! Component: 200 - total flux cleaned = 0.0482643 Jy ! Total flux subtracted in 200 components = 0.0482643 Jy ! Clean residual min=-0.004244 max=0.007949 Jy/beam ! Clean residual mean=0.000964 rms=0.001474 Jy/beam ! Combined flux in latest and established models = 0.468612 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 59 components and 0.468612 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.188419Jy sigma=0.001518 ! Fit after self-cal, rms=0.188215Jy sigma=0.001515 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 369 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.62 FD 1.72 FT 1.00* HH 1.06 ! HN 1.00* KK 1.79 KP 0.53 MC 0.91 ! MK 1.46 NL 1.00* NY 1.00* OV 1.63 ! PT 0.52 SC 1.00* SH 1.68 TS 1.00* ! WF 1.00* WZ 1.01 ! ! A total of 117 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 66 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.92 FD 1.02 FT 1.00* HH 1.00* ! HN 1.00* KK 1.02 KP 1.01 MC 1.00* ! MK 1.04 NL 1.01 NY 1.00* OV 0.91 ! PT 1.02 SC 1.00* SH 1.00* TS 1.00* ! WF 1.00* WZ 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.94 FD 1.01 FT 1.00* HH 1.01 ! HN 1.02 KK 0.98 KP 0.97 MC 0.97 ! MK 1.02 NL 1.00 NY 1.00* OV 0.95 ! PT 1.00 SC 1.02 SH 1.00* TS 1.00* ! WF 1.05 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! 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 369 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.79 FD 1.72 FT 1.00* HH 1.12 ! HN 1.00* KK 1.55 KP 0.52 MC 0.85 ! MK 1.50 NL 1.00* NY 1.00* OV 1.73 ! PT 0.50 SC 1.00* SH 1.47 TS 1.00* ! WF 1.00* WZ 0.92 ! ! A total of 117 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 66 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.97 FD 0.96 FT 1.00* HH 1.00* ! HN 1.00* KK 1.03 KP 1.02 MC 1.00* ! MK 1.00 NL 1.00 NY 1.00* OV 0.97 ! PT 1.00 SC 1.00* SH 1.00* TS 1.00* ! WF 1.00* WZ 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.98 FD 0.99 FT 1.00* HH 1.01 ! HN 1.02 KK 0.97 KP 1.01 MC 0.93 ! MK 1.00 NL 0.97 NY 1.00* OV 1.00 ! PT 0.98 SC 1.00 SH 1.00* TS 1.00* ! WF 1.00 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! 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 369 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.38 FD 1.29 FT 1.00* HH 1.10 ! HN 1.00* KK 1.22 KP 0.65 MC 1.00 ! MK 1.18 NL 1.00* NY 1.00* OV 1.34 ! PT 0.66 SC 1.00* SH 1.27 TS 1.14 ! WF 1.00* WZ 0.88 ! ! A total of 117 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 1.02 FT 1.00* HH 1.00* ! HN 1.00* KK 1.10 KP 1.01 MC 1.00* ! MK 0.97 NL 1.07 NY 1.00* OV 0.97 ! PT 1.04 SC 1.00* SH 0.98 TS 1.03 ! WF 1.00* WZ 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 FD 0.99 FT 1.00* HH 1.03 ! HN 1.03 KK 0.98 KP 1.02 MC 0.97 ! MK 0.99 NL 1.03 NY 1.00* OV 1.02 ! PT 1.02 SC 1.01 SH 1.00* TS 1.00* ! WF 1.04 WZ 1.00* ! ! A total of 45 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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 369 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.20 FD 1.24 FT 1.00* HH 1.11 ! HN 1.00* KK 1.08 KP 0.74 MC 0.83 ! MK 1.14 NL 1.00* NY 1.00* OV 1.19 ! PT 0.69 SC 1.00* SH 1.08 TS 1.24 ! WF 1.00* WZ 0.95 ! ! A total of 117 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 24 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 1.00 FT 1.00* HH 1.00* ! HN 1.00* KK 1.00 KP 1.03 MC 1.00* ! MK 1.00 NL 1.06 NY 1.00* OV 0.96 ! PT 1.03 SC 1.00* SH 0.92 TS 0.92 ! WF 1.00* WZ 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.96 FD 1.00 FT 1.00* HH 1.01 ! HN 1.01 KK 0.98 KP 1.00 MC 0.85 ! MK 1.01 NL 1.00 NY 1.00* OV 0.99 ! PT 1.00 SC 1.02 SH 1.00* TS 1.00* ! WF 1.02 WZ 1.00* ! ! A total of 45 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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.188215Jy sigma=0.001515 ! Fit after self-cal, rms=0.182220Jy sigma=0.001472 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.246 mas, bmaj=5.001 mas, bpa=8.584 degrees ! Estimated noise=350.85 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 (18, 51.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0051757 Jy ! Component: 100 - total flux cleaned = 0.00830487 Jy ! Component: 150 - total flux cleaned = 0.00947117 Jy ! Component: 200 - total flux cleaned = 0.0101094 Jy ! Total flux subtracted in 200 components = 0.0101094 Jy ! Clean residual min=-0.002940 max=0.004457 Jy/beam ! Clean residual mean=0.000399 rms=0.000992 Jy/beam ! Combined flux in latest and established models = 0.478721 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 80 components and 0.478721 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181894Jy sigma=0.001470 ! Fit after self-cal, rms=0.182776Jy sigma=0.001462 ! 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.182776Jy sigma=0.001462 ! Fit after self-cal, rms=0.192552Jy sigma=0.001460 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.247 mas, bmaj=4.999 mas, bpa=8.59 degrees ! Estimated noise=351.851 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00302656 Jy ! Component: 100 - total flux cleaned = 0.00468931 Jy ! Component: 150 - total flux cleaned = 0.00556641 Jy ! Component: 200 - total flux cleaned = 0.00574349 Jy ! Total flux subtracted in 200 components = 0.00574349 Jy ! Clean residual min=-0.003152 max=0.004668 Jy/beam ! Clean residual mean=0.000348 rms=0.001008 Jy/beam ! Combined flux in latest and established models = 0.484465 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 92 components and 0.484465 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192368Jy sigma=0.001460 ! Fit after self-cal, rms=0.192473Jy sigma=0.001459 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 (-1, -33.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00402478 Jy ! Component: 100 - total flux cleaned = 0.00470061 Jy ! Component: 150 - total flux cleaned = 0.00487173 Jy ! Component: 200 - total flux cleaned = 0.00511158 Jy ! Total flux subtracted in 200 components = 0.00511158 Jy ! Clean residual min=-0.002943 max=0.004761 Jy/beam ! Clean residual mean=0.000322 rms=0.001004 Jy/beam ! Combined flux in latest and established models = 0.489576 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 108 components and 0.489576 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192196Jy sigma=0.001458 ! Fit after self-cal, rms=0.192176Jy sigma=0.001458 ! Inverting map ! Added new window around map position (1, 34). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00376964 Jy ! Component: 100 - total flux cleaned = 0.00479588 Jy ! Component: 150 - total flux cleaned = 0.0052176 Jy ! Component: 200 - total flux cleaned = 0.0053754 Jy ! Total flux subtracted in 200 components = 0.0053754 Jy ! Clean residual min=-0.003141 max=0.004067 Jy/beam ! Clean residual mean=0.000289 rms=0.001012 Jy/beam ! Combined flux in latest and established models = 0.494952 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 123 components and 0.494952 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.191942Jy sigma=0.001457 ! Fit after self-cal, rms=0.191931Jy sigma=0.001457 ! 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.191931Jy sigma=0.001457 ! Fit after self-cal, rms=0.187261Jy sigma=0.001457 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.245 mas, bmaj=4.993 mas, bpa=8.628 degrees ! Estimated noise=352.505 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00192395 Jy ! Component: 100 - total flux cleaned = 0.00232662 Jy ! Component: 150 - total flux cleaned = 0.00261986 Jy ! Component: 200 - total flux cleaned = 0.00276085 Jy ! Total flux subtracted in 200 components = 0.00276085 Jy ! Clean residual min=-0.003278 max=0.004063 Jy/beam ! Clean residual mean=0.000265 rms=0.001011 Jy/beam ! Combined flux in latest and established models = 0.497713 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 132 components and 0.497713 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187113Jy sigma=0.001456 ! Fit after self-cal, rms=0.187073Jy sigma=0.001456 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.000577699 Jy ! Component: 100 - total flux cleaned = 0.000577465 Jy ! Component: 150 - total flux cleaned = 0.000638242 Jy ! Component: 200 - total flux cleaned = 0.000577383 Jy ! Total flux subtracted in 200 components = 0.000577383 Jy ! Clean residual min=-0.003312 max=0.004090 Jy/beam ! Clean residual mean=0.000262 rms=0.001016 Jy/beam ! Combined flux in latest and established models = 0.49829 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 140 components and 0.49829 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.186984Jy sigma=0.001456 ! Fit after self-cal, rms=0.186983Jy sigma=0.001456 ! 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.186983Jy sigma=0.001456 ! Fit after self-cal, rms=0.185893Jy sigma=0.001456 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.141 mas, bmaj=14.3 mas, bpa=-2.536 degrees ! Estimated noise=447.193 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000405579 Jy ! Component: 100 - total flux cleaned = -0.000181331 Jy ! Component: 150 - total flux cleaned = 8.55627e-05 Jy ! Component: 200 - total flux cleaned = 0.00020387 Jy ! Total flux subtracted in 200 components = 0.00020387 Jy ! Clean residual min=-0.004756 max=0.006129 Jy/beam ! Clean residual mean=0.000639 rms=0.001710 Jy/beam ! Combined flux in latest and established models = 0.498494 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 153 components and 0.498494 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185843Jy sigma=0.001456 ! Fit after self-cal, rms=0.185824Jy sigma=0.001456 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.000140508 Jy ! Component: 100 - total flux cleaned = 0.000140854 Jy ! Component: 150 - total flux cleaned = 0.000141709 Jy ! Component: 200 - total flux cleaned = 0.000141658 Jy ! Total flux subtracted in 200 components = 0.000141658 Jy ! Clean residual min=-0.004695 max=0.006141 Jy/beam ! Clean residual mean=0.000637 rms=0.001704 Jy/beam ! Combined flux in latest and established models = 0.498635 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 162 components and 0.498635 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185810Jy sigma=0.001456 ! Fit after self-cal, rms=0.185811Jy sigma=0.001456 ! 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.243 mas, bmaj=4.987 mas, bpa=8.697 degrees ! Estimated noise=353.47 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.243 x 4.987 at 8.697 degrees (North through East) ! Clean map min=-0.0043297 max=0.36959 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.76 mas, bmaj=4.228 mas, bpa=10.95 degrees ! Estimated noise=867.113 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.293413 Jy ! Component: 100 - total flux cleaned = 0.365778 Jy ! Component: 150 - total flux cleaned = 0.401031 Jy ! Component: 200 - total flux cleaned = 0.42495 Jy ! Total flux subtracted in 200 components = 0.42495 Jy ! Clean residual min=-0.009883 max=0.016512 Jy/beam ! Clean residual mean=0.000055 rms=0.003474 Jy/beam ! Combined flux in latest and established models = 0.42495 Jy uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 25 model components to the UV plane model. ! The established model now contains 25 components and 0.42495 Jy ! Inverting map and beam ! Estimated beam: bmin=2.243 mas, bmaj=4.987 mas, bpa=8.697 degrees ! Estimated noise=353.47 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.0240486 Jy ! Component: 100 - total flux cleaned = 0.0383151 Jy ! Component: 150 - total flux cleaned = 0.0473727 Jy ! Component: 200 - total flux cleaned = 0.0534872 Jy ! Total flux subtracted in 200 components = 0.0534872 Jy ! Clean residual min=-0.003668 max=0.004783 Jy/beam ! Clean residual mean=0.000383 rms=0.001161 Jy/beam ! Combined flux in latest and established models = 0.478437 Jy ! Adding 49 model components to the UV plane model. ! The established model now contains 73 components and 0.478437 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.986 mas, bmaj=4.26 mas, bpa=12.09 degrees ! Estimated noise=454.511 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00927335 Jy ! Component: 100 - total flux cleaned = 0.0153894 Jy ! Component: 150 - total flux cleaned = 0.0186677 Jy ! Component: 200 - total flux cleaned = 0.0203442 Jy ! Total flux subtracted in 200 components = 0.0203442 Jy ! Clean residual min=-0.005626 max=0.007319 Jy/beam ! Clean residual mean=0.000012 rms=0.001778 Jy/beam ! Combined flux in latest and established models = 0.498781 Jy clean ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00153715 Jy ! Component: 100 - total flux cleaned = 0.00242659 Jy ! Component: 150 - total flux cleaned = 0.00342991 Jy ! Component: 200 - total flux cleaned = 0.00437123 Jy ! Total flux subtracted in 200 components = 0.00437123 Jy ! Clean residual min=-0.005849 max=0.007041 Jy/beam ! Clean residual mean=-0.000041 rms=0.001762 Jy/beam ! Combined flux in latest and established models = 0.503152 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 122 components and 0.503152 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185504Jy sigma=0.001459 ! Fit after self-cal, rms=0.185472Jy sigma=0.001459 wmodel J2146-1525_S_2011_06_28_pus_map.mod ! Writing 122 model components to file: J2146-1525_S_2011_06_28_pus_map.mod wobs J2146-1525_S_2011_06_28_pus_uvs.fits ! Writing UV FITS file: J2146-1525_S_2011_06_28_pus_uvs.fits wwins J2146-1525_S_2011_06_28_pus_map.win ! wwins: Wrote 6 windows to J2146-1525_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.00326175 Jy ! Component: 100 - total flux cleaned = 0.0062633 Jy ! Component: 150 - total flux cleaned = 0.00886653 Jy ! Component: 200 - total flux cleaned = 0.0111442 Jy ! Component: 250 - total flux cleaned = 0.0133521 Jy ! Component: 300 - total flux cleaned = 0.0154931 Jy ! Component: 350 - total flux cleaned = 0.0175713 Jy ! Component: 400 - total flux cleaned = 0.0193049 Jy ! Component: 450 - total flux cleaned = 0.0206208 Jy ! Component: 500 - total flux cleaned = 0.0220014 Jy ! Component: 550 - total flux cleaned = 0.0231731 Jy ! Component: 600 - total flux cleaned = 0.0242339 Jy ! Component: 650 - total flux cleaned = 0.025361 Jy ! Component: 700 - total flux cleaned = 0.0262968 Jy ! Component: 750 - total flux cleaned = 0.0272171 Jy ! Component: 800 - total flux cleaned = 0.0280399 Jy ! Component: 850 - total flux cleaned = 0.0289307 Jy ! Component: 900 - total flux cleaned = 0.029489 Jy ! Component: 950 - total flux cleaned = 0.0302758 Jy ! Component: 1000 - total flux cleaned = 0.0308969 Jy ! Component: 1050 - total flux cleaned = 0.0313578 Jy ! Component: 1100 - total flux cleaned = 0.0317383 Jy ! Component: 1150 - total flux cleaned = 0.0322644 Jy ! Component: 1200 - total flux cleaned = 0.0324876 Jy ! Component: 1250 - total flux cleaned = 0.0327831 Jy ! Component: 1300 - total flux cleaned = 0.0332211 Jy ! Component: 1350 - total flux cleaned = 0.03351 Jy ! Component: 1400 - total flux cleaned = 0.033725 Jy ! Component: 1450 - total flux cleaned = 0.0339385 Jy ! Component: 1500 - total flux cleaned = 0.0340796 Jy ! Component: 1550 - total flux cleaned = 0.0343597 Jy ! Component: 1600 - total flux cleaned = 0.0343598 Jy ! Component: 1650 - total flux cleaned = 0.0345667 Jy ! Component: 1700 - total flux cleaned = 0.0347034 Jy ! Component: 1750 - total flux cleaned = 0.0346352 Jy ! Component: 1800 - total flux cleaned = 0.0347031 Jy ! Component: 1850 - total flux cleaned = 0.0346358 Jy ! Component: 1900 - total flux cleaned = 0.0346358 Jy ! Component: 1950 - total flux cleaned = 0.0347023 Jy ! Component: 2000 - total flux cleaned = 0.0346364 Jy ! Total flux subtracted in 2048 components = 0.0347021 Jy ! Clean residual min=-0.003959 max=0.004517 Jy/beam ! Clean residual mean=-0.000605 rms=0.001403 Jy/beam ! Combined flux in latest and established models = 0.537854 Jy keep ! Adding 175 model components to the UV plane model. ! The established model now contains 297 components and 0.537854 Jy ![Exited script file: muppet] wmap "J2146-1525_S_2011_06_28_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.986 x 4.26 at 12.09 degrees (North through East) ! Clean map min=-0.0052243 max=0.37174 Jy/beam ! Writing clean map to FITS file: J2146-1525_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.371736 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00140035 device "J2146-1525_S_2011_06_28_pus\_map.eps/VPS" ! Attempting to open device: 'J2146-1525_S_2011_06_28_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.13011 1.13011 2.26023 4.52045 9.0409 18.0818 36.1636 72.3272 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0042 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.371 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.986 x 4.26 at 12.09 degrees (North through East) ! Clean map min=-0.0052243 max=0.37174 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J2146-1525_S_2011_06_28_pus\_map.gif/vgif" ! Attempting to open device: 'J2146-1525_S_2011_06_28_pus_map.gif/vgif' mapplot cln