! Started logfile: J1924-2914_X_2003_09_17_pus_uvs.log on Mon May 14 14:43:19 2007 obs J1924-2914_X_2003_09_17_pus_uva.fits ! Reading UV FITS file: J1924-2914_X_2003_09_17_pus_uva.fits ! AN table 1: 31 integrations on 120 of 120 possible baselines. ! AN table 2: 89 integrations on 120 of 120 possible baselines. ! AN table 3: 97 integrations on 120 of 120 possible baselines. ! AN table 4: 27 integrations on 120 of 120 possible baselines. ! AN table 5: 82 integrations on 120 of 120 possible baselines. ! Apparent sampling: 0.0669223 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1924-2914 ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! ! 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 1489 lines of history. ! ! Reading 10472 visibilities. select I ! Polarization I is unavailable. ! 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 float field_size field_size = 512 float field_cell field_cell = 0.2 float taper_size taper_size = 70. integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 512x512 pixels with 0.200x0.200 milli-arcsec cellsize. 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 496 telescope corrections were flagged in sub-array 1. ! A total of 1327 telescope corrections were flagged in sub-array 2. ! A total of 1133 telescope corrections were flagged in sub-array 3. ! A total of 302 telescope corrections were flagged in sub-array 4. ! A total of 1184 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 496 telescope corrections were flagged in sub-array 1. ! A total of 1327 telescope corrections were flagged in sub-array 2. ! A total of 1133 telescope corrections were flagged in sub-array 3. ! A total of 302 telescope corrections were flagged in sub-array 4. ! A total of 1184 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 496 telescope corrections were flagged in sub-array 1. ! A total of 1327 telescope corrections were flagged in sub-array 2. ! A total of 1133 telescope corrections were flagged in sub-array 3. ! A total of 302 telescope corrections were flagged in sub-array 4. ! A total of 1184 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 496 telescope corrections were flagged in sub-array 1. ! A total of 1327 telescope corrections were flagged in sub-array 2. ! A total of 1133 telescope corrections were flagged in sub-array 3. ! A total of 302 telescope corrections were flagged in sub-array 4. ! A total of 1184 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=4.446073Jy sigma=0.039329 ! Fit after self-cal, rms=4.444833Jy sigma=0.039310 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvw 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. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8668 mas, bmaj=2.817 mas, bpa=12.91 degrees ! Estimated noise=1479.21 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.16364 Jy ! Component: 100 - total flux cleaned = 5.07439 Jy ! Total flux subtracted in 100 components = 5.07439 Jy ! Clean residual min=-0.122568 max=0.293181 Jy/beam ! Clean residual mean=-0.000963 rms=0.047480 Jy/beam ! Combined flux in latest and established models = 5.07439 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 5.07439 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.816463Jy sigma=0.007772 ! Fit after self-cal, rms=0.816422Jy sigma=0.007770 ! Inverting map ! Added new window around map position (3.2, 7). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.338592 Jy ! Component: 100 - total flux cleaned = 0.561947 Jy ! Total flux subtracted in 100 components = 0.561947 Jy ! Clean residual min=-0.091874 max=0.187802 Jy/beam ! Clean residual mean=-0.000522 rms=0.031505 Jy/beam ! Combined flux in latest and established models = 5.63633 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 38 components and 5.63633 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.610262Jy sigma=0.005688 ! Fit after self-cal, rms=0.604122Jy sigma=0.005653 ! Inverting map ! Added new window around map position (1.6, 7.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.200161 Jy ! Component: 100 - total flux cleaned = 0.342753 Jy ! Total flux subtracted in 100 components = 0.342753 Jy ! Clean residual min=-0.071308 max=0.121031 Jy/beam ! Clean residual mean=-0.000374 rms=0.023442 Jy/beam ! Combined flux in latest and established models = 5.97909 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 62 components and 5.97909 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.478000Jy sigma=0.004373 ! Fit after self-cal, rms=0.468726Jy sigma=0.004319 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 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 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.168 mas, bmaj=3.139 mas, bpa=2.485 degrees ! Estimated noise=898.729 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.254442 Jy ! Component: 100 - total flux cleaned = 0.405171 Jy ! Component: 150 - total flux cleaned = 0.512748 Jy ! Component: 200 - total flux cleaned = 0.592931 Jy ! Total flux subtracted in 200 components = 0.592931 Jy ! Clean residual min=-0.060633 max=0.137713 Jy/beam ! Clean residual mean=0.000829 rms=0.024926 Jy/beam ! Combined flux in latest and established models = 6.57202 Jy ! Performing phase self-cal ! Adding 71 model components to the UV plane model. ! The established model now contains 126 components and 6.57202 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.305493Jy sigma=0.002596 ! Fit after self-cal, rms=0.278418Jy sigma=0.002326 ! Inverting map ! Added new window around map position (5, 5.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.130384 Jy ! Component: 100 - total flux cleaned = 0.20784 Jy ! Component: 150 - total flux cleaned = 0.264325 Jy ! Component: 200 - total flux cleaned = 0.307911 Jy ! Total flux subtracted in 200 components = 0.307911 Jy ! Clean residual min=-0.030772 max=0.043910 Jy/beam ! Clean residual mean=-0.000528 rms=0.011026 Jy/beam ! Combined flux in latest and established models = 6.87993 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 175 components and 6.87993 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.202603Jy sigma=0.001527 ! Fit after self-cal, rms=0.197146Jy sigma=0.001470 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! FD 0.00* GC 0.00* HN 0.00* KP 0.00* ! LA 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* KK 0.00* WZ 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00* GC 1.00* HN 1.00* KP 1.00 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.99 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.01 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 0.98 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.99 GC 1.00* HN 0.99 KP 1.00 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 0.99 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 0.00* GC 0.00* HN 0.00* KP 0.00* ! LA 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* KK 0.00* WZ 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00* GC 1.00* HN 1.00* KP 1.00 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.99 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.01 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 0.99 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.99 GC 1.00* HN 0.99 KP 1.00 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 0.99 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 0.00* GC 0.00* HN 0.00* KP 0.00* ! LA 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* KK 0.00* WZ 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00* GC 1.00* HN 1.00* KP 1.00 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.99 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 GC 1.00* HN 1.00 KP 0.99 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 0.99 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.00 GC 1.00* HN 0.99 KP 1.00 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 0.99 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 0.00* GC 0.00* HN 0.00* KP 0.00* ! LA 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* KK 0.00* WZ 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00* GC 1.00* HN 1.00* KP 1.00 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.99 GC 1.00* HN 1.01 KP 0.99 ! LA 0.99 MK 1.00* NL 1.01 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.01 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.00 GC 1.00* HN 1.00 KP 0.99 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 1.00 PT 0.99 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! FD 1.00 GC 1.00* HN 0.99 KP 1.00 ! LA 0.99 MK 1.00* NL 1.00 NY 1.00* ! ON 1.00* OV 0.99 PT 1.00 SC 1.00 ! TS 1.00* KK 1.00* WZ 1.00* WF 0.99 ! ! ! Fit before self-cal, rms=0.197146Jy sigma=0.001470 ! Fit after self-cal, rms=0.181014Jy sigma=0.001359 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.171 mas, bmaj=3.149 mas, bpa=2.523 degrees ! Estimated noise=894.814 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0329372 Jy ! Component: 100 - total flux cleaned = 0.0469403 Jy ! Component: 150 - total flux cleaned = 0.0568578 Jy ! Component: 200 - total flux cleaned = 0.0608231 Jy ! Total flux subtracted in 200 components = 0.0608231 Jy ! Clean residual min=-0.022737 max=0.024606 Jy/beam ! Clean residual mean=-0.000435 rms=0.005972 Jy/beam ! Combined flux in latest and established models = 6.94075 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 192 components and 6.94075 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.172498Jy sigma=0.001274 ! Fit after self-cal, rms=0.171200Jy sigma=0.001264 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171200Jy sigma=0.001264 ! Fit after self-cal, rms=0.170753Jy sigma=0.001260 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.171 mas, bmaj=3.15 mas, bpa=2.517 degrees ! Estimated noise=894.992 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.008406 Jy ! Component: 100 - total flux cleaned = 0.0132463 Jy ! Component: 150 - total flux cleaned = 0.0154679 Jy ! Component: 200 - total flux cleaned = 0.0158557 Jy ! Total flux subtracted in 200 components = 0.0158557 Jy ! Clean residual min=-0.021245 max=0.021097 Jy/beam ! Clean residual mean=-0.000354 rms=0.005257 Jy/beam ! Combined flux in latest and established models = 6.95661 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 199 components and 6.95661 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167580Jy sigma=0.001232 ! Fit after self-cal, rms=0.166921Jy sigma=0.001228 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00323285 Jy ! Component: 100 - total flux cleaned = 0.00359849 Jy ! Component: 150 - total flux cleaned = 0.00360803 Jy ! Component: 200 - total flux cleaned = 0.00392363 Jy ! Total flux subtracted in 200 components = 0.00392363 Jy ! Clean residual min=-0.020029 max=0.017981 Jy/beam ! Clean residual mean=-0.000372 rms=0.004912 Jy/beam ! Combined flux in latest and established models = 6.96053 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 208 components and 6.96053 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.165318Jy sigma=0.001215 ! Fit after self-cal, rms=0.165006Jy sigma=0.001214 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.165006Jy sigma=0.001214 ! Fit after self-cal, rms=0.164598Jy sigma=0.001211 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.171 mas, bmaj=3.152 mas, bpa=2.512 degrees ! Estimated noise=895.227 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00194092 Jy ! Component: 100 - total flux cleaned = 0.00286666 Jy ! Component: 150 - total flux cleaned = 0.00313591 Jy ! Component: 200 - total flux cleaned = 0.00337525 Jy ! Total flux subtracted in 200 components = 0.00337525 Jy ! Clean residual min=-0.019049 max=0.018007 Jy/beam ! Clean residual mean=-0.000287 rms=0.004699 Jy/beam ! Combined flux in latest and established models = 6.96391 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 214 components and 6.96391 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163398Jy sigma=0.001200 ! Fit after self-cal, rms=0.163090Jy sigma=0.001199 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00182822 Jy ! Component: 100 - total flux cleaned = 0.00158215 Jy ! Component: 150 - total flux cleaned = 0.00202559 Jy ! Component: 200 - total flux cleaned = 0.00243343 Jy ! Total flux subtracted in 200 components = 0.00243343 Jy ! Clean residual min=-0.017934 max=0.016391 Jy/beam ! Clean residual mean=-0.000297 rms=0.004549 Jy/beam ! Combined flux in latest and established models = 6.96634 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 220 components and 6.96634 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162389Jy sigma=0.001192 ! Fit after self-cal, rms=0.162222Jy sigma=0.001192 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162222Jy sigma=0.001192 ! Fit after self-cal, rms=0.162067Jy sigma=0.001190 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.349 mas, bmaj=6.074 mas, bpa=-8.312 degrees ! Estimated noise=1006.75 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00459123 Jy ! Component: 100 - total flux cleaned = 0.00720683 Jy ! Component: 150 - total flux cleaned = 0.00791697 Jy ! Component: 200 - total flux cleaned = 0.00813908 Jy ! Total flux subtracted in 200 components = 0.00813908 Jy ! Clean residual min=-0.018473 max=0.016129 Jy/beam ! Clean residual mean=-0.000438 rms=0.005541 Jy/beam ! Combined flux in latest and established models = 6.97448 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 226 components and 6.97448 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161591Jy sigma=0.001186 ! Fit after self-cal, rms=0.161363Jy sigma=0.001184 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00126651 Jy ! Component: 100 - total flux cleaned = 0.00150692 Jy ! Component: 150 - total flux cleaned = 0.00171293 Jy ! Component: 200 - total flux cleaned = 0.00189516 Jy ! Total flux subtracted in 200 components = 0.00189516 Jy ! Clean residual min=-0.018066 max=0.016529 Jy/beam ! Clean residual mean=-0.000451 rms=0.005409 Jy/beam ! Combined flux in latest and established models = 6.97638 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 230 components and 6.97638 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161150Jy sigma=0.001182 ! Fit after self-cal, rms=0.160997Jy sigma=0.001181 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=1.172 mas, bmaj=3.154 mas, bpa=2.51 degrees ! Estimated noise=895.379 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.172 x 3.154 at 2.51 degrees (North through East) ! Clean map min=-0.015275 max=5.3035 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal 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 new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160997Jy sigma=0.001181 ! Fit after self-cal, rms=0.160802Jy sigma=0.001179 ! Inverting map and beam ! Estimated beam: bmin=1.172 mas, bmaj=3.155 mas, bpa=2.519 degrees ! Estimated noise=895.675 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00165816 Jy ! Component: 100 - total flux cleaned = -0.00306602 Jy ! Component: 150 - total flux cleaned = -0.00409675 Jy ! Component: 200 - total flux cleaned = -0.00466723 Jy ! Total flux subtracted in 200 components = -0.00466723 Jy ! Clean residual min=-0.015744 max=0.014941 Jy/beam ! Clean residual mean=-0.000233 rms=0.004224 Jy/beam ! Combined flux in latest and established models = 6.97171 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 241 components and 6.97171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160189Jy sigma=0.001173 ! Fit after self-cal, rms=0.160118Jy sigma=0.001173 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000361083 Jy ! Component: 100 - total flux cleaned = 0.000182223 Jy ! Component: 150 - total flux cleaned = 1.12928e-05 Jy ! Component: 200 - total flux cleaned = -0.000148377 Jy ! Total flux subtracted in 200 components = -0.000148377 Jy ! Clean residual min=-0.015132 max=0.014324 Jy/beam ! Clean residual mean=-0.000231 rms=0.004175 Jy/beam ! Combined flux in latest and established models = 6.97156 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 245 components and 6.97156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.159770Jy sigma=0.001170 ! Fit after self-cal, rms=0.159716Jy sigma=0.001170 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.159716Jy sigma=0.001170 ! Fit after self-cal, rms=0.159741Jy sigma=0.001169 print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. uvw 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. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8685 mas, bmaj=2.829 mas, bpa=12.87 degrees ! Estimated noise=1475.27 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.1498 Jy ! Component: 100 - total flux cleaned = 5.06397 Jy ! Component: 150 - total flux cleaned = 5.42583 Jy ! Component: 200 - total flux cleaned = 5.68309 Jy ! Total flux subtracted in 200 components = 5.68308 Jy ! Clean residual min=-0.071183 max=0.148292 Jy/beam ! Clean residual mean=-0.000326 rms=0.030252 Jy/beam ! Combined flux in latest and established models = 5.68308 Jy uvw 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. flux_cutoff = imstat(rms) * dynam ! Adding 45 model components to the UV plane model. ! The established model now contains 45 components and 5.68309 Jy ! Inverting map and beam ! Estimated beam: bmin=1.172 mas, bmaj=3.156 mas, bpa=2.519 degrees ! Estimated noise=895.939 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.33689 Jy ! Component: 100 - total flux cleaned = 0.560179 Jy ! Component: 150 - total flux cleaned = 0.721253 Jy ! Component: 200 - total flux cleaned = 0.841264 Jy ! Total flux subtracted in 200 components = 0.841264 Jy ! Clean residual min=-0.042422 max=0.069459 Jy/beam ! Clean residual mean=0.001526 rms=0.020860 Jy/beam ! Combined flux in latest and established models = 6.52435 Jy ! Adding 100 model components to the UV plane model. ! The established model now contains 135 components and 6.52435 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.35 mas, bmaj=6.075 mas, bpa=-8.297 degrees ! Estimated noise=1006.95 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.138379 Jy ! Component: 100 - total flux cleaned = 0.229583 Jy ! Component: 150 - total flux cleaned = 0.292008 Jy ! Component: 200 - total flux cleaned = 0.335471 Jy ! Total flux subtracted in 200 components = 0.335471 Jy ! Clean residual min=-0.023751 max=0.030256 Jy/beam ! Clean residual mean=0.000387 rms=0.009707 Jy/beam ! Combined flux in latest and established models = 6.85982 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0306346 Jy ! Component: 100 - total flux cleaned = 0.0525595 Jy ! Component: 150 - total flux cleaned = 0.0683351 Jy ! Component: 200 - total flux cleaned = 0.0797029 Jy ! Total flux subtracted in 200 components = 0.0797029 Jy ! Clean residual min=-0.015971 max=0.017965 Jy/beam ! Clean residual mean=-0.000176 rms=0.005669 Jy/beam ! Combined flux in latest and established models = 6.93952 Jy ! Adding 124 model components to the UV plane model. ! The established model now contains 245 components and 6.93952 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.172 mas, bmaj=3.156 mas, bpa=2.519 degrees ! Estimated noise=895.939 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00256107 Jy ! Component: 100 - total flux cleaned = 0.00425044 Jy ! Component: 150 - total flux cleaned = 0.00706298 Jy ! Component: 200 - total flux cleaned = 0.00872102 Jy ! Total flux subtracted in 200 components = 0.00872102 Jy ! Clean residual min=-0.015703 max=0.016419 Jy/beam ! Clean residual mean=-0.000126 rms=0.004236 Jy/beam ! Combined flux in latest and established models = 6.94825 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 256 components and 6.94825 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162770Jy sigma=0.001191 ! Fit after self-cal, rms=0.161796Jy sigma=0.001183 wmodel J1924-2914_X_2003_09_17_pus_map.mod ! Writing 256 model components to file: J1924-2914_X_2003_09_17_pus_map.mod wobs J1924-2914_X_2003_09_17_pus_uvs.fits ! Writing UV FITS file: J1924-2914_X_2003_09_17_pus_uvs.fits wwins J1924-2914_X_2003_09_17_pus_map.win ! wwins: Wrote 4 windows to J1924-2914_X_2003_09_17_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00523307 Jy ! Component: 100 - total flux cleaned = 0.00891036 Jy ! Component: 150 - total flux cleaned = 0.0125849 Jy ! Component: 200 - total flux cleaned = 0.0162303 Jy ! Component: 250 - total flux cleaned = 0.0185787 Jy ! Component: 300 - total flux cleaned = 0.0204038 Jy ! Component: 350 - total flux cleaned = 0.0225418 Jy ! Component: 400 - total flux cleaned = 0.0238448 Jy ! Component: 450 - total flux cleaned = 0.0252777 Jy ! Component: 500 - total flux cleaned = 0.0261369 Jy ! Component: 550 - total flux cleaned = 0.0272982 Jy ! Component: 600 - total flux cleaned = 0.0280984 Jy ! Component: 650 - total flux cleaned = 0.02887 Jy ! Component: 700 - total flux cleaned = 0.0294681 Jy ! Component: 750 - total flux cleaned = 0.029615 Jy ! Component: 800 - total flux cleaned = 0.0300391 Jy ! Component: 850 - total flux cleaned = 0.0299007 Jy ! Component: 900 - total flux cleaned = 0.0303083 Jy ! Component: 950 - total flux cleaned = 0.0303103 Jy ! Component: 1000 - total flux cleaned = 0.0307005 Jy ! Component: 1050 - total flux cleaned = 0.0306996 Jy ! Component: 1100 - total flux cleaned = 0.0306989 Jy ! Component: 1150 - total flux cleaned = 0.0304562 Jy ! Component: 1200 - total flux cleaned = 0.0304581 Jy ! Component: 1250 - total flux cleaned = 0.0304564 Jy ! Component: 1300 - total flux cleaned = 0.0302259 Jy ! Component: 1350 - total flux cleaned = 0.0300011 Jy ! Component: 1400 - total flux cleaned = 0.0296688 Jy ! Component: 1450 - total flux cleaned = 0.0294512 Jy ! Component: 1500 - total flux cleaned = 0.0292347 Jy ! Component: 1550 - total flux cleaned = 0.0288106 Jy ! Component: 1600 - total flux cleaned = 0.0283948 Jy ! Component: 1650 - total flux cleaned = 0.0280894 Jy ! Component: 1700 - total flux cleaned = 0.0276885 Jy ! Component: 1750 - total flux cleaned = 0.0271937 Jy ! Component: 1800 - total flux cleaned = 0.0269986 Jy ! Component: 1850 - total flux cleaned = 0.0266118 Jy ! Component: 1900 - total flux cleaned = 0.0262316 Jy ! Component: 1950 - total flux cleaned = 0.0254843 Jy ! Component: 2000 - total flux cleaned = 0.0252095 Jy ! Total flux subtracted in 2048 components = 0.0250269 Jy ! Clean residual min=-0.005806 max=0.004756 Jy/beam ! Clean residual mean=-0.000180 rms=0.001882 Jy/beam ! Combined flux in latest and established models = 6.97327 Jy keep ! Adding 133 model components to the UV plane model. ! The established model now contains 389 components and 6.97327 Jy wmap "J1924-2914_X_2003_09_17_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.172 x 3.156 at 2.519 degrees (North through East) ! Clean map min=-0.011437 max=5.3038 Jy/beam ! Writing clean map to FITS file: J1924-2914_X_2003_09_17_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 5.30375 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00187097 device "J1924-2914_X_2003_09_17_pus\_map.eps/VPS" ! Attempting to open device: 'J1924-2914_X_2003_09_17_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.105829 0.105829 0.211658 0.423316 0.846632 1.69326 3.38653 6.77306 13.5461 27.0922 54.1845 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0056 print "peak=",int(1000*peak_flux)/1000 ! peak= 5.303 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.172 x 3.156 at 2.519 degrees (North through East) ! Clean map min=-0.011437 max=5.3038 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1924-2914_X_2003_09_17_pus\_map.gif/vgif" ! Attempting to open device: 'J1924-2914_X_2003_09_17_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J1924-2914_X_2003_09_17_pus_uvs.log closed on Mon May 14 14:44:54 2007