! Started logfile: VR422201_X.log on Tue Dec 12 21:30:38 2006 observe VR422201_X.edt ! Reading UV FITS file: VR422201_X.edt ! AN table 1: 18 integrations on 153 of 153 possible baselines. ! AN table 2: 217 integrations on 153 of 153 possible baselines. ! AN table 3: 32 integrations on 153 of 153 possible baselines. ! AN table 4: 15 integrations on 153 of 153 possible baselines. ! Apparent sampling: 0.182103 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: VR422201 ! ! 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 1304 lines of history. ! ! Reading 31428 visibilities. ![@gonzo6-0.1 VR422201_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 = 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.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 302 telescope corrections were flagged in sub-array 1. ! A total of 2680 telescope corrections were flagged in sub-array 2. ! A total of 407 telescope corrections were flagged in sub-array 3. ! A total of 255 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 299 telescope corrections were flagged in sub-array 1. ! A total of 2683 telescope corrections were flagged in sub-array 2. ! A total of 428 telescope corrections were flagged in sub-array 3. ! A total of 261 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 308 telescope corrections were flagged in sub-array 1. ! A total of 2693 telescope corrections were flagged in sub-array 2. ! A total of 427 telescope corrections were flagged in sub-array 3. ! A total of 249 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 305 telescope corrections were flagged in sub-array 1. ! A total of 2698 telescope corrections were flagged in sub-array 2. ! A total of 423 telescope corrections were flagged in sub-array 3. ! A total of 252 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.357095Jy sigma=0.005467 ! Fit after self-cal, rms=0.356549Jy sigma=0.005473 ! 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.4893 mas, bmaj=0.759 mas, bpa=-13.32 degrees ! Estimated noise=695.032 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.6696 Jy ! Total flux subtracted in 50 components = 0.6696 Jy ! Clean residual min=-0.032386 max=0.187297 Jy/beam ! Clean residual mean=0.000006 rms=0.008935 Jy/beam ! Combined flux in latest and established models = 0.6696 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 2 model components to the UV plane model. ! The established model now contains 2 components and 0.6696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.460707Jy sigma=0.007322 ! Fit after self-cal, rms=0.460697Jy sigma=0.007322 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.178409 Jy ! Total flux subtracted in 50 components = 0.178409 Jy ! Clean residual min=-0.025426 max=0.116740 Jy/beam ! Clean residual mean=0.000007 rms=0.006804 Jy/beam ! Combined flux in latest and established models = 0.848009 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.848009 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.367857Jy sigma=0.005855 ! Fit after self-cal, rms=0.367764Jy sigma=0.005853 ! Inverting map ! Added new window around map position (-0.4, -1.4). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.113452 Jy ! Total flux subtracted in 50 components = 0.113452 Jy ! Clean residual min=-0.022317 max=0.057181 Jy/beam ! Clean residual mean=0.000011 rms=0.005782 Jy/beam ! Combined flux in latest and established models = 0.961461 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 16 components and 0.961461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317727Jy sigma=0.005037 ! Fit after self-cal, rms=0.314990Jy sigma=0.005018 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0769466 Jy ! Total flux subtracted in 50 components = 0.0769466 Jy ! Clean residual min=-0.020320 max=0.043493 Jy/beam ! Clean residual mean=0.000014 rms=0.005128 Jy/beam ! Combined flux in latest and established models = 1.03841 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 26 components and 1.03841 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.285297Jy sigma=0.004520 ! Fit after self-cal, rms=0.284111Jy sigma=0.004513 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0570332 Jy ! Total flux subtracted in 50 components = 0.0570332 Jy ! Clean residual min=-0.019154 max=0.042989 Jy/beam ! Clean residual mean=0.000017 rms=0.004758 Jy/beam ! Combined flux in latest and established models = 1.09544 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 22 model components to the UV plane model. ! The established model now contains 41 components and 1.09544 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.264714Jy sigma=0.004177 ! Fit after self-cal, rms=0.264000Jy sigma=0.004173 ! Inverting map ! Added new window around map position (-0.9, -3.5). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0521687 Jy ! Total flux subtracted in 50 components = 0.0521687 Jy ! Clean residual min=-0.017696 max=0.031501 Jy/beam ! Clean residual mean=0.000017 rms=0.004405 Jy/beam ! Combined flux in latest and established models = 1.14761 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.14761 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.247808Jy sigma=0.003882 ! Fit after self-cal, rms=0.246930Jy sigma=0.003878 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0439287 Jy ! Total flux subtracted in 50 components = 0.0439287 Jy ! Clean residual min=-0.016201 max=0.028443 Jy/beam ! Clean residual mean=0.000019 rms=0.004121 Jy/beam ! Combined flux in latest and established models = 1.19154 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 59 components and 1.19154 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.233922Jy sigma=0.003640 ! Fit after self-cal, rms=0.233391Jy sigma=0.003638 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0386123 Jy ! Total flux subtracted in 50 components = 0.0386123 Jy ! Clean residual min=-0.015040 max=0.026127 Jy/beam ! Clean residual mean=0.000020 rms=0.003902 Jy/beam ! Combined flux in latest and established models = 1.23015 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 64 components and 1.23015 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.222576Jy sigma=0.003437 ! Fit after self-cal, rms=0.222147Jy sigma=0.003436 ! Inverting map ! Added new window around map position (-1, -2.6). ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0348838 Jy ! Total flux subtracted in 50 components = 0.0348838 Jy ! Clean residual min=-0.014160 max=0.021976 Jy/beam ! Clean residual mean=0.000021 rms=0.003713 Jy/beam ! Combined flux in latest and established models = 1.26503 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 76 components and 1.26503 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.212886Jy sigma=0.003260 ! Fit after self-cal, rms=0.212504Jy sigma=0.003259 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0310747 Jy ! Total flux subtracted in 50 components = 0.0310747 Jy ! Clean residual min=-0.013427 max=0.020617 Jy/beam ! Clean residual mean=0.000021 rms=0.003554 Jy/beam ! Combined flux in latest and established models = 1.29611 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 81 components and 1.29611 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.204768Jy sigma=0.003109 ! Fit after self-cal, rms=0.204440Jy sigma=0.003108 ! 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.6004 mas, bmaj=0.9845 mas, bpa=-13.78 degrees ! Estimated noise=341.879 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.0797293 Jy ! Component: 100 - total flux cleaned = 0.13818 Jy ! Total flux subtracted in 100 components = 0.13818 Jy ! Clean residual min=-0.016855 max=0.040341 Jy/beam ! Clean residual mean=0.001369 rms=0.006509 Jy/beam ! Combined flux in latest and established models = 1.43429 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 63 model components to the UV plane model. ! The established model now contains 117 components and 1.43429 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176855Jy sigma=0.002537 ! Fit after self-cal, rms=0.174812Jy sigma=0.002517 ! Inverting map ! Added new window around map position (-0.1, -4.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0519795 Jy ! Component: 100 - total flux cleaned = 0.0920773 Jy ! Total flux subtracted in 100 components = 0.0920773 Jy ! Clean residual min=-0.014324 max=0.029260 Jy/beam ! Clean residual mean=0.001035 rms=0.005236 Jy/beam ! Combined flux in latest and established models = 1.52637 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 62 model components to the UV plane model. ! The established model now contains 161 components and 1.52637 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161345Jy sigma=0.002214 ! Fit after self-cal, rms=0.160044Jy sigma=0.002200 ! Inverting map ! Added new window around map position (0.5, -1.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0368493 Jy ! Component: 100 - total flux cleaned = 0.0659441 Jy ! Total flux subtracted in 100 components = 0.0659441 Jy ! Clean residual min=-0.012368 max=0.025449 Jy/beam ! Clean residual mean=0.000797 rms=0.004482 Jy/beam ! Combined flux in latest and established models = 1.59231 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 62 model components to the UV plane model. ! The established model now contains 206 components and 1.59231 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.152792Jy sigma=0.002024 ! Fit after self-cal, rms=0.152152Jy sigma=0.002017 ! Inverting map ! Added new window around map position (2.2, -6.1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0288219 Jy ! Component: 100 - total flux cleaned = 0.0522445 Jy ! Total flux subtracted in 100 components = 0.0522445 Jy ! Clean residual min=-0.012572 max=0.020083 Jy/beam ! Clean residual mean=0.000613 rms=0.003952 Jy/beam ! Combined flux in latest and established models = 1.64455 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 52 model components to the UV plane model. ! The established model now contains 235 components and 1.64455 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147210Jy sigma=0.001902 ! Fit after self-cal, rms=0.146876Jy sigma=0.001898 ! 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 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00* GC 1.01 HH 1.00* HN 1.00* ! LA 1.00* MC 0.99 MK 1.00* NY 1.00* ! ON 1.00* OV 1.00* PT 1.00* SC 1.00* ! TS 0.98 WF 1.12 BR 1.00* KK 1.00* ! KP 1.00* WZ 1.00* ! ! A total of 120 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! FD 0.99 GC 0.97 HH 1.00* HN 1.00 ! LA 0.95 MC 1.13 MK 0.99 NY 0.98 ! ON 1.01 OV 0.98 PT 0.95 SC 1.00* ! TS 1.03 WF 0.96 BR 1.01 KK 1.13 ! KP 0.96 WZ 0.93 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.00* GC 1.00 HH 1.23 HN 0.96 ! LA 1.06 MC 0.97 MK 1.00* NY 0.99 ! ON 0.93 OV 1.00* PT 1.00* SC 1.00* ! TS 1.00* WF 0.92 BR 1.00 KK 1.00* ! KP 1.00* WZ 0.94 ! ! A total of 15 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* GC 0.00* HH 0.00* HN 0.00* ! LA 0.00* MC 0.00* MK 0.00* NY 0.00* ! ON 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* BR 0.00* KK 0.00* ! KP 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00* GC 1.03 HH 1.00* HN 1.00* ! LA 1.00* MC 0.95 MK 1.00* NY 1.00* ! ON 1.00* OV 1.00* PT 1.00* SC 1.00* ! TS 1.00 WF 1.16 BR 1.00* KK 1.00* ! KP 1.00* WZ 1.00* ! ! A total of 117 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! FD 0.99 GC 0.99 HH 1.00* HN 0.99 ! LA 0.95 MC 1.13 MK 0.99 NY 0.98 ! ON 1.02 OV 0.97 PT 0.98 SC 1.00* ! TS 1.03 WF 0.98 BR 1.02 KK 1.14 ! KP 0.92 WZ 0.92 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.00* GC 1.02 HH 1.00* HN 0.96 ! LA 1.07 MC 0.99 MK 1.00* NY 0.97 ! ON 0.91 OV 1.00* PT 1.00* SC 1.00* ! TS 1.00* WF 0.92 BR 1.01 KK 1.00* ! KP 1.00* WZ 0.91 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* GC 0.00* HH 0.00* HN 0.00* ! LA 0.00* MC 0.00* MK 0.00* NY 0.00* ! ON 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* BR 0.00* KK 0.00* ! KP 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00* GC 1.05 HH 1.00* HN 1.00* ! LA 1.00* MC 0.98 MK 1.00* NY 1.00* ! ON 1.00* OV 1.00* PT 1.00* SC 1.00* ! TS 1.00 WF 1.23 BR 1.00* KK 1.00* ! KP 1.00* WZ 1.00* ! ! A total of 111 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! FD 0.97 GC 1.01 HH 1.00* HN 1.00 ! LA 0.96 MC 1.16 MK 1.01 NY 1.00 ! ON 1.02 OV 1.00 PT 1.01 SC 1.00* ! TS 1.02 WF 1.11 BR 1.05 KK 1.13 ! KP 0.92 WZ 0.93 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.00* GC 1.06 HH 1.29 HN 0.97 ! LA 1.05 MC 1.01 MK 1.00* NY 1.00 ! ON 0.95 OV 1.00* PT 1.00* SC 1.00* ! TS 1.00* WF 1.06 BR 1.01 KK 1.00* ! KP 1.00* WZ 0.92 ! ! A total of 21 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* GC 0.00* HH 0.00* HN 0.00* ! LA 0.00* MC 0.00* MK 0.00* NY 0.00* ! ON 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* BR 0.00* KK 0.00* ! KP 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00* GC 1.11 HH 1.00* HN 1.00* ! LA 1.00* MC 0.97 MK 1.00* NY 1.00* ! ON 1.00* OV 1.00* PT 1.00* SC 1.00* ! TS 1.03 WF 1.26 BR 1.00* KK 1.00* ! KP 1.00* WZ 1.00* ! ! A total of 132 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! FD 0.96 GC 1.05 HH 1.00* HN 1.00 ! LA 0.94 MC 1.24 MK 1.04 NY 1.04 ! ON 1.02 OV 0.99 PT 1.00 SC 1.00* ! TS 1.03 WF 1.10 BR 1.02 KK 1.14 ! KP 0.94 WZ 0.92 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.00* GC 1.08 HH 1.43 HN 0.99 ! LA 1.04 MC 1.08 MK 1.00* NY 1.01 ! ON 0.94 OV 1.00* PT 1.00* SC 1.00* ! TS 1.00* WF 1.04 BR 1.00 KK 1.00* ! KP 1.00* WZ 0.95 ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! FD 0.00* GC 0.00* HH 0.00* HN 0.00* ! LA 0.00* MC 0.00* MK 0.00* NY 0.00* ! ON 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* BR 0.00* KK 0.00* ! KP 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.146876Jy sigma=0.001898 ! Fit after self-cal, rms=0.126358Jy sigma=0.001524 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6031 mas, bmaj=0.9944 mas, bpa=-14.55 degrees ! Estimated noise=340.64 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 (1.1, -5.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0189631 Jy ! Component: 100 - total flux cleaned = 0.0323604 Jy ! Total flux subtracted in 100 components = 0.0323604 Jy ! Clean residual min=-0.006469 max=0.011091 Jy/beam ! Clean residual mean=0.000649 rms=0.002066 Jy/beam ! Combined flux in latest and established models = 1.67691 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 32 model components to the UV plane model. ! The established model now contains 264 components and 1.67692 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123987Jy sigma=0.001462 ! Fit after self-cal, rms=0.123751Jy sigma=0.001459 ! Inverting map ! Added new window around map position (-1.7, -3.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0128842 Jy ! Component: 100 - total flux cleaned = 0.0222952 Jy ! Total flux subtracted in 100 components = 0.0222952 Jy ! Clean residual min=-0.006279 max=0.009025 Jy/beam ! Clean residual mean=0.000573 rms=0.001865 Jy/beam ! Combined flux in latest and established models = 1.69921 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 34 model components to the UV plane model. ! The established model now contains 287 components and 1.69921 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122523Jy sigma=0.001428 ! Fit after self-cal, rms=0.122368Jy sigma=0.001427 ! 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 (4.6, -9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00838908 Jy ! Component: 100 - total flux cleaned = 0.0133337 Jy ! Total flux subtracted in 100 components = 0.0133337 Jy ! Clean residual min=-0.006292 max=0.007397 Jy/beam ! Clean residual mean=0.000525 rms=0.001736 Jy/beam ! Combined flux in latest and established models = 1.71254 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 39 model components to the UV plane model. ! The established model now contains 306 components and 1.71254 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121634Jy sigma=0.001407 ! Fit after self-cal, rms=0.112231Jy sigma=0.001271 ! Inverting map and beam ! Estimated beam: bmin=0.6009 mas, bmaj=0.9966 mas, bpa=-14.19 degrees ! Estimated noise=340.271 mJy/beam. ! Added new window around map position (-0.7, -6.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00927161 Jy ! Component: 100 - total flux cleaned = 0.0162323 Jy ! Total flux subtracted in 100 components = 0.0162323 Jy ! Clean residual min=-0.004167 max=0.006363 Jy/beam ! Clean residual mean=0.000419 rms=0.001180 Jy/beam ! Combined flux in latest and established models = 1.72878 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 329 components and 1.72878 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.111517Jy sigma=0.001253 ! Fit after self-cal, rms=0.111199Jy sigma=0.001249 ! Inverting map and beam ! Estimated beam: bmin=0.6008 mas, bmaj=0.9965 mas, bpa=-14.17 degrees ! Estimated noise=340.305 mJy/beam. ! Added new window around map position (0.9, -9.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00762572 Jy ! Component: 100 - total flux cleaned = 0.0130581 Jy ! Total flux subtracted in 100 components = 0.0130581 Jy ! Clean residual min=-0.003825 max=0.005304 Jy/beam ! Clean residual mean=0.000376 rms=0.001071 Jy/beam ! Combined flux in latest and established models = 1.74183 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 39 model components to the UV plane model. ! The established model now contains 355 components and 1.74183 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110670Jy sigma=0.001237 ! Fit after self-cal, rms=0.110582Jy sigma=0.001235 ! Inverting map and beam ! Estimated beam: bmin=0.6007 mas, bmaj=0.9965 mas, bpa=-14.16 degrees ! Estimated noise=340.395 mJy/beam. ! Added new window around map position (-2.9, -5.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00645342 Jy ! Component: 100 - total flux cleaned = 0.0112393 Jy ! Total flux subtracted in 100 components = 0.0112393 Jy ! Clean residual min=-0.003510 max=0.004981 Jy/beam ! Clean residual mean=0.000350 rms=0.000990 Jy/beam ! Combined flux in latest and established models = 1.75307 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 46 model components to the UV plane model. ! The established model now contains 381 components and 1.75307 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110191Jy sigma=0.001226 ! Fit after self-cal, rms=0.110182Jy sigma=0.001225 ! Inverting map and beam ! Estimated beam: bmin=0.6007 mas, bmaj=0.9965 mas, bpa=-14.15 degrees ! Estimated noise=340.515 mJy/beam. ! Added new window around map position (3.2, -9.9). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00593118 Jy ! Component: 100 - total flux cleaned = 0.00967638 Jy ! Total flux subtracted in 100 components = 0.00967638 Jy ! Clean residual min=-0.003308 max=0.004039 Jy/beam ! Clean residual mean=0.000329 rms=0.000919 Jy/beam ! Combined flux in latest and established models = 1.76275 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 38 model components to the UV plane model. ! The established model now contains 397 components and 1.76275 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109856Jy sigma=0.001218 ! Fit after self-cal, rms=0.109840Jy sigma=0.001217 ! Inverting map and beam ! Estimated beam: bmin=0.6005 mas, bmaj=0.9964 mas, bpa=-14.15 degrees ! Estimated noise=340.657 mJy/beam. ! Added new window around map position (3.5, -6.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00513809 Jy ! Component: 100 - total flux cleaned = 0.00834211 Jy ! Total flux subtracted in 100 components = 0.00834211 Jy ! Clean residual min=-0.003062 max=0.003881 Jy/beam ! Clean residual mean=0.000312 rms=0.000866 Jy/beam ! Combined flux in latest and established models = 1.77109 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 51 model components to the UV plane model. ! The established model now contains 423 components and 1.77109 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109575Jy sigma=0.001212 ! Fit after self-cal, rms=0.109563Jy sigma=0.001211 ! Inverting map and beam ! Estimated beam: bmin=0.6005 mas, bmaj=0.9964 mas, bpa=-14.14 degrees ! Estimated noise=340.76 mJy/beam. ! Added new window around map position (6.2, -13). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00465942 Jy ! Component: 100 - total flux cleaned = 0.00760455 Jy ! Total flux subtracted in 100 components = 0.00760455 Jy ! Clean residual min=-0.002902 max=0.003699 Jy/beam ! Clean residual mean=0.000293 rms=0.000823 Jy/beam ! Combined flux in latest and established models = 1.7787 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 50 model components to the UV plane model. ! The established model now contains 437 components and 1.7787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109339Jy sigma=0.001207 ! Fit after self-cal, rms=0.109342Jy sigma=0.001206 ! Inverting map and beam ! Estimated beam: bmin=0.6004 mas, bmaj=0.9963 mas, bpa=-14.14 degrees ! Estimated noise=340.859 mJy/beam. ! Added new window around map position (3.5, -19.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00413988 Jy ! Component: 100 - total flux cleaned = 0.00702688 Jy ! Total flux subtracted in 100 components = 0.00702688 Jy ! Clean residual min=-0.002899 max=0.003507 Jy/beam ! Clean residual mean=0.000275 rms=0.000788 Jy/beam ! Combined flux in latest and established models = 1.78572 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 59 model components to the UV plane model. ! The established model now contains 457 components and 1.78572 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109146Jy sigma=0.001202 ! Fit after self-cal, rms=0.109153Jy sigma=0.001202 ! Inverting map and beam ! Estimated beam: bmin=0.6003 mas, bmaj=0.9963 mas, bpa=-14.14 degrees ! Estimated noise=340.942 mJy/beam. ! Added new window around map position (8.6, -6.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0040077 Jy ! Component: 100 - total flux cleaned = 0.00643111 Jy ! Total flux subtracted in 100 components = 0.00643111 Jy ! Clean residual min=-0.002831 max=0.003159 Jy/beam ! Clean residual mean=0.000262 rms=0.000755 Jy/beam ! Combined flux in latest and established models = 1.79216 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 52 model components to the UV plane model. ! The established model now contains 467 components and 1.79216 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108978Jy sigma=0.001199 ! Fit after self-cal, rms=0.108995Jy sigma=0.001198 ! Inverting map and beam ! Estimated beam: bmin=0.6003 mas, bmaj=0.9962 mas, bpa=-14.14 degrees ! Estimated noise=341.027 mJy/beam. ! Added new window around map position (3.7, -2.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00358418 Jy ! Component: 100 - total flux cleaned = 0.00585195 Jy ! Total flux subtracted in 100 components = 0.00585195 Jy ! Clean residual min=-0.002746 max=0.002971 Jy/beam ! Clean residual mean=0.000250 rms=0.000723 Jy/beam ! Combined flux in latest and established models = 1.79801 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 481 components and 1.79801 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108840Jy sigma=0.001195 ! Fit after self-cal, rms=0.108863Jy sigma=0.001195 ! Inverting map and beam ! Estimated beam: bmin=0.6002 mas, bmaj=0.9962 mas, bpa=-14.14 degrees ! Estimated noise=341.111 mJy/beam. ! Added new window around map position (1.7, -13.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00340854 Jy ! Component: 100 - total flux cleaned = 0.0057973 Jy ! Total flux subtracted in 100 components = 0.0057973 Jy ! Clean residual min=-0.002662 max=0.002927 Jy/beam ! Clean residual mean=0.000235 rms=0.000692 Jy/beam ! Combined flux in latest and established models = 1.80381 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 501 components and 1.80381 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108721Jy sigma=0.001193 ! Fit after self-cal, rms=0.108741Jy sigma=0.001193 ! Inverting map and beam ! Estimated beam: bmin=0.6002 mas, bmaj=0.9962 mas, bpa=-14.14 degrees ! Estimated noise=341.191 mJy/beam. ! Added new window around map position (12.3, -17.1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00343606 Jy ! Component: 100 - total flux cleaned = 0.00560263 Jy ! Total flux subtracted in 100 components = 0.00560263 Jy ! Clean residual min=-0.002647 max=0.002813 Jy/beam ! Clean residual mean=0.000225 rms=0.000663 Jy/beam ! Combined flux in latest and established models = 1.80941 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 55 model components to the UV plane model. ! The established model now contains 515 components and 1.80941 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108610Jy sigma=0.001190 ! Fit after self-cal, rms=0.108630Jy sigma=0.001190 ! Inverting map and beam ! Estimated beam: bmin=0.6001 mas, bmaj=0.9962 mas, bpa=-14.15 degrees ! Estimated noise=341.255 mJy/beam. ! Added new window around map position (-6.7, 0.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.003053 Jy ! Component: 100 - total flux cleaned = 0.00486254 Jy ! Total flux subtracted in 100 components = 0.00486254 Jy ! Clean residual min=-0.002640 max=0.002575 Jy/beam ! Clean residual mean=0.000216 rms=0.000638 Jy/beam ! Combined flux in latest and established models = 1.81427 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 523 components and 1.81427 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108513Jy sigma=0.001188 ! Fit after self-cal, rms=0.108537Jy sigma=0.001188 ! Inverting map and beam ! Estimated beam: bmin=0.6001 mas, bmaj=0.9963 mas, bpa=-14.15 degrees ! Estimated noise=341.303 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.00227693 Jy ! Component: 100 - total flux cleaned = 0.00406976 Jy ! Total flux subtracted in 100 components = 0.00406976 Jy ! Clean residual min=-0.002610 max=0.002509 Jy/beam ! Clean residual mean=0.000205 rms=0.000618 Jy/beam ! Combined flux in latest and established models = 1.81834 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 537 components and 1.81834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108443Jy sigma=0.001186 ! Fit after self-cal, rms=0.108396Jy sigma=0.001186 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save VR422201_X ! Writing UV FITS file: VR422201_X.uvf ! Writing 537 model components to file: VR422201_X.mod ! wwins: Wrote 22 windows to VR422201_X.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6001 x 0.9963 at -14.15 degrees (North through East) ! Clean map min=-0.0025374 max=0.89475 Jy/beam ! Writing clean map to FITS file: VR422201_X.fits ! Writing difmap environment to: VR422201_X.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot VR422201_X] device VR422201_Xuv.ps/vps,1,2 ! Attempting to open device: 'VR422201_Xuv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device VR422201_Xim.ps/vps ! Attempting to open device: 'VR422201_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.6001 mas, bmaj=0.9963 mas, bpa=-14.15 degrees ! Estimated noise=341.303 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6001 x 0.9963 at -14.15 degrees (North through East) ! Clean map min=-0.0025374 max=0.89475 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo6-0.1] quit ! Quitting program ! Log file VR422201_X.log closed on Tue Dec 12 21:31:57 2006