! Started logfile: 0014+813_S.log on Mon Aug 18 11:16:51 2008 observe 0014+813_S.edt2 ! Reading UV FITS file: 0014+813_S.edt2 ! AN table 1: 102 integrations on 171 of 171 possible baselines. ! AN table 2: 71 integrations on 171 of 171 possible baselines. ! AN table 3: 51 integrations on 171 of 171 possible baselines. ! AN table 4: 106 integrations on 171 of 171 possible baselines. ! AN table 5: 60 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.258599 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 0014+813 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.22497e+09 8e+06 1 8e+06 ! 02 2 2.24497e+09 8e+06 1 8e+06 ! 03 3 2.33497e+09 8e+06 1 8e+06 ! 04 4 2.36497e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1619 lines of history. ! ! Reading 68984 visibilities. ![@gonzo6-0.4 0014+813_S] print "***********************************" ! *********************************** print "*** Using gonzo instruction set ***" ! *** Using gonzo instruction set *** print "***********************************" ! *********************************** float field_size field_size = 1024 float field_cell field_cell = 0.4 integer clean_niter clean_niter = 50 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float win_mult win_mult = 1.8 float solint solint = 3.6/6.0 logical doflt doflt = true logical doamp float new_peak float flux_cutoff mapsize field_size,field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 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 1214 telescope corrections were flagged in sub-array 1. ! A total of 890 telescope corrections were flagged in sub-array 2. ! A total of 641 telescope corrections were flagged in sub-array 3. ! A total of 1178 telescope corrections were flagged in sub-array 4. ! A total of 839 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 1217 telescope corrections were flagged in sub-array 1. ! A total of 908 telescope corrections were flagged in sub-array 2. ! A total of 648 telescope corrections were flagged in sub-array 3. ! A total of 1218 telescope corrections were flagged in sub-array 4. ! A total of 838 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 1215 telescope corrections were flagged in sub-array 1. ! A total of 900 telescope corrections were flagged in sub-array 2. ! A total of 677 telescope corrections were flagged in sub-array 3. ! A total of 1158 telescope corrections were flagged in sub-array 4. ! A total of 837 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 1194 telescope corrections were flagged in sub-array 1. ! A total of 898 telescope corrections were flagged in sub-array 2. ! A total of 665 telescope corrections were flagged in sub-array 3. ! A total of 1164 telescope corrections were flagged in sub-array 4. ! A total of 837 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.572222Jy sigma=8.694957 ! Fit after self-cal, rms=0.570250Jy sigma=8.702277 ! 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=1.917 mas, bmaj=2.11 mas, bpa=-40.29 degrees ! Estimated noise=0.437331 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.331162 Jy ! Total flux subtracted in 50 components = 0.331162 Jy ! Clean residual min=-0.008641 max=0.092519 Jy/beam ! Clean residual mean=-0.000128 rms=0.002444 Jy/beam ! Combined flux in latest and established models = 0.331162 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 0.331162 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207091Jy sigma=3.024578 ! Fit after self-cal, rms=0.207092Jy sigma=3.024579 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0766336 Jy ! Total flux subtracted in 50 components = 0.0766336 Jy ! Clean residual min=-0.007481 max=0.030582 Jy/beam ! Clean residual mean=-0.000100 rms=0.001706 Jy/beam ! Combined flux in latest and established models = 0.407796 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 0.407796 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.172529Jy sigma=2.310243 ! Fit after self-cal, rms=0.172538Jy sigma=2.310101 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0320914 Jy ! Total flux subtracted in 50 components = 0.0320914 Jy ! Clean residual min=-0.007550 max=0.028939 Jy/beam ! Clean residual mean=-0.000089 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 0.439887 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 6 model components to the UV plane model. ! The established model now contains 9 components and 0.439887 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.164267Jy sigma=2.110289 ! Fit after self-cal, rms=0.164269Jy sigma=2.110059 ! Inverting map ! Added new window around map position (-1.6, -8.8). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.026127 Jy ! Total flux subtracted in 50 components = 0.026127 Jy ! Clean residual min=-0.007462 max=0.024567 Jy/beam ! Clean residual mean=-0.000080 rms=0.001419 Jy/beam ! Combined flux in latest and established models = 0.466014 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 20 components and 0.466014 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158418Jy sigma=1.979157 ! Fit after self-cal, rms=0.158019Jy sigma=1.974621 ! Inverting map ! Added new window around map position (-1.2, -4.8). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0241608 Jy ! Total flux subtracted in 50 components = 0.0241608 Jy ! Clean residual min=-0.006833 max=0.012085 Jy/beam ! Clean residual mean=-0.000072 rms=0.001295 Jy/beam ! Combined flux in latest and established models = 0.490175 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 31 components and 0.490175 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.152579Jy sigma=1.850874 ! Fit after self-cal, rms=0.152146Jy sigma=1.845460 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0171897 Jy ! Total flux subtracted in 50 components = 0.0171897 Jy ! Clean residual min=-0.006336 max=0.010065 Jy/beam ! Clean residual mean=-0.000066 rms=0.001219 Jy/beam ! Combined flux in latest and established models = 0.507365 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 37 components and 0.507365 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.149092Jy sigma=1.773129 ! Fit after self-cal, rms=0.148941Jy sigma=1.771517 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0143009 Jy ! Total flux subtracted in 50 components = 0.0143009 Jy ! Clean residual min=-0.006044 max=0.008662 Jy/beam ! Clean residual mean=-0.000061 rms=0.001170 Jy/beam ! Combined flux in latest and established models = 0.521665 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 45 components and 0.521666 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146722Jy sigma=1.718789 ! Fit after self-cal, rms=0.146638Jy sigma=1.717942 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.01237 Jy ! Total flux subtracted in 50 components = 0.01237 Jy ! Clean residual min=-0.005833 max=0.007560 Jy/beam ! Clean residual mean=-0.000057 rms=0.001134 Jy/beam ! Combined flux in latest and established models = 0.534036 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 51 components and 0.534036 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144937Jy sigma=1.677668 ! Fit after self-cal, rms=0.144880Jy sigma=1.677098 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0108206 Jy ! Total flux subtracted in 50 components = 0.0108206 Jy ! Clean residual min=-0.005648 max=0.007040 Jy/beam ! Clean residual mean=-0.000053 rms=0.001105 Jy/beam ! Combined flux in latest and established models = 0.544856 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 54 components and 0.544856 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143538Jy sigma=1.645673 ! Fit after self-cal, rms=0.143500Jy sigma=1.645310 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00950775 Jy ! Total flux subtracted in 50 components = 0.00950775 Jy ! Clean residual min=-0.005500 max=0.006592 Jy/beam ! Clean residual mean=-0.000050 rms=0.001082 Jy/beam ! Combined flux in latest and established models = 0.554364 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 56 components and 0.554364 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.142434Jy sigma=1.620910 ! Fit after self-cal, rms=0.142405Jy sigma=1.620616 ! 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=2.439 mas, bmaj=2.654 mas, bpa=-30.76 degrees ! Estimated noise=0.229453 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.0213396 Jy ! Component: 100 - total flux cleaned = 0.0359498 Jy ! Total flux subtracted in 100 components = 0.0359498 Jy ! Clean residual min=-0.004240 max=0.009022 Jy/beam ! Clean residual mean=-0.000043 rms=0.001244 Jy/beam ! Combined flux in latest and established models = 0.590314 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 42 model components to the UV plane model. ! The established model now contains 84 components and 0.590314 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139385Jy sigma=1.558761 ! Fit after self-cal, rms=0.139272Jy sigma=1.555761 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0114084 Jy ! Component: 100 - total flux cleaned = 0.0193957 Jy ! Total flux subtracted in 100 components = 0.0193957 Jy ! Clean residual min=-0.003915 max=0.006604 Jy/beam ! Clean residual mean=-0.000082 rms=0.001077 Jy/beam ! Combined flux in latest and established models = 0.609709 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 44 model components to the UV plane model. ! The established model now contains 100 components and 0.609709 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138297Jy sigma=1.542031 ! Fit after self-cal, rms=0.138269Jy sigma=1.540773 ! Inverting map ! Added new window around map position (2.4, -5.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00757791 Jy ! Component: 100 - total flux cleaned = 0.012748 Jy ! Total flux subtracted in 100 components = 0.0127479 Jy ! Clean residual min=-0.003821 max=0.005600 Jy/beam ! Clean residual mean=-0.000107 rms=0.001000 Jy/beam ! Combined flux in latest and established models = 0.622457 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 111 components and 0.622457 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137782Jy sigma=1.536955 ! Fit after self-cal, rms=0.137783Jy sigma=1.536208 ! 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 18 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.00 GC 1.03 HH 1.00* ! HN 1.00 KK 0.88 KP 0.99 LA 1.01 ! MA 0.94 MK 1.02 NL 1.02 NY 1.01 ! OV 1.03 PT 0.99 SC 1.02 TS 1.06 ! WF 1.01 WZ 0.92 GG 1.00* ! ! A total of 21 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.90 FD 1.12 GC 0.97 HH 1.00* ! HN 0.95 KK 0.79 KP 0.93 LA 1.09 ! MA 1.07 MK 0.91 NL 0.85 NY 1.14 ! OV 1.00 PT 1.08 SC 0.81 TS 1.02 ! WF 1.04 WZ 1.17 GG 1.00* ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 1.09 GC 1.07 HH 1.00* ! HN 0.93 KK 1.00* KP 0.91 LA 1.09 ! MA 1.16 MK 0.94 NL 0.95 NY 1.17 ! OV 0.99 PT 1.06 SC 0.88 TS 1.00* ! WF 1.08 WZ 1.00* GG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.92 FD 1.11 GC 1.48 HH 1.00* ! HN 0.93 KK 1.00* KP 0.90 LA 1.10 ! MA 1.39 MK 0.91 NL 0.92 NY 1.10 ! OV 1.01 PT 1.04 SC 0.84 TS 1.00* ! WF 0.98 WZ 1.06 GG 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.07 GC 1.00* HH 1.00* ! HN 0.98 KK 0.78 KP 0.96 LA 1.10 ! MA 1.20 MK 0.97 NL 0.99 NY 1.19 ! OV 1.03 PT 1.08 SC 0.86 TS 1.05 ! WF 1.09 WZ 1.01 GG 1.00* ! ! ! Correcting IF 2. ! A total of 15 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.02 GC 1.02 HH 1.00* ! HN 1.02 KK 0.95 KP 1.01 LA 1.01 ! MA 1.02 MK 1.00 NL 1.02 NY 1.04 ! OV 1.03 PT 1.00 SC 1.01 TS 1.09 ! WF 0.97 WZ 0.92 GG 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.89 FD 1.01 GC 0.91 HH 1.00* ! HN 1.00 KK 0.81 KP 1.04 LA 1.11 ! MA 1.12 MK 0.99 NL 0.90 NY 1.11 ! OV 0.95 PT 0.92 SC 0.96 TS 1.01 ! WF 1.09 WZ 1.05 GG 1.00* ! ! A total of 27 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 1.00 GC 1.01 HH 1.00* ! HN 0.98 KK 1.00* KP 1.02 LA 1.08 ! MA 1.10 MK 0.98 NL 1.03 NY 1.08 ! OV 0.95 PT 0.90 SC 0.99 TS 1.00* ! WF 1.07 WZ 1.00* GG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.89 FD 1.04 GC 1.51 HH 1.00* ! HN 0.98 KK 1.00* KP 1.05 LA 1.08 ! MA 1.27 MK 0.96 NL 1.07 NY 1.08 ! OV 0.96 PT 0.86 SC 1.00 TS 1.00* ! WF 0.93 WZ 0.99 GG 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.09 GC 1.00* HH 1.00* ! HN 1.04 KK 0.81 KP 1.10 LA 1.07 ! MA 1.07 MK 1.05 NL 1.11 NY 1.14 ! OV 1.01 PT 0.92 SC 1.01 TS 1.06 ! WF 1.03 WZ 1.01 GG 1.00* ! ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 1.00 GC 1.02 HH 1.00* ! HN 1.03 KK 0.90 KP 1.02 LA 1.02 ! MA 0.80 MK 0.99 NL 1.01 NY 1.01 ! OV 1.01 PT 0.99 SC 1.02 TS 1.10 ! WF 0.99 WZ 0.92 GG 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 0.99 GC 0.97 HH 1.00* ! HN 1.22 KK 0.79 KP 1.10 LA 0.92 ! MA 0.85 MK 0.94 NL 0.96 NY 1.08 ! OV 1.05 PT 0.97 SC 1.15 TS 1.01 ! WF 1.03 WZ 1.07 GG 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.96 FD 0.99 GC 0.99 HH 1.00* ! HN 1.19 KK 1.00* KP 1.12 LA 0.94 ! MA 0.84 MK 0.93 NL 0.97 NY 1.10 ! OV 1.08 PT 0.96 SC 1.03 TS 1.00* ! WF 1.13 WZ 1.00* GG 1.00* ! ! A total of 45 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.91 FD 0.97 GC 0.79 HH 1.00* ! HN 1.17 KK 1.00* KP 1.09 LA 0.90 ! MA 0.98 MK 0.93 NL 0.95 NY 1.12 ! OV 1.03 PT 0.88 SC 1.01 TS 1.00* ! WF 1.00 WZ 0.98 GG 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 0.95 GC 1.00* HH 1.00* ! HN 1.20 KK 0.86 KP 1.11 LA 0.96 ! MA 0.98 MK 0.99 NL 0.99 NY 1.07 ! OV 1.11 PT 0.98 SC 1.19 TS 1.03 ! WF 1.11 WZ 0.92 GG 1.00* ! ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00 GC 1.00 HH 1.00* ! HN 0.99 KK 1.01 KP 1.03 LA 1.02 ! MA 0.83 MK 0.98 NL 1.02 NY 1.00 ! OV 1.01 PT 0.99 SC 1.00 TS 1.09 ! WF 0.96 WZ 0.91 GG 1.00* ! ! A total of 18 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 0.89 GC 0.99 HH 1.00* ! HN 1.13 KK 0.84 KP 1.14 LA 0.77 ! MA 0.97 MK 0.83 NL 0.91 NY 1.11 ! OV 0.95 PT 0.98 SC 1.04 TS 1.00 ! WF 1.04 WZ 1.05 GG 1.00* ! ! A total of 15 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.92 FD 0.89 GC 1.02 HH 1.00* ! HN 1.13 KK 1.00* KP 1.12 LA 0.80 ! MA 0.95 MK 0.88 NL 0.97 NY 1.08 ! OV 0.95 PT 0.93 SC 1.05 TS 1.00* ! WF 1.09 WZ 1.00* GG 1.00* ! ! A total of 24 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.91 FD 0.89 GC 1.37 HH 1.00* ! HN 1.06 KK 1.00* KP 1.07 LA 0.78 ! MA 1.00 MK 0.88 NL 0.99 NY 1.06 ! OV 0.95 PT 0.91 SC 1.04 TS 1.00* ! WF 1.00 WZ 0.98 GG 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 0.90 GC 0.81 HH 1.00* ! HN 1.14 KK 0.89 KP 1.12 LA 0.79 ! MA 1.04 MK 0.95 NL 1.04 NY 1.06 ! OV 0.95 PT 0.96 SC 1.03 TS 1.06 ! WF 1.14 WZ 0.93 GG 1.00* ! ! ! Fit before self-cal, rms=0.137783Jy sigma=1.536208 ! Fit after self-cal, rms=0.128667Jy sigma=1.297426 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.446 mas, bmaj=2.649 mas, bpa=-34 degrees ! Estimated noise=0.227384 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.00386062 Jy ! Component: 100 - total flux cleaned = 0.00444544 Jy ! Total flux subtracted in 100 components = 0.00444544 Jy ! Clean residual min=-0.002276 max=0.002870 Jy/beam ! Clean residual mean=0.000015 rms=0.000395 Jy/beam ! Combined flux in latest and established models = 0.626903 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 117 components and 0.626903 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128484Jy sigma=1.295401 ! Fit after self-cal, rms=0.128507Jy sigma=1.294990 ! Inverting map ! Added new window around map position (-0.4, -11.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00198853 Jy ! Component: 100 - total flux cleaned = 0.00236552 Jy ! Total flux subtracted in 100 components = 0.00236552 Jy ! Clean residual min=-0.002209 max=0.001762 Jy/beam ! Clean residual mean=0.000011 rms=0.000380 Jy/beam ! Combined flux in latest and established models = 0.629268 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 127 components and 0.629268 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128403Jy sigma=1.293776 ! Fit after self-cal, rms=0.128415Jy sigma=1.293720 ! 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.000632299 Jy ! Component: 100 - total flux cleaned = 0.000928781 Jy ! Total flux subtracted in 100 components = 0.000928781 Jy ! Clean residual min=-0.002178 max=0.001700 Jy/beam ! Clean residual mean=0.000009 rms=0.000374 Jy/beam ! Combined flux in latest and established models = 0.630197 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 132 components and 0.630197 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128358Jy sigma=1.293168 ! Fit after self-cal, rms=0.133309Jy sigma=1.249534 ! Inverting map and beam ! Estimated beam: bmin=2.444 mas, bmaj=2.651 mas, bpa=-34.25 degrees ! Estimated noise=0.22743 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00168442 Jy ! Component: 100 - total flux cleaned = 0.00189068 Jy ! Total flux subtracted in 100 components = 0.00189068 Jy ! Clean residual min=-0.001478 max=0.001603 Jy/beam ! Clean residual mean=0.000003 rms=0.000331 Jy/beam ! Combined flux in latest and established models = 0.632088 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 141 components and 0.632088 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133255Jy sigma=1.248997 ! Fit after self-cal, rms=0.133625Jy sigma=1.245220 ! Inverting map and beam ! Estimated beam: bmin=2.444 mas, bmaj=2.65 mas, bpa=-34.45 degrees ! Estimated noise=0.227312 mJy/beam. ! Added new window around map position (-13.6, -28). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00125504 Jy ! Component: 100 - total flux cleaned = 0.00175648 Jy ! Total flux subtracted in 100 components = 0.00175648 Jy ! Clean residual min=-0.001422 max=0.001427 Jy/beam ! Clean residual mean=-0.000004 rms=0.000321 Jy/beam ! Combined flux in latest and established models = 0.633844 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 150 components and 0.633844 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133584Jy sigma=1.244757 ! Fit after self-cal, rms=0.134973Jy sigma=1.244039 ! Inverting map and beam ! Estimated beam: bmin=2.443 mas, bmaj=2.65 mas, bpa=-34.61 degrees ! Estimated noise=0.227192 mJy/beam. ! Added new window around map position (-3.2, 26.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00055503 Jy ! Component: 100 - total flux cleaned = 0.000840492 Jy ! Total flux subtracted in 100 components = 0.000840492 Jy ! Clean residual min=-0.001393 max=0.001286 Jy/beam ! Clean residual mean=-0.000006 rms=0.000317 Jy/beam ! Combined flux in latest and established models = 0.634685 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 157 components and 0.634685 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134942Jy sigma=1.243700 ! Fit after self-cal, rms=0.136605Jy sigma=1.243586 ! Inverting map and beam ! Estimated beam: bmin=2.443 mas, bmaj=2.649 mas, bpa=-34.72 degrees ! Estimated noise=0.227062 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 = 4.03748e-05 Jy ! Component: 100 - total flux cleaned = 0.000349188 Jy ! Total flux subtracted in 100 components = 0.000349188 Jy ! Clean residual min=-0.001374 max=0.001274 Jy/beam ! Clean residual mean=-0.000008 rms=0.000313 Jy/beam ! Combined flux in latest and established models = 0.635034 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 31 model components to the UV plane model. ! The established model now contains 161 components and 0.635034 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136580Jy sigma=1.243326 ! Fit after self-cal, rms=0.136511Jy sigma=1.242517 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 0014+813_S ! Writing UV FITS file: 0014+813_S.uvf ! Writing 161 model components to file: 0014+813_S.mod ! wwins: Wrote 7 windows to 0014+813_S.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.443 x 2.649 at -34.72 degrees (North through East) ! Clean map min=-0.0012542 max=0.43289 Jy/beam ! Writing clean map to FITS file: 0014+813_S.fits ! Writing difmap environment to: 0014+813_S.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 0014+813_S] device 0014+813_Suv.ps/vps,1,2 ! Attempting to open device: '0014+813_Suv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 0014+813_Sim.ps/vps ! Attempting to open device: '0014+813_Sim.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=2.443 mas, bmaj=2.649 mas, bpa=-34.72 degrees ! Estimated noise=0.227062 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.443 x 2.649 at -34.72 degrees (North through East) ! Clean map min=-0.0012542 max=0.43289 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo6-0.4] quit ! Quitting program ! Log file 0014+813_S.log closed on Mon Aug 18 11:17:23 2008