! Started logfile: J1408-0752_X_2001_03_12_pus_uvs.log on Wed Jan 2 19:58:39 2013 obs J1408-0752_X_2001_03_12_pus_uva.fits ! Reading UV FITS file: J1408-0752_X_2001_03_12_pus_uva.fits ! AN table 1: 51 integrations on 136 of 136 possible baselines. ! AN table 2: 18 integrations on 136 of 136 possible baselines. ! AN table 3: 11 integrations on 136 of 136 possible baselines. ! AN table 4: 28 integrations on 136 of 136 possible baselines. ! Apparent sampling: 0.203772 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1408-0752 ! ! 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 1749 lines of history. ! ! Reading 11972 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 J1408-0752_X_2001_03_12_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 459 telescope corrections were flagged in sub-array 1. ! A total of 169 telescope corrections were flagged in sub-array 2. ! A total of 187 telescope corrections were flagged in sub-array 3. ! A total of 252 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 459 telescope corrections were flagged in sub-array 1. ! A total of 169 telescope corrections were flagged in sub-array 2. ! A total of 187 telescope corrections were flagged in sub-array 3. ! A total of 252 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 459 telescope corrections were flagged in sub-array 1. ! A total of 169 telescope corrections were flagged in sub-array 2. ! A total of 187 telescope corrections were flagged in sub-array 3. ! A total of 252 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 459 telescope corrections were flagged in sub-array 1. ! A total of 169 telescope corrections were flagged in sub-array 2. ! A total of 187 telescope corrections were flagged in sub-array 3. ! A total of 252 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.736624Jy sigma=0.006423 ! Fit after self-cal, rms=0.243973Jy sigma=0.002148 ! 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 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.5807 mas, bmaj=2.101 mas, bpa=-2.424 degrees ! Estimated noise=1056.58 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.881783 Jy ! Component: 100 - total flux cleaned = 1.07792 Jy ! Total flux subtracted in 100 components = 1.07792 Jy ! Clean residual min=-0.017645 max=0.066316 Jy/beam ! Clean residual mean=-0.000020 rms=0.005562 Jy/beam ! Combined flux in latest and established models = 1.07792 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 1.07792 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.200175Jy sigma=0.001720 ! Fit after self-cal, rms=0.200167Jy sigma=0.001720 ! Inverting map ! Your choice of large map pixels excluded 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0701088 Jy ! Component: 100 - total flux cleaned = 0.104189 Jy ! Total flux subtracted in 100 components = 0.104189 Jy ! Clean residual min=-0.012219 max=0.018415 Jy/beam ! Clean residual mean=-0.000018 rms=0.003289 Jy/beam ! Combined flux in latest and established models = 1.18211 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 9 components and 1.18211 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167192Jy sigma=0.001388 ! Fit after self-cal, rms=0.167074Jy sigma=0.001387 ! Inverting map ! Your choice of large map pixels excluded 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 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 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7907 mas, bmaj=2.573 mas, bpa=-5.093 degrees ! Estimated noise=736.663 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 (-1.2, -1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0288253 Jy ! Component: 100 - total flux cleaned = 0.0433205 Jy ! Component: 150 - total flux cleaned = 0.0496798 Jy ! Component: 200 - total flux cleaned = 0.0513729 Jy ! Total flux subtracted in 200 components = 0.0513729 Jy ! Clean residual min=-0.011835 max=0.020367 Jy/beam ! Clean residual mean=-0.000050 rms=0.003329 Jy/beam ! Combined flux in latest and established models = 1.23348 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 31 components and 1.23348 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162068Jy sigma=0.001340 ! Fit after self-cal, rms=0.161765Jy sigma=0.001337 ! Inverting map ! Your choice of large map pixels excluded 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Added new window around map position (-6, -2.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0198345 Jy ! Component: 100 - total flux cleaned = 0.0299852 Jy ! Component: 150 - total flux cleaned = 0.0336255 Jy ! Component: 200 - total flux cleaned = 0.035175 Jy ! Total flux subtracted in 200 components = 0.035175 Jy ! Clean residual min=-0.009289 max=0.012305 Jy/beam ! Clean residual mean=-0.000070 rms=0.002634 Jy/beam ! Combined flux in latest and established models = 1.26866 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 67 components and 1.26866 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.159334Jy sigma=0.001312 ! Fit after self-cal, rms=0.159199Jy sigma=0.001311 ! Inverting map ! Your choice of large map pixels excluded 1% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.99 HH 1.00* HN 1.17 ! KK 1.00 KP 1.03 LA 0.97 MA 1.00* ! MC 1.00* MK 0.95 NL 1.00 ON 1.00* ! OV 0.98 PT 0.99 SC 1.00 TS 1.03 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 1.05 HH 1.00* HN 1.00* ! KK 0.94 KP 1.03 LA 1.03 MA 1.00* ! MC 1.00* MK 0.97 NL 1.00* ON 1.00* ! OV 1.02 PT 1.00 SC 1.00* TS 0.99 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* HH 0.00* HN 0.00* ! KK 0.00* KP 0.00* LA 0.00* MA 0.00* ! MC 0.00* MK 0.00* NL 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TS 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.99 FD 1.03 HH 1.00* HN 1.00 ! KK 1.00* KP 1.02 LA 1.00 MA 1.00* ! MC 1.00* MK 1.00* NL 1.01 ON 1.00* ! OV 0.99 PT 0.98 SC 0.99 TS 1.00* ! WZ 1.00* ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.97 HH 1.00* HN 1.11 ! KK 1.00 KP 0.97 LA 0.98 MA 1.00* ! MC 1.00* MK 0.96 NL 0.97 ON 1.00* ! OV 0.97 PT 1.03 SC 1.00 TS 1.02 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.04 HH 1.00* HN 1.00* ! KK 0.96 KP 0.99 LA 1.05 MA 1.00* ! MC 1.00* MK 0.95 NL 1.00* ON 1.00* ! OV 1.01 PT 1.02 SC 1.00* TS 0.98 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* HH 0.00* HN 0.00* ! KK 0.00* KP 0.00* LA 0.00* MA 0.00* ! MC 0.00* MK 0.00* NL 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TS 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.01 FD 1.02 HH 1.00* HN 0.99 ! KK 1.00* KP 0.96 LA 1.00 MA 1.00* ! MC 1.00* MK 1.00* NL 0.97 ON 1.00* ! OV 1.00 PT 1.00 SC 0.98 TS 1.00* ! WZ 1.00* ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 0.96 HH 1.00* HN 1.13 ! KK 1.09 KP 0.95 LA 0.99 MA 1.00* ! MC 1.00* MK 0.97 NL 1.00 ON 1.00* ! OV 1.01 PT 1.04 SC 1.01 TS 1.01 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 1.04 HH 1.00* HN 1.00* ! KK 1.03 KP 0.96 LA 1.02 MA 1.00* ! MC 1.00* MK 0.96 NL 1.00* ON 1.00* ! OV 1.03 PT 1.06 SC 1.00* TS 0.95 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* HH 0.00* HN 0.00* ! KK 0.00* KP 0.00* LA 0.00* MA 0.00* ! MC 0.00* MK 0.00* NL 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TS 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.04 FD 0.99 HH 1.00* HN 1.00 ! KK 1.00* KP 0.95 LA 1.02 MA 1.00* ! MC 1.00* MK 1.00* NL 0.99 ON 1.00* ! OV 1.02 PT 1.03 SC 1.00 TS 1.00* ! WZ 1.00* ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.95 HH 1.00* HN 1.15 ! KK 1.01 KP 0.96 LA 0.96 MA 1.00* ! MC 1.00* MK 1.00 NL 1.01 ON 1.00* ! OV 0.99 PT 1.04 SC 1.08 TS 0.98 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 1.01 HH 1.00* HN 1.00* ! KK 0.96 KP 0.99 LA 1.01 MA 1.00* ! MC 1.00* MK 1.01 NL 1.00* ON 1.00* ! OV 1.03 PT 1.06 SC 1.00* TS 0.96 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* HH 0.00* HN 0.00* ! KK 0.00* KP 0.00* LA 0.00* MA 0.00* ! MC 0.00* MK 0.00* NL 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TS 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 FD 0.99 HH 1.00* HN 1.00 ! KK 1.00* KP 0.96 LA 0.99 MA 1.00* ! MC 1.00* MK 1.00* NL 1.00 ON 1.00* ! OV 1.00 PT 1.04 SC 1.05 TS 1.00* ! WZ 1.00* ! ! ! Fit before self-cal, rms=0.159199Jy sigma=0.001311 ! Fit after self-cal, rms=0.149800Jy sigma=0.001221 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7929 mas, bmaj=2.576 mas, bpa=-5.077 degrees ! Estimated noise=735.63 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.4, -2.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00842908 Jy ! Component: 100 - total flux cleaned = 0.0123126 Jy ! Component: 150 - total flux cleaned = 0.014549 Jy ! Component: 200 - total flux cleaned = 0.0155697 Jy ! Total flux subtracted in 200 components = 0.0155697 Jy ! Clean residual min=-0.004421 max=0.005471 Jy/beam ! Clean residual mean=0.000046 rms=0.001409 Jy/beam ! Combined flux in latest and established models = 1.28423 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 96 components and 1.28423 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.149275Jy sigma=0.001215 ! Fit after self-cal, rms=0.149300Jy sigma=0.001215 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.149300Jy sigma=0.001215 ! Fit after self-cal, rms=0.147672Jy sigma=0.001203 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7932 mas, bmaj=2.577 mas, bpa=-5.033 degrees ! Estimated noise=735.002 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00265681 Jy ! Component: 100 - total flux cleaned = 0.0043164 Jy ! Component: 150 - total flux cleaned = 0.00511721 Jy ! Component: 200 - total flux cleaned = 0.00592069 Jy ! Total flux subtracted in 200 components = 0.00592069 Jy ! Clean residual min=-0.003441 max=0.004507 Jy/beam ! Clean residual mean=0.000046 rms=0.001032 Jy/beam ! Combined flux in latest and established models = 1.29015 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 125 components and 1.29015 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147499Jy sigma=0.001201 ! Fit after self-cal, rms=0.147497Jy sigma=0.001201 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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, -15.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00366114 Jy ! Component: 100 - total flux cleaned = 0.00479054 Jy ! Component: 150 - total flux cleaned = 0.00541946 Jy ! Component: 200 - total flux cleaned = 0.00571106 Jy ! Total flux subtracted in 200 components = 0.00571106 Jy ! Clean residual min=-0.003576 max=0.004197 Jy/beam ! Clean residual mean=0.000045 rms=0.000994 Jy/beam ! Combined flux in latest and established models = 1.29586 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 139 components and 1.29586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147359Jy sigma=0.001200 ! Fit after self-cal, rms=0.147361Jy sigma=0.001200 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Added new window around map position (-3, -0.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00334533 Jy ! Component: 100 - total flux cleaned = 0.0035996 Jy ! Component: 150 - total flux cleaned = 0.0037422 Jy ! Component: 200 - total flux cleaned = 0.00400411 Jy ! Total flux subtracted in 200 components = 0.00400411 Jy ! Clean residual min=-0.003241 max=0.003983 Jy/beam ! Clean residual mean=0.000045 rms=0.000951 Jy/beam ! Combined flux in latest and established models = 1.29986 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 158 components and 1.29986 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147246Jy sigma=0.001199 ! Fit after self-cal, rms=0.147241Jy sigma=0.001199 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.147241Jy sigma=0.001199 ! Fit after self-cal, rms=0.147225Jy sigma=0.001199 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.793 mas, bmaj=2.577 mas, bpa=-5.026 degrees ! Estimated noise=734.77 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000904807 Jy ! Component: 100 - total flux cleaned = 0.00131612 Jy ! Component: 150 - total flux cleaned = 0.00156911 Jy ! Component: 200 - total flux cleaned = 0.00174744 Jy ! Total flux subtracted in 200 components = 0.00174744 Jy ! Clean residual min=-0.003322 max=0.003940 Jy/beam ! Clean residual mean=0.000046 rms=0.000924 Jy/beam ! Combined flux in latest and established models = 1.30161 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 174 components and 1.30161 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147166Jy sigma=0.001198 ! Fit after self-cal, rms=0.147168Jy sigma=0.001198 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.000200571 Jy ! Component: 100 - total flux cleaned = 0.000314036 Jy ! Component: 150 - total flux cleaned = 0.000366177 Jy ! Component: 200 - total flux cleaned = 0.000468072 Jy ! Total flux subtracted in 200 components = 0.000468072 Jy ! Clean residual min=-0.003249 max=0.003831 Jy/beam ! Clean residual mean=0.000045 rms=0.000911 Jy/beam ! Combined flux in latest and established models = 1.30208 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 183 components and 1.30208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147131Jy sigma=0.001198 ! Fit after self-cal, rms=0.147131Jy sigma=0.001198 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.147131Jy sigma=0.001198 ! Fit after self-cal, rms=0.147144Jy sigma=0.001198 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=2.337 mas, bmaj=3.807 mas, bpa=-4.737 degrees ! Estimated noise=912.364 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00234244 Jy ! Component: 100 - total flux cleaned = 0.00351182 Jy ! Component: 150 - total flux cleaned = 0.00367209 Jy ! Component: 200 - total flux cleaned = 0.00372109 Jy ! Total flux subtracted in 200 components = 0.00372109 Jy ! Clean residual min=-0.003341 max=0.005481 Jy/beam ! Clean residual mean=0.000083 rms=0.001241 Jy/beam ! Combined flux in latest and established models = 1.3058 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 200 components and 1.3058 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147127Jy sigma=0.001198 ! Fit after self-cal, rms=0.147128Jy sigma=0.001198 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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 (-12.2, -5.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00497572 Jy ! Component: 100 - total flux cleaned = 0.00763052 Jy ! Component: 150 - total flux cleaned = 0.00901105 Jy ! Component: 200 - total flux cleaned = 0.00941941 Jy ! Total flux subtracted in 200 components = 0.00941941 Jy ! Clean residual min=-0.002853 max=0.003548 Jy/beam ! Clean residual mean=0.000055 rms=0.001014 Jy/beam ! Combined flux in latest and established models = 1.31522 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 255 components and 1.31522 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147032Jy sigma=0.001197 ! Fit after self-cal, rms=0.147027Jy sigma=0.001197 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.793 mas, bmaj=2.577 mas, bpa=-5.022 degrees ! Estimated noise=734.708 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.793 x 2.577 at -5.022 degrees (North through East) ! Clean map min=-0.0027515 max=1.1532 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.147027Jy sigma=0.001197 ! Fit after self-cal, rms=0.136878Jy sigma=0.001100 ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7934 mas, bmaj=2.576 mas, bpa=-5.021 degrees ! Estimated noise=734.432 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00169056 Jy ! Component: 100 - total flux cleaned = 0.00234393 Jy ! Component: 150 - total flux cleaned = 0.00240756 Jy ! Component: 200 - total flux cleaned = 0.00265075 Jy ! Total flux subtracted in 200 components = 0.00265075 Jy ! Clean residual min=-0.002932 max=0.002963 Jy/beam ! Clean residual mean=0.000020 rms=0.000754 Jy/beam ! Combined flux in latest and established models = 1.31787 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 287 components and 1.31787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136807Jy sigma=0.001099 ! Fit after self-cal, rms=0.136809Jy sigma=0.001099 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000450808 Jy ! Component: 100 - total flux cleaned = 0.000624522 Jy ! Component: 150 - total flux cleaned = 0.000735611 Jy ! Component: 200 - total flux cleaned = 0.000894724 Jy ! Total flux subtracted in 200 components = 0.000894724 Jy ! Clean residual min=-0.002831 max=0.002911 Jy/beam ! Clean residual mean=0.000019 rms=0.000741 Jy/beam ! Combined flux in latest and established models = 1.31876 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 300 components and 1.31876 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136765Jy sigma=0.001099 ! Fit after self-cal, rms=0.136766Jy sigma=0.001099 ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.136766Jy sigma=0.001099 ! Fit after self-cal, rms=0.136767Jy sigma=0.001099 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.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.5811 mas, bmaj=2.105 mas, bpa=-2.378 degrees ! Estimated noise=1053.39 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.878012 Jy ! Component: 100 - total flux cleaned = 1.07342 Jy ! Component: 150 - total flux cleaned = 1.14404 Jy ! Component: 200 - total flux cleaned = 1.17735 Jy ! Total flux subtracted in 200 components = 1.17735 Jy ! Clean residual min=-0.006850 max=0.015213 Jy/beam ! Clean residual mean=0.000007 rms=0.002061 Jy/beam ! Combined flux in latest and established models = 1.17735 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 8 model components to the UV plane model. ! The established model now contains 8 components and 1.17735 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7934 mas, bmaj=2.577 mas, bpa=-5.02 degrees ! Estimated noise=734.342 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.0295697 Jy ! Component: 100 - total flux cleaned = 0.0495424 Jy ! Component: 150 - total flux cleaned = 0.0643134 Jy ! Component: 200 - total flux cleaned = 0.0755677 Jy ! Total flux subtracted in 200 components = 0.0755677 Jy ! Clean residual min=-0.005662 max=0.007589 Jy/beam ! Clean residual mean=0.000060 rms=0.001558 Jy/beam ! Combined flux in latest and established models = 1.25291 Jy ! Adding 55 model components to the UV plane model. ! The established model now contains 58 components and 1.25291 Jy ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Added new window around map position (-4.4, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00900261 Jy ! Component: 100 - total flux cleaned = 0.0162262 Jy ! Component: 150 - total flux cleaned = 0.0222882 Jy ! Component: 200 - total flux cleaned = 0.0275583 Jy ! Total flux subtracted in 200 components = 0.0275583 Jy ! Clean residual min=-0.003871 max=0.003806 Jy/beam ! Clean residual mean=0.000042 rms=0.001063 Jy/beam ! Combined flux in latest and established models = 1.28047 Jy ! Adding 70 model components to the UV plane model. ! The established model now contains 117 components and 1.28047 Jy ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 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.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=2.338 mas, bmaj=3.807 mas, bpa=-4.724 degrees ! Estimated noise=911.569 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00776003 Jy ! Component: 100 - total flux cleaned = 0.0135224 Jy ! Component: 150 - total flux cleaned = 0.0180113 Jy ! Component: 200 - total flux cleaned = 0.02172 Jy ! Total flux subtracted in 200 components = 0.02172 Jy ! Clean residual min=-0.003961 max=0.003855 Jy/beam ! Clean residual mean=0.000055 rms=0.001088 Jy/beam ! Combined flux in latest and established models = 1.30219 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.00313127 Jy ! Component: 100 - total flux cleaned = 0.00579047 Jy ! Component: 150 - total flux cleaned = 0.00806199 Jy ! Component: 200 - total flux cleaned = 0.0100174 Jy ! Total flux subtracted in 200 components = 0.0100174 Jy ! Clean residual min=-0.003244 max=0.003543 Jy/beam ! Clean residual mean=0.000042 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 1.31221 Jy ! Adding 140 model components to the UV plane model. ! The established model now contains 236 components and 1.31221 Jy ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! Estimated beam: bmin=0.7934 mas, bmaj=2.577 mas, bpa=-5.02 degrees ! Estimated noise=734.342 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00208861 Jy ! Component: 100 - total flux cleaned = 0.00299311 Jy ! Component: 150 - total flux cleaned = 0.003507 Jy ! Component: 200 - total flux cleaned = 0.00399121 Jy ! Total flux subtracted in 200 components = 0.00399121 Jy ! Clean residual min=-0.002909 max=0.002962 Jy/beam ! Clean residual mean=0.000019 rms=0.000748 Jy/beam ! Combined flux in latest and established models = 1.3162 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 266 components and 1.3162 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136808Jy sigma=0.001099 ! Fit after self-cal, rms=0.136763Jy sigma=0.001099 wmodel J1408-0752_X_2001_03_12_pus_map.mod ! Writing 266 model components to file: J1408-0752_X_2001_03_12_pus_map.mod wobs J1408-0752_X_2001_03_12_pus_uvs.fits ! Writing UV FITS file: J1408-0752_X_2001_03_12_pus_uvs.fits wwins J1408-0752_X_2001_03_12_pus_map.win ! wwins: Wrote 8 windows to J1408-0752_X_2001_03_12_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.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000328856 Jy ! Component: 100 - total flux cleaned = 0.000379572 Jy ! Component: 150 - total flux cleaned = 0.00033283 Jy ! Component: 200 - total flux cleaned = 0.000422712 Jy ! Component: 250 - total flux cleaned = 0.000509279 Jy ! Component: 300 - total flux cleaned = 0.000593977 Jy ! Component: 350 - total flux cleaned = 0.000677745 Jy ! Component: 400 - total flux cleaned = 0.000759624 Jy ! Component: 450 - total flux cleaned = 0.00083974 Jy ! Component: 500 - total flux cleaned = 0.000997335 Jy ! Component: 550 - total flux cleaned = 0.00103582 Jy ! Component: 600 - total flux cleaned = 0.00118879 Jy ! Component: 650 - total flux cleaned = 0.00126368 Jy ! Component: 700 - total flux cleaned = 0.00133751 Jy ! Component: 750 - total flux cleaned = 0.00148277 Jy ! Component: 800 - total flux cleaned = 0.00155403 Jy ! Component: 850 - total flux cleaned = 0.00158913 Jy ! Component: 900 - total flux cleaned = 0.00162398 Jy ! Component: 950 - total flux cleaned = 0.00158995 Jy ! Component: 1000 - total flux cleaned = 0.00162403 Jy ! Component: 1050 - total flux cleaned = 0.00155644 Jy ! Component: 1100 - total flux cleaned = 0.00159011 Jy ! Component: 1150 - total flux cleaned = 0.00152445 Jy ! Component: 1200 - total flux cleaned = 0.00145948 Jy ! Component: 1250 - total flux cleaned = 0.00155657 Jy ! Component: 1300 - total flux cleaned = 0.00139614 Jy ! Component: 1350 - total flux cleaned = 0.00133295 Jy ! Component: 1400 - total flux cleaned = 0.00139581 Jy ! Component: 1450 - total flux cleaned = 0.00124002 Jy ! Component: 1500 - total flux cleaned = 0.0012709 Jy ! Component: 1550 - total flux cleaned = 0.00117878 Jy ! Component: 1600 - total flux cleaned = 0.00105695 Jy ! Component: 1650 - total flux cleaned = 0.000905995 Jy ! Component: 1700 - total flux cleaned = 0.000846526 Jy ! Component: 1750 - total flux cleaned = 0.000727731 Jy ! Component: 1800 - total flux cleaned = 0.000698241 Jy ! Component: 1850 - total flux cleaned = 0.000580935 Jy ! Component: 1900 - total flux cleaned = 0.000493492 Jy ! Component: 1950 - total flux cleaned = 0.000406786 Jy ! Component: 2000 - total flux cleaned = 0.00040637 Jy ! Total flux subtracted in 2048 components = 0.00032071 Jy ! Clean residual min=-0.001671 max=0.001708 Jy/beam ! Clean residual mean=0.000016 rms=0.000551 Jy/beam ! Combined flux in latest and established models = 1.31652 Jy keep ! Adding 148 model components to the UV plane model. ! The established model now contains 414 components and 1.31652 Jy ![Exited script file: muppet] wmap "J1408-0752_X_2001_03_12_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7934 x 2.577 at -5.02 degrees (North through East) ! Clean map min=-0.0024448 max=1.1538 Jy/beam ! Writing clean map to FITS file: J1408-0752_X_2001_03_12_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.1538 invert ! Inverting map ! Your choice of large map pixels excluded 0.993% of the data. ! The x-axis pixel size should ideally be below 0.1912 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000549599 device "J1408-0752_X_2001_03_12_pus\_map.eps/VPS" ! Attempting to open device: 'J1408-0752_X_2001_03_12_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.142901 0.142901 0.285802 0.571603 1.14321 2.28641 4.57283 9.14566 18.2913 36.5826 73.1652 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0016 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.153 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.7934 x 2.577 at -5.02 degrees (North through East) ! Clean map min=-0.0024448 max=1.1538 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J1408-0752_X_2001_03_12_pus_uvs.log closed on Wed Jan 2 19:58:45 2013