! Started logfile: 0552+398_X.log on Fri Jun 27 12:14:55 2008 observe 0552+398_X.edt3 ! Reading UV FITS file: 0552+398_X.edt3 ! AN table 1: 46 integrations on 190 of 190 possible baselines. ! AN table 2: 14 integrations on 190 of 190 possible baselines. ! AN table 3: 118 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.25958 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 0552+398 ! ! 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 1480 lines of history. ! ! Reading 35116 visibilities. ![@gonzo5-0.1 0552+398_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 696 telescope corrections were flagged in sub-array 1. ! A total of 268 telescope corrections were flagged in sub-array 2. ! A total of 1654 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 695 telescope corrections were flagged in sub-array 1. ! A total of 268 telescope corrections were flagged in sub-array 2. ! A total of 1647 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 696 telescope corrections were flagged in sub-array 1. ! A total of 268 telescope corrections were flagged in sub-array 2. ! A total of 1651 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 694 telescope corrections were flagged in sub-array 1. ! A total of 268 telescope corrections were flagged in sub-array 2. ! A total of 1663 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=2.558586Jy sigma=35.129851 ! Fit after self-cal, rms=2.560107Jy sigma=35.205420 ! 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.5545 mas, bmaj=0.8207 mas, bpa=0.1233 degrees ! Estimated noise=0.855636 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 = 2.39045 Jy ! Total flux subtracted in 50 components = 2.39045 Jy ! Clean residual min=-0.087949 max=0.674439 Jy/beam ! Clean residual mean=0.000101 rms=0.026754 Jy/beam ! Combined flux in latest and established models = 2.39045 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 2 model components to the UV plane model. ! The established model now contains 2 components and 2.39045 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.354087Jy sigma=18.399407 ! Fit after self-cal, rms=1.353845Jy sigma=18.398482 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.700109 Jy ! Total flux subtracted in 50 components = 0.700109 Jy ! Clean residual min=-0.057710 max=0.310861 Jy/beam ! Clean residual mean=0.000105 rms=0.016656 Jy/beam ! Combined flux in latest and established models = 3.09056 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 8 components and 3.09056 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.908677Jy sigma=11.592188 ! Fit after self-cal, rms=0.905329Jy sigma=11.569913 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.33935 Jy ! Total flux subtracted in 50 components = 0.33935 Jy ! Clean residual min=-0.052942 max=0.209783 Jy/beam ! Clean residual mean=0.000105 rms=0.013249 Jy/beam ! Combined flux in latest and established models = 3.42991 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 9 model components to the UV plane model. ! The established model now contains 13 components and 3.42991 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.739596Jy sigma=8.842403 ! Fit after self-cal, rms=0.737650Jy sigma=8.828328 ! Inverting map ! Added new window around map position (-0.7, 0.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.218133 Jy ! Total flux subtracted in 50 components = 0.218133 Jy ! Clean residual min=-0.051018 max=0.113769 Jy/beam ! Clean residual mean=0.000098 rms=0.011549 Jy/beam ! Combined flux in latest and established models = 3.64804 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 10 model components to the UV plane model. ! The established model now contains 21 components and 3.64804 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.657109Jy sigma=7.399625 ! Fit after self-cal, rms=0.653530Jy sigma=7.367416 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.150408 Jy ! Total flux subtracted in 50 components = 0.150408 Jy ! Clean residual min=-0.047412 max=0.083594 Jy/beam ! Clean residual mean=0.000094 rms=0.010627 Jy/beam ! Combined flux in latest and established models = 3.79845 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 32 components and 3.79845 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.608524Jy sigma=6.560317 ! Fit after self-cal, rms=0.607425Jy sigma=6.550014 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.111642 Jy ! Total flux subtracted in 50 components = 0.111642 Jy ! Clean residual min=-0.046462 max=0.063177 Jy/beam ! Clean residual mean=0.000090 rms=0.010093 Jy/beam ! Combined flux in latest and established models = 3.91009 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 14 model components to the UV plane model. ! The established model now contains 40 components and 3.91009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.580730Jy sigma=6.070157 ! Fit after self-cal, rms=0.580152Jy sigma=6.063505 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0853516 Jy ! Total flux subtracted in 50 components = 0.0853516 Jy ! Clean residual min=-0.045804 max=0.048894 Jy/beam ! Clean residual mean=0.000087 rms=0.009759 Jy/beam ! Combined flux in latest and established models = 3.99544 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 12 model components to the UV plane model. ! The established model now contains 43 components and 3.99545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.563703Jy sigma=5.776880 ! Fit after self-cal, rms=0.563409Jy sigma=5.772881 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0659253 Jy ! Total flux subtracted in 50 components = 0.0659253 Jy ! Clean residual min=-0.045341 max=0.043546 Jy/beam ! Clean residual mean=0.000084 rms=0.009554 Jy/beam ! Combined flux in latest and established models = 4.06137 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 49 components and 4.06137 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.553132Jy sigma=5.605299 ! Fit after self-cal, rms=0.552970Jy sigma=5.602599 ! 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.6689 mas, bmaj=0.9872 mas, bpa=-1.132 degrees ! Estimated noise=0.52506 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.0806026 Jy ! Component: 100 - total flux cleaned = 0.128235 Jy ! Total flux subtracted in 100 components = 0.128235 Jy ! Clean residual min=-0.054962 max=0.038692 Jy/beam ! Clean residual mean=0.000145 rms=0.009871 Jy/beam ! Combined flux in latest and established models = 4.18961 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 18 model components to the UV plane model. ! The established model now contains 58 components and 4.18961 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.538590Jy sigma=5.426859 ! Fit after self-cal, rms=0.538302Jy sigma=5.419010 ! 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 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.12 FD 1.07 GC 1.06 HH 1.00* ! HN 1.07 KK 0.94 KP 1.03 LA 1.07 ! MA 1.00* MC 1.00* MK 1.02 NL 0.98 ! NY 1.11 ON 1.00* OV 0.99 PT 1.03 ! SC 1.00* TS 1.03 WF 1.04 WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.03 GC 1.07 HH 1.00* ! HN 0.99 KK 0.97 KP 1.04 LA 1.04 ! MA 1.00* MC 1.00* MK 1.04 NL 0.99 ! NY 0.81 ON 1.00* OV 0.96 PT 1.00 ! SC 1.00 TS 0.94 WF 0.84 WZ 1.00* ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.99 GC 1.03 HH 1.00* ! HN 1.00 KK 1.07 KP 1.00 LA 0.88 ! MA 1.00* MC 1.00* MK 0.87 NL 1.07 ! NY 1.04 ON 1.00* OV 0.95 PT 1.08 ! SC 1.00* TS 1.01 WF 0.86 WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.91 FD 1.07 GC 1.09 HH 1.00* ! HN 0.93 KK 0.71 KP 1.06 LA 1.02 ! MA 1.00* MC 1.00* MK 0.96 NL 0.96 ! NY 1.03 ON 1.00* OV 0.95 PT 0.97 ! SC 1.12 TS 1.01 WF 0.76 WZ 1.00* ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.86 FD 1.05 GC 0.92 HH 1.00* ! HN 1.09 KK 1.02 KP 0.91 LA 1.06 ! MA 1.00* MC 1.00* MK 0.88 NL 1.02 ! NY 1.08 ON 1.00* OV 0.94 PT 1.05 ! SC 1.00* TS 0.98 WF 0.89 WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.02 FD 1.01 GC 1.12 HH 1.00* ! HN 0.95 KK 0.80 KP 1.09 LA 0.96 ! MA 1.00* MC 1.00* MK 1.05 NL 0.97 ! NY 0.99 ON 1.00* OV 0.93 PT 1.04 ! SC 0.94 TS 0.94 WF 0.79 WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 FD 1.03 GC 1.05 HH 1.00* ! HN 1.04 KK 1.10 KP 0.99 LA 1.03 ! MA 1.00* MC 1.00* MK 0.92 NL 1.04 ! NY 0.85 ON 1.00* OV 0.91 PT 1.08 ! SC 1.00* TS 1.01 WF 0.92 WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! 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* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.01 FD 1.02 GC 1.09 HH 1.00* ! HN 0.97 KK 0.82 KP 1.03 LA 1.02 ! MA 1.00* MC 1.00* MK 1.03 NL 1.01 ! NY 0.77 ON 1.00* OV 0.84 PT 0.92 ! SC 1.09 TS 0.95 WF 0.82 WZ 1.00* ! ! ! Fit before self-cal, rms=0.538302Jy sigma=5.419010 ! Fit after self-cal, rms=0.330124Jy sigma=3.208144 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6648 mas, bmaj=0.9888 mas, bpa=-1.582 degrees ! Estimated noise=0.519608 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.0412879 Jy ! Component: 100 - total flux cleaned = 0.0492936 Jy ! Total flux subtracted in 100 components = 0.0492936 Jy ! Clean residual min=-0.028802 max=0.028806 Jy/beam ! Clean residual mean=0.000134 rms=0.005204 Jy/beam ! Combined flux in latest and established models = 4.2389 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 64 components and 4.2389 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323630Jy sigma=3.116376 ! Fit after self-cal, rms=0.323438Jy sigma=3.110864 ! Inverting map ! Added new window around map position (1.7, 4.1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0193992 Jy ! Component: 100 - total flux cleaned = 0.0194027 Jy ! Total flux subtracted in 100 components = 0.0194027 Jy ! Clean residual min=-0.027473 max=0.024564 Jy/beam ! Clean residual mean=0.000125 rms=0.005122 Jy/beam ! Combined flux in latest and established models = 4.2583 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 14 model components to the UV plane model. ! The established model now contains 73 components and 4.2583 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.320332Jy sigma=3.073403 ! Fit after self-cal, rms=0.320295Jy sigma=3.071039 ! 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.00332725 Jy ! Component: 100 - total flux cleaned = 0.00201818 Jy ! Total flux subtracted in 100 components = 0.00201818 Jy ! Clean residual min=-0.027844 max=0.024658 Jy/beam ! Clean residual mean=0.000123 rms=0.005077 Jy/beam ! Combined flux in latest and established models = 4.26032 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 81 components and 4.26032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.318450Jy sigma=3.055860 ! Fit after self-cal, rms=0.230780Jy sigma=1.949401 ! Inverting map and beam ! Estimated beam: bmin=0.6614 mas, bmaj=0.9856 mas, bpa=-1.631 degrees ! Estimated noise=0.519345 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0225195 Jy ! Component: 100 - total flux cleaned = 0.0245156 Jy ! Total flux subtracted in 100 components = 0.0245156 Jy ! Clean residual min=-0.014109 max=0.018688 Jy/beam ! Clean residual mean=0.000054 rms=0.002726 Jy/beam ! Combined flux in latest and established models = 4.28484 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 11 model components to the UV plane model. ! The established model now contains 87 components and 4.28484 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227332Jy sigma=1.913008 ! Fit after self-cal, rms=0.224515Jy sigma=1.874514 ! Inverting map and beam ! Estimated beam: bmin=0.66 mas, bmaj=0.9845 mas, bpa=-1.723 degrees ! Estimated noise=0.519563 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0206215 Jy ! Component: 100 - total flux cleaned = 0.0230257 Jy ! Total flux subtracted in 100 components = 0.0230257 Jy ! Clean residual min=-0.013912 max=0.016954 Jy/beam ! Clean residual mean=0.000054 rms=0.002543 Jy/beam ! Combined flux in latest and established models = 4.30786 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 16 model components to the UV plane model. ! The established model now contains 93 components and 4.30786 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.221930Jy sigma=1.845406 ! Fit after self-cal, rms=0.220157Jy sigma=1.825934 ! Inverting map and beam ! Estimated beam: bmin=0.6587 mas, bmaj=0.9836 mas, bpa=-1.799 degrees ! Estimated noise=0.519733 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0178215 Jy ! Component: 100 - total flux cleaned = 0.0212804 Jy ! Total flux subtracted in 100 components = 0.0212804 Jy ! Clean residual min=-0.013710 max=0.015539 Jy/beam ! Clean residual mean=0.000052 rms=0.002416 Jy/beam ! Combined flux in latest and established models = 4.32914 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 16 model components to the UV plane model. ! The established model now contains 94 components and 4.32914 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.218046Jy sigma=1.800515 ! Fit after self-cal, rms=0.216614Jy sigma=1.787219 ! Inverting map and beam ! Estimated beam: bmin=0.6575 mas, bmaj=0.9826 mas, bpa=-1.852 degrees ! Estimated noise=0.519867 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0158095 Jy ! Component: 100 - total flux cleaned = 0.0201654 Jy ! Total flux subtracted in 100 components = 0.0201654 Jy ! Clean residual min=-0.013686 max=0.014395 Jy/beam ! Clean residual mean=0.000049 rms=0.002311 Jy/beam ! Combined flux in latest and established models = 4.34931 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 16 model components to the UV plane model. ! The established model now contains 97 components and 4.34931 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.214801Jy sigma=1.764673 ! Fit after self-cal, rms=0.213515Jy sigma=1.753919 ! Inverting map and beam ! Estimated beam: bmin=0.6563 mas, bmaj=0.9816 mas, bpa=-1.889 degrees ! Estimated noise=0.520009 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.01475 Jy ! Component: 100 - total flux cleaned = 0.0188791 Jy ! Total flux subtracted in 100 components = 0.0188791 Jy ! Clean residual min=-0.013501 max=0.013385 Jy/beam ! Clean residual mean=0.000045 rms=0.002219 Jy/beam ! Combined flux in latest and established models = 4.36819 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 15 model components to the UV plane model. ! The established model now contains 97 components and 4.36819 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.211920Jy sigma=1.733413 ! Fit after self-cal, rms=0.210719Jy sigma=1.724465 ! Inverting map and beam ! Estimated beam: bmin=0.6553 mas, bmaj=0.9806 mas, bpa=-1.915 degrees ! Estimated noise=0.520101 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0138254 Jy ! Component: 100 - total flux cleaned = 0.01775 Jy ! Total flux subtracted in 100 components = 0.01775 Jy ! Clean residual min=-0.013408 max=0.012596 Jy/beam ! Clean residual mean=0.000042 rms=0.002134 Jy/beam ! Combined flux in latest and established models = 4.38594 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 19 model components to the UV plane model. ! The established model now contains 100 components and 4.38594 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.209307Jy sigma=1.705922 ! Fit after self-cal, rms=0.208202Jy sigma=1.698233 ! Inverting map and beam ! Estimated beam: bmin=0.6544 mas, bmaj=0.9796 mas, bpa=-1.934 degrees ! Estimated noise=0.520213 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0130382 Jy ! Component: 100 - total flux cleaned = 0.0163372 Jy ! Total flux subtracted in 100 components = 0.0163372 Jy ! Clean residual min=-0.013189 max=0.011927 Jy/beam ! Clean residual mean=0.000038 rms=0.002057 Jy/beam ! Combined flux in latest and established models = 4.40227 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 102 components and 4.40227 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206927Jy sigma=1.681353 ! Fit after self-cal, rms=0.205919Jy sigma=1.674633 ! Inverting map and beam ! Estimated beam: bmin=0.6536 mas, bmaj=0.9786 mas, bpa=-1.952 degrees ! Estimated noise=0.520331 mJy/beam. ! Added new window around map position (-1.5, 2.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0145388 Jy ! Component: 100 - total flux cleaned = 0.0182394 Jy ! Total flux subtracted in 100 components = 0.0182394 Jy ! Clean residual min=-0.012446 max=0.009998 Jy/beam ! Clean residual mean=0.000034 rms=0.001975 Jy/beam ! Combined flux in latest and established models = 4.42051 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 106 components and 4.42051 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.204679Jy sigma=1.658245 ! Fit after self-cal, rms=0.203768Jy sigma=1.652121 ! Inverting map and beam ! Estimated beam: bmin=0.6529 mas, bmaj=0.9777 mas, bpa=-1.964 degrees ! Estimated noise=0.520448 mJy/beam. ! Added new window around map position (-1.4, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00809248 Jy ! Component: 100 - total flux cleaned = 0.0119958 Jy ! Total flux subtracted in 100 components = 0.0119958 Jy ! Clean residual min=-0.009211 max=0.009446 Jy/beam ! Clean residual mean=0.000035 rms=0.001892 Jy/beam ! Combined flux in latest and established models = 4.43251 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 111 components and 4.43251 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.202418Jy sigma=1.633724 ! Fit after self-cal, rms=0.201723Jy sigma=1.629626 ! Inverting map and beam ! Estimated beam: bmin=0.6524 mas, bmaj=0.9771 mas, bpa=-1.975 degrees ! Estimated noise=0.520692 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00989922 Jy ! Component: 100 - total flux cleaned = 0.0116496 Jy ! Total flux subtracted in 100 components = 0.0116496 Jy ! Clean residual min=-0.009052 max=0.009168 Jy/beam ! Clean residual mean=0.000032 rms=0.001832 Jy/beam ! Combined flux in latest and established models = 4.44416 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 112 components and 4.44416 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.200658Jy sigma=1.615414 ! Fit after self-cal, rms=0.199980Jy sigma=1.611676 ! Inverting map and beam ! Estimated beam: bmin=0.6518 mas, bmaj=0.9765 mas, bpa=-1.983 degrees ! Estimated noise=0.520846 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00983356 Jy ! Component: 100 - total flux cleaned = 0.0110712 Jy ! Total flux subtracted in 100 components = 0.0110712 Jy ! Clean residual min=-0.008836 max=0.008845 Jy/beam ! Clean residual mean=0.000029 rms=0.001779 Jy/beam ! Combined flux in latest and established models = 4.45523 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 19 model components to the UV plane model. ! The established model now contains 113 components and 4.45523 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199037Jy sigma=1.599142 ! Fit after self-cal, rms=0.198368Jy sigma=1.595647 ! Inverting map and beam ! Estimated beam: bmin=0.6513 mas, bmaj=0.9758 mas, bpa=-1.986 degrees ! Estimated noise=0.520963 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00976695 Jy ! Component: 100 - total flux cleaned = 0.0105466 Jy ! Total flux subtracted in 100 components = 0.0105466 Jy ! Clean residual min=-0.008629 max=0.008541 Jy/beam ! Clean residual mean=0.000027 rms=0.001727 Jy/beam ! Combined flux in latest and established models = 4.46578 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 115 components and 4.46578 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.197525Jy sigma=1.584090 ! Fit after self-cal, rms=0.196879Jy sigma=1.580875 ! Inverting map and beam ! Estimated beam: bmin=0.6507 mas, bmaj=0.9751 mas, bpa=-1.985 degrees ! Estimated noise=0.521054 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00893418 Jy ! Component: 100 - total flux cleaned = 0.0100705 Jy ! Total flux subtracted in 100 components = 0.0100705 Jy ! Clean residual min=-0.008501 max=0.008268 Jy/beam ! Clean residual mean=0.000024 rms=0.001680 Jy/beam ! Combined flux in latest and established models = 4.47585 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 22 model components to the UV plane model. ! The established model now contains 117 components and 4.47585 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.196112Jy sigma=1.570166 ! Fit after self-cal, rms=0.195501Jy sigma=1.567201 ! Inverting map and beam ! Estimated beam: bmin=0.6502 mas, bmaj=0.9744 mas, bpa=-1.986 degrees ! Estimated noise=0.521109 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00853477 Jy ! Component: 100 - total flux cleaned = 0.00963176 Jy ! Total flux subtracted in 100 components = 0.00963176 Jy ! Clean residual min=-0.008417 max=0.007998 Jy/beam ! Clean residual mean=0.000022 rms=0.001635 Jy/beam ! Combined flux in latest and established models = 4.48548 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 25 model components to the UV plane model. ! The established model now contains 118 components and 4.48548 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.194802Jy sigma=1.557205 ! Fit after self-cal, rms=0.194218Jy sigma=1.554445 ! Inverting map and beam ! Estimated beam: bmin=0.6497 mas, bmaj=0.9738 mas, bpa=-1.986 degrees ! Estimated noise=0.521141 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00778523 Jy ! Component: 100 - total flux cleaned = 0.00919536 Jy ! Total flux subtracted in 100 components = 0.00919536 Jy ! Clean residual min=-0.008266 max=0.007742 Jy/beam ! Clean residual mean=0.000020 rms=0.001591 Jy/beam ! Combined flux in latest and established models = 4.49467 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 121 components and 4.49467 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.193572Jy sigma=1.545288 ! Fit after self-cal, rms=0.193018Jy sigma=1.542678 ! Inverting map and beam ! Estimated beam: bmin=0.6492 mas, bmaj=0.9731 mas, bpa=-1.986 degrees ! Estimated noise=0.521168 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00745228 Jy ! Component: 100 - total flux cleaned = 0.00880682 Jy ! Total flux subtracted in 100 components = 0.00880682 Jy ! Clean residual min=-0.008095 max=0.007524 Jy/beam ! Clean residual mean=0.000018 rms=0.001550 Jy/beam ! Combined flux in latest and established models = 4.50348 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 123 components and 4.50348 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192417Jy sigma=1.534138 ! Fit after self-cal, rms=0.191914Jy sigma=1.531841 ! Inverting map and beam ! Estimated beam: bmin=0.6488 mas, bmaj=0.9725 mas, bpa=-1.987 degrees ! Estimated noise=0.521179 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00679052 Jy ! Component: 100 - total flux cleaned = 0.00841984 Jy ! Total flux subtracted in 100 components = 0.00841984 Jy ! Clean residual min=-0.007962 max=0.007309 Jy/beam ! Clean residual mean=0.000016 rms=0.001512 Jy/beam ! Combined flux in latest and established models = 4.5119 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 24 model components to the UV plane model. ! The established model now contains 124 components and 4.5119 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.191359Jy sigma=1.523923 ! Fit after self-cal, rms=0.190880Jy sigma=1.521732 ! Inverting map and beam ! Estimated beam: bmin=0.6484 mas, bmaj=0.9719 mas, bpa=-1.988 degrees ! Estimated noise=0.521177 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00650586 Jy ! Component: 100 - total flux cleaned = 0.00806805 Jy ! Total flux subtracted in 100 components = 0.00806805 Jy ! Clean residual min=-0.007866 max=0.007095 Jy/beam ! Clean residual mean=0.000015 rms=0.001476 Jy/beam ! Combined flux in latest and established models = 4.51997 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 25 model components to the UV plane model. ! The established model now contains 125 components and 4.51997 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190365Jy sigma=1.514409 ! Fit after self-cal, rms=0.189918Jy sigma=1.512394 ! Inverting map and beam ! Estimated beam: bmin=0.648 mas, bmaj=0.9713 mas, bpa=-1.988 degrees ! Estimated noise=0.521164 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00623663 Jy ! Component: 100 - total flux cleaned = 0.00772318 Jy ! Total flux subtracted in 100 components = 0.00772318 Jy ! Clean residual min=-0.007726 max=0.006884 Jy/beam ! Clean residual mean=0.000013 rms=0.001441 Jy/beam ! Combined flux in latest and established models = 4.52769 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 22 model components to the UV plane model. ! The established model now contains 126 components and 4.52769 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.189435Jy sigma=1.505577 ! Fit after self-cal, rms=0.189012Jy sigma=1.503713 ! Inverting map and beam ! Estimated beam: bmin=0.6476 mas, bmaj=0.9708 mas, bpa=-1.989 degrees ! Estimated noise=0.521126 mJy/beam. ! Added new window around map position (1.6, -2.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00798012 Jy ! Component: 100 - total flux cleaned = 0.0106532 Jy ! Total flux subtracted in 100 components = 0.0106532 Jy ! Clean residual min=-0.007535 max=0.006004 Jy/beam ! Clean residual mean=0.000009 rms=0.001404 Jy/beam ! Combined flux in latest and established models = 4.53834 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 130 components and 4.53834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.188545Jy sigma=1.496307 ! Fit after self-cal, rms=0.188166Jy sigma=1.494546 ! Inverting map and beam ! Estimated beam: bmin=0.6473 mas, bmaj=0.9703 mas, bpa=-1.992 degrees ! Estimated noise=0.521101 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.00576772 Jy ! Component: 100 - total flux cleaned = 0.00775498 Jy ! Total flux subtracted in 100 components = 0.00775498 Jy ! Clean residual min=-0.007454 max=0.005814 Jy/beam ! Clean residual mean=0.000008 rms=0.001372 Jy/beam ! Combined flux in latest and established models = 4.5461 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 25 model components to the UV plane model. ! The established model now contains 131 components and 4.5461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187744Jy sigma=1.488498 ! Fit after self-cal, rms=0.186976Jy sigma=1.475528 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 0552+398_X ! Writing UV FITS file: 0552+398_X.uvf ! Writing 131 model components to file: 0552+398_X.mod ! wwins: Wrote 6 windows to 0552+398_X.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6473 x 0.9703 at -1.992 degrees (North through East) ! Clean map min=-0.0079312 max=3.1789 Jy/beam ! Writing clean map to FITS file: 0552+398_X.fits ! Writing difmap environment to: 0552+398_X.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 0552+398_X] device 0552+398_Xuv.ps/vps,1,2 ! Attempting to open device: '0552+398_Xuv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 0552+398_Xim.ps/vps ! Attempting to open device: '0552+398_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.6473 mas, bmaj=0.9703 mas, bpa=-1.992 degrees ! Estimated noise=0.521101 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6473 x 0.9703 at -1.992 degrees (North through East) ! Clean map min=-0.0079312 max=3.1789 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo5-0.1] quit ! Quitting program ! Log file 0552+398_X.log closed on Fri Jun 27 12:15:33 2008