! Started logfile: 0742+103_S.log on Tue Jun 17 12:22:57 2008 observe 0742+103_S.edt ! Reading UV FITS file: 0742+103_S.edt ! AN table 1: 11 integrations on 190 of 190 possible baselines. ! AN table 2: 83 integrations on 190 of 190 possible baselines. ! AN table 3: 94 integrations on 190 of 190 possible baselines. ! AN table 4: 12 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.146711 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: 0742+103 ! ! 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 1901 lines of history. ! ! Reading 22300 visibilities. ![@gonzo7-0.4 0742+103_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 = 7.0 float win_mult win_mult = 1.8 float solint solint = 3.6/6.0 logical doflt doflt = true logical doamp float new_peak float flux_cutoff mapsize field_size,field_cell ! Map grid = 1024x1024 pixels with 0.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 220 telescope corrections were flagged in sub-array 1. ! A total of 1211 telescope corrections were flagged in sub-array 2. ! A total of 1455 telescope corrections were flagged in sub-array 3. ! A total of 219 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 220 telescope corrections were flagged in sub-array 1. ! A total of 1211 telescope corrections were flagged in sub-array 2. ! A total of 1455 telescope corrections were flagged in sub-array 3. ! A total of 219 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 220 telescope corrections were flagged in sub-array 1. ! A total of 1211 telescope corrections were flagged in sub-array 2. ! A total of 1461 telescope corrections were flagged in sub-array 3. ! A total of 219 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 220 telescope corrections were flagged in sub-array 1. ! A total of 1211 telescope corrections were flagged in sub-array 2. ! A total of 1463 telescope corrections were flagged in sub-array 3. ! A total of 219 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.843654Jy sigma=28.314109 ! Fit after self-cal, rms=1.843675Jy sigma=28.389257 ! 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.945 mas, bmaj=5.198 mas, bpa=-5.369 degrees ! Estimated noise=0.740086 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 = 1.88926 Jy ! Total flux subtracted in 50 components = 1.88926 Jy ! Clean residual min=-0.080377 max=0.529960 Jy/beam ! Clean residual mean=-0.000766 rms=0.026202 Jy/beam ! Combined flux in latest and established models = 1.88926 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 1.88926 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.148616Jy sigma=16.587437 ! Fit after self-cal, rms=1.148657Jy sigma=16.587179 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.532899 Jy ! Total flux subtracted in 50 components = 0.532899 Jy ! Clean residual min=-0.068111 max=0.238340 Jy/beam ! Clean residual mean=-0.000478 rms=0.018509 Jy/beam ! Combined flux in latest and established models = 2.42216 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 15 components and 2.42216 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.887495Jy sigma=11.348472 ! Fit after self-cal, rms=0.887333Jy sigma=11.347239 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.284606 Jy ! Total flux subtracted in 50 components = 0.284606 Jy ! Clean residual min=-0.064249 max=0.149776 Jy/beam ! Clean residual mean=-0.000327 rms=0.015749 Jy/beam ! Combined flux in latest and established models = 2.70676 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 15 model components to the UV plane model. ! The established model now contains 27 components and 2.70676 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.795287Jy sigma=9.285598 ! Fit after self-cal, rms=0.795104Jy sigma=9.284684 ! Inverting map ! clean: niter=50 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.18525 Jy ! Total flux subtracted in 50 components = 0.18525 Jy ! Clean residual min=-0.064217 max=0.100232 Jy/beam ! Clean residual mean=-0.000228 rms=0.014468 Jy/beam ! Combined flux in latest and established models = 2.89201 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 9 model components to the UV plane model. ! The established model now contains 32 components and 2.89201 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.751734Jy sigma=8.295347 ! Fit after self-cal, rms=0.751589Jy sigma=8.294727 ! 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.423 mas, bmaj=5.907 mas, bpa=-4.367 degrees ! Estimated noise=0.456079 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.172313 Jy ! Component: 100 - total flux cleaned = 0.269027 Jy ! Total flux subtracted in 100 components = 0.269027 Jy ! Clean residual min=-0.047250 max=0.113806 Jy/beam ! Clean residual mean=0.005764 rms=0.013946 Jy/beam ! Combined flux in latest and established models = 3.16104 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 53 components and 3.16104 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.711529Jy sigma=7.522862 ! Fit after self-cal, rms=0.710919Jy sigma=7.517678 ! Inverting map ! Added new window around map position (5.6, -5.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.122541 Jy ! Component: 100 - total flux cleaned = 0.196133 Jy ! Total flux subtracted in 100 components = 0.196133 Jy ! Clean residual min=-0.049169 max=0.087466 Jy/beam ! Clean residual mean=0.005038 rms=0.012445 Jy/beam ! Combined flux in latest and established models = 3.35717 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 79 components and 3.35717 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.683523Jy sigma=7.187935 ! Fit after self-cal, rms=0.682087Jy sigma=7.132140 ! Inverting map print "****************************************" ! **************************************** print "*** Finished Natural Weighting CLEAN ***" ! *** Finished Natural Weighting CLEAN *** print "****************************************" ! **************************************** gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 27 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.17 FD 1.04 GC 1.10 HH 1.00* ! HN 0.63 KK 0.84 KP 1.01 LA 0.97 ! MA 1.00* MC 1.00* MK 0.83 NL 0.94 ! NY 1.00* ON 1.00* OV 1.04 PT 0.96 ! SC 0.84 TS 1.01 WF 0.77 WZ 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.96 FD 1.01 GC 0.94 HH 1.00* ! HN 0.97 KK 1.23 KP 0.93 LA 1.02 ! MA 6.59 MC 1.26 MK 1.13 NL 1.01 ! NY 1.16 ON 1.08 OV 1.00 PT 1.02 ! SC 0.80 TS 1.00* WF 0.86 WZ 1.19 ! ! A total of 21 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 27 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 1.02 GC 0.99 HH 1.00* ! HN 0.79 KK 0.93 KP 0.97 LA 0.92 ! MA 1.00* MC 1.00* MK 0.98 NL 0.93 ! NY 1.00* ON 1.00* OV 1.05 PT 0.99 ! SC 0.74 TS 0.95 WF 0.87 WZ 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.98 FD 1.00 GC 1.14 HH 1.00* ! HN 1.07 KK 1.17 KP 1.02 LA 0.98 ! MA 2.91 MC 1.18 MK 1.06 NL 1.06 ! NY 1.16 ON 1.12 OV 1.01 PT 0.92 ! SC 1.09 TS 1.00* WF 0.99 WZ 1.09 ! ! A total of 21 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 27 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.14 FD 1.04 GC 1.11 HH 1.00* ! HN 0.77 KK 0.82 KP 1.04 LA 0.83 ! MA 1.00* MC 1.00* MK 0.65 NL 0.92 ! NY 1.00* ON 1.00* OV 1.09 PT 0.99 ! SC 0.99 TS 0.97 WF 0.82 WZ 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 1.00 GC 1.14 HH 1.00* ! HN 1.07 KK 1.20 KP 1.02 LA 1.03 ! MA 2.29 MC 2.45 MK 1.15 NL 1.03 ! NY 1.20 ON 1.21 OV 1.02 PT 1.02 ! SC 0.98 TS 1.00* WF 1.00 WZ 1.13 ! ! A total of 21 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 27 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.10 FD 1.03 GC 1.09 HH 1.00* ! HN 0.88 KK 0.91 KP 1.05 LA 0.80 ! MA 1.00* MC 1.00* MK 0.77 NL 0.94 ! NY 1.00* ON 1.00* OV 1.03 PT 1.04 ! SC 0.97 TS 0.96 WF 0.85 WZ 1.00* ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.97 FD 0.98 GC 1.14 HH 1.00* ! HN 1.09 KK 1.21 KP 1.02 LA 1.02 ! MA 2.75 MC 3.66 MK 1.01 NL 1.05 ! NY 1.17 ON 1.23 OV 1.04 PT 0.99 ! SC 1.06 TS 1.00* WF 1.04 WZ 1.20 ! ! A total of 21 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* HH 0.00* ! HN 0.00* KK 0.00* KP 0.00* LA 0.00* ! MA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* ON 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.682087Jy sigma=7.132140 ! Fit after self-cal, rms=0.309866Jy sigma=2.684272 dynam = dynam - 0.5 doamp = false flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.465 mas, bmaj=5.944 mas, bpa=-4.401 degrees ! Estimated noise=0.456711 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.0221408 Jy ! Component: 100 - total flux cleaned = 0.0219644 Jy ! Total flux subtracted in 100 components = 0.0219644 Jy ! Clean residual min=-0.022323 max=0.024964 Jy/beam ! Clean residual mean=0.001064 rms=0.003750 Jy/beam ! Combined flux in latest and established models = 3.37913 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 90 components and 3.37913 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.303053Jy sigma=2.549530 ! Fit after self-cal, rms=0.302035Jy sigma=2.529991 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0084104 Jy ! Component: 100 - total flux cleaned = 0.00841402 Jy ! Total flux subtracted in 100 components = 0.00841402 Jy ! Clean residual min=-0.021111 max=0.018858 Jy/beam ! Clean residual mean=0.001044 rms=0.003558 Jy/beam ! Combined flux in latest and established models = 3.38755 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 93 components and 3.38755 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.299135Jy sigma=2.478459 ! Fit after self-cal, rms=0.298879Jy sigma=2.474689 ! 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.00360123 Jy ! Component: 100 - total flux cleaned = 0.00362428 Jy ! Total flux subtracted in 100 components = 0.00362428 Jy ! Clean residual min=-0.019675 max=0.018705 Jy/beam ! Clean residual mean=0.001035 rms=0.003455 Jy/beam ! Combined flux in latest and established models = 3.39117 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 15 model components to the UV plane model. ! The established model now contains 95 components and 3.39117 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.296948Jy sigma=2.443516 ! Fit after self-cal, rms=0.332553Jy sigma=1.785750 ! Inverting map and beam ! Estimated beam: bmin=2.462 mas, bmaj=5.948 mas, bpa=-4.621 degrees ! Estimated noise=0.457527 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0113061 Jy ! Component: 100 - total flux cleaned = 0.0126471 Jy ! Total flux subtracted in 100 components = 0.0126471 Jy ! Clean residual min=-0.014811 max=0.014489 Jy/beam ! Clean residual mean=0.000690 rms=0.002302 Jy/beam ! Combined flux in latest and established models = 3.40382 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 104 components and 3.40382 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.330243Jy sigma=1.737390 ! Fit after self-cal, rms=0.317758Jy sigma=1.693516 ! Inverting map and beam ! Estimated beam: bmin=2.462 mas, bmaj=5.948 mas, bpa=-4.701 degrees ! Estimated noise=0.458072 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00861896 Jy ! Component: 100 - total flux cleaned = 0.0108755 Jy ! Total flux subtracted in 100 components = 0.0108755 Jy ! Clean residual min=-0.013443 max=0.013900 Jy/beam ! Clean residual mean=0.000568 rms=0.002109 Jy/beam ! Combined flux in latest and established models = 3.4147 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 17 model components to the UV plane model. ! The established model now contains 107 components and 3.4147 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316235Jy sigma=1.661603 ! Fit after self-cal, rms=0.315661Jy sigma=1.642103 ! Inverting map and beam ! Estimated beam: bmin=2.462 mas, bmaj=5.948 mas, bpa=-4.751 degrees ! Estimated noise=0.458589 mJy/beam. ! Added new window around map position (2, 2.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00713042 Jy ! Component: 100 - total flux cleaned = 0.00861795 Jy ! Total flux subtracted in 100 components = 0.00861795 Jy ! Clean residual min=-0.012553 max=0.013517 Jy/beam ! Clean residual mean=0.000494 rms=0.001979 Jy/beam ! Combined flux in latest and established models = 3.42331 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 22 model components to the UV plane model. ! The established model now contains 119 components and 3.42331 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.314516Jy sigma=1.617102 ! Fit after self-cal, rms=0.317045Jy sigma=1.606383 ! Inverting map and beam ! Estimated beam: bmin=2.462 mas, bmaj=5.947 mas, bpa=-4.785 degrees ! Estimated noise=0.459079 mJy/beam. ! Added new window around map position (11.2, -5.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0128208 Jy ! Component: 100 - total flux cleaned = 0.0153429 Jy ! Total flux subtracted in 100 components = 0.0153429 Jy ! Clean residual min=-0.011940 max=0.010028 Jy/beam ! Clean residual mean=0.000417 rms=0.001818 Jy/beam ! Combined flux in latest and established models = 3.43866 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 24 model components to the UV plane model. ! The established model now contains 127 components and 3.43866 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316048Jy sigma=1.578889 ! Fit after self-cal, rms=0.318211Jy sigma=1.572007 ! Inverting map and beam ! Estimated beam: bmin=2.461 mas, bmaj=5.946 mas, bpa=-4.8 degrees ! Estimated noise=0.459446 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00781539 Jy ! Component: 100 - total flux cleaned = 0.00871712 Jy ! Total flux subtracted in 100 components = 0.00871712 Jy ! Clean residual min=-0.011344 max=0.009575 Jy/beam ! Clean residual mean=0.000378 rms=0.001724 Jy/beam ! Combined flux in latest and established models = 3.44737 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 28 model components to the UV plane model. ! The established model now contains 132 components and 3.44737 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317414Jy sigma=1.551810 ! Fit after self-cal, rms=0.319313Jy sigma=1.547668 ! Inverting map and beam ! Estimated beam: bmin=2.461 mas, bmaj=5.944 mas, bpa=-4.813 degrees ! Estimated noise=0.459813 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00535814 Jy ! Component: 100 - total flux cleaned = 0.00621819 Jy ! Total flux subtracted in 100 components = 0.00621819 Jy ! Clean residual min=-0.010529 max=0.008972 Jy/beam ! Clean residual mean=0.000348 rms=0.001643 Jy/beam ! Combined flux in latest and established models = 3.45359 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 24 model components to the UV plane model. ! The established model now contains 134 components and 3.45359 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.318617Jy sigma=1.530430 ! Fit after self-cal, rms=0.320284Jy sigma=1.527490 ! Inverting map and beam ! Estimated beam: bmin=2.461 mas, bmaj=5.943 mas, bpa=-4.824 degrees ! Estimated noise=0.460159 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00456166 Jy ! Component: 100 - total flux cleaned = 0.00575413 Jy ! Total flux subtracted in 100 components = 0.00575413 Jy ! Clean residual min=-0.010011 max=0.008396 Jy/beam ! Clean residual mean=0.000324 rms=0.001568 Jy/beam ! Combined flux in latest and established models = 3.45934 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 29 model components to the UV plane model. ! The established model now contains 138 components and 3.45934 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.319663Jy sigma=1.512537 ! Fit after self-cal, rms=0.321168Jy sigma=1.510308 ! Inverting map and beam ! Estimated beam: bmin=2.461 mas, bmaj=5.942 mas, bpa=-4.832 degrees ! Estimated noise=0.460489 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00425924 Jy ! Component: 100 - total flux cleaned = 0.00499413 Jy ! Total flux subtracted in 100 components = 0.00499413 Jy ! Clean residual min=-0.009472 max=0.008035 Jy/beam ! Clean residual mean=0.000305 rms=0.001503 Jy/beam ! Combined flux in latest and established models = 3.46434 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 28 model components to the UV plane model. ! The established model now contains 142 components and 3.46434 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.320617Jy sigma=1.497192 ! Fit after self-cal, rms=0.321954Jy sigma=1.495456 ! Inverting map and beam ! Estimated beam: bmin=2.46 mas, bmaj=5.94 mas, bpa=-4.839 degrees ! Estimated noise=0.460789 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00399307 Jy ! Component: 100 - total flux cleaned = 0.00468983 Jy ! Total flux subtracted in 100 components = 0.00468983 Jy ! Clean residual min=-0.009044 max=0.007641 Jy/beam ! Clean residual mean=0.000286 rms=0.001444 Jy/beam ! Combined flux in latest and established models = 3.46903 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 29 model components to the UV plane model. ! The established model now contains 145 components and 3.46903 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.321467Jy sigma=1.483880 ! Fit after self-cal, rms=0.322680Jy sigma=1.482483 ! Inverting map and beam ! Estimated beam: bmin=2.46 mas, bmaj=5.939 mas, bpa=-4.846 degrees ! Estimated noise=0.461094 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00342973 Jy ! Component: 100 - total flux cleaned = 0.00409622 Jy ! Total flux subtracted in 100 components = 0.00409622 Jy ! Clean residual min=-0.008608 max=0.007162 Jy/beam ! Clean residual mean=0.000271 rms=0.001392 Jy/beam ! Combined flux in latest and established models = 3.47312 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 31 model components to the UV plane model. ! The established model now contains 151 components and 3.47313 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.322247Jy sigma=1.472206 ! Fit after self-cal, rms=0.323325Jy sigma=1.471123 ! Inverting map and beam ! Estimated beam: bmin=2.46 mas, bmaj=5.937 mas, bpa=-4.851 degrees ! Estimated noise=0.461384 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00357842 Jy ! Component: 100 - total flux cleaned = 0.00420126 Jy ! Total flux subtracted in 100 components = 0.00420126 Jy ! Clean residual min=-0.008255 max=0.006783 Jy/beam ! Clean residual mean=0.000255 rms=0.001344 Jy/beam ! Combined flux in latest and established models = 3.47733 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 30 model components to the UV plane model. ! The established model now contains 155 components and 3.47733 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.322937Jy sigma=1.461928 ! Fit after self-cal, rms=0.323904Jy sigma=1.460971 ! Inverting map and beam ! Estimated beam: bmin=2.46 mas, bmaj=5.936 mas, bpa=-4.855 degrees ! Estimated noise=0.46168 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00339717 Jy ! Component: 100 - total flux cleaned = 0.00399081 Jy ! Total flux subtracted in 100 components = 0.00399081 Jy ! Clean residual min=-0.007928 max=0.006525 Jy/beam ! Clean residual mean=0.000240 rms=0.001300 Jy/beam ! Combined flux in latest and established models = 3.48132 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 27 model components to the UV plane model. ! The established model now contains 157 components and 3.48132 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.323551Jy sigma=1.452768 ! Fit after self-cal, rms=0.324447Jy sigma=1.451942 ! Inverting map and beam ! Estimated beam: bmin=2.459 mas, bmaj=5.934 mas, bpa=-4.86 degrees ! Estimated noise=0.461977 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00323527 Jy ! Component: 100 - total flux cleaned = 0.00380167 Jy ! Total flux subtracted in 100 components = 0.00380167 Jy ! Clean residual min=-0.007578 max=0.006416 Jy/beam ! Clean residual mean=0.000225 rms=0.001260 Jy/beam ! Combined flux in latest and established models = 3.48512 Jy ! Performing amp+phase self-cal over 1 minute time intervals ! Adding 30 model components to the UV plane model. ! The established model now contains 159 components and 3.48512 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.324127Jy sigma=1.444599 ! Fit after self-cal, rms=0.324921Jy sigma=1.443977 ! Inverting map and beam ! Estimated beam: bmin=2.459 mas, bmaj=5.933 mas, bpa=-4.863 degrees ! Estimated noise=0.462236 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.00335278 Jy ! Component: 100 - total flux cleaned = 0.00361839 Jy ! Total flux subtracted in 100 components = 0.00361839 Jy ! Clean residual min=-0.007304 max=0.006186 Jy/beam ! Clean residual mean=0.000211 rms=0.001223 Jy/beam ! Combined flux in latest and established models = 3.48874 Jy ! Performing phase self-cal over 0.6 minute time intervals ! Adding 28 model components to the UV plane model. ! The established model now contains 162 components and 3.48874 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.324632Jy sigma=1.437402 ! Fit after self-cal, rms=0.324510Jy sigma=1.425006 ! Inverting map print "***********************************" ! *********************************** print "*** Finished Amplitude Self-Cal ***" ! *** Finished Amplitude Self-Cal *** print "***********************************" ! *********************************** save 0742+103_S ! Writing UV FITS file: 0742+103_S.uvf ! Writing 162 model components to file: 0742+103_S.mod ! wwins: Wrote 4 windows to 0742+103_S.win ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.459 x 5.933 at -4.863 degrees (North through East) ! Clean map min=-0.0047512 max=2.6009 Jy/beam ! Writing clean map to FITS file: 0742+103_S.fits ! Writing difmap environment to: 0742+103_S.par print "****************************" ! **************************** print "*** Now, make some plots ***" ! *** Now, make some plots *** print "****************************" ! **************************** ![@gonzo_plot 0742+103_S] device 0742+103_Suv.ps/vps,1,2 ! Attempting to open device: '0742+103_Suv.ps/vps' radplot ! Using default options string "m1" uvplot ! Using default options string "" device 0742+103_Sim.ps/vps ! Attempting to open device: '0742+103_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.459 mas, bmaj=5.933 mas, bpa=-4.863 degrees ! Estimated noise=0.462236 mJy/beam. mappl cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.459 x 5.933 at -4.863 degrees (North through East) ! Clean map min=-0.0047512 max=2.6009 Jy/beam ![Exited script file: gonzo_plot] ![Exited script file: gonzo7-0.4] quit ! Quitting program ! Log file 0742+103_S.log closed on Tue Jun 17 12:23:18 2008