! Started logfile: J0241-0815_X_2001_07_05_pus_uvs.log on Thu Jul 5 16:00:18 2007 obs J0241-0815_X_2001_07_05_pus_uva.fits ! Reading UV FITS file: J0241-0815_X_2001_07_05_pus_uva.fits ! AN table 1: 11 integrations on 171 of 171 possible baselines. ! AN table 2: 177 integrations on 171 of 171 possible baselines. ! AN table 3: 220 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.125688 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0241-0815 ! 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 2374 lines of history. ! ! Reading 35076 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 print "field_size =",field_size ! field_size = 512 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",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 209 telescope corrections were flagged in sub-array 1. ! A total of 2361 telescope corrections were flagged in sub-array 2. ! A total of 2975 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 209 telescope corrections were flagged in sub-array 1. ! A total of 2363 telescope corrections were flagged in sub-array 2. ! A total of 2975 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 209 telescope corrections were flagged in sub-array 1. ! A total of 2363 telescope corrections were flagged in sub-array 2. ! A total of 2975 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 209 telescope corrections were flagged in sub-array 1. ! A total of 2363 telescope corrections were flagged in sub-array 2. ! A total of 2975 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.861145Jy sigma=0.006988 ! Fit after self-cal, rms=0.565121Jy sigma=0.004340 ! 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 ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.5733 mas, bmaj=1.694 mas, bpa=-2.059 degrees ! Estimated noise=878.119 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 = 0.379194 Jy ! Component: 100 - total flux cleaned = 0.531266 Jy ! Total flux subtracted in 100 components = 0.531266 Jy ! Clean residual min=-0.028409 max=0.070873 Jy/beam ! Clean residual mean=0.000460 rms=0.009671 Jy/beam ! Combined flux in latest and established models = 0.531266 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.531266 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.307942Jy sigma=0.002579 ! Fit after self-cal, rms=0.307537Jy sigma=0.002576 ! Inverting map ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0791336 Jy ! Component: 100 - total flux cleaned = 0.12626 Jy ! Total flux subtracted in 100 components = 0.12626 Jy ! Clean residual min=-0.026480 max=0.047181 Jy/beam ! Clean residual mean=0.000362 rms=0.008116 Jy/beam ! Combined flux in latest and established models = 0.657525 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 31 components and 0.657525 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272018Jy sigma=0.002190 ! Fit after self-cal, rms=0.271983Jy sigma=0.002190 ! Inverting map ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7164 mas, bmaj=2.092 mas, bpa=-3.254 degrees ! Estimated noise=472.329 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 (4.8, 1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0964902 Jy ! Component: 100 - total flux cleaned = 0.160061 Jy ! Component: 150 - total flux cleaned = 0.20684 Jy ! Component: 200 - total flux cleaned = 0.242949 Jy ! Total flux subtracted in 200 components = 0.242949 Jy ! Clean residual min=-0.020059 max=0.054663 Jy/beam ! Clean residual mean=0.002412 rms=0.008758 Jy/beam ! Combined flux in latest and established models = 0.900474 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 85 components and 0.900474 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.231117Jy sigma=0.001752 ! Fit after self-cal, rms=0.227351Jy sigma=0.001711 ! Inverting map ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-3.2, -1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0494297 Jy ! Component: 100 - total flux cleaned = 0.0825177 Jy ! Component: 150 - total flux cleaned = 0.108329 Jy ! Component: 200 - total flux cleaned = 0.129316 Jy ! Total flux subtracted in 200 components = 0.129316 Jy ! Clean residual min=-0.015259 max=0.029909 Jy/beam ! Clean residual mean=0.001284 rms=0.005845 Jy/beam ! Combined flux in latest and established models = 1.02979 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 113 components and 1.02979 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.214192Jy sigma=0.001570 ! Fit after self-cal, rms=0.213556Jy sigma=0.001564 ! Inverting map ! Your choice of large map pixels excluded 0.374% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 51 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.94 GC 1.00* GG 1.00* ! HH 1.00* HN 1.00 KK 1.00* KP 1.09 ! LA 1.06 MC 1.05 MK 1.00* NL 1.03 ! NY 1.00* OV 0.99 PT 0.96 SC 0.98 ! TS 1.00* WF 0.93 WZ 1.02 ! ! A total of 156 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.89 FD 0.89 GC 1.61 GG 1.00* ! HH 1.00* HN 1.11 KK 0.92 KP 0.91 ! LA 0.85 MC 1.00* MK 0.98 NL 1.03 ! NY 1.00* OV 0.94 PT 0.96 SC 1.01 ! TS 1.00* WF 1.17 WZ 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 51 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 0.92 GC 1.00* GG 1.00* ! HH 1.00* HN 0.99 KK 1.00* KP 1.09 ! LA 1.05 MC 1.10 MK 1.00* NL 1.00 ! NY 1.00* OV 0.99 PT 0.93 SC 0.94 ! TS 1.00* WF 0.91 WZ 0.98 ! ! A total of 156 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.91 FD 0.88 GC 1.51 GG 1.00* ! HH 1.00* HN 1.09 KK 0.93 KP 0.89 ! LA 0.84 MC 1.00* MK 0.98 NL 1.00 ! NY 1.00* OV 0.94 PT 0.94 SC 0.98 ! TS 1.00* WF 1.13 WZ 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 51 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 0.91 GC 1.00* GG 1.00* ! HH 1.00* HN 1.01 KK 1.00* KP 1.14 ! LA 1.08 MC 1.15 MK 1.00* NL 1.02 ! NY 1.00* OV 1.03 PT 0.96 SC 0.95 ! TS 1.00* WF 0.92 WZ 0.92 ! ! A total of 156 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 0.87 GC 1.57 GG 1.00* ! HH 1.00* HN 1.11 KK 0.88 KP 0.93 ! LA 0.86 MC 1.00* MK 1.02 NL 1.01 ! NY 1.00* OV 0.97 PT 0.96 SC 1.01 ! TS 1.00* WF 1.11 WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MC 0.00* MK 0.00* NL 0.00* ! NY 0.00* OV 0.00* PT 0.00* SC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 51 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 0.90 GC 1.00* GG 1.00* ! HH 1.00* HN 1.04 KK 1.00* KP 1.10 ! LA 1.06 MC 1.09 MK 1.00* NL 1.03 ! NY 1.00* OV 1.01 PT 0.96 SC 0.99 ! TS 1.00* WF 0.94 WZ 0.98 ! ! A total of 156 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.91 FD 0.87 GC 1.58 GG 1.00* ! HH 1.00* HN 1.13 KK 0.87 KP 0.89 ! LA 0.85 MC 1.00* MK 1.03 NL 1.02 ! NY 1.00* OV 0.96 PT 0.97 SC 1.05 ! TS 1.00* WF 1.11 WZ 1.00* ! ! ! Fit before self-cal, rms=0.213556Jy sigma=0.001564 ! Fit after self-cal, rms=0.217656Jy sigma=0.001409 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7249 mas, bmaj=2.184 mas, bpa=-3.81 degrees ! Estimated noise=467.793 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 (-7.2, -4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.02215 Jy ! Component: 100 - total flux cleaned = 0.0331482 Jy ! Component: 150 - total flux cleaned = 0.0400566 Jy ! Component: 200 - total flux cleaned = 0.0448688 Jy ! Total flux subtracted in 200 components = 0.0448688 Jy ! Clean residual min=-0.008210 max=0.019995 Jy/beam ! Clean residual mean=0.001800 rms=0.003028 Jy/beam ! Combined flux in latest and established models = 1.07466 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 148 components and 1.07466 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.213929Jy sigma=0.001377 ! Fit after self-cal, rms=0.214214Jy sigma=0.001375 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (10.2, 5.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0188743 Jy ! Component: 100 - total flux cleaned = 0.0307465 Jy ! Component: 150 - total flux cleaned = 0.0382634 Jy ! Component: 200 - total flux cleaned = 0.0428563 Jy ! Total flux subtracted in 200 components = 0.0428563 Jy ! Clean residual min=-0.008158 max=0.013604 Jy/beam ! Clean residual mean=0.001556 rms=0.002571 Jy/beam ! Combined flux in latest and established models = 1.11751 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 178 components and 1.11751 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.211631Jy sigma=0.001352 ! Fit after self-cal, rms=0.211546Jy sigma=0.001351 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (3.2, 1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0107077 Jy ! Component: 100 - total flux cleaned = 0.014201 Jy ! Component: 150 - total flux cleaned = 0.0170828 Jy ! Component: 200 - total flux cleaned = 0.0194871 Jy ! Total flux subtracted in 200 components = 0.0194871 Jy ! Clean residual min=-0.007915 max=0.012350 Jy/beam ! Clean residual mean=0.001421 rms=0.002453 Jy/beam ! Combined flux in latest and established models = 1.137 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 199 components and 1.137 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.209936Jy sigma=0.001342 ! Fit after self-cal, rms=0.209894Jy sigma=0.001341 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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.209894Jy sigma=0.001341 ! Fit after self-cal, rms=0.298512Jy sigma=0.001272 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7262 mas, bmaj=2.191 mas, bpa=-3.698 degrees ! Estimated noise=468.197 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00755232 Jy ! Component: 100 - total flux cleaned = 0.0128797 Jy ! Component: 150 - total flux cleaned = 0.0178231 Jy ! Component: 200 - total flux cleaned = 0.021315 Jy ! Total flux subtracted in 200 components = 0.021315 Jy ! Clean residual min=-0.005827 max=0.010313 Jy/beam ! Clean residual mean=0.001252 rms=0.001893 Jy/beam ! Combined flux in latest and established models = 1.15832 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 225 components and 1.15832 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297778Jy sigma=0.001265 ! Fit after self-cal, rms=0.298710Jy sigma=0.001265 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 (-13.6, -7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00969013 Jy ! Component: 100 - total flux cleaned = 0.0155577 Jy ! Component: 150 - total flux cleaned = 0.0192216 Jy ! Component: 200 - total flux cleaned = 0.022408 Jy ! Total flux subtracted in 200 components = 0.022408 Jy ! Clean residual min=-0.005637 max=0.009799 Jy/beam ! Clean residual mean=0.001126 rms=0.001693 Jy/beam ! Combined flux in latest and established models = 1.18072 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 258 components and 1.18072 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.298167Jy sigma=0.001258 ! Fit after self-cal, rms=0.298179Jy sigma=0.001257 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-0.8, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0076939 Jy ! Component: 100 - total flux cleaned = 0.0107453 Jy ! Component: 150 - total flux cleaned = 0.013501 Jy ! Component: 200 - total flux cleaned = 0.0159012 Jy ! Total flux subtracted in 200 components = 0.0159012 Jy ! Clean residual min=-0.004962 max=0.007633 Jy/beam ! Clean residual mean=0.001014 rms=0.001600 Jy/beam ! Combined flux in latest and established models = 1.19663 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 276 components and 1.19663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297768Jy sigma=0.001253 ! Fit after self-cal, rms=0.297749Jy sigma=0.001252 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (13, -7.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00723781 Jy ! Component: 100 - total flux cleaned = 0.0103213 Jy ! Component: 150 - total flux cleaned = 0.0123062 Jy ! Component: 200 - total flux cleaned = 0.0143259 Jy ! Total flux subtracted in 200 components = 0.0143259 Jy ! Clean residual min=-0.004553 max=0.006797 Jy/beam ! Clean residual mean=0.000929 rms=0.001518 Jy/beam ! Combined flux in latest and established models = 1.21095 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 291 components and 1.21095 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297448Jy sigma=0.001249 ! Fit after self-cal, rms=0.297433Jy sigma=0.001248 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-4.6, 22.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00615388 Jy ! Component: 100 - total flux cleaned = 0.00810158 Jy ! Component: 150 - total flux cleaned = 0.00898682 Jy ! Component: 200 - total flux cleaned = 0.00983872 Jy ! Total flux subtracted in 200 components = 0.00983872 Jy ! Clean residual min=-0.004941 max=0.006594 Jy/beam ! Clean residual mean=0.000870 rms=0.001441 Jy/beam ! Combined flux in latest and established models = 1.22079 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 313 components and 1.22079 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297221Jy sigma=0.001246 ! Fit after self-cal, rms=0.297258Jy sigma=0.001245 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (5.4, -10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0060727 Jy ! Component: 100 - total flux cleaned = 0.00868947 Jy ! Component: 150 - total flux cleaned = 0.0103916 Jy ! Component: 200 - total flux cleaned = 0.0118714 Jy ! Total flux subtracted in 200 components = 0.0118714 Jy ! Clean residual min=-0.004827 max=0.005958 Jy/beam ! Clean residual mean=0.000803 rms=0.001343 Jy/beam ! Combined flux in latest and established models = 1.23266 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 330 components and 1.23266 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.297000Jy sigma=0.001243 ! Fit after self-cal, rms=0.296944Jy sigma=0.001243 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (0, -12). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00562176 Jy ! Component: 100 - total flux cleaned = 0.00880227 Jy ! Component: 150 - total flux cleaned = 0.0105029 Jy ! Component: 200 - total flux cleaned = 0.0119989 Jy ! Total flux subtracted in 200 components = 0.0119989 Jy ! Clean residual min=-0.004536 max=0.005574 Jy/beam ! Clean residual mean=0.000721 rms=0.001282 Jy/beam ! Combined flux in latest and established models = 1.24466 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 345 components and 1.24466 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.296730Jy sigma=0.001240 ! Fit after self-cal, rms=0.296716Jy sigma=0.001240 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-7.6, 4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00504929 Jy ! Component: 100 - total flux cleaned = 0.00702633 Jy ! Component: 150 - total flux cleaned = 0.00863398 Jy ! Component: 200 - total flux cleaned = 0.00970752 Jy ! Total flux subtracted in 200 components = 0.00970752 Jy ! Clean residual min=-0.004335 max=0.005155 Jy/beam ! Clean residual mean=0.000664 rms=0.001222 Jy/beam ! Combined flux in latest and established models = 1.25437 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 366 components and 1.25437 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.296534Jy sigma=0.001239 ! Fit after self-cal, rms=0.296534Jy sigma=0.001238 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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.296534Jy sigma=0.001238 ! Fit after self-cal, rms=0.314326Jy sigma=0.001235 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7275 mas, bmaj=2.183 mas, bpa=-3.664 degrees ! Estimated noise=467.349 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00451276 Jy ! Component: 100 - total flux cleaned = 0.00760249 Jy ! Component: 150 - total flux cleaned = 0.0103805 Jy ! Component: 200 - total flux cleaned = 0.0123365 Jy ! Total flux subtracted in 200 components = 0.0123365 Jy ! Clean residual min=-0.004008 max=0.006008 Jy/beam ! Clean residual mean=0.000564 rms=0.001119 Jy/beam ! Combined flux in latest and established models = 1.2667 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 393 components and 1.2667 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.314191Jy sigma=0.001233 ! Fit after self-cal, rms=0.314266Jy sigma=0.001233 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 (7, 2.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00572871 Jy ! Component: 100 - total flux cleaned = 0.00805285 Jy ! Component: 150 - total flux cleaned = 0.00949291 Jy ! Component: 200 - total flux cleaned = 0.0102723 Jy ! Total flux subtracted in 200 components = 0.0102723 Jy ! Clean residual min=-0.003624 max=0.004515 Jy/beam ! Clean residual mean=0.000499 rms=0.001064 Jy/beam ! Combined flux in latest and established models = 1.27698 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 420 components and 1.27698 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.314086Jy sigma=0.001230 ! Fit after self-cal, rms=0.314076Jy sigma=0.001230 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-9, -21.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00430588 Jy ! Component: 100 - total flux cleaned = 0.00558423 Jy ! Component: 150 - total flux cleaned = 0.00632615 Jy ! Component: 200 - total flux cleaned = 0.00652409 Jy ! Total flux subtracted in 200 components = 0.00652409 Jy ! Clean residual min=-0.003653 max=0.004069 Jy/beam ! Clean residual mean=0.000466 rms=0.001009 Jy/beam ! Combined flux in latest and established models = 1.2835 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 444 components and 1.2835 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.313958Jy sigma=0.001229 ! Fit after self-cal, rms=0.313966Jy sigma=0.001229 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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.313966Jy sigma=0.001229 ! Fit after self-cal, rms=0.315231Jy sigma=0.001228 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=2.16 mas, bmaj=3.476 mas, bpa=-3.875 degrees ! Estimated noise=588.985 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0057027 Jy ! Component: 100 - total flux cleaned = 0.00989559 Jy ! Component: 150 - total flux cleaned = 0.0129322 Jy ! Component: 200 - total flux cleaned = 0.0153344 Jy ! Total flux subtracted in 200 components = 0.0153344 Jy ! Clean residual min=-0.003669 max=0.005778 Jy/beam ! Clean residual mean=0.000834 rms=0.001260 Jy/beam ! Combined flux in latest and established models = 1.29883 Jy selfcal ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 479 components and 1.29883 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.315175Jy sigma=0.001227 ! Fit after self-cal, rms=0.315181Jy sigma=0.001227 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 (-15.2, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00531341 Jy ! Component: 100 - total flux cleaned = 0.00878144 Jy ! Component: 150 - total flux cleaned = 0.0113125 Jy ! Component: 200 - total flux cleaned = 0.0130397 Jy ! Total flux subtracted in 200 components = 0.0130397 Jy ! Clean residual min=-0.003448 max=0.004682 Jy/beam ! Clean residual mean=0.000653 rms=0.001141 Jy/beam ! Combined flux in latest and established models = 1.31187 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 514 components and 1.31187 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.315080Jy sigma=0.001226 ! Fit after self-cal, rms=0.315088Jy sigma=0.001226 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7276 mas, bmaj=2.181 mas, bpa=-3.652 degrees ! Estimated noise=467.698 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7276 x 2.181 at -3.652 degrees (North through East) ! Clean map min=-0.0053987 max=0.52462 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.315088Jy sigma=0.001226 ! Fit after self-cal, rms=0.551201Jy sigma=0.001133 ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7223 mas, bmaj=2.177 mas, bpa=-3.796 degrees ! Estimated noise=465.975 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00257727 Jy ! Component: 100 - total flux cleaned = 0.00396497 Jy ! Component: 150 - total flux cleaned = 0.00515117 Jy ! Component: 200 - total flux cleaned = 0.0056099 Jy ! Total flux subtracted in 200 components = 0.0056099 Jy ! Clean residual min=-0.003149 max=0.003431 Jy/beam ! Clean residual mean=0.000191 rms=0.000811 Jy/beam ! Combined flux in latest and established models = 1.31748 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 537 components and 1.31748 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.551147Jy sigma=0.001132 ! Fit after self-cal, rms=0.551140Jy sigma=0.001132 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000471626 Jy ! Component: 100 - total flux cleaned = 0.000560833 Jy ! Component: 150 - total flux cleaned = 0.00064618 Jy ! Component: 200 - total flux cleaned = 0.000728295 Jy ! Total flux subtracted in 200 components = 0.000728295 Jy ! Clean residual min=-0.003075 max=0.003374 Jy/beam ! Clean residual mean=0.000187 rms=0.000789 Jy/beam ! Combined flux in latest and established models = 1.31821 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 552 components and 1.31821 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.551108Jy sigma=0.001131 ! Fit after self-cal, rms=0.551100Jy sigma=0.001131 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (24.8, 14.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00272377 Jy ! Component: 100 - total flux cleaned = 0.00314303 Jy ! Component: 150 - total flux cleaned = 0.00322271 Jy ! Component: 200 - total flux cleaned = 0.00314756 Jy ! Total flux subtracted in 200 components = 0.00314756 Jy ! Clean residual min=-0.002912 max=0.003003 Jy/beam ! Clean residual mean=0.000169 rms=0.000762 Jy/beam ! Combined flux in latest and established models = 1.32136 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 573 components and 1.32136 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.551071Jy sigma=0.001130 ! Fit after self-cal, rms=0.551066Jy sigma=0.001130 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.551066Jy sigma=0.001130 ! Fit after self-cal, rms=0.544684Jy sigma=0.001130 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 ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.5768 mas, bmaj=1.76 mas, bpa=-3.417 degrees ! Estimated noise=876.037 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.381335 Jy ! Component: 100 - total flux cleaned = 0.535209 Jy ! Component: 150 - total flux cleaned = 0.617659 Jy ! Component: 200 - total flux cleaned = 0.676758 Jy ! Total flux subtracted in 200 components = 0.676758 Jy ! Clean residual min=-0.015158 max=0.035486 Jy/beam ! Clean residual mean=0.000412 rms=0.005760 Jy/beam ! Combined flux in latest and established models = 0.676758 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 27 model components to the UV plane model. ! The established model now contains 27 components and 0.676758 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7214 mas, bmaj=2.177 mas, bpa=-3.809 degrees ! Estimated noise=465.809 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.0885847 Jy ! Component: 100 - total flux cleaned = 0.1532 Jy ! Component: 150 - total flux cleaned = 0.202985 Jy ! Component: 200 - total flux cleaned = 0.242521 Jy ! Total flux subtracted in 200 components = 0.242521 Jy ! Clean residual min=-0.009499 max=0.031576 Jy/beam ! Clean residual mean=0.002941 rms=0.005580 Jy/beam ! Combined flux in latest and established models = 0.919279 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 81 components and 0.919279 Jy ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (0.8, 0). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0346239 Jy ! Component: 100 - total flux cleaned = 0.0624462 Jy ! Component: 150 - total flux cleaned = 0.0861247 Jy ! Component: 200 - total flux cleaned = 0.106751 Jy ! Total flux subtracted in 200 components = 0.106751 Jy ! Clean residual min=-0.006304 max=0.013493 Jy/beam ! Clean residual mean=0.002124 rms=0.003387 Jy/beam ! Combined flux in latest and established models = 1.02603 Jy ! Adding 78 model components to the UV plane model. ! The established model now contains 141 components and 1.02603 Jy ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec 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 ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=2.16 mas, bmaj=3.472 mas, bpa=-3.853 degrees ! Estimated noise=589.514 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0362117 Jy ! Component: 100 - total flux cleaned = 0.0649067 Jy ! Component: 150 - total flux cleaned = 0.0883853 Jy ! Component: 200 - total flux cleaned = 0.108098 Jy ! Total flux subtracted in 200 components = 0.108098 Jy ! Clean residual min=-0.005769 max=0.014605 Jy/beam ! Clean residual mean=0.003179 rms=0.003510 Jy/beam ! Combined flux in latest and established models = 1.13413 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.0167157 Jy ! Component: 100 - total flux cleaned = 0.0310162 Jy ! Component: 150 - total flux cleaned = 0.0433443 Jy ! Component: 200 - total flux cleaned = 0.0541366 Jy ! Total flux subtracted in 200 components = 0.0541366 Jy ! Clean residual min=-0.003626 max=0.009521 Jy/beam ! Clean residual mean=0.002351 rms=0.002143 Jy/beam ! Combined flux in latest and established models = 1.18827 Jy ! Adding 170 model components to the UV plane model. ! The established model now contains 283 components and 1.18827 Jy ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (-9.6, -4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0103195 Jy ! Component: 100 - total flux cleaned = 0.0189967 Jy ! Component: 150 - total flux cleaned = 0.0268109 Jy ! Component: 200 - total flux cleaned = 0.0339187 Jy ! Total flux subtracted in 200 components = 0.0339187 Jy ! Clean residual min=-0.002524 max=0.007256 Jy/beam ! Clean residual mean=0.001845 rms=0.001499 Jy/beam ! Combined flux in latest and established models = 1.22218 Jy ! Adding 76 model components to the UV plane model. ! The established model now contains 343 components and 1.22218 Jy ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Added new window around map position (12.8, 4.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00740043 Jy ! Component: 100 - total flux cleaned = 0.0135646 Jy ! Component: 150 - total flux cleaned = 0.0191929 Jy ! Component: 200 - total flux cleaned = 0.0243815 Jy ! Total flux subtracted in 200 components = 0.0243815 Jy ! Clean residual min=-0.002794 max=0.004588 Jy/beam ! Clean residual mean=0.001486 rms=0.001144 Jy/beam ! Combined flux in latest and established models = 1.24657 Jy ! Adding 63 model components to the UV plane model. ! The established model now contains 390 components and 1.24657 Jy ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! Estimated beam: bmin=0.7214 mas, bmaj=2.177 mas, bpa=-3.809 degrees ! Estimated noise=465.809 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00551523 Jy ! Component: 100 - total flux cleaned = 0.00945817 Jy ! Component: 150 - total flux cleaned = 0.0125065 Jy ! Component: 200 - total flux cleaned = 0.0158293 Jy ! Total flux subtracted in 200 components = 0.0158293 Jy ! Clean residual min=-0.002544 max=0.004171 Jy/beam ! Clean residual mean=0.000539 rms=0.000868 Jy/beam ! Combined flux in latest and established models = 1.2624 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 427 components and 1.26239 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.545143Jy sigma=0.001134 ! Fit after self-cal, rms=0.545027Jy sigma=0.001132 wmodel J0241-0815_X_2001_07_05_pus_map.mod ! Writing 427 model components to file: J0241-0815_X_2001_07_05_pus_map.mod wobs J0241-0815_X_2001_07_05_pus_uvs.fits ! Writing UV FITS file: J0241-0815_X_2001_07_05_pus_uvs.fits wwins J0241-0815_X_2001_07_05_pus_map.win ! wwins: Wrote 20 windows to J0241-0815_X_2001_07_05_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00174328 Jy ! Component: 100 - total flux cleaned = 0.00329691 Jy ! Component: 150 - total flux cleaned = 0.00478337 Jy ! Component: 200 - total flux cleaned = 0.00621917 Jy ! Component: 250 - total flux cleaned = 0.0076142 Jy ! Component: 300 - total flux cleaned = 0.0089738 Jy ! Component: 350 - total flux cleaned = 0.0103022 Jy ! Component: 400 - total flux cleaned = 0.0116008 Jy ! Component: 450 - total flux cleaned = 0.0128715 Jy ! Component: 500 - total flux cleaned = 0.014117 Jy ! Component: 550 - total flux cleaned = 0.0153388 Jy ! Component: 600 - total flux cleaned = 0.016396 Jy ! Component: 650 - total flux cleaned = 0.017483 Jy ! Component: 700 - total flux cleaned = 0.0184606 Jy ! Component: 750 - total flux cleaned = 0.0194699 Jy ! Component: 800 - total flux cleaned = 0.0205098 Jy ! Component: 850 - total flux cleaned = 0.0214907 Jy ! Component: 900 - total flux cleaned = 0.0224146 Jy ! Component: 950 - total flux cleaned = 0.0233268 Jy ! Component: 1000 - total flux cleaned = 0.0242276 Jy ! Component: 1050 - total flux cleaned = 0.025117 Jy ! Component: 1100 - total flux cleaned = 0.0259957 Jy ! Component: 1150 - total flux cleaned = 0.0268639 Jy ! Component: 1200 - total flux cleaned = 0.0276815 Jy ! Component: 1250 - total flux cleaned = 0.0285305 Jy ! Component: 1300 - total flux cleaned = 0.0292902 Jy ! Component: 1350 - total flux cleaned = 0.0301605 Jy ! Component: 1400 - total flux cleaned = 0.0308651 Jy ! Component: 1450 - total flux cleaned = 0.0316015 Jy ! Component: 1500 - total flux cleaned = 0.0322538 Jy ! Component: 1550 - total flux cleaned = 0.0329005 Jy ! Component: 1600 - total flux cleaned = 0.0336163 Jy ! Component: 1650 - total flux cleaned = 0.0343628 Jy ! Component: 1700 - total flux cleaned = 0.0349914 Jy ! Component: 1750 - total flux cleaned = 0.0355412 Jy ! Component: 1800 - total flux cleaned = 0.036268 Jy ! Component: 1850 - total flux cleaned = 0.0368444 Jy ! Component: 1900 - total flux cleaned = 0.0374873 Jy ! Component: 1950 - total flux cleaned = 0.0380895 Jy ! Component: 2000 - total flux cleaned = 0.0386864 Jy ! Total flux subtracted in 2048 components = 0.0392089 Jy ! Clean residual min=-0.002087 max=0.001755 Jy/beam ! Clean residual mean=0.000299 rms=0.000605 Jy/beam ! Combined flux in latest and established models = 1.3016 Jy keep ! Adding 150 model components to the UV plane model. ! The established model now contains 567 components and 1.3016 Jy wmap "J0241-0815_X_2001_07_05_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7214 x 2.177 at -3.809 degrees (North through East) ! Clean map min=-0.0022839 max=0.51765 Jy/beam ! Writing clean map to FITS file: J0241-0815_X_2001_07_05_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.517651 invert ! Inverting map ! Your choice of large map pixels excluded 0.383% of the data. ! The x-axis pixel size should ideally be below 0.1947 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000604112 device "J0241-0815_X_2001_07_05_pus\_map.eps/VPS" ! Attempting to open device: 'J0241-0815_X_2001_07_05_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.350108 0.350108 0.700215 1.40043 2.80086 5.60172 11.2034 22.4069 44.8138 89.6275 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0018 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.517 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7214 x 2.177 at -3.809 degrees (North through East) ! Clean map min=-0.0022839 max=0.51765 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J0241-0815_X_2001_07_05_pus\_map.gif/vgif" ! Attempting to open device: 'J0241-0815_X_2001_07_05_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J0241-0815_X_2001_07_05_pus_uvs.log closed on Thu Jul 5 16:03:46 2007