! Started logfile: 3C274_X.log on Fri Jul 11 14:29:49 2008 observe 3C274_X.edt ! Reading UV FITS file: 3C274_X.edt ! AN table 1: 126 integrations on 190 of 190 possible baselines. ! AN table 2: 107 integrations on 190 of 190 possible baselines. ! AN table 3: 34 integrations on 190 of 190 possible baselines. ! AN table 4: 38 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.148041 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 3C274 ! ! 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 34316 visibilities. ![@gonzo5-0.1 3C274_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 = 5.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 1790 telescope corrections were flagged in sub-array 1. ! A total of 1826 telescope corrections were flagged in sub-array 2. ! A total of 586 telescope corrections were flagged in sub-array 3. ! A total of 530 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 1787 telescope corrections were flagged in sub-array 1. ! A total of 1851 telescope corrections were flagged in sub-array 2. ! A total of 588 telescope corrections were flagged in sub-array 3. ! A total of 530 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 1790 telescope corrections were flagged in sub-array 1. ! A total of 1826 telescope corrections were flagged in sub-array 2. ! A total of 586 telescope corrections were flagged in sub-array 3. ! A total of 531 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 1791 telescope corrections were flagged in sub-array 1. ! A total of 1829 telescope corrections were flagged in sub-array 2. ! A total of 587 telescope corrections were flagged in sub-array 3. ! A total of 530 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.488883Jy sigma=7.633554 ! Fit after self-cal, rms=0.480190Jy sigma=7.625748 ! 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.5838 mas, bmaj=1.218 mas, bpa=-4.861 degrees ! Estimated noise=0.686611 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 = 0.545794 Jy ! Total flux subtracted in 50 components = 0.545794 Jy ! Clean residual min=-0.035219 max=0.197946 Jy/beam ! Clean residual mean=-0.000735 rms=0.010871 Jy/beam ! Combined flux in latest and established models = 0.545794 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.545794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.569969Jy sigma=8.870706 ! Fit after self-cal, rms=0.569893Jy sigma=8.870168 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.210917 Jy ! Total flux subtracted in 50 components = 0.210917 Jy ! Clean residual min=-0.026497 max=0.114621 Jy/beam ! Clean residual mean=-0.000492 rms=0.007792 Jy/beam ! Combined flux in latest and established models = 0.756712 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 13 model components to the UV plane model. ! The established model now contains 16 components and 0.756712 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.450062Jy sigma=6.605268 ! Fit after self-cal, rms=0.449588Jy sigma=6.601259 ! Inverting map ! Added new window around map position (-0.9, 0.1). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.124713 Jy ! Total flux subtracted in 50 components = 0.124713 Jy ! Clean residual min=-0.022227 max=0.067112 Jy/beam ! Clean residual mean=-0.000353 rms=0.006461 Jy/beam ! Combined flux in latest and established models = 0.881425 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 15 model components to the UV plane model. ! The established model now contains 24 components and 0.881425 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.392486Jy sigma=5.481682 ! Fit after self-cal, rms=0.390771Jy sigma=5.459609 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0830428 Jy ! Total flux subtracted in 50 components = 0.0830428 Jy ! Clean residual min=-0.020595 max=0.051679 Jy/beam ! Clean residual mean=-0.000259 rms=0.005735 Jy/beam ! Combined flux in latest and established models = 0.964468 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 17 model components to the UV plane model. ! The established model now contains 36 components and 0.964468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.358898Jy sigma=4.814529 ! Fit after self-cal, rms=0.358114Jy sigma=4.805741 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0615835 Jy ! Total flux subtracted in 50 components = 0.0615835 Jy ! Clean residual min=-0.018685 max=0.041928 Jy/beam ! Clean residual mean=-0.000191 rms=0.005294 Jy/beam ! Combined flux in latest and established models = 1.02605 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 16 model components to the UV plane model. ! The established model now contains 50 components and 1.02605 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.338510Jy sigma=4.394482 ! Fit after self-cal, rms=0.337993Jy sigma=4.388554 ! Inverting map ! Added new window around map position (-1.6, -0.6). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0494871 Jy ! Total flux subtracted in 50 components = 0.0494871 Jy ! Clean residual min=-0.018188 max=0.035162 Jy/beam ! Clean residual mean=-0.000139 rms=0.004989 Jy/beam ! Combined flux in latest and established models = 1.07554 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 16 model components to the UV plane model. ! The established model now contains 62 components and 1.07554 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.324502Jy sigma=4.104143 ! Fit after self-cal, rms=0.324078Jy sigma=4.099272 ! Inverting map ! Added new window around map position (-2.2, 1.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0436692 Jy ! Total flux subtracted in 50 components = 0.0436692 Jy ! Clean residual min=-0.017378 max=0.030537 Jy/beam ! Clean residual mean=-0.000100 rms=0.004731 Jy/beam ! Combined flux in latest and established models = 1.11921 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 17 model components to the UV plane model. ! The established model now contains 70 components and 1.11921 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.313786Jy sigma=3.883566 ! Fit after self-cal, rms=0.313418Jy sigma=3.878723 ! Inverting map ! Added new window around map position (0.9, -0.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0370638 Jy ! Total flux subtracted in 50 components = 0.0370638 Jy ! Clean residual min=-0.017010 max=0.024639 Jy/beam ! Clean residual mean=-0.000064 rms=0.004544 Jy/beam ! Combined flux in latest and established models = 1.15627 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 26 model components to the UV plane model. ! The established model now contains 86 components and 1.15627 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.305488Jy sigma=3.709082 ! Fit after self-cal, rms=0.305291Jy sigma=3.707249 ! Inverting map ! Added new window around map position (-3, 0.5). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0321519 Jy ! Total flux subtracted in 50 components = 0.0321519 Jy ! Clean residual min=-0.017030 max=0.020809 Jy/beam ! Clean residual mean=-0.000036 rms=0.004402 Jy/beam ! Combined flux in latest and established models = 1.18842 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 97 components and 1.18842 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.299216Jy sigma=3.581502 ! Fit after self-cal, rms=0.299070Jy sigma=3.579798 ! 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.7579 mas, bmaj=1.576 mas, bpa=-0.897 degrees ! Estimated noise=0.388524 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) ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0621709 Jy ! Component: 100 - total flux cleaned = 0.103836 Jy ! Total flux subtracted in 100 components = 0.103836 Jy ! Clean residual min=-0.021955 max=0.030554 Jy/beam ! Clean residual mean=0.002111 rms=0.006688 Jy/beam ! Combined flux in latest and established models = 1.29226 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 64 model components to the UV plane model. ! The established model now contains 152 components and 1.29226 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.282675Jy sigma=3.280926 ! Fit after self-cal, rms=0.282391Jy sigma=3.272278 ! Inverting map ! Added new window around map position (-7.2, 3.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0355797 Jy ! Component: 100 - total flux cleaned = 0.0627606 Jy ! Total flux subtracted in 100 components = 0.0627606 Jy ! Clean residual min=-0.021436 max=0.026694 Jy/beam ! Clean residual mean=0.002020 rms=0.006144 Jy/beam ! Combined flux in latest and established models = 1.35502 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 55 model components to the UV plane model. ! The established model now contains 189 components and 1.35502 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.275876Jy sigma=3.177059 ! Fit after self-cal, rms=0.275686Jy sigma=3.169539 ! 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 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.88 FD 1.04 GC 1.14 HH 1.00* ! HN 0.99 KK 0.90 KP 0.99 LA 0.99 ! MA 1.00* MC 1.00* MK 1.11 NL 1.10 ! NY 1.14 ON 1.00* OV 1.05 PT 1.00 ! SC 1.06 TS 0.93 WF 0.99 WZ 1.00* ! ! A total of 222 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 0.92 GC 1.01 HH 1.00* ! HN 1.00* KK 1.00 KP 1.02 LA 1.10 ! MA 1.00* MC 1.00* MK 1.03 NL 1.00* ! NY 1.00* ON 1.00* OV 1.12 PT 0.94 ! SC 1.00* TS 0.96 WF 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 0.84 MC 0.79 MK 1.00* NL 1.00* ! NY 1.27 ON 0.92 OV 1.00* PT 1.00* ! SC 1.00* TS 0.95 WF 1.00* WZ 0.91 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.90 FD 1.02 GC 1.10 HH 1.00* ! HN 0.92 KK 1.16 KP 1.02 LA 0.89 ! MA 1.00* MC 1.00* MK 1.01 NL 1.01 ! NY 0.62 ON 1.00* OV 1.07 PT 0.97 ! SC 1.03 TS 1.00* WF 0.86 WZ 1.00* ! ! ! Correcting IF 2. ! A total of 21 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.74 FD 0.97 GC 1.07 HH 1.00* ! HN 0.91 KK 0.96 KP 0.97 LA 0.81 ! MA 1.00* MC 1.00* MK 0.90 NL 1.19 ! NY 1.06 ON 1.00* OV 1.00 PT 1.06 ! SC 1.15 TS 0.98 WF 0.89 WZ 1.00* ! ! A total of 213 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.05 GC 1.08 HH 1.00* ! HN 1.00* KK 1.03 KP 1.05 LA 1.10 ! MA 1.00* MC 1.00* MK 1.00 NL 1.00* ! NY 1.00* ON 1.00* OV 1.29 PT 0.95 ! SC 1.00* TS 0.96 WF 1.00* WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 0.87 MC 0.50 MK 1.00* NL 1.00* ! NY 1.25 ON 0.86 OV 1.00* PT 1.00* ! SC 1.00* TS 0.96 WF 1.00* WZ 0.99 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.81 FD 1.03 GC 0.91 HH 1.00* ! HN 0.97 KK 1.21 KP 0.98 LA 0.87 ! MA 1.00* MC 1.00* MK 0.94 NL 0.99 ! NY 0.52 ON 1.00* OV 1.06 PT 0.96 ! SC 1.02 TS 1.00* WF 0.80 WZ 1.00* ! ! ! Correcting IF 3. ! A total of 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.62 FD 1.03 GC 0.94 HH 1.00* ! HN 0.95 KK 1.03 KP 0.86 LA 1.01 ! MA 1.00* MC 1.00* MK 0.95 NL 1.15 ! NY 1.15 ON 1.00* OV 0.98 PT 1.01 ! SC 0.98 TS 0.92 WF 1.07 WZ 1.00* ! ! A total of 222 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.13 GC 0.93 HH 1.00* ! HN 1.00* KK 1.08 KP 0.99 LA 1.11 ! MA 1.00* MC 1.00* MK 1.06 NL 1.00* ! NY 1.00* ON 1.00* OV 1.26 PT 0.94 ! SC 1.00* TS 0.92 WF 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 1.05 MC 0.67 MK 1.00* NL 1.00* ! NY 1.24 ON 0.96 OV 1.00* PT 1.00* ! SC 1.00* TS 0.95 WF 1.00* WZ 0.83 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.87 FD 1.05 GC 1.06 HH 1.00* ! HN 1.00 KK 1.12 KP 0.98 LA 0.88 ! MA 1.00* MC 1.00* MK 1.04 NL 0.98 ! NY 0.59 ON 1.00* OV 1.05 PT 0.81 ! SC 1.27 TS 1.00* WF 0.76 WZ 1.00* ! ! ! Correcting IF 4. ! A total of 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.75 FD 0.93 GC 1.09 HH 1.00* ! HN 0.97 KK 1.01 KP 0.98 LA 1.00 ! MA 1.00* MC 1.00* MK 0.97 NL 1.16 ! NY 0.86 ON 1.00* OV 0.94 PT 0.97 ! SC 0.94 TS 1.00 WF 0.87 WZ 1.00* ! ! A total of 219 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.11 GC 1.00 HH 1.00* ! HN 1.00* KK 1.07 KP 1.04 LA 1.07 ! MA 1.00* MC 1.00* MK 1.02 NL 1.00* ! NY 1.00* ON 1.00* OV 1.22 PT 0.94 ! SC 1.00* TS 0.95 WF 1.00* WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 1.05 MC 0.58 MK 1.00* NL 1.00* ! NY 1.22 ON 0.86 OV 1.00* PT 1.00* ! SC 1.00* TS 0.88 WF 1.00* WZ 0.95 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.77 FD 1.04 GC 1.07 HH 1.00* ! HN 0.94 KK 1.11 KP 0.96 LA 0.87 ! MA 1.00* MC 1.00* MK 0.68 NL 1.03 ! NY 0.78 ON 1.00* OV 1.07 PT 0.97 ! SC 1.14 TS 1.00* WF 0.78 WZ 1.00* ! ! ! Fit before self-cal, rms=0.275686Jy sigma=3.169539 ! Fit after self-cal, rms=0.199420Jy sigma=2.444055 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.7793 mas, bmaj=1.564 mas, bpa=-3.105 degrees ! Estimated noise=0.394981 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 (-11, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0250868 Jy ! Component: 100 - total flux cleaned = 0.0429434 Jy ! Total flux subtracted in 100 components = 0.0429434 Jy ! Clean residual min=-0.017608 max=0.022672 Jy/beam ! Clean residual mean=0.001620 rms=0.005000 Jy/beam ! Combined flux in latest and established models = 1.39796 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 45 model components to the UV plane model. ! The established model now contains 234 components and 1.39796 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.196136Jy sigma=2.377061 ! Fit after self-cal, rms=0.195781Jy sigma=2.369320 ! Inverting map ! Added new window around map position (-14.7, 6.7). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0234922 Jy ! Component: 100 - total flux cleaned = 0.0400278 Jy ! Total flux subtracted in 100 components = 0.0400278 Jy ! Clean residual min=-0.016072 max=0.016144 Jy/beam ! Clean residual mean=0.001512 rms=0.004629 Jy/beam ! Combined flux in latest and established models = 1.43799 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 50 model components to the UV plane model. ! The established model now contains 275 components and 1.43799 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192692Jy sigma=2.305384 ! Fit after self-cal, rms=0.192580Jy sigma=2.299725 ! 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) ! Added new window around map position (-1.6, 3.1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0172443 Jy ! Component: 100 - total flux cleaned = 0.0288089 Jy ! Total flux subtracted in 100 components = 0.0288089 Jy ! Clean residual min=-0.015207 max=0.015554 Jy/beam ! Clean residual mean=0.001444 rms=0.004481 Jy/beam ! Combined flux in latest and established models = 1.4668 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 47 model components to the UV plane model. ! The established model now contains 302 components and 1.4668 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190635Jy sigma=2.264554 ! Fit after self-cal, rms=0.175256Jy sigma=1.626666 ! Inverting map and beam ! Estimated beam: bmin=0.7812 mas, bmaj=1.565 mas, bpa=-2.845 degrees ! Estimated noise=0.396227 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0166638 Jy ! Component: 100 - total flux cleaned = 0.0256036 Jy ! Total flux subtracted in 100 components = 0.0256036 Jy ! Clean residual min=-0.006469 max=0.011234 Jy/beam ! Clean residual mean=0.001141 rms=0.002187 Jy/beam ! Combined flux in latest and established models = 1.4924 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 21 model components to the UV plane model. ! The established model now contains 320 components and 1.4924 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.173177Jy sigma=1.599347 ! Fit after self-cal, rms=0.171613Jy sigma=1.576388 ! Inverting map and beam ! Estimated beam: bmin=0.7793 mas, bmaj=1.564 mas, bpa=-2.779 degrees ! Estimated noise=0.396367 mJy/beam. ! Added new window around map position (-1.6, 0.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0125129 Jy ! Component: 100 - total flux cleaned = 0.0184098 Jy ! Total flux subtracted in 100 components = 0.0184098 Jy ! Clean residual min=-0.006099 max=0.011210 Jy/beam ! Clean residual mean=0.000967 rms=0.002070 Jy/beam ! Combined flux in latest and established models = 1.51081 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 32 model components to the UV plane model. ! The established model now contains 335 components and 1.51081 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170336Jy sigma=1.560896 ! Fit after self-cal, rms=0.171335Jy sigma=1.551569 ! Inverting map and beam ! Estimated beam: bmin=0.7784 mas, bmaj=1.563 mas, bpa=-2.753 degrees ! Estimated noise=0.396572 mJy/beam. ! Added new window around map position (-5, 2.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0126911 Jy ! Component: 100 - total flux cleaned = 0.0211981 Jy ! Total flux subtracted in 100 components = 0.0211981 Jy ! Clean residual min=-0.006322 max=0.008905 Jy/beam ! Clean residual mean=0.000830 rms=0.001930 Jy/beam ! Combined flux in latest and established models = 1.53201 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 356 components and 1.53201 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170136Jy sigma=1.534692 ! Fit after self-cal, rms=0.171901Jy sigma=1.528038 ! Inverting map and beam ! Estimated beam: bmin=0.7786 mas, bmaj=1.562 mas, bpa=-2.749 degrees ! Estimated noise=0.397093 mJy/beam. ! Added new window around map position (-10.2, 4.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0102776 Jy ! Component: 100 - total flux cleaned = 0.0180806 Jy ! Total flux subtracted in 100 components = 0.0180806 Jy ! Clean residual min=-0.005795 max=0.007913 Jy/beam ! Clean residual mean=0.000760 rms=0.001796 Jy/beam ! Combined flux in latest and established models = 1.55009 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 38 model components to the UV plane model. ! The established model now contains 379 components and 1.55009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171048Jy sigma=1.514847 ! Fit after self-cal, rms=0.172387Jy sigma=1.511585 ! Inverting map and beam ! Estimated beam: bmin=0.7787 mas, bmaj=1.562 mas, bpa=-2.748 degrees ! Estimated noise=0.397544 mJy/beam. ! Added new window around map position (-6.6, 0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00893288 Jy ! Component: 100 - total flux cleaned = 0.014989 Jy ! Total flux subtracted in 100 components = 0.014989 Jy ! Clean residual min=-0.005794 max=0.007260 Jy/beam ! Clean residual mean=0.000708 rms=0.001718 Jy/beam ! Combined flux in latest and established models = 1.56508 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 42 model components to the UV plane model. ! The established model now contains 396 components and 1.56508 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171715Jy sigma=1.501768 ! Fit after self-cal, rms=0.172917Jy sigma=1.499707 ! Inverting map and beam ! Estimated beam: bmin=0.7785 mas, bmaj=1.561 mas, bpa=-2.736 degrees ! Estimated noise=0.397822 mJy/beam. ! Added new window around map position (-17.8, 8.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00804131 Jy ! Component: 100 - total flux cleaned = 0.0134395 Jy ! Total flux subtracted in 100 components = 0.0134395 Jy ! Clean residual min=-0.005753 max=0.006248 Jy/beam ! Clean residual mean=0.000650 rms=0.001646 Jy/beam ! Combined flux in latest and established models = 1.57852 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 43 model components to the UV plane model. ! The established model now contains 418 components and 1.57852 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.172381Jy sigma=1.491403 ! Fit after self-cal, rms=0.173403Jy sigma=1.490041 ! Inverting map and beam ! Estimated beam: bmin=0.7781 mas, bmaj=1.561 mas, bpa=-2.74 degrees ! Estimated noise=0.398051 mJy/beam. ! Added new window around map position (-20.4, 8.7). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00752889 Jy ! Component: 100 - total flux cleaned = 0.0131845 Jy ! Total flux subtracted in 100 components = 0.0131845 Jy ! Clean residual min=-0.005699 max=0.006173 Jy/beam ! Clean residual mean=0.000599 rms=0.001574 Jy/beam ! Combined flux in latest and established models = 1.5917 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 50 model components to the UV plane model. ! The established model now contains 444 components and 1.5917 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.172932Jy sigma=1.482791 ! Fit after self-cal, rms=0.173991Jy sigma=1.481522 ! Inverting map and beam ! Estimated beam: bmin=0.7779 mas, bmaj=1.56 mas, bpa=-2.738 degrees ! Estimated noise=0.398224 mJy/beam. ! Added new window around map position (15.1, 10.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00688417 Jy ! Component: 100 - total flux cleaned = 0.0119606 Jy ! Total flux subtracted in 100 components = 0.0119606 Jy ! Clean residual min=-0.005574 max=0.005691 Jy/beam ! Clean residual mean=0.000560 rms=0.001513 Jy/beam ! Combined flux in latest and established models = 1.60366 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 47 model components to the UV plane model. ! The established model now contains 455 components and 1.60366 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.173584Jy sigma=1.474881 ! Fit after self-cal, rms=0.174359Jy sigma=1.473946 ! Inverting map and beam ! Estimated beam: bmin=0.7778 mas, bmaj=1.56 mas, bpa=-2.742 degrees ! Estimated noise=0.398375 mJy/beam. ! Added new window around map position (-4.2, -0.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00587752 Jy ! Component: 100 - total flux cleaned = 0.0095317 Jy ! Total flux subtracted in 100 components = 0.0095317 Jy ! Clean residual min=-0.005475 max=0.005521 Jy/beam ! Clean residual mean=0.000526 rms=0.001466 Jy/beam ! Combined flux in latest and established models = 1.6132 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 46 model components to the UV plane model. ! The established model now contains 469 components and 1.6132 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174025Jy sigma=1.468203 ! Fit after self-cal, rms=0.174662Jy sigma=1.467393 ! Inverting map and beam ! Estimated beam: bmin=0.7778 mas, bmaj=1.559 mas, bpa=-2.749 degrees ! Estimated noise=0.398548 mJy/beam. ! Added new window around map position (-20.7, 4.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00588056 Jy ! Component: 100 - total flux cleaned = 0.00912725 Jy ! Total flux subtracted in 100 components = 0.00912725 Jy ! Clean residual min=-0.005445 max=0.005314 Jy/beam ! Clean residual mean=0.000496 rms=0.001408 Jy/beam ! Combined flux in latest and established models = 1.62232 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 42 model components to the UV plane model. ! The established model now contains 485 components and 1.62232 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174360Jy sigma=1.462133 ! Fit after self-cal, rms=0.174910Jy sigma=1.461421 ! Inverting map and beam ! Estimated beam: bmin=0.7777 mas, bmaj=1.559 mas, bpa=-2.749 degrees ! Estimated noise=0.398669 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.00445947 Jy ! Component: 100 - total flux cleaned = 0.00710766 Jy ! Total flux subtracted in 100 components = 0.00710766 Jy ! Clean residual min=-0.005374 max=0.005249 Jy/beam ! Clean residual mean=0.000475 rms=0.001373 Jy/beam ! Combined flux in latest and established models = 1.62943 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 56 model components to the UV plane model. ! The established model now contains 502 components and 1.62943 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174673Jy sigma=1.457149 ! Fit after self-cal, rms=0.174050Jy sigma=1.449127 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 3C274_X ! Writing UV FITS file: 3C274_X.uvf ! Writing 502 model components to file: 3C274_X.mod ! wwins: Wrote 19 windows to 3C274_X.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7777 x 1.559 at -2.749 degrees (North through East) ! Clean map min=-0.0053403 max=0.75324 Jy/beam ! Writing clean map to FITS file: 3C274_X.fits ! Writing difmap environment to: 3C274_X.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 3C274_X] device 3C274_Xuv.ps/vps,1,2 ! Attempting to open device: '3C274_Xuv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 3C274_Xim.ps/vps ! Attempting to open device: '3C274_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.7777 mas, bmaj=1.559 mas, bpa=-2.749 degrees ! Estimated noise=0.398669 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7777 x 1.559 at -2.749 degrees (North through East) ! Clean map min=-0.0053403 max=0.75324 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo5-0.1] quit ! Quitting program ! Log file 3C274_X.log closed on Fri Jul 11 14:30:26 2008