! Started logfile: 3C274_S.log on Mon Aug 25 16:13:13 2008 observe 3C274_S.edt ! Reading UV FITS file: 3C274_S.edt ! AN table 1: 75 integrations on 171 of 171 possible baselines. ! AN table 2: 90 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.119121 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 3C274 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 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 1615 lines of history. ! ! Reading 13444 visibilities. ![@gonzo5-0.4 3C274_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 = 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.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 998 telescope corrections were flagged in sub-array 1. ! A total of 1478 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 986 telescope corrections were flagged in sub-array 1. ! A total of 1486 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 989 telescope corrections were flagged in sub-array 1. ! A total of 1480 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 986 telescope corrections were flagged in sub-array 1. ! A total of 1484 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.585752Jy sigma=7.250388 ! Fit after self-cal, rms=0.582643Jy sigma=7.263975 ! 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=2.169 mas, bmaj=3.82 mas, bpa=2.918 degrees ! Estimated noise=0.967009 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.679743 Jy ! Total flux subtracted in 50 components = 0.679743 Jy ! Clean residual min=-0.053219 max=0.257136 Jy/beam ! Clean residual mean=-0.000301 rms=0.015078 Jy/beam ! Combined flux in latest and established models = 0.679743 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 10 components and 0.679743 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.598861Jy sigma=7.688690 ! Fit after self-cal, rms=0.598399Jy sigma=7.683820 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.271396 Jy ! Total flux subtracted in 50 components = 0.271396 Jy ! Clean residual min=-0.034371 max=0.144167 Jy/beam ! Clean residual mean=-0.000240 rms=0.010765 Jy/beam ! Combined flux in latest and established models = 0.951139 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 11 model components to the UV plane model. ! The established model now contains 17 components and 0.951139 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.467306Jy sigma=5.772116 ! Fit after self-cal, rms=0.466784Jy sigma=5.767122 ! Inverting map ! Added new window around map position (-2.8, 0.4). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.152801 Jy ! Total flux subtracted in 50 components = 0.152801 Jy ! Clean residual min=-0.028585 max=0.091172 Jy/beam ! Clean residual mean=-0.000205 rms=0.009073 Jy/beam ! Combined flux in latest and established models = 1.10394 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 26 components and 1.10394 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.410003Jy sigma=4.948411 ! Fit after self-cal, rms=0.408962Jy sigma=4.934918 ! Inverting map ! Added new window around map position (-6, 1.6). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.105817 Jy ! Total flux subtracted in 50 components = 0.105817 Jy ! Clean residual min=-0.026093 max=0.074983 Jy/beam ! Clean residual mean=-0.000178 rms=0.008116 Jy/beam ! Combined flux in latest and established models = 1.20976 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 40 components and 1.20976 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.375586Jy sigma=4.466266 ! Fit after self-cal, rms=0.374787Jy sigma=4.451908 ! Inverting map ! Added new window around map position (-12, 3.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0862349 Jy ! Total flux subtracted in 50 components = 0.0862349 Jy ! Clean residual min=-0.023540 max=0.059046 Jy/beam ! Clean residual mean=-0.000158 rms=0.007395 Jy/beam ! Combined flux in latest and established models = 1.29599 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 1.29599 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.350177Jy sigma=4.119457 ! Fit after self-cal, rms=0.349568Jy sigma=4.110054 ! Inverting map ! Added new window around map position (-8.8, 2.8). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.070622 Jy ! Total flux subtracted in 50 components = 0.070622 Jy ! Clean residual min=-0.022754 max=0.054495 Jy/beam ! Clean residual mean=-0.000140 rms=0.006888 Jy/beam ! Combined flux in latest and established models = 1.36661 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 65 components and 1.36661 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.331426Jy sigma=3.865494 ! Fit after self-cal, rms=0.331119Jy sigma=3.859916 ! Inverting map ! Added new window around map position (-18.8, 7.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0627162 Jy ! Total flux subtracted in 50 components = 0.0627162 Jy ! Clean residual min=-0.021535 max=0.044760 Jy/beam ! Clean residual mean=-0.000131 rms=0.006463 Jy/beam ! Combined flux in latest and established models = 1.42933 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 80 components and 1.42933 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316132Jy sigma=3.665661 ! Fit after self-cal, rms=0.315740Jy sigma=3.660251 ! Inverting map ! Added new window around map position (-15.2, 4.4). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.054413 Jy ! Total flux subtracted in 50 components = 0.054413 Jy ! Clean residual min=-0.020606 max=0.040730 Jy/beam ! Clean residual mean=-0.000122 rms=0.006114 Jy/beam ! Combined flux in latest and established models = 1.48374 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 26 model components to the UV plane model. ! The established model now contains 91 components and 1.48374 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.303555Jy sigma=3.504187 ! Fit after self-cal, rms=0.303311Jy sigma=3.500512 ! Inverting map ! Added new window around map position (3.2, -1.2). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0486306 Jy ! Total flux subtracted in 50 components = 0.0486306 Jy ! Clean residual min=-0.019749 max=0.029802 Jy/beam ! Clean residual mean=-0.000112 rms=0.005831 Jy/beam ! Combined flux in latest and established models = 1.53237 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 29 model components to the UV plane model. ! The established model now contains 105 components and 1.53237 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.293085Jy sigma=3.367324 ! Fit after self-cal, rms=0.292932Jy sigma=3.365704 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0425971 Jy ! Total flux subtracted in 50 components = 0.0425971 Jy ! Clean residual min=-0.019265 max=0.028216 Jy/beam ! Clean residual mean=-0.000105 rms=0.005607 Jy/beam ! Combined flux in latest and established models = 1.57497 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 38 model components to the UV plane model. ! The established model now contains 124 components and 1.57497 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.284696Jy sigma=3.259888 ! Fit after self-cal, rms=0.284589Jy sigma=3.258551 ! Inverting map ! Added new window around map position (-23.6, 8). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0386294 Jy ! Total flux subtracted in 50 components = 0.0386294 Jy ! Clean residual min=-0.018806 max=0.024492 Jy/beam ! Clean residual mean=-0.000100 rms=0.005415 Jy/beam ! Combined flux in latest and established models = 1.6136 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 37 model components to the UV plane model. ! The established model now contains 147 components and 1.6136 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.277654Jy sigma=3.170575 ! Fit after self-cal, rms=0.277545Jy sigma=3.169306 ! 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.394 mas, bmaj=4.411 mas, bpa=-1.208 degrees ! Estimated noise=0.641615 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.0672212 Jy ! Component: 100 - total flux cleaned = 0.120415 Jy ! Total flux subtracted in 100 components = 0.120415 Jy ! Clean residual min=-0.026064 max=0.032412 Jy/beam ! Clean residual mean=0.000654 rms=0.007386 Jy/beam ! Combined flux in latest and established models = 1.73401 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 59 model components to the UV plane model. ! The established model now contains 195 components and 1.73401 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.258316Jy sigma=2.940226 ! Fit after self-cal, rms=0.257962Jy sigma=2.930288 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0471797 Jy ! Component: 100 - total flux cleaned = 0.0860641 Jy ! Total flux subtracted in 100 components = 0.0860641 Jy ! Clean residual min=-0.023039 max=0.027825 Jy/beam ! Clean residual mean=0.000374 rms=0.006494 Jy/beam ! Combined flux in latest and established models = 1.82008 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 69 model components to the UV plane model. ! The established model now contains 241 components and 1.82008 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.247376Jy sigma=2.812342 ! Fit after self-cal, rms=0.247246Jy sigma=2.805765 ! Inverting map ! Added new window around map position (-27.2, 10.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0362749 Jy ! Component: 100 - total flux cleaned = 0.0662598 Jy ! Total flux subtracted in 100 components = 0.0662598 Jy ! Clean residual min=-0.021470 max=0.025688 Jy/beam ! Clean residual mean=0.000146 rms=0.005922 Jy/beam ! Combined flux in latest and established models = 1.88634 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 64 model components to the UV plane model. ! The established model now contains 278 components and 1.88634 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.240587Jy sigma=2.735808 ! Fit after self-cal, rms=0.240553Jy sigma=2.731130 ! Inverting map ! Added new window around map position (-33.6, 15.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0302693 Jy ! Component: 100 - total flux cleaned = 0.0551898 Jy ! Total flux subtracted in 100 components = 0.0551898 Jy ! Clean residual min=-0.021072 max=0.021040 Jy/beam ! Clean residual mean=-0.000065 rms=0.005496 Jy/beam ! Combined flux in latest and established models = 1.94153 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 64 model components to the UV plane model. ! The established model now contains 313 components and 1.94153 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.235804Jy sigma=2.678882 ! Fit after self-cal, rms=0.235828Jy sigma=2.675938 ! 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 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 1.14 GC 1.12 HH 1.00* ! HN 0.89 KK 0.85 KP 0.91 LA 1.09 ! MA 1.00* MK 0.84 NL 0.99 NY 1.25 ! OV 1.03 PT 1.03 SC 0.78 TS 1.14 ! WF 0.96 WZ 1.00* GG 1.00* ! ! A total of 42 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.05 GC 1.00* HH 0.93 ! HN 0.93 KK 0.91 KP 0.90 LA 1.00 ! MA 1.00* MK 0.95 NL 0.95 NY 1.28 ! OV 1.05 PT 1.00 SC 0.75 TS 1.00* ! WF 1.08 WZ 1.04 GG 1.00* ! ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.88 FD 1.02 GC 0.99 HH 1.00* ! HN 0.96 KK 0.86 KP 1.02 LA 1.08 ! MA 0.70 MK 0.92 NL 1.12 NY 1.29 ! OV 0.95 PT 0.89 SC 1.00 TS 1.17 ! WF 0.99 WZ 1.00* GG 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.92 FD 0.95 GC 1.00* HH 0.95 ! HN 0.97 KK 0.93 KP 1.03 LA 1.03 ! MA 1.00* MK 1.04 NL 1.08 NY 1.21 ! OV 0.97 PT 0.88 SC 0.96 TS 1.00* ! WF 1.07 WZ 1.07 GG 1.00* ! ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.00 GC 1.08 HH 1.00* ! HN 1.19 KK 0.87 KP 1.14 LA 0.90 ! MA 1.48 MK 0.92 NL 1.08 NY 1.31 ! OV 1.02 PT 0.95 SC 1.07 TS 1.06 ! WF 1.00 WZ 1.00* GG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 0.98 GC 1.00* HH 0.91 ! HN 1.21 KK 0.94 KP 1.13 LA 0.85 ! MA 1.00* MK 1.04 NL 1.03 NY 1.27 ! OV 1.09 PT 0.92 SC 1.00 TS 1.00* ! WF 1.13 WZ 0.94 GG 1.00* ! ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 0.89 GC 1.06 HH 1.00* ! HN 1.15 KK 0.90 KP 1.12 LA 0.75 ! MA 1.37 MK 0.86 NL 1.04 NY 1.30 ! OV 0.95 PT 0.90 SC 1.02 TS 1.05 ! WF 0.97 WZ 1.00* GG 1.00* ! ! A total of 36 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 0.86 GC 1.00* HH 0.88 ! HN 1.14 KK 0.93 KP 1.11 LA 0.80 ! MA 1.00* MK 0.99 NL 0.99 NY 1.32 ! OV 0.96 PT 0.90 SC 0.98 TS 1.00* ! WF 1.11 WZ 0.96 GG 1.00* ! ! ! Fit before self-cal, rms=0.235828Jy sigma=2.675938 ! Fit after self-cal, rms=0.169514Jy sigma=1.768892 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.405 mas, bmaj=4.521 mas, bpa=-0.3966 degrees ! Estimated noise=0.640345 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 (-32, 8.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0232748 Jy ! Component: 100 - total flux cleaned = 0.0391553 Jy ! Total flux subtracted in 100 components = 0.0391553 Jy ! Clean residual min=-0.011039 max=0.013308 Jy/beam ! Clean residual mean=0.000757 rms=0.002951 Jy/beam ! Combined flux in latest and established models = 1.98068 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 35 model components to the UV plane model. ! The established model now contains 345 components and 1.98068 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166010Jy sigma=1.729647 ! Fit after self-cal, rms=0.166101Jy sigma=1.724096 ! Inverting map ! Added new window around map position (-58.8, 24). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0161239 Jy ! Component: 100 - total flux cleaned = 0.0290585 Jy ! Total flux subtracted in 100 components = 0.0290585 Jy ! Clean residual min=-0.010693 max=0.012245 Jy/beam ! Clean residual mean=0.000636 rms=0.002720 Jy/beam ! Combined flux in latest and established models = 2.00974 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 43 model components to the UV plane model. ! The established model now contains 372 components and 2.00974 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.164076Jy sigma=1.702888 ! Fit after self-cal, rms=0.164148Jy sigma=1.701671 ! Inverting map ! Added new window around map position (-42.4, 18.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0139586 Jy ! Component: 100 - total flux cleaned = 0.0241799 Jy ! Total flux subtracted in 100 components = 0.0241799 Jy ! Clean residual min=-0.010202 max=0.012283 Jy/beam ! Clean residual mean=0.000533 rms=0.002538 Jy/beam ! Combined flux in latest and established models = 2.03392 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 399 components and 2.03392 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162575Jy sigma=1.685241 ! Fit after self-cal, rms=0.162625Jy sigma=1.684499 ! Inverting map ! Added new window around map position (-74, 32). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0129055 Jy ! Component: 100 - total flux cleaned = 0.0227626 Jy ! Total flux subtracted in 100 components = 0.0227626 Jy ! Clean residual min=-0.010224 max=0.010217 Jy/beam ! Clean residual mean=0.000436 rms=0.002381 Jy/beam ! Combined flux in latest and established models = 2.05668 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 49 model components to the UV plane model. ! The established model now contains 427 components and 2.05668 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161276Jy sigma=1.667378 ! Fit after self-cal, rms=0.161303Jy sigma=1.666729 ! Inverting map ! Added new window around map position (-65.2, 29.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.010296 Jy ! Component: 100 - total flux cleaned = 0.0187729 Jy ! Total flux subtracted in 100 components = 0.0187729 Jy ! Clean residual min=-0.009824 max=0.008908 Jy/beam ! Clean residual mean=0.000354 rms=0.002261 Jy/beam ! Combined flux in latest and established models = 2.07546 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 54 model components to the UV plane model. ! The established model now contains 460 components and 2.07546 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160199Jy sigma=1.655395 ! Fit after self-cal, rms=0.160243Jy sigma=1.654977 ! 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.00774991 Jy ! Component: 100 - total flux cleaned = 0.0145941 Jy ! Total flux subtracted in 100 components = 0.0145941 Jy ! Clean residual min=-0.009634 max=0.008360 Jy/beam ! Clean residual mean=0.000293 rms=0.002177 Jy/beam ! Combined flux in latest and established models = 2.09005 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 60 model components to the UV plane model. ! The established model now contains 477 components and 2.09005 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.159455Jy sigma=1.647296 ! Fit after self-cal, rms=0.156406Jy sigma=1.555067 ! Inverting map and beam ! Estimated beam: bmin=2.409 mas, bmaj=4.522 mas, bpa=-0.5474 degrees ! Estimated noise=0.639375 mJy/beam. ! Added new window around map position (-37.6, 16.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0103933 Jy ! Component: 100 - total flux cleaned = 0.0179766 Jy ! Total flux subtracted in 100 components = 0.0179766 Jy ! Clean residual min=-0.006284 max=0.007817 Jy/beam ! Clean residual mean=0.000450 rms=0.001810 Jy/beam ! Combined flux in latest and established models = 2.10803 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 45 model components to the UV plane model. ! The established model now contains 506 components and 2.10803 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155462Jy sigma=1.546041 ! Fit after self-cal, rms=0.156237Jy sigma=1.539657 ! Inverting map and beam ! Estimated beam: bmin=2.409 mas, bmaj=4.516 mas, bpa=-0.5384 degrees ! Estimated noise=0.639219 mJy/beam. ! Added new window around map position (-47.6, 24). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00912541 Jy ! Component: 100 - total flux cleaned = 0.0156832 Jy ! Total flux subtracted in 100 components = 0.0156832 Jy ! Clean residual min=-0.006079 max=0.006915 Jy/beam ! Clean residual mean=0.000404 rms=0.001706 Jy/beam ! Combined flux in latest and established models = 2.12371 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 42 model components to the UV plane model. ! The established model now contains 530 components and 2.12371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155470Jy sigma=1.532318 ! Fit after self-cal, rms=0.156178Jy sigma=1.530803 ! Inverting map and beam ! Estimated beam: bmin=2.41 mas, bmaj=4.514 mas, bpa=-0.5307 degrees ! Estimated noise=0.639272 mJy/beam. ! Added new window around map position (-66.4, 11.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00820605 Jy ! Component: 100 - total flux cleaned = 0.0142088 Jy ! Total flux subtracted in 100 components = 0.0142088 Jy ! Clean residual min=-0.005645 max=0.006768 Jy/beam ! Clean residual mean=0.000363 rms=0.001628 Jy/beam ! Combined flux in latest and established models = 2.13792 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 48 model components to the UV plane model. ! The established model now contains 549 components and 2.13792 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155570Jy sigma=1.524291 ! Fit after self-cal, rms=0.156035Jy sigma=1.523453 ! Inverting map and beam ! Estimated beam: bmin=2.41 mas, bmaj=4.511 mas, bpa=-0.5068 degrees ! Estimated noise=0.639373 mJy/beam. ! Added new window around map position (-80, 26.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0075438 Jy ! Component: 100 - total flux cleaned = 0.0133833 Jy ! Total flux subtracted in 100 components = 0.0133833 Jy ! Clean residual min=-0.005522 max=0.005963 Jy/beam ! Clean residual mean=0.000327 rms=0.001565 Jy/beam ! Combined flux in latest and established models = 2.1513 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 56 model components to the UV plane model. ! The established model now contains 565 components and 2.1513 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155492Jy sigma=1.518303 ! Fit after self-cal, rms=0.155771Jy sigma=1.517521 ! Inverting map and beam ! Estimated beam: bmin=2.41 mas, bmaj=4.509 mas, bpa=-0.4873 degrees ! Estimated noise=0.639384 mJy/beam. ! Added new window around map position (101.6, 65.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00655878 Jy ! Component: 100 - total flux cleaned = 0.0119618 Jy ! Total flux subtracted in 100 components = 0.0119618 Jy ! Clean residual min=-0.005403 max=0.005822 Jy/beam ! Clean residual mean=0.000303 rms=0.001510 Jy/beam ! Combined flux in latest and established models = 2.16326 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 56 model components to the UV plane model. ! The established model now contains 581 components and 2.16326 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155320Jy sigma=1.512926 ! Fit after self-cal, rms=0.155476Jy sigma=1.512156 ! Inverting map and beam ! Estimated beam: bmin=2.41 mas, bmaj=4.507 mas, bpa=-0.4659 degrees ! Estimated noise=0.639346 mJy/beam. ! Added new window around map position (90.8, 18). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00659036 Jy ! Component: 100 - total flux cleaned = 0.0114681 Jy ! Total flux subtracted in 100 components = 0.0114681 Jy ! Clean residual min=-0.005383 max=0.005316 Jy/beam ! Clean residual mean=0.000280 rms=0.001463 Jy/beam ! Combined flux in latest and established models = 2.17473 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 57 model components to the UV plane model. ! The established model now contains 599 components and 2.17473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.155085Jy sigma=1.508054 ! Fit after self-cal, rms=0.155149Jy sigma=1.507153 ! Inverting map and beam ! Estimated beam: bmin=2.41 mas, bmaj=4.505 mas, bpa=-0.4449 degrees ! Estimated noise=0.639253 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.00527836 Jy ! Component: 100 - total flux cleaned = 0.00980868 Jy ! Total flux subtracted in 100 components = 0.00980868 Jy ! Clean residual min=-0.005321 max=0.005232 Jy/beam ! Clean residual mean=0.000265 rms=0.001419 Jy/beam ! Combined flux in latest and established models = 2.18454 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 69 model components to the UV plane model. ! The established model now contains 615 components and 2.18454 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154818Jy sigma=1.504071 ! Fit after self-cal, rms=0.154304Jy sigma=1.498381 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 3C274_S ! Writing UV FITS file: 3C274_S.uvf ! Writing 615 model components to file: 3C274_S.mod ! wwins: Wrote 22 windows to 3C274_S.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.41 x 4.505 at -0.4449 degrees (North through East) ! Clean map min=-0.0049635 max=0.79436 Jy/beam ! Writing clean map to FITS file: 3C274_S.fits ! Writing difmap environment to: 3C274_S.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 3C274_S] device 3C274_Suv.ps/vps,1,2 ! Attempting to open device: '3C274_Suv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 3C274_Sim.ps/vps ! Attempting to open device: '3C274_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.41 mas, bmaj=4.505 mas, bpa=-0.4449 degrees ! Estimated noise=0.639253 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.41 x 4.505 at -0.4449 degrees (North through East) ! Clean map min=-0.0049635 max=0.79436 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo5-0.4] quit ! Quitting program ! Log file 3C274_S.log closed on Mon Aug 25 16:13:36 2008