! Started logfile: 0234+285_X.log on Thu Jun 26 15:49:31 2008 observe 0234+285_X.edt ! Reading UV FITS file: 0234+285_X.edt ! AN table 1: 114 integrations on 190 of 190 possible baselines. ! AN table 2: 20 integrations on 190 of 190 possible baselines. ! AN table 3: 116 integrations on 190 of 190 possible baselines. ! AN table 4: 14 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.174601 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 0234+285 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40997e+09 8e+06 1 8e+06 ! 02 2 8.47997e+09 8e+06 1 8e+06 ! 03 3 8.79497e+09 8e+06 1 8e+06 ! 04 4 8.89997e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1472 lines of history. ! ! Reading 35032 visibilities. ![@gonzo6-0.1 0234+285_X] print "***********************************" ! *********************************** print "*** Using gonzo instruction set ***" ! *** Using gonzo instruction set *** print "***********************************" ! *********************************** float field_size field_size = 1024 float field_cell field_cell = 0.1 integer clean_niter clean_niter = 50 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float win_mult win_mult = 1.8 float solint solint = 3.6/6.0 logical doflt doflt = true logical doamp float new_peak float flux_cutoff mapsize field_size,field_cell ! Map grid = 1024x1024 pixels with 0.100x0.100 milli-arcsec cellsize. select ! 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 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 1749 telescope corrections were flagged in sub-array 1. ! A total of 247 telescope corrections were flagged in sub-array 2. ! A total of 1860 telescope corrections were flagged in sub-array 3. ! A total of 247 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 1749 telescope corrections were flagged in sub-array 1. ! A total of 247 telescope corrections were flagged in sub-array 2. ! A total of 1803 telescope corrections were flagged in sub-array 3. ! A total of 242 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 1746 telescope corrections were flagged in sub-array 1. ! A total of 247 telescope corrections were flagged in sub-array 2. ! A total of 1803 telescope corrections were flagged in sub-array 3. ! A total of 247 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 1753 telescope corrections were flagged in sub-array 1. ! A total of 247 telescope corrections were flagged in sub-array 2. ! A total of 1804 telescope corrections were flagged in sub-array 3. ! A total of 242 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.105918Jy sigma=19.401350 ! Fit after self-cal, rms=1.105118Jy sigma=19.415702 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvweight 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. doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.5071 mas, bmaj=0.9422 mas, bpa=-6.331 degrees ! Estimated noise=0.649697 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal doamp, doflt, solint new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.27163 Jy ! Total flux subtracted in 50 components = 1.27163 Jy ! Clean residual min=-0.047430 max=0.354882 Jy/beam ! Clean residual mean=0.000224 rms=0.013539 Jy/beam ! Combined flux in latest and established models = 1.27163 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1.27163 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.930508Jy sigma=13.993909 ! Fit after self-cal, rms=0.930508Jy sigma=13.993914 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.293779 Jy ! Total flux subtracted in 50 components = 0.293779 Jy ! Clean residual min=-0.061516 max=0.114884 Jy/beam ! Clean residual mean=0.000241 rms=0.010724 Jy/beam ! Combined flux in latest and established models = 1.56541 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.56541 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.806849Jy sigma=9.247208 ! Fit after self-cal, rms=0.806817Jy sigma=9.246867 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.124892 Jy ! Total flux subtracted in 50 components = 0.124892 Jy ! Clean residual min=-0.069825 max=0.057190 Jy/beam ! Clean residual mean=0.000248 rms=0.010274 Jy/beam ! Combined flux in latest and established models = 1.6903 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 4 model components to the UV plane model. ! The established model now contains 9 components and 1.6903 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.778954Jy sigma=7.854984 ! Fit after self-cal, rms=0.778914Jy sigma=7.854484 ! Inverting map print "****************************************" ! **************************************** print "*** Finished Uniform Weighting CLEAN ***" ! *** Finished Uniform Weighting CLEAN *** print "****************************************" ! **************************************** uvweight 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 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6449 mas, bmaj=1.142 mas, bpa=-9.394 degrees ! Estimated noise=0.382405 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal doamp, doflt, solint new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-1, 4.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0840042 Jy ! Component: 100 - total flux cleaned = 0.137287 Jy ! Total flux subtracted in 100 components = 0.137287 Jy ! Clean residual min=-0.056147 max=0.050007 Jy/beam ! Clean residual mean=0.001622 rms=0.010639 Jy/beam ! Combined flux in latest and established models = 1.82759 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 21 model components to the UV plane model. ! The established model now contains 29 components and 1.82759 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.765967Jy sigma=6.689073 ! Fit after self-cal, rms=0.764874Jy sigma=6.651994 ! Inverting map print "****************************************" ! **************************************** print "*** Finished Natural Weighting CLEAN ***" ! *** Finished Natural Weighting CLEAN *** print "****************************************" ! **************************************** gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 162 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.72 FD 1.15 GC 1.04 HH 1.00* ! HN 1.07 KK 0.95 KP 1.08 LA 1.09 ! MA 1.00* MC 1.00* MK 1.01 NL 1.00 ! NY 1.18 ON 1.00* OV 1.06 PT 1.00 ! SC 1.17 TS 0.93 WF 1.01 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.20 FD 1.11 GC 1.05 HH 1.00* ! HN 1.01 KK 1.00 KP 1.17 LA 1.22 ! MA 1.00* MC 1.00* MK 0.79 NL 0.94 ! NY 1.27 ON 1.00* OV 0.98 PT 0.96 ! SC 1.00* TS 0.88 WF 0.65 WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.54 FD 1.12 GC 1.12 HH 1.00* ! HN 1.00 KK 1.22 KP 1.06 LA 0.95 ! MA 0.88 MC 1.16 MK 1.05 NL 1.04 ! NY 0.50 ON 0.39 OV 1.05 PT 1.00 ! SC 1.05 TS 1.00* WF 0.96 WZ 1.13 ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! A total of 162 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.59 FD 1.04 GC 1.04 HH 1.00* ! HN 0.96 KK 1.07 KP 1.07 LA 0.87 ! MA 1.00* MC 1.00* MK 0.87 NL 1.08 ! NY 1.09 ON 1.00* OV 1.03 PT 1.09 ! SC 1.15 TS 0.98 WF 0.83 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.42 FD 1.10 GC 1.09 HH 1.00* ! HN 1.15 KK 1.00 KP 1.15 LA 1.32 ! MA 1.00* MC 1.00* MK 0.75 NL 0.92 ! NY 1.22 ON 1.00* OV 0.99 PT 0.89 ! SC 1.00* TS 0.88 WF 0.75 WZ 1.00* ! ! A total of 54 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.47 FD 1.12 GC 0.91 HH 1.00* ! HN 1.08 KK 1.23 KP 1.02 LA 0.94 ! MA 0.87 MC 1.16 MK 0.99 NL 1.02 ! NY 0.40 ON 1.17 OV 1.09 PT 0.97 ! SC 1.04 TS 1.00* WF 0.85 WZ 1.17 ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 1.33 KK 1.00* KP 1.00* LA 1.00* ! MA 1.00* MC 1.00* MK 1.00* NL 1.00* ! NY 1.16 ON 1.00* OV 1.00* PT 1.00* ! SC 1.18 TS 1.00* WF 0.94 WZ 1.00* ! ! ! Correcting IF 3. ! A total of 162 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.53 FD 1.12 GC 0.89 HH 1.00* ! HN 1.02 KK 1.00 KP 0.92 LA 1.06 ! MA 1.00* MC 1.00* MK 0.84 NL 1.02 ! NY 1.17 ON 1.00* OV 1.00 PT 1.01 ! SC 1.06 TS 0.89 WF 0.82 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.15 FD 1.15 GC 0.98 HH 1.00* ! HN 0.91 KK 0.95 KP 1.16 LA 1.32 ! MA 1.00* MC 1.00* MK 0.94 NL 0.90 ! NY 1.28 ON 1.00* OV 0.87 PT 1.08 ! SC 1.00* TS 0.85 WF 0.66 WZ 1.00* ! ! A total of 54 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.49 FD 1.11 GC 1.07 HH 1.00* ! HN 1.00 KK 1.17 KP 1.01 LA 0.93 ! MA 1.02 MC 1.21 MK 1.06 NL 1.04 ! NY 0.44 ON 1.18 OV 1.09 PT 0.84 ! SC 1.04 TS 1.00* WF 0.81 WZ 1.14 ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! A total of 159 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.59 FD 1.07 GC 1.04 HH 1.00* ! HN 0.99 KK 1.02 KP 1.01 LA 1.01 ! MA 1.00* MC 1.00* MK 0.89 NL 1.06 ! NY 0.88 ON 1.00* OV 0.97 PT 1.03 ! SC 0.96 TS 0.92 WF 0.88 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.31 FD 1.20 GC 1.07 HH 1.00* ! HN 0.77 KK 0.94 KP 1.11 LA 1.21 ! MA 1.00* MC 1.00* MK 0.78 NL 0.87 ! NY 1.29 ON 1.00* OV 0.89 PT 1.04 ! SC 1.00* TS 0.86 WF 0.80 WZ 1.00* ! ! A total of 54 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.36 FD 1.10 GC 1.10 HH 1.00* ! HN 1.04 KK 1.14 KP 1.00 LA 0.91 ! MA 1.03 MC 1.28 MK 0.75 NL 1.08 ! NY 0.63 ON 1.07 OV 1.11 PT 0.98 ! SC 1.10 TS 1.00* WF 0.82 WZ 1.20 ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 4.04 KK 1.00* KP 1.00* LA 1.00* ! MA 1.00* MC 1.00* MK 1.00* NL 1.00* ! NY 0.66 ON 1.00* OV 1.00* PT 1.00* ! SC 1.16 TS 1.00* WF 1.27 WZ 1.00* ! ! ! Fit before self-cal, rms=0.764874Jy sigma=6.651994 ! Fit after self-cal, rms=0.423013Jy sigma=2.536628 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6429 mas, bmaj=1.172 mas, bpa=-11.3 degrees ! Estimated noise=0.381888 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal doamp, doflt, solint new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-0.8, 2.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0331877 Jy ! Component: 100 - total flux cleaned = 0.0500154 Jy ! Total flux subtracted in 100 components = 0.0500154 Jy ! Clean residual min=-0.015856 max=0.021954 Jy/beam ! Clean residual mean=0.000315 rms=0.003693 Jy/beam ! Combined flux in latest and established models = 1.8776 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 23 model components to the UV plane model. ! The established model now contains 51 components and 1.8776 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.420222Jy sigma=2.580484 ! Fit after self-cal, rms=0.419768Jy sigma=2.562940 ! Inverting map ! Added new window around map position (-0.9, 7.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0204002 Jy ! Component: 100 - total flux cleaned = 0.0284914 Jy ! Total flux subtracted in 100 components = 0.0284914 Jy ! Clean residual min=-0.015950 max=0.013210 Jy/beam ! Clean residual mean=0.000208 rms=0.003485 Jy/beam ! Combined flux in latest and established models = 1.9061 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 24 model components to the UV plane model. ! The established model now contains 67 components and 1.9061 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.418168Jy sigma=2.718064 ! Fit after self-cal, rms=0.418089Jy sigma=2.713225 ! Inverting map dynam = dynam - 1.0 doamp = true solint = solint * (6.0/3.6) 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 doamp, doflt, solint new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00687548 Jy ! Component: 100 - total flux cleaned = 0.0107255 Jy ! Total flux subtracted in 100 components = 0.0107255 Jy ! Clean residual min=-0.016152 max=0.012073 Jy/beam ! Clean residual mean=0.000170 rms=0.003432 Jy/beam ! Combined flux in latest and established models = 1.91682 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 31 model components to the UV plane model. ! The established model now contains 81 components and 1.91682 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.417353Jy sigma=2.789347 ! Fit after self-cal, rms=0.439108Jy sigma=1.510855 ! Inverting map and beam ! Estimated beam: bmin=0.6399 mas, bmaj=1.167 mas, bpa=-11.24 degrees ! Estimated noise=0.383076 mJy/beam. ! Added new window around map position (-0.3, 0.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0161381 Jy ! Component: 100 - total flux cleaned = 0.0266898 Jy ! Total flux subtracted in 100 components = 0.0266898 Jy ! Clean residual min=-0.005609 max=0.009015 Jy/beam ! Clean residual mean=0.000223 rms=0.001439 Jy/beam ! Combined flux in latest and established models = 1.94351 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 28 model components to the UV plane model. ! The established model now contains 108 components and 1.94351 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.438432Jy sigma=1.538999 ! Fit after self-cal, rms=0.422671Jy sigma=1.455763 ! Inverting map and beam ! Estimated beam: bmin=0.6395 mas, bmaj=1.164 mas, bpa=-11.18 degrees ! Estimated noise=0.383706 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0125189 Jy ! Component: 100 - total flux cleaned = 0.0195342 Jy ! Total flux subtracted in 100 components = 0.0195342 Jy ! Clean residual min=-0.005300 max=0.008534 Jy/beam ! Clean residual mean=0.000219 rms=0.001316 Jy/beam ! Combined flux in latest and established models = 1.96305 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 26 model components to the UV plane model. ! The established model now contains 118 components and 1.96305 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.422196Jy sigma=1.469315 ! Fit after self-cal, rms=0.433202Jy sigma=1.425475 ! Inverting map and beam ! Estimated beam: bmin=0.639 mas, bmaj=1.162 mas, bpa=-11.14 degrees ! Estimated noise=0.384301 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0106005 Jy ! Component: 100 - total flux cleaned = 0.0152308 Jy ! Total flux subtracted in 100 components = 0.0152308 Jy ! Clean residual min=-0.004740 max=0.008276 Jy/beam ! Clean residual mean=0.000213 rms=0.001232 Jy/beam ! Combined flux in latest and established models = 1.97828 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 27 model components to the UV plane model. ! The established model now contains 127 components and 1.97828 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.432859Jy sigma=1.431767 ! Fit after self-cal, rms=0.443268Jy sigma=1.405113 ! Inverting map and beam ! Estimated beam: bmin=0.6385 mas, bmaj=1.161 mas, bpa=-11.12 degrees ! Estimated noise=0.38477 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00922729 Jy ! Component: 100 - total flux cleaned = 0.012546 Jy ! Total flux subtracted in 100 components = 0.012546 Jy ! Clean residual min=-0.004291 max=0.008112 Jy/beam ! Clean residual mean=0.000207 rms=0.001166 Jy/beam ! Combined flux in latest and established models = 1.99082 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 31 model components to the UV plane model. ! The established model now contains 138 components and 1.99082 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.443003Jy sigma=1.408191 ! Fit after self-cal, rms=0.452105Jy sigma=1.390297 ! Inverting map and beam ! Estimated beam: bmin=0.6381 mas, bmaj=1.159 mas, bpa=-11.11 degrees ! Estimated noise=0.385171 mJy/beam. ! Added new window around map position (-2.5, 5.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00944998 Jy ! Component: 100 - total flux cleaned = 0.0135814 Jy ! Total flux subtracted in 100 components = 0.0135814 Jy ! Clean residual min=-0.004056 max=0.006132 Jy/beam ! Clean residual mean=0.000191 rms=0.001082 Jy/beam ! Combined flux in latest and established models = 2.0044 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 20 model components to the UV plane model. ! The established model now contains 147 components and 2.0044 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.451879Jy sigma=1.394146 ! Fit after self-cal, rms=0.460750Jy sigma=1.374898 ! Inverting map and beam ! Estimated beam: bmin=0.6377 mas, bmaj=1.158 mas, bpa=-11.1 degrees ! Estimated noise=0.385537 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00765079 Jy ! Component: 100 - total flux cleaned = 0.0105294 Jy ! Total flux subtracted in 100 components = 0.0105294 Jy ! Clean residual min=-0.003799 max=0.005855 Jy/beam ! Clean residual mean=0.000184 rms=0.001031 Jy/beam ! Combined flux in latest and established models = 2.01493 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 29 model components to the UV plane model. ! The established model now contains 157 components and 2.01493 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.460582Jy sigma=1.376732 ! Fit after self-cal, rms=0.468022Jy sigma=1.364982 ! Inverting map and beam ! Estimated beam: bmin=0.6374 mas, bmaj=1.157 mas, bpa=-11.09 degrees ! Estimated noise=0.385838 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00661641 Jy ! Component: 100 - total flux cleaned = 0.00880621 Jy ! Total flux subtracted in 100 components = 0.00880621 Jy ! Clean residual min=-0.003730 max=0.005597 Jy/beam ! Clean residual mean=0.000178 rms=0.000990 Jy/beam ! Combined flux in latest and established models = 2.02374 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 28 model components to the UV plane model. ! The established model now contains 163 components and 2.02374 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.467879Jy sigma=1.365498 ! Fit after self-cal, rms=0.474431Jy sigma=1.357280 ! Inverting map and beam ! Estimated beam: bmin=0.6371 mas, bmaj=1.156 mas, bpa=-11.08 degrees ! Estimated noise=0.386077 mJy/beam. ! Added new window around map position (0.4, 3.7). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00683011 Jy ! Component: 100 - total flux cleaned = 0.0087226 Jy ! Total flux subtracted in 100 components = 0.0087226 Jy ! Clean residual min=-0.003647 max=0.005102 Jy/beam ! Clean residual mean=0.000169 rms=0.000942 Jy/beam ! Combined flux in latest and established models = 2.03246 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 23 model components to the UV plane model. ! The established model now contains 170 components and 2.03246 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.474283Jy sigma=1.357980 ! Fit after self-cal, rms=0.480215Jy sigma=1.350125 ! Inverting map and beam ! Estimated beam: bmin=0.6368 mas, bmaj=1.156 mas, bpa=-11.08 degrees ! Estimated noise=0.386301 mJy/beam. ! Added new window around map position (2, 5.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0059068 Jy ! Component: 100 - total flux cleaned = 0.00788064 Jy ! Total flux subtracted in 100 components = 0.00788064 Jy ! Clean residual min=-0.003585 max=0.004682 Jy/beam ! Clean residual mean=0.000160 rms=0.000901 Jy/beam ! Combined flux in latest and established models = 2.04034 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 29 model components to the UV plane model. ! The established model now contains 175 components and 2.04034 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.480085Jy sigma=1.350380 ! Fit after self-cal, rms=0.485708Jy sigma=1.344246 ! Inverting map and beam ! Estimated beam: bmin=0.6366 mas, bmaj=1.155 mas, bpa=-11.07 degrees ! Estimated noise=0.386476 mJy/beam. ! Added new window around map position (0.9, 6.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00585533 Jy ! Component: 100 - total flux cleaned = 0.00772767 Jy ! Total flux subtracted in 100 components = 0.00772767 Jy ! Clean residual min=-0.003460 max=0.004351 Jy/beam ! Clean residual mean=0.000149 rms=0.000862 Jy/beam ! Combined flux in latest and established models = 2.04807 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 26 model components to the UV plane model. ! The established model now contains 183 components and 2.04807 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.485596Jy sigma=1.345016 ! Fit after self-cal, rms=0.490322Jy sigma=1.339401 ! Inverting map and beam ! Estimated beam: bmin=0.6364 mas, bmaj=1.155 mas, bpa=-11.07 degrees ! Estimated noise=0.386626 mJy/beam. ! Added new window around map position (-2.8, 3.7). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0047257 Jy ! Component: 100 - total flux cleaned = 0.00629899 Jy ! Total flux subtracted in 100 components = 0.00629899 Jy ! Clean residual min=-0.003323 max=0.003742 Jy/beam ! Clean residual mean=0.000145 rms=0.000833 Jy/beam ! Combined flux in latest and established models = 2.05437 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 26 model components to the UV plane model. ! The established model now contains 191 components and 2.05437 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.490226Jy sigma=1.339179 ! Fit after self-cal, rms=0.494794Jy sigma=1.335425 ! Inverting map and beam ! Estimated beam: bmin=0.6363 mas, bmaj=1.154 mas, bpa=-11.07 degrees ! Estimated noise=0.386738 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00362037 Jy ! Component: 100 - total flux cleaned = 0.00475628 Jy ! Total flux subtracted in 100 components = 0.00475628 Jy ! Clean residual min=-0.003220 max=0.003676 Jy/beam ! Clean residual mean=0.000142 rms=0.000812 Jy/beam ! Combined flux in latest and established models = 2.05913 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 35 model components to the UV plane model. ! The established model now contains 202 components and 2.05913 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.494721Jy sigma=1.334375 ! Fit after self-cal, rms=0.498945Jy sigma=1.332162 ! Inverting map and beam ! Estimated beam: bmin=0.6361 mas, bmaj=1.154 mas, bpa=-11.07 degrees ! Estimated noise=0.386824 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00287757 Jy ! Component: 100 - total flux cleaned = 0.00395213 Jy ! Total flux subtracted in 100 components = 0.00395213 Jy ! Clean residual min=-0.003160 max=0.003701 Jy/beam ! Clean residual mean=0.000139 rms=0.000795 Jy/beam ! Combined flux in latest and established models = 2.06308 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 36 model components to the UV plane model. ! The established model now contains 210 components and 2.06308 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.498881Jy sigma=1.330837 ! Fit after self-cal, rms=0.502730Jy sigma=1.329267 ! Inverting map and beam ! Estimated beam: bmin=0.636 mas, bmaj=1.153 mas, bpa=-11.07 degrees ! Estimated noise=0.386881 mJy/beam. ! Added new window around map position (-2.4, 8.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00343774 Jy ! Component: 100 - total flux cleaned = 0.00447979 Jy ! Total flux subtracted in 100 components = 0.00447978 Jy ! Clean residual min=-0.003035 max=0.003362 Jy/beam ! Clean residual mean=0.000132 rms=0.000773 Jy/beam ! Combined flux in latest and established models = 2.06756 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 35 model components to the UV plane model. ! The established model now contains 219 components and 2.06756 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.502668Jy sigma=1.328326 ! Fit after self-cal, rms=0.506155Jy sigma=1.326450 ! Inverting map and beam ! Estimated beam: bmin=0.6359 mas, bmaj=1.153 mas, bpa=-11.06 degrees ! Estimated noise=0.386926 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00241718 Jy ! Component: 100 - total flux cleaned = 0.00339924 Jy ! Total flux subtracted in 100 components = 0.00339924 Jy ! Clean residual min=-0.002974 max=0.003300 Jy/beam ! Clean residual mean=0.000130 rms=0.000758 Jy/beam ! Combined flux in latest and established models = 2.07096 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 35 model components to the UV plane model. ! The established model now contains 222 components and 2.07096 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.506102Jy sigma=1.325254 ! Fit after self-cal, rms=0.509414Jy sigma=1.324184 ! Inverting map and beam ! Estimated beam: bmin=0.6357 mas, bmaj=1.152 mas, bpa=-11.06 degrees ! Estimated noise=0.38694 mJy/beam. ! Added new window around map position (12.6, 21.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00285498 Jy ! Component: 100 - total flux cleaned = 0.00393704 Jy ! Total flux subtracted in 100 components = 0.00393704 Jy ! Clean residual min=-0.002923 max=0.003212 Jy/beam ! Clean residual mean=0.000125 rms=0.000744 Jy/beam ! Combined flux in latest and established models = 2.07489 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 34 model components to the UV plane model. ! The established model now contains 226 components and 2.07489 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.509358Jy sigma=1.322784 ! Fit after self-cal, rms=0.512466Jy sigma=1.322078 ! Inverting map and beam ! Estimated beam: bmin=0.6356 mas, bmaj=1.152 mas, bpa=-11.06 degrees ! Estimated noise=0.386948 mJy/beam. ! Added new window around map position (2.7, 12.1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00272831 Jy ! Component: 100 - total flux cleaned = 0.00377141 Jy ! Total flux subtracted in 100 components = 0.00377141 Jy ! Clean residual min=-0.002867 max=0.003014 Jy/beam ! Clean residual mean=0.000117 rms=0.000726 Jy/beam ! Combined flux in latest and established models = 2.07866 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 35 model components to the UV plane model. ! The established model now contains 233 components and 2.07867 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.512414Jy sigma=1.321242 ! Fit after self-cal, rms=0.514870Jy sigma=1.320255 ! Inverting map and beam ! Estimated beam: bmin=0.6355 mas, bmaj=1.152 mas, bpa=-11.06 degrees ! Estimated noise=0.386937 mJy/beam. ! Added new window around map position (-1.9, -3.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00259154 Jy ! Component: 100 - total flux cleaned = 0.0038459 Jy ! Total flux subtracted in 100 components = 0.0038459 Jy ! Clean residual min=-0.002850 max=0.002939 Jy/beam ! Clean residual mean=0.000112 rms=0.000712 Jy/beam ! Combined flux in latest and established models = 2.08251 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 45 model components to the UV plane model. ! The established model now contains 241 components and 2.08251 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.514822Jy sigma=1.319532 ! Fit after self-cal, rms=0.517181Jy sigma=1.318470 ! Inverting map and beam ! Estimated beam: bmin=0.6354 mas, bmaj=1.151 mas, bpa=-11.05 degrees ! Estimated noise=0.386911 mJy/beam. ! Added new window around map position (0.2, 19.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00251365 Jy ! Component: 100 - total flux cleaned = 0.00384769 Jy ! Total flux subtracted in 100 components = 0.00384769 Jy ! Clean residual min=-0.002798 max=0.002808 Jy/beam ! Clean residual mean=0.000104 rms=0.000698 Jy/beam ! Combined flux in latest and established models = 2.08636 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 40 model components to the UV plane model. ! The established model now contains 247 components and 2.08636 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.517135Jy sigma=1.317875 ! Fit after self-cal, rms=0.519495Jy sigma=1.316858 ! Inverting map and beam ! Estimated beam: bmin=0.6353 mas, bmaj=1.151 mas, bpa=-11.05 degrees ! Estimated noise=0.386881 mJy/beam. ! Added new window around map position (-3.8, 13.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0020082 Jy ! Component: 100 - total flux cleaned = 0.00329226 Jy ! Total flux subtracted in 100 components = 0.00329226 Jy ! Clean residual min=-0.002764 max=0.002710 Jy/beam ! Clean residual mean=0.000098 rms=0.000686 Jy/beam ! Combined flux in latest and established models = 2.08965 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 40 model components to the UV plane model. ! The established model now contains 254 components and 2.08965 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.519451Jy sigma=1.316307 ! Fit after self-cal, rms=0.521630Jy sigma=1.315541 ! Inverting map and beam ! Estimated beam: bmin=0.6352 mas, bmaj=1.151 mas, bpa=-11.05 degrees ! Estimated noise=0.386829 mJy/beam. doamp = false solint = solint / (6.0/3.6) 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 doamp, doflt, solint new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00143713 Jy ! Component: 100 - total flux cleaned = 0.00244047 Jy ! Total flux subtracted in 100 components = 0.00244047 Jy ! Clean residual min=-0.002717 max=0.002688 Jy/beam ! Clean residual mean=0.000096 rms=0.000676 Jy/beam ! Combined flux in latest and established models = 2.09209 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 49 model components to the UV plane model. ! The established model now contains 256 components and 2.09209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.521590Jy sigma=1.314824 ! Fit after self-cal, rms=0.521529Jy sigma=1.306636 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 0234+285_X ! Writing UV FITS file: 0234+285_X.uvf ! Writing 256 model components to file: 0234+285_X.mod ! wwins: Wrote 16 windows to 0234+285_X.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6352 x 1.151 at -11.05 degrees (North through East) ! Clean map min=-0.0026463 max=1.6321 Jy/beam ! Writing clean map to FITS file: 0234+285_X.fits ! Writing difmap environment to: 0234+285_X.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 0234+285_X] device 0234+285_Xuv.ps/vps,1,2 ! Attempting to open device: '0234+285_Xuv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 0234+285_Xim.ps/vps ! Attempting to open device: '0234+285_Xim.ps/vps' uvweight 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. loglevs -1,2048,2 ! The new contour levels are: ! -1 1 2 4 8 16 32 64 128 256 512 1024 2048 cmul = imstat(rms)*3 ! Inverting map and beam ! Estimated beam: bmin=0.6352 mas, bmaj=1.151 mas, bpa=-11.05 degrees ! Estimated noise=0.386829 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6352 x 1.151 at -11.05 degrees (North through East) ! Clean map min=-0.0026463 max=1.6321 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo6-0.1] quit ! Quitting program ! Log file 0234+285_X.log closed on Thu Jun 26 15:50:09 2008