! Started logfile: 0552+398_X.log on Thu Aug 7 14:18:35 2008 observe 0552+398_X.edt ! Reading UV FITS file: 0552+398_X.edt ! AN table 1: 148 integrations on 171 of 171 possible baselines. ! AN table 2: 33 integrations on 171 of 171 possible baselines. ! AN table 3: 74 integrations on 171 of 171 possible baselines. ! AN table 4: 23 integrations on 171 of 171 possible baselines. ! AN table 5: 24 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.147748 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 1526 lines of history. ! ! Reading 30520 visibilities. ![@gonzo7-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 = 7.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 1892 telescope corrections were flagged in sub-array 1. ! A total of 382 telescope corrections were flagged in sub-array 2. ! A total of 1315 telescope corrections were flagged in sub-array 3. ! A total of 352 telescope corrections were flagged in sub-array 4. ! A total of 417 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 1892 telescope corrections were flagged in sub-array 1. ! A total of 379 telescope corrections were flagged in sub-array 2. ! A total of 1315 telescope corrections were flagged in sub-array 3. ! A total of 352 telescope corrections were flagged in sub-array 4. ! A total of 410 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 1892 telescope corrections were flagged in sub-array 1. ! A total of 379 telescope corrections were flagged in sub-array 2. ! A total of 1314 telescope corrections were flagged in sub-array 3. ! A total of 352 telescope corrections were flagged in sub-array 4. ! A total of 410 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 1892 telescope corrections were flagged in sub-array 1. ! A total of 381 telescope corrections were flagged in sub-array 2. ! A total of 1314 telescope corrections were flagged in sub-array 3. ! A total of 352 telescope corrections were flagged in sub-array 4. ! A total of 410 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=2.834150Jy sigma=33.947752 ! Fit after self-cal, rms=2.836992Jy sigma=33.991391 ! 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.5364 mas, bmaj=0.8058 mas, bpa=-4.342 degrees ! Estimated noise=0.865954 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.39795 Jy ! Total flux subtracted in 50 components = 2.39795 Jy ! Clean residual min=-0.109758 max=0.721876 Jy/beam ! Clean residual mean=-0.000388 rms=0.029648 Jy/beam ! Combined flux in latest and established models = 2.39795 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 4 components and 2.39795 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.641754Jy sigma=20.058059 ! Fit after self-cal, rms=1.641682Jy sigma=20.057970 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.768944 Jy ! Total flux subtracted in 50 components = 0.768944 Jy ! Clean residual min=-0.077862 max=0.358478 Jy/beam ! Clean residual mean=-0.000250 rms=0.019903 Jy/beam ! Combined flux in latest and established models = 3.16689 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 10 components and 3.16689 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.120228Jy sigma=13.814326 ! Fit after self-cal, rms=1.118745Jy sigma=13.802806 ! Inverting map ! Added new window around map position (-0.5, 0.1). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.440906 Jy ! Total flux subtracted in 50 components = 0.440906 Jy ! Clean residual min=-0.059579 max=0.236770 Jy/beam ! Clean residual mean=-0.000173 rms=0.015299 Jy/beam ! Combined flux in latest and established models = 3.6078 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 19 components and 3.6078 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.863049Jy sigma=10.661566 ! Fit after self-cal, rms=0.861349Jy sigma=10.646199 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.302519 Jy ! Total flux subtracted in 50 components = 0.302519 Jy ! Clean residual min=-0.045405 max=0.167686 Jy/beam ! Clean residual mean=-0.000121 rms=0.012462 Jy/beam ! Combined flux in latest and established models = 3.91032 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 30 components and 3.91032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.702778Jy sigma=8.648333 ! Fit after self-cal, rms=0.701565Jy sigma=8.637651 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.219948 Jy ! Total flux subtracted in 50 components = 0.219948 Jy ! Clean residual min=-0.038496 max=0.124142 Jy/beam ! Clean residual mean=-0.000086 rms=0.010652 Jy/beam ! Combined flux in latest and established models = 4.13026 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 36 components and 4.13026 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.597073Jy sigma=7.297755 ! Fit after self-cal, rms=0.596005Jy sigma=7.286482 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.164078 Jy ! Total flux subtracted in 50 components = 0.164078 Jy ! Clean residual min=-0.035075 max=0.093109 Jy/beam ! Clean residual mean=-0.000060 rms=0.009480 Jy/beam ! Combined flux in latest and established models = 4.29434 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 43 components and 4.29434 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.525520Jy sigma=6.365242 ! Fit after self-cal, rms=0.524948Jy sigma=6.357386 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.124943 Jy ! Total flux subtracted in 50 components = 0.124943 Jy ! Clean residual min=-0.033544 max=0.071673 Jy/beam ! Clean residual mean=-0.000040 rms=0.008741 Jy/beam ! Combined flux in latest and established models = 4.41928 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 19 model components to the UV plane model. ! The established model now contains 52 components and 4.41928 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.476823Jy sigma=5.717949 ! Fit after self-cal, rms=0.476395Jy sigma=5.711674 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0972876 Jy ! Total flux subtracted in 50 components = 0.0972876 Jy ! Clean residual min=-0.034286 max=0.056939 Jy/beam ! Clean residual mean=-0.000026 rms=0.008261 Jy/beam ! Combined flux in latest and established models = 4.51657 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 57 components and 4.51657 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.443044Jy sigma=5.258900 ! Fit after self-cal, rms=0.442811Jy sigma=5.254490 ! Inverting map ! Added new window around map position (-1, 0.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0769758 Jy ! Total flux subtracted in 50 components = 0.0769758 Jy ! Clean residual min=-0.034544 max=0.046840 Jy/beam ! Clean residual mean=-0.000015 rms=0.007952 Jy/beam ! Combined flux in latest and established models = 4.59355 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 61 components and 4.59355 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.419598Jy sigma=4.935018 ! Fit after self-cal, rms=0.419520Jy sigma=4.932319 ! 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.6801 mas, bmaj=0.9948 mas, bpa=-1.547 degrees ! Estimated noise=0.515871 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.140978 Jy ! Component: 100 - total flux cleaned = 0.227035 Jy ! Total flux subtracted in 100 components = 0.227035 Jy ! Clean residual min=-0.031925 max=0.053047 Jy/beam ! Clean residual mean=0.001182 rms=0.010591 Jy/beam ! Combined flux in latest and established models = 4.82058 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 84 components and 4.82058 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.370910Jy sigma=4.234906 ! Fit after self-cal, rms=0.370100Jy sigma=4.217593 ! 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 24 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.96 GC 1.01 HH 1.00* ! HN 0.99 KK 0.90 KP 0.98 LA 0.96 ! MA 1.04 MK 1.02 NL 0.98 NY 0.99 ! OV 0.97 PT 0.96 SC 1.02 TS 1.01 ! WF 1.00 WZ 1.01 GG 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.11 FD 0.91 GC 0.99 HH 1.00* ! HN 1.11 KK 1.00* KP 1.13 LA 0.93 ! MA 1.00 MK 1.00* NL 1.14 NY 0.94 ! OV 1.01 PT 0.94 SC 1.17 TS 1.00* ! WF 0.91 WZ 0.90 GG 1.00* ! ! A total of 39 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.07 FD 0.93 GC 0.76 HH 1.00* ! HN 1.11 KK 1.00* KP 1.13 LA 0.97 ! MA 1.00* MK 1.15 NL 1.10 NY 0.93 ! OV 0.98 PT 1.00* SC 1.00* TS 1.00* ! WF 1.10 WZ 1.00* GG 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 1.05 GC 0.99 HH 1.00* ! HN 1.00* KK 1.22 KP 1.13 LA 0.99 ! MA 1.00* MK 1.32 NL 1.14 NY 0.90 ! OV 1.00 PT 1.04 SC 1.00* TS 1.00 ! WF 1.00* WZ 0.93 GG 1.00* ! ! A total of 27 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.07 HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 1.01 MK 1.00* NL 1.00* NY 0.90 ! OV 1.00* PT 1.00* SC 1.00* TS 1.00* ! WF 1.00* WZ 0.86 GG 1.00* ! ! ! Correcting IF 2. ! A total of 24 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.95 GC 1.01 HH 1.00* ! HN 1.00 KK 0.89 KP 0.98 LA 0.96 ! MA 1.05 MK 1.02 NL 0.97 NY 0.99 ! OV 0.97 PT 0.96 SC 1.01 TS 1.00 ! WF 0.99 WZ 1.00 GG 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.12 FD 1.01 GC 1.02 HH 1.00* ! HN 1.02 KK 1.00* KP 0.98 LA 0.93 ! MA 0.99 MK 1.00* NL 0.95 NY 0.97 ! OV 1.08 PT 1.13 SC 1.00 TS 1.00* ! WF 0.90 WZ 0.96 GG 1.00* ! ! A total of 39 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.11 FD 1.01 GC 1.09 HH 1.00* ! HN 0.95 KK 1.00* KP 0.96 LA 0.97 ! MA 1.00* MK 1.11 NL 0.97 NY 0.92 ! OV 0.99 PT 1.00* SC 1.00* TS 1.00* ! WF 1.11 WZ 1.00* GG 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.07 FD 1.14 GC 1.07 HH 1.00* ! HN 1.00* KK 1.22 KP 0.97 LA 0.97 ! MA 1.00* MK 1.23 NL 1.02 NY 0.90 ! OV 1.06 PT 1.21 SC 1.00* TS 0.99 ! WF 1.00* WZ 0.97 GG 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.12 HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 0.94 MK 1.00* NL 1.00* NY 0.93 ! OV 1.00* PT 1.00* SC 1.00* TS 1.00* ! WF 1.00* WZ 0.89 GG 1.00* ! ! ! Correcting IF 3. ! A total of 24 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 0.96 GC 1.02 HH 1.00* ! HN 1.01 KK 0.99 KP 0.98 LA 0.96 ! MA 0.88 MK 1.01 NL 1.02 NY 1.00 ! OV 0.98 PT 0.97 SC 1.01 TS 1.02 ! WF 1.00 WZ 0.99 GG 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 1.06 GC 1.05 HH 1.00* ! HN 0.86 KK 1.00* KP 0.89 LA 1.10 ! MA 1.00 MK 1.00* NL 1.08 NY 0.97 ! OV 0.99 PT 1.05 SC 0.87 TS 1.00* ! WF 0.90 WZ 0.92 GG 1.00* ! ! A total of 39 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.10 FD 1.05 GC 1.03 HH 1.00* ! HN 0.86 KK 1.00* KP 0.91 LA 1.06 ! MA 1.00* MK 1.10 NL 1.07 NY 0.92 ! OV 0.95 PT 1.21 SC 1.00* TS 1.00* ! WF 0.99 WZ 1.00* GG 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 1.26 GC 1.04 HH 1.00* ! HN 1.00* KK 1.20 KP 0.94 LA 1.23 ! MA 1.00* MK 1.24 NL 1.09 NY 0.96 ! OV 0.98 PT 1.12 SC 1.00* TS 0.99 ! WF 1.00* WZ 0.98 GG 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.13 HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 0.99 MK 1.00* NL 1.00* NY 0.91 ! OV 1.00* PT 1.00* SC 1.00* TS 1.00* ! WF 1.00* WZ 0.88 GG 1.00* ! ! ! Correcting IF 4. ! A total of 24 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.97 GC 1.01 HH 1.00* ! HN 1.00 KK 0.95 KP 0.98 LA 0.96 ! MA 0.98 MK 1.01 NL 0.99 NY 1.00 ! OV 0.98 PT 0.97 SC 1.02 TS 1.00 ! WF 0.99 WZ 1.02 GG 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 FD 1.18 GC 1.00 HH 1.00* ! HN 0.92 KK 1.00* KP 0.91 LA 1.36 ! MA 0.93 MK 1.00* NL 1.09 NY 0.92 ! OV 1.08 PT 1.08 SC 0.92 TS 1.00* ! WF 0.89 WZ 0.96 GG 1.00* ! ! A total of 39 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.08 FD 1.17 GC 0.72 HH 1.00* ! HN 0.91 KK 1.00* KP 0.94 LA 1.14 ! MA 1.00* MK 1.22 NL 1.07 NY 0.93 ! OV 1.04 PT 0.91 SC 1.00* TS 1.00* ! WF 1.10 WZ 1.00* GG 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.05 FD 1.32 GC 1.01 HH 1.00* ! HN 1.00* KK 1.21 KP 0.95 LA 1.39 ! MA 1.00* MK 1.45 NL 1.10 NY 0.92 ! OV 1.09 PT 1.18 SC 1.00* TS 1.01 ! WF 1.00* WZ 0.93 GG 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.08 HH 1.00* ! HN 1.00* KK 1.00* KP 1.00* LA 1.00* ! MA 0.93 MK 1.00* NL 1.00* NY 0.91 ! OV 1.00* PT 1.00* SC 1.00* TS 1.00* ! WF 1.00* WZ 0.90 GG 1.00* ! ! ! Fit before self-cal, rms=0.370100Jy sigma=4.217593 ! Fit after self-cal, rms=0.236196Jy sigma=1.914963 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6817 mas, bmaj=1.011 mas, bpa=-1.91 degrees ! Estimated noise=0.513286 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.0297514 Jy ! Component: 100 - total flux cleaned = 0.0291408 Jy ! Total flux subtracted in 100 components = 0.0291408 Jy ! Clean residual min=-0.013576 max=0.015737 Jy/beam ! Clean residual mean=0.000897 rms=0.002725 Jy/beam ! Combined flux in latest and established models = 4.84972 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 91 components and 4.84972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.228491Jy sigma=1.871082 ! Fit after self-cal, rms=0.228342Jy sigma=1.854604 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00638775 Jy ! Component: 100 - total flux cleaned = 0.00536026 Jy ! Total flux subtracted in 100 components = 0.00536026 Jy ! Clean residual min=-0.012466 max=0.014260 Jy/beam ! Clean residual mean=0.000865 rms=0.002609 Jy/beam ! Combined flux in latest and established models = 4.85508 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 96 components and 4.85508 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225906Jy sigma=1.847198 ! Fit after self-cal, rms=0.226313Jy sigma=1.843655 ! 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 (0.1, -10.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00949365 Jy ! Component: 100 - total flux cleaned = 0.0110984 Jy ! Total flux subtracted in 100 components = 0.0110984 Jy ! Clean residual min=-0.012588 max=0.013799 Jy/beam ! Clean residual mean=0.000797 rms=0.002529 Jy/beam ! Combined flux in latest and established models = 4.86618 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 105 components and 4.86618 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.224844Jy sigma=1.835881 ! Fit after self-cal, rms=0.190028Jy sigma=1.376679 ! Inverting map and beam ! Estimated beam: bmin=0.6808 mas, bmaj=1.012 mas, bpa=-1.836 degrees ! Estimated noise=0.513012 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0159133 Jy ! Component: 100 - total flux cleaned = 0.017617 Jy ! Total flux subtracted in 100 components = 0.017617 Jy ! Clean residual min=-0.007705 max=0.007703 Jy/beam ! Clean residual mean=0.000489 rms=0.001326 Jy/beam ! Combined flux in latest and established models = 4.8838 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 14 model components to the UV plane model. ! The established model now contains 113 components and 4.8838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187370Jy sigma=1.367933 ! Fit after self-cal, rms=0.186391Jy sigma=1.346460 ! Inverting map and beam ! Estimated beam: bmin=0.6801 mas, bmaj=1.012 mas, bpa=-1.815 degrees ! Estimated noise=0.513012 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0108163 Jy ! Component: 100 - total flux cleaned = 0.0126281 Jy ! Total flux subtracted in 100 components = 0.0126281 Jy ! Clean residual min=-0.006548 max=0.007016 Jy/beam ! Clean residual mean=0.000399 rms=0.001167 Jy/beam ! Combined flux in latest and established models = 4.89643 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 114 components and 4.89643 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.184725Jy sigma=1.341656 ! Fit after self-cal, rms=0.184487Jy sigma=1.333251 ! Inverting map and beam ! Estimated beam: bmin=0.6795 mas, bmaj=1.012 mas, bpa=-1.803 degrees ! Estimated noise=0.513002 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00843509 Jy ! Component: 100 - total flux cleaned = 0.00969761 Jy ! Total flux subtracted in 100 components = 0.00969761 Jy ! Clean residual min=-0.005753 max=0.006555 Jy/beam ! Clean residual mean=0.000346 rms=0.001070 Jy/beam ! Combined flux in latest and established models = 4.90613 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 118 components and 4.90613 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183285Jy sigma=1.329841 ! Fit after self-cal, rms=0.183226Jy sigma=1.325588 ! Inverting map and beam ! Estimated beam: bmin=0.6788 mas, bmaj=1.012 mas, bpa=-1.798 degrees ! Estimated noise=0.512942 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00644644 Jy ! Component: 100 - total flux cleaned = 0.00786957 Jy ! Total flux subtracted in 100 components = 0.00786957 Jy ! Clean residual min=-0.005357 max=0.006103 Jy/beam ! Clean residual mean=0.000306 rms=0.001005 Jy/beam ! Combined flux in latest and established models = 4.91399 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 123 components and 4.91399 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182323Jy sigma=1.322894 ! Fit after self-cal, rms=0.182317Jy sigma=1.320372 ! Inverting map and beam ! Estimated beam: bmin=0.6782 mas, bmaj=1.011 mas, bpa=-1.796 degrees ! Estimated noise=0.512847 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00519829 Jy ! Component: 100 - total flux cleaned = 0.00649474 Jy ! Total flux subtracted in 100 components = 0.00649474 Jy ! Clean residual min=-0.005046 max=0.005730 Jy/beam ! Clean residual mean=0.000275 rms=0.000958 Jy/beam ! Combined flux in latest and established models = 4.92049 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 17 model components to the UV plane model. ! The established model now contains 125 components and 4.92049 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181607Jy sigma=1.318002 ! Fit after self-cal, rms=0.181619Jy sigma=1.316363 ! Inverting map and beam ! Estimated beam: bmin=0.6777 mas, bmaj=1.011 mas, bpa=-1.795 degrees ! Estimated noise=0.512742 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00426271 Jy ! Component: 100 - total flux cleaned = 0.00543623 Jy ! Total flux subtracted in 100 components = 0.00543623 Jy ! Clean residual min=-0.004812 max=0.005471 Jy/beam ! Clean residual mean=0.000249 rms=0.000920 Jy/beam ! Combined flux in latest and established models = 4.92593 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 14 model components to the UV plane model. ! The established model now contains 127 components and 4.92593 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181033Jy sigma=1.314321 ! Fit after self-cal, rms=0.181052Jy sigma=1.313147 ! Inverting map and beam ! Estimated beam: bmin=0.6773 mas, bmaj=1.011 mas, bpa=-1.795 degrees ! Estimated noise=0.512613 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00326842 Jy ! Component: 100 - total flux cleaned = 0.00456499 Jy ! Total flux subtracted in 100 components = 0.00456499 Jy ! Clean residual min=-0.004542 max=0.005302 Jy/beam ! Clean residual mean=0.000226 rms=0.000891 Jy/beam ! Combined flux in latest and established models = 4.93049 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 130 components and 4.93049 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180553Jy sigma=1.311522 ! Fit after self-cal, rms=0.180587Jy sigma=1.310606 ! Inverting map and beam ! Estimated beam: bmin=0.6768 mas, bmaj=1.011 mas, bpa=-1.796 degrees ! Estimated noise=0.512475 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00265841 Jy ! Component: 100 - total flux cleaned = 0.00364708 Jy ! Total flux subtracted in 100 components = 0.00364708 Jy ! Clean residual min=-0.004305 max=0.005109 Jy/beam ! Clean residual mean=0.000205 rms=0.000868 Jy/beam ! Combined flux in latest and established models = 4.93414 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 13 model components to the UV plane model. ! The established model now contains 130 components and 4.93414 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180155Jy sigma=1.309324 ! Fit after self-cal, rms=0.180189Jy sigma=1.308585 ! Inverting map and beam ! Estimated beam: bmin=0.6765 mas, bmaj=1.011 mas, bpa=-1.798 degrees ! Estimated noise=0.512307 mJy/beam. ! Added new window around map position (1.2, -2.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00308923 Jy ! Component: 100 - total flux cleaned = 0.00422587 Jy ! Total flux subtracted in 100 components = 0.00422587 Jy ! Clean residual min=-0.004318 max=0.005093 Jy/beam ! Clean residual mean=0.000179 rms=0.000845 Jy/beam ! Combined flux in latest and established models = 4.93836 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 134 components and 4.93836 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179789Jy sigma=1.307433 ! Fit after self-cal, rms=0.179849Jy sigma=1.306755 ! Inverting map and beam ! Estimated beam: bmin=0.6762 mas, bmaj=1.011 mas, bpa=-1.797 degrees ! Estimated noise=0.512142 mJy/beam. ! Added new window around map position (0.6, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00129294 Jy ! Component: 100 - total flux cleaned = 0.00207714 Jy ! Total flux subtracted in 100 components = 0.00207714 Jy ! Clean residual min=-0.003957 max=0.004165 Jy/beam ! Clean residual mean=0.000169 rms=0.000831 Jy/beam ! Combined flux in latest and established models = 4.94044 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 17 model components to the UV plane model. ! The established model now contains 141 components and 4.94044 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179553Jy sigma=1.303597 ! Fit after self-cal, rms=0.179482Jy sigma=1.303018 ! Inverting map and beam ! Estimated beam: bmin=0.6758 mas, bmaj=1.011 mas, bpa=-1.791 degrees ! Estimated noise=0.511962 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00112431 Jy ! Component: 100 - total flux cleaned = 0.00149859 Jy ! Total flux subtracted in 100 components = 0.00149859 Jy ! Clean residual min=-0.003771 max=0.003911 Jy/beam ! Clean residual mean=0.000160 rms=0.000809 Jy/beam ! Combined flux in latest and established models = 4.94194 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 18 model components to the UV plane model. ! The established model now contains 143 components and 4.94194 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179205Jy sigma=1.300919 ! Fit after self-cal, rms=0.179181Jy sigma=1.300445 ! Inverting map and beam ! Estimated beam: bmin=0.6756 mas, bmaj=1.011 mas, bpa=-1.786 degrees ! Estimated noise=0.511769 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000979819 Jy ! Component: 100 - total flux cleaned = 0.000981929 Jy ! Total flux subtracted in 100 components = 0.000981929 Jy ! Clean residual min=-0.003613 max=0.003785 Jy/beam ! Clean residual mean=0.000151 rms=0.000792 Jy/beam ! Combined flux in latest and established models = 4.94292 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 146 components and 4.94292 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.178915Jy sigma=1.298933 ! Fit after self-cal, rms=0.178921Jy sigma=1.298527 ! Inverting map and beam ! Estimated beam: bmin=0.6753 mas, bmaj=1.011 mas, bpa=-1.781 degrees ! Estimated noise=0.511566 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.00068372 Jy ! Component: 100 - total flux cleaned = 0.000682778 Jy ! Total flux subtracted in 100 components = 0.000682778 Jy ! Clean residual min=-0.003455 max=0.003576 Jy/beam ! Clean residual mean=0.000143 rms=0.000777 Jy/beam ! Combined flux in latest and established models = 4.94361 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 150 components and 4.9436 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.178671Jy sigma=1.297323 ! Fit after self-cal, rms=0.177882Jy sigma=1.288034 ! Inverting map ! Added new window around map position (-1.6, -0.3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00103818 Jy ! Component: 100 - total flux cleaned = 0.000883711 Jy ! Total flux subtracted in 100 components = 0.000883711 Jy ! Clean residual min=-0.003289 max=0.003430 Jy/beam ! Clean residual mean=0.000138 rms=0.000754 Jy/beam ! Combined flux in latest and established models = 4.94449 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 158 components and 4.94449 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177657Jy sigma=1.287505 ! Fit after self-cal, rms=0.177860Jy sigma=1.287206 ! 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 158 model components to file: 0552+398_X.mod ! wwins: Wrote 7 windows to 0552+398_X.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6753 x 1.011 at -1.781 degrees (North through East) ! Clean map min=-0.003364 max=3.3658 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.6753 mas, bmaj=1.011 mas, bpa=-1.781 degrees ! Estimated noise=0.511566 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6753 x 1.011 at -1.781 degrees (North through East) ! Clean map min=-0.003364 max=3.3658 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo7-0.1] quit ! Quitting program ! Log file 0552+398_X.log closed on Thu Aug 7 14:19:05 2008