! Started logfile: J2038+5119_X_2010_06_23_pus_uvs.log on Thu Jan 27 22:06:13 2011 obs J2038+5119_X_2010_06_23_pus_uva.fits ! Reading UV FITS file: J2038+5119_X_2010_06_23_pus_uva.fits ! AN table 1: 32 integrations on 105 of 105 possible baselines. ! AN table 2: 23 integrations on 105 of 105 possible baselines. ! AN table 3: 63 integrations on 105 of 105 possible baselines. ! AN table 4: 90 integrations on 105 of 105 possible baselines. ! Apparent sampling: 0.231319 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2038+5119 ! ! 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.40996e+09 8e+06 1 8e+06 ! 02 2 8.47996e+09 8e+06 1 8e+06 ! 03 3 8.79496e+09 8e+06 1 8e+06 ! 04 4 8.89996e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 3588 lines of history. ! ! Reading 20208 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 ![@muppet J2038+5119_X_2010_06_23_pus] 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 394 telescope corrections were flagged in sub-array 1. ! A total of 270 telescope corrections were flagged in sub-array 2. ! A total of 622 telescope corrections were flagged in sub-array 3. ! A total of 611 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 408 telescope corrections were flagged in sub-array 1. ! A total of 302 telescope corrections were flagged in sub-array 2. ! A total of 671 telescope corrections were flagged in sub-array 3. ! A total of 635 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 394 telescope corrections were flagged in sub-array 1. ! A total of 270 telescope corrections were flagged in sub-array 2. ! A total of 622 telescope corrections were flagged in sub-array 3. ! A total of 611 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 394 telescope corrections were flagged in sub-array 1. ! A total of 270 telescope corrections were flagged in sub-array 2. ! A total of 622 telescope corrections were flagged in sub-array 3. ! A total of 611 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.567108Jy sigma=0.006674 ! Fit after self-cal, rms=0.505131Jy sigma=0.006041 ! 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 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.5041 mas, bmaj=0.6978 mas, bpa=12.54 degrees ! Estimated noise=727.031 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.92478 Jy ! Component: 100 - total flux cleaned = 1.18846 Jy ! Total flux subtracted in 100 components = 1.18846 Jy ! Clean residual min=-0.024849 max=0.122389 Jy/beam ! Clean residual mean=0.000210 rms=0.007951 Jy/beam ! Combined flux in latest and established models = 1.18846 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 1.18846 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.401269Jy sigma=0.004783 ! Fit after self-cal, rms=0.400614Jy sigma=0.004777 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.141009 Jy ! Component: 100 - total flux cleaned = 0.233429 Jy ! Total flux subtracted in 100 components = 0.233429 Jy ! Clean residual min=-0.018681 max=0.051164 Jy/beam ! Clean residual mean=0.000198 rms=0.005530 Jy/beam ! Combined flux in latest and established models = 1.42189 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 16 components and 1.42189 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.306317Jy sigma=0.003531 ! Fit after self-cal, rms=0.305729Jy sigma=0.003526 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0664515 Jy ! Component: 100 - total flux cleaned = 0.11326 Jy ! Total flux subtracted in 100 components = 0.11326 Jy ! Clean residual min=-0.017272 max=0.036217 Jy/beam ! Clean residual mean=0.000192 rms=0.004735 Jy/beam ! Combined flux in latest and established models = 1.53515 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 18 components and 1.53515 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.277087Jy sigma=0.003093 ! Fit after self-cal, rms=0.276802Jy sigma=0.003090 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-4.2, -5.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0397242 Jy ! Component: 100 - total flux cleaned = 0.0718918 Jy ! Total flux subtracted in 100 components = 0.0718918 Jy ! Clean residual min=-0.015926 max=0.031231 Jy/beam ! Clean residual mean=0.000180 rms=0.004384 Jy/beam ! Combined flux in latest and established models = 1.60704 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 29 components and 1.60704 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262748Jy sigma=0.002872 ! Fit after self-cal, rms=0.262459Jy sigma=0.002869 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-5.2, -5.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0353495 Jy ! Component: 100 - total flux cleaned = 0.062792 Jy ! Total flux subtracted in 100 components = 0.062792 Jy ! Clean residual min=-0.015407 max=0.028264 Jy/beam ! Clean residual mean=0.000163 rms=0.004096 Jy/beam ! Combined flux in latest and established models = 1.66983 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 44 components and 1.66983 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.250530Jy sigma=0.002684 ! Fit after self-cal, rms=0.250267Jy sigma=0.002682 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-0.4, 1.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0281899 Jy ! Component: 100 - total flux cleaned = 0.0509411 Jy ! Total flux subtracted in 100 components = 0.0509411 Jy ! Clean residual min=-0.014435 max=0.025721 Jy/beam ! Clean residual mean=0.000153 rms=0.003874 Jy/beam ! Combined flux in latest and established models = 1.72078 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 51 components and 1.72078 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.242709Jy sigma=0.002557 ! Fit after self-cal, rms=0.242567Jy sigma=0.002555 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-1.6, -1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0255609 Jy ! Component: 100 - total flux cleaned = 0.0455261 Jy ! Total flux subtracted in 100 components = 0.0455261 Jy ! Clean residual min=-0.014218 max=0.023618 Jy/beam ! Clean residual mean=0.000145 rms=0.003710 Jy/beam ! Combined flux in latest and established models = 1.7663 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 58 components and 1.7663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.236522Jy sigma=0.002457 ! Fit after self-cal, rms=0.236383Jy sigma=0.002456 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (1, 3.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0222542 Jy ! Component: 100 - total flux cleaned = 0.0396927 Jy ! Total flux subtracted in 100 components = 0.0396927 Jy ! Clean residual min=-0.013561 max=0.021824 Jy/beam ! Clean residual mean=0.000137 rms=0.003563 Jy/beam ! Combined flux in latest and established models = 1.80599 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 64 components and 1.80599 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.232360Jy sigma=0.002388 ! Fit after self-cal, rms=0.232331Jy sigma=0.002387 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-6.2, -6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0227669 Jy ! Component: 100 - total flux cleaned = 0.0393894 Jy ! Total flux subtracted in 100 components = 0.0393894 Jy ! Clean residual min=-0.013297 max=0.020042 Jy/beam ! Clean residual mean=0.000126 rms=0.003439 Jy/beam ! Combined flux in latest and established models = 1.84538 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 80 components and 1.84538 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227131Jy sigma=0.002306 ! Fit after self-cal, rms=0.226938Jy sigma=0.002304 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 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 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6407 mas, bmaj=0.8606 mas, bpa=10.89 degrees ! Estimated noise=486.975 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 (-5.6, -6.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0539949 Jy ! Component: 100 - total flux cleaned = 0.0951868 Jy ! Component: 150 - total flux cleaned = 0.129708 Jy ! Component: 200 - total flux cleaned = 0.159194 Jy ! Total flux subtracted in 200 components = 0.159194 Jy ! Clean residual min=-0.017532 max=0.029289 Jy/beam ! Clean residual mean=0.000754 rms=0.005083 Jy/beam ! Combined flux in latest and established models = 2.00458 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 114 components and 2.00458 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.211759Jy sigma=0.002028 ! Fit after self-cal, rms=0.210931Jy sigma=0.002018 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-2.4, -3.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.030902 Jy ! Component: 100 - total flux cleaned = 0.0556643 Jy ! Component: 150 - total flux cleaned = 0.0773054 Jy ! Component: 200 - total flux cleaned = 0.0963797 Jy ! Total flux subtracted in 200 components = 0.0963797 Jy ! Clean residual min=-0.016408 max=0.022018 Jy/beam ! Clean residual mean=0.000375 rms=0.004322 Jy/beam ! Combined flux in latest and established models = 2.10096 Jy ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 142 components and 2.10096 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.205181Jy sigma=0.001908 ! Fit after self-cal, rms=0.204938Jy sigma=0.001905 ! Inverting map ! Your choice of large map pixels excluded 4.55% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.00* HN 1.00* KK 1.17 ! KP 1.00* LA 1.00* MA 1.05 MK 1.00* ! NL 1.00* NY 1.09 OV 1.00* PT 1.00* ! SC 1.00* TS 0.97 ZC 1.03 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.00* HN 1.24 KK 1.00* ! KP 1.00* LA 1.00* MA 1.01 MK 1.00* ! NL 1.00* NY 1.06 OV 1.00* PT 1.00* ! SC 1.00* TS 0.99 ZC 1.07 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.95 FD 1.03 HN 1.15 KK 1.00* ! KP 1.00* LA 0.97 MA 1.10 MK 1.00* ! NL 0.99 NY 1.06 OV 1.00* PT 1.06 ! SC 1.07 TS 1.03 ZC 1.05 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.96 FD 0.98 HN 1.07 KK 0.98 ! KP 0.98 LA 0.95 MA 1.04 MK 1.00 ! NL 0.99 NY 1.04 OV 0.98 PT 0.96 ! SC 1.01 TS 0.98 ZC 1.01 ! ! ! Correcting IF 2. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00* HN 1.00* KK 1.16 ! KP 1.00* LA 1.00* MA 1.00* MK 1.00* ! NL 1.00* NY 1.11 OV 1.00* PT 1.00* ! SC 1.00* TS 0.98 ZC 1.03 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 1.00* HN 1.29 KK 1.00* ! KP 1.00* LA 1.00* MA 1.00* MK 1.00* ! NL 1.00* NY 1.07 OV 1.00* PT 1.00* ! SC 1.00* TS 0.99 ZC 1.08 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 1.03 HN 1.16 KK 1.00* ! KP 1.00* LA 1.01 MA 1.00* MK 1.00* ! NL 1.01 NY 1.06 OV 1.00* PT 1.05 ! SC 1.09 TS 1.06 ZC 1.05 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 FD 0.97 HN 1.09 KK 1.02 ! KP 0.98 LA 0.98 MA 1.00* MK 0.99 ! NL 1.01 NY 1.04 OV 1.00 PT 0.95 ! SC 1.03 TS 0.98 ZC 1.01 ! ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.00* HN 1.00* KK 1.12 ! KP 1.00* LA 1.00* MA 1.02 MK 1.00* ! NL 1.00* NY 1.09 OV 1.00* PT 1.00* ! SC 1.00* TS 0.99 ZC 1.00 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 1.00* HN 1.33 KK 1.00* ! KP 1.00* LA 1.00* MA 1.04 MK 1.00* ! NL 1.00* NY 1.08 OV 1.00* PT 1.00* ! SC 1.00* TS 0.99 ZC 1.02 ! ! 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 HN 1.15 KK 1.00* ! KP 1.00* LA 0.99 MA 0.97 MK 1.00* ! NL 0.95 NY 1.07 OV 1.00* PT 1.06 ! SC 1.05 TS 1.02 ZC 1.05 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.97 FD 0.96 HN 1.10 KK 0.95 ! KP 0.97 LA 0.97 MA 1.03 MK 0.98 ! NL 0.96 NY 1.04 OV 0.99 PT 0.94 ! SC 0.99 TS 0.96 ZC 0.98 ! ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.00* HN 1.00* KK 1.09 ! KP 1.00* LA 1.00* MA 1.02 MK 1.00* ! NL 1.00* NY 1.08 OV 1.00* PT 1.00* ! SC 1.00* TS 0.99 ZC 1.06 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 1.00* HN 1.29 KK 1.00* ! KP 1.00* LA 1.00* MA 1.02 MK 1.00* ! NL 1.00* NY 1.12 OV 1.00* PT 1.00* ! SC 1.00* TS 0.96 ZC 1.10 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.04 FD 1.01 HN 1.15 KK 1.00* ! KP 1.00* LA 0.98 MA 1.04 MK 1.00* ! NL 0.98 NY 1.05 OV 1.00* PT 1.08 ! SC 1.09 TS 1.02 ZC 1.09 ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 0.97 HN 1.09 KK 0.94 ! KP 0.98 LA 0.96 MA 1.04 MK 0.98 ! NL 0.99 NY 1.04 OV 0.97 PT 0.95 ! SC 1.00 TS 0.96 ZC 1.02 ! ! ! Fit before self-cal, rms=0.204938Jy sigma=0.001905 ! Fit after self-cal, rms=0.188697Jy sigma=0.001646 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6451 mas, bmaj=0.8676 mas, bpa=10.76 degrees ! Estimated noise=487.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 (-3.2, -6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.014807 Jy ! Component: 100 - total flux cleaned = 0.0224181 Jy ! Component: 150 - total flux cleaned = 0.028223 Jy ! Component: 200 - total flux cleaned = 0.033566 Jy ! Total flux subtracted in 200 components = 0.033566 Jy ! Clean residual min=-0.009637 max=0.016775 Jy/beam ! Clean residual mean=0.000491 rms=0.002808 Jy/beam ! Combined flux in latest and established models = 2.13452 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 168 components and 2.13452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185879Jy sigma=0.001594 ! Fit after self-cal, rms=0.185817Jy sigma=0.001592 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-5.2, -8.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0148758 Jy ! Component: 100 - total flux cleaned = 0.0196766 Jy ! Component: 150 - total flux cleaned = 0.0240625 Jy ! Component: 200 - total flux cleaned = 0.0273542 Jy ! Total flux subtracted in 200 components = 0.0273542 Jy ! Clean residual min=-0.008887 max=0.013503 Jy/beam ! Clean residual mean=0.000377 rms=0.002607 Jy/beam ! Combined flux in latest and established models = 2.16188 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 189 components and 2.16188 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.184228Jy sigma=0.001565 ! Fit after self-cal, rms=0.184056Jy sigma=0.001563 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-9.2, -8.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0120617 Jy ! Component: 100 - total flux cleaned = 0.0155294 Jy ! Component: 150 - total flux cleaned = 0.018745 Jy ! Component: 200 - total flux cleaned = 0.0211836 Jy ! Total flux subtracted in 200 components = 0.0211836 Jy ! Clean residual min=-0.008651 max=0.013013 Jy/beam ! Clean residual mean=0.000289 rms=0.002448 Jy/beam ! Combined flux in latest and established models = 2.18306 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 207 components and 2.18306 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183044Jy sigma=0.001547 ! Fit after self-cal, rms=0.182952Jy sigma=0.001546 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-0.4, -2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0104716 Jy ! Component: 100 - total flux cleaned = 0.0125925 Jy ! Component: 150 - total flux cleaned = 0.0146997 Jy ! Component: 200 - total flux cleaned = 0.0168548 Jy ! Total flux subtracted in 200 components = 0.0168548 Jy ! Clean residual min=-0.008368 max=0.011672 Jy/beam ! Clean residual mean=0.000224 rms=0.002342 Jy/beam ! Combined flux in latest and established models = 2.19992 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 221 components and 2.19992 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182140Jy sigma=0.001532 ! Fit after self-cal, rms=0.182092Jy sigma=0.001531 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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.182092Jy sigma=0.001531 ! Fit after self-cal, rms=0.177818Jy sigma=0.001474 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6459 mas, bmaj=0.8689 mas, bpa=10.59 degrees ! Estimated noise=486.993 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00334641 Jy ! Component: 100 - total flux cleaned = -0.00358785 Jy ! Component: 150 - total flux cleaned = -0.00380257 Jy ! Component: 200 - total flux cleaned = -0.00341441 Jy ! Total flux subtracted in 200 components = -0.00341441 Jy ! Clean residual min=-0.007223 max=0.012530 Jy/beam ! Clean residual mean=0.000180 rms=0.001990 Jy/beam ! Combined flux in latest and established models = 2.1965 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 234 components and 2.1965 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177236Jy sigma=0.001466 ! Fit after self-cal, rms=0.177191Jy sigma=0.001466 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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 (-6, -4.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0130266 Jy ! Component: 100 - total flux cleaned = 0.0156377 Jy ! Component: 150 - total flux cleaned = 0.0154049 Jy ! Component: 200 - total flux cleaned = 0.0160698 Jy ! Total flux subtracted in 200 components = 0.0160698 Jy ! Clean residual min=-0.006575 max=0.008238 Jy/beam ! Clean residual mean=0.000112 rms=0.001825 Jy/beam ! Combined flux in latest and established models = 2.21257 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 253 components and 2.21257 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176085Jy sigma=0.001448 ! Fit after self-cal, rms=0.176056Jy sigma=0.001447 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-4.6, -4.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00529941 Jy ! Component: 100 - total flux cleaned = 0.00505574 Jy ! Component: 150 - total flux cleaned = 0.00524503 Jy ! Component: 200 - total flux cleaned = 0.00524608 Jy ! Total flux subtracted in 200 components = 0.00524608 Jy ! Clean residual min=-0.006261 max=0.007575 Jy/beam ! Clean residual mean=0.000092 rms=0.001775 Jy/beam ! Combined flux in latest and established models = 2.21782 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 264 components and 2.21782 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.175488Jy sigma=0.001441 ! Fit after self-cal, rms=0.175468Jy sigma=0.001441 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-13.6, -16.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00740621 Jy ! Component: 100 - total flux cleaned = 0.00948975 Jy ! Component: 150 - total flux cleaned = 0.0103245 Jy ! Component: 200 - total flux cleaned = 0.010795 Jy ! Total flux subtracted in 200 components = 0.010795 Jy ! Clean residual min=-0.006431 max=0.006955 Jy/beam ! Clean residual mean=0.000059 rms=0.001690 Jy/beam ! Combined flux in latest and established models = 2.22861 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 277 components and 2.22861 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.175047Jy sigma=0.001435 ! Fit after self-cal, rms=0.175018Jy sigma=0.001434 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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.175018Jy sigma=0.001434 ! Fit after self-cal, rms=0.171658Jy sigma=0.001397 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6463 mas, bmaj=0.8699 mas, bpa=10.54 degrees ! Estimated noise=486.949 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000379357 Jy ! Component: 100 - total flux cleaned = -0.000206296 Jy ! Component: 150 - total flux cleaned = -0.000206284 Jy ! Component: 200 - total flux cleaned = -3.69814e-05 Jy ! Total flux subtracted in 200 components = -3.69814e-05 Jy ! Clean residual min=-0.004764 max=0.006025 Jy/beam ! Clean residual mean=0.000021 rms=0.001291 Jy/beam ! Combined flux in latest and established models = 2.22857 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 291 components and 2.22857 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171241Jy sigma=0.001392 ! Fit after self-cal, rms=0.171235Jy sigma=0.001392 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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 (-10.2, -16.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00626859 Jy ! Component: 100 - total flux cleaned = 0.00750709 Jy ! Component: 150 - total flux cleaned = 0.00782524 Jy ! Component: 200 - total flux cleaned = 0.00812582 Jy ! Total flux subtracted in 200 components = 0.00812582 Jy ! Clean residual min=-0.004723 max=0.005858 Jy/beam ! Clean residual mean=-0.000004 rms=0.001222 Jy/beam ! Combined flux in latest and established models = 2.2367 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 307 components and 2.2367 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170894Jy sigma=0.001387 ! Fit after self-cal, rms=0.170891Jy sigma=0.001386 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-2.8, -8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00324504 Jy ! Component: 100 - total flux cleaned = 0.00338605 Jy ! Component: 150 - total flux cleaned = 0.00381041 Jy ! Component: 200 - total flux cleaned = 0.00394963 Jy ! Total flux subtracted in 200 components = 0.00394963 Jy ! Clean residual min=-0.004762 max=0.005152 Jy/beam ! Clean residual mean=-0.000021 rms=0.001180 Jy/beam ! Combined flux in latest and established models = 2.24065 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 321 components and 2.24065 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170602Jy sigma=0.001383 ! Fit after self-cal, rms=0.170595Jy sigma=0.001383 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-10.8, -21). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00556967 Jy ! Component: 100 - total flux cleaned = 0.00725633 Jy ! Component: 150 - total flux cleaned = 0.00752236 Jy ! Component: 200 - total flux cleaned = 0.00752079 Jy ! Total flux subtracted in 200 components = 0.00752079 Jy ! Clean residual min=-0.004743 max=0.004714 Jy/beam ! Clean residual mean=-0.000035 rms=0.001109 Jy/beam ! Combined flux in latest and established models = 2.24817 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 341 components and 2.24817 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170329Jy sigma=0.001379 ! Fit after self-cal, rms=0.170316Jy sigma=0.001379 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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.170316Jy sigma=0.001379 ! Fit after self-cal, rms=0.170429Jy sigma=0.001377 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=2.456 mas, bmaj=2.825 mas, bpa=-67.49 degrees ! Estimated noise=683.043 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00365348 Jy ! Component: 100 - total flux cleaned = 0.00621518 Jy ! Component: 150 - total flux cleaned = 0.00821254 Jy ! Component: 200 - total flux cleaned = 0.00944978 Jy ! Total flux subtracted in 200 components = 0.00944978 Jy ! Clean residual min=-0.006581 max=0.006766 Jy/beam ! Clean residual mean=-0.000132 rms=0.001979 Jy/beam ! Combined flux in latest and established models = 2.25762 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 367 components and 2.25762 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170332Jy sigma=0.001375 ! Fit after self-cal, rms=0.170324Jy sigma=0.001375 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00145384 Jy ! Component: 100 - total flux cleaned = 0.00239578 Jy ! Component: 150 - total flux cleaned = 0.00296832 Jy ! Component: 200 - total flux cleaned = 0.00336663 Jy ! Total flux subtracted in 200 components = 0.00336663 Jy ! Clean residual min=-0.006612 max=0.006572 Jy/beam ! Clean residual mean=-0.000168 rms=0.001889 Jy/beam ! Combined flux in latest and established models = 2.26099 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 383 components and 2.26099 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.170277Jy sigma=0.001374 ! Fit after self-cal, rms=0.170268Jy sigma=0.001374 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6467 mas, bmaj=0.8707 mas, bpa=10.51 degrees ! Estimated noise=486.959 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6467 x 0.8707 at 10.51 degrees (North through East) ! Clean map min=-0.0039496 max=1.2874 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.170268Jy sigma=0.001374 ! Fit after self-cal, rms=0.161314Jy sigma=0.001268 ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6468 mas, bmaj=0.8709 mas, bpa=10.51 degrees ! Estimated noise=486.548 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00314267 Jy ! Component: 100 - total flux cleaned = -0.00437445 Jy ! Component: 150 - total flux cleaned = -0.00476118 Jy ! Component: 200 - total flux cleaned = -0.00512793 Jy ! Total flux subtracted in 200 components = -0.00512793 Jy ! Clean residual min=-0.004392 max=0.004029 Jy/beam ! Clean residual mean=0.000002 rms=0.000991 Jy/beam ! Combined flux in latest and established models = 2.25586 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 393 components and 2.25586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161086Jy sigma=0.001266 ! Fit after self-cal, rms=0.161083Jy sigma=0.001266 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000124456 Jy ! Component: 100 - total flux cleaned = -0.000468707 Jy ! Component: 150 - total flux cleaned = -0.000576732 Jy ! Component: 200 - total flux cleaned = -0.000575544 Jy ! Total flux subtracted in 200 components = -0.000575544 Jy ! Clean residual min=-0.004265 max=0.003995 Jy/beam ! Clean residual mean=0.000005 rms=0.000977 Jy/beam ! Combined flux in latest and established models = 2.25528 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 403 components and 2.25528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160927Jy sigma=0.001264 ! Fit after self-cal, rms=0.160927Jy sigma=0.001264 ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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.160927Jy sigma=0.001264 ! Fit after self-cal, rms=0.160894Jy sigma=0.001264 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 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.5035 mas, bmaj=0.7027 mas, bpa=12.04 degrees ! Estimated noise=736.766 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.936835 Jy ! Component: 100 - total flux cleaned = 1.2049 Jy ! Component: 150 - total flux cleaned = 1.34109 Jy ! Component: 200 - total flux cleaned = 1.42991 Jy ! Total flux subtracted in 200 components = 1.42991 Jy ! Clean residual min=-0.013781 max=0.048422 Jy/beam ! Clean residual mean=0.000224 rms=0.004393 Jy/beam ! Combined flux in latest and established models = 1.42992 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 13 model components to the UV plane model. ! The established model now contains 13 components and 1.42992 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6472 mas, bmaj=0.8715 mas, bpa=10.5 degrees ! Estimated noise=486.593 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.122257 Jy ! Component: 100 - total flux cleaned = 0.213639 Jy ! Component: 150 - total flux cleaned = 0.285753 Jy ! Component: 200 - total flux cleaned = 0.344559 Jy ! Total flux subtracted in 200 components = 0.344559 Jy ! Clean residual min=-0.011369 max=0.037716 Jy/beam ! Clean residual mean=0.001779 rms=0.005639 Jy/beam ! Combined flux in latest and established models = 1.77447 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 83 components and 1.77447 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-4.4, -6.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0494838 Jy ! Component: 100 - total flux cleaned = 0.0913751 Jy ! Component: 150 - total flux cleaned = 0.127294 Jy ! Component: 200 - total flux cleaned = 0.158359 Jy ! Total flux subtracted in 200 components = 0.158359 Jy ! Clean residual min=-0.007174 max=0.020735 Jy/beam ! Clean residual mean=0.001175 rms=0.003538 Jy/beam ! Combined flux in latest and established models = 1.93283 Jy ! Adding 96 model components to the UV plane model. ! The established model now contains 153 components and 1.93283 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-6.8, -7.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0272404 Jy ! Component: 100 - total flux cleaned = 0.0510815 Jy ! Component: 150 - total flux cleaned = 0.072179 Jy ! Component: 200 - total flux cleaned = 0.0909824 Jy ! Total flux subtracted in 200 components = 0.0909824 Jy ! Clean residual min=-0.005695 max=0.013828 Jy/beam ! Clean residual mean=0.000828 rms=0.002500 Jy/beam ! Combined flux in latest and established models = 2.02381 Jy ! Adding 105 model components to the UV plane model. ! The established model now contains 216 components and 2.02382 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-0.6, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0171088 Jy ! Component: 100 - total flux cleaned = 0.032245 Jy ! Component: 150 - total flux cleaned = 0.0458893 Jy ! Component: 200 - total flux cleaned = 0.0582683 Jy ! Total flux subtracted in 200 components = 0.0582683 Jy ! Clean residual min=-0.004825 max=0.008618 Jy/beam ! Clean residual mean=0.000607 rms=0.001928 Jy/beam ! Combined flux in latest and established models = 2.08208 Jy ! Adding 108 model components to the UV plane model. ! The established model now contains 272 components and 2.08208 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-3.4, -3.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0113997 Jy ! Component: 100 - total flux cleaned = 0.0217944 Jy ! Component: 150 - total flux cleaned = 0.0314124 Jy ! Component: 200 - total flux cleaned = 0.0403768 Jy ! Total flux subtracted in 200 components = 0.0403768 Jy ! Clean residual min=-0.004326 max=0.006851 Jy/beam ! Clean residual mean=0.000461 rms=0.001586 Jy/beam ! Combined flux in latest and established models = 2.12246 Jy ! Adding 115 model components to the UV plane model. ! The established model now contains 322 components and 2.12246 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Added new window around map position (-7.4, -9.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00856965 Jy ! Component: 100 - total flux cleaned = 0.0164033 Jy ! Component: 150 - total flux cleaned = 0.0237163 Jy ! Component: 200 - total flux cleaned = 0.0305664 Jy ! Total flux subtracted in 200 components = 0.0305664 Jy ! Clean residual min=-0.004001 max=0.005635 Jy/beam ! Clean residual mean=0.000352 rms=0.001350 Jy/beam ! Combined flux in latest and established models = 2.15303 Jy ! Adding 106 model components to the UV plane model. ! The established model now contains 362 components and 2.15303 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 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 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=2.458 mas, bmaj=2.826 mas, bpa=-67.36 degrees ! Estimated noise=681.304 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0167118 Jy ! Component: 100 - total flux cleaned = 0.0302208 Jy ! Component: 150 - total flux cleaned = 0.0414745 Jy ! Component: 200 - total flux cleaned = 0.050951 Jy ! Total flux subtracted in 200 components = 0.050951 Jy ! Clean residual min=-0.005867 max=0.007876 Jy/beam ! Clean residual mean=0.000549 rms=0.002293 Jy/beam ! Combined flux in latest and established models = 2.20398 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.00803809 Jy ! Component: 100 - total flux cleaned = 0.014905 Jy ! Component: 150 - total flux cleaned = 0.0208031 Jy ! Component: 200 - total flux cleaned = 0.0258839 Jy ! Total flux subtracted in 200 components = 0.0258839 Jy ! Clean residual min=-0.005541 max=0.006963 Jy/beam ! Clean residual mean=0.000285 rms=0.001850 Jy/beam ! Combined flux in latest and established models = 2.22986 Jy ! Adding 168 model components to the UV plane model. ! The established model now contains 476 components and 2.22986 Jy ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! Estimated beam: bmin=0.6472 mas, bmaj=0.8715 mas, bpa=10.5 degrees ! Estimated noise=486.593 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00166117 Jy ! Component: 100 - total flux cleaned = 0.00361753 Jy ! Component: 150 - total flux cleaned = 0.00518057 Jy ! Component: 200 - total flux cleaned = 0.00622482 Jy ! Total flux subtracted in 200 components = 0.00622482 Jy ! Clean residual min=-0.003904 max=0.004690 Jy/beam ! Clean residual mean=0.000060 rms=0.000915 Jy/beam ! Combined flux in latest and established models = 2.23609 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 490 components and 2.23609 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160070Jy sigma=0.001255 ! Fit after self-cal, rms=0.159802Jy sigma=0.001252 wmodel J2038+5119_X_2010_06_23_pus_map.mod ! Writing 490 model components to file: J2038+5119_X_2010_06_23_pus_map.mod wobs J2038+5119_X_2010_06_23_pus_uvs.fits ! Writing UV FITS file: J2038+5119_X_2010_06_23_pus_uvs.fits wwins J2038+5119_X_2010_06_23_pus_map.win ! wwins: Wrote 24 windows to J2038+5119_X_2010_06_23_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 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000836455 Jy ! Component: 100 - total flux cleaned = 0.000833241 Jy ! Component: 150 - total flux cleaned = 0.000957184 Jy ! Component: 200 - total flux cleaned = 0.00125867 Jy ! Component: 250 - total flux cleaned = 0.00184071 Jy ! Component: 300 - total flux cleaned = 0.00229318 Jy ! Component: 350 - total flux cleaned = 0.00262406 Jy ! Component: 400 - total flux cleaned = 0.00289322 Jy ! Component: 450 - total flux cleaned = 0.0031578 Jy ! Component: 500 - total flux cleaned = 0.00357153 Jy ! Component: 550 - total flux cleaned = 0.00387646 Jy ! Component: 600 - total flux cleaned = 0.00417702 Jy ! Component: 650 - total flux cleaned = 0.00427593 Jy ! Component: 700 - total flux cleaned = 0.00456895 Jy ! Component: 750 - total flux cleaned = 0.00471334 Jy ! Component: 800 - total flux cleaned = 0.00499952 Jy ! Component: 850 - total flux cleaned = 0.00523509 Jy ! Component: 900 - total flux cleaned = 0.00542165 Jy ! Component: 950 - total flux cleaned = 0.0056987 Jy ! Component: 1000 - total flux cleaned = 0.00592742 Jy ! Component: 1050 - total flux cleaned = 0.00610888 Jy ! Component: 1100 - total flux cleaned = 0.00651201 Jy ! Component: 1150 - total flux cleaned = 0.00664508 Jy ! Component: 1200 - total flux cleaned = 0.00690918 Jy ! Component: 1250 - total flux cleaned = 0.00695319 Jy ! Component: 1300 - total flux cleaned = 0.00725585 Jy ! Component: 1350 - total flux cleaned = 0.00729857 Jy ! Component: 1400 - total flux cleaned = 0.00742662 Jy ! Component: 1450 - total flux cleaned = 0.00759586 Jy ! Component: 1500 - total flux cleaned = 0.00763816 Jy ! Component: 1550 - total flux cleaned = 0.00788811 Jy ! Component: 1600 - total flux cleaned = 0.00797079 Jy ! Component: 1650 - total flux cleaned = 0.00813539 Jy ! Component: 1700 - total flux cleaned = 0.00829844 Jy ! Component: 1750 - total flux cleaned = 0.00842021 Jy ! Component: 1800 - total flux cleaned = 0.0085811 Jy ! Component: 1850 - total flux cleaned = 0.00854092 Jy ! Component: 1900 - total flux cleaned = 0.00870028 Jy ! Component: 1950 - total flux cleaned = 0.00877958 Jy ! Component: 2000 - total flux cleaned = 0.00885801 Jy ! Total flux subtracted in 2048 components = 0.00897541 Jy ! Clean residual min=-0.002372 max=0.002289 Jy/beam ! Clean residual mean=0.000042 rms=0.000718 Jy/beam ! Combined flux in latest and established models = 2.24506 Jy keep ! Adding 177 model components to the UV plane model. ! The established model now contains 666 components and 2.24506 Jy ![Exited script file: muppet] wmap "J2038+5119_X_2010_06_23_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6472 x 0.8715 at 10.5 degrees (North through East) ! Clean map min=-0.0032314 max=1.2938 Jy/beam ! Writing clean map to FITS file: J2038+5119_X_2010_06_23_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.29381 invert ! Inverting map ! Your choice of large map pixels excluded 4.62% of the data. ! The x-axis pixel size should ideally be below 0.1584 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000718867 device "J2038+5119_X_2010_06_23_pus\_map.eps/VPS" ! Attempting to open device: 'J2038+5119_X_2010_06_23_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.166686 0.166686 0.333372 0.666745 1.33349 2.66698 5.33396 10.6679 21.3358 42.6717 85.3433 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0021 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.293 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.6472 x 0.8715 at 10.5 degrees (North through East) ! Clean map min=-0.0032314 max=1.2938 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J2038+5119_X_2010_06_23_pus\_map.gif/vgif" ! Attempting to open device: 'J2038+5119_X_2010_06_23_pus_map.gif/vgif' mapplot cln