! Started logfile: J1740+5211_S_2003_07_09_pus_uvs.log on Mon Apr 8 11:53:35 2013 obs J1740+5211_S_2003_07_09_pus_uva.fits ! Reading UV FITS file: J1740+5211_S_2003_07_09_pus_uva.fits ! AN table 1: 95 integrations on 171 of 171 possible baselines. ! AN table 2: 266 integrations on 171 of 171 possible baselines. ! AN table 3: 21 integrations on 171 of 171 possible baselines. ! AN table 4: 52 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.217196 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1740+5211 ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23697e+09 8e+06 1 8e+06 ! 02 2 2.26697e+09 8e+06 1 8e+06 ! 03 3 2.35697e+09 8e+06 1 8e+06 ! 04 4 2.37697e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1520 lines of history. ! ! Reading 64476 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.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J1740+5211_S_2003_07_09_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.500x0.500 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 1071 telescope corrections were flagged in sub-array 1. ! A total of 2847 telescope corrections were flagged in sub-array 2. ! A total of 393 telescope corrections were flagged in sub-array 3. ! A total of 816 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 1117 telescope corrections were flagged in sub-array 1. ! A total of 2885 telescope corrections were flagged in sub-array 2. ! A total of 396 telescope corrections were flagged in sub-array 3. ! A total of 816 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 1070 telescope corrections were flagged in sub-array 1. ! A total of 2836 telescope corrections were flagged in sub-array 2. ! A total of 393 telescope corrections were flagged in sub-array 3. ! A total of 816 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 1070 telescope corrections were flagged in sub-array 1. ! A total of 2872 telescope corrections were flagged in sub-array 2. ! A total of 396 telescope corrections were flagged in sub-array 3. ! A total of 822 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.096545Jy sigma=0.007710 ! Fit after self-cal, rms=0.357315Jy sigma=0.002354 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.05 mas, bmaj=2.381 mas, bpa=11.61 degrees ! Estimated noise=829.581 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.729937 Jy ! Component: 100 - total flux cleaned = 0.892842 Jy ! Total flux subtracted in 100 components = 0.892842 Jy ! Clean residual min=-0.019386 max=0.057941 Jy/beam ! Clean residual mean=0.000066 rms=0.005868 Jy/beam ! Combined flux in latest and established models = 0.892842 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.892842 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.370968Jy sigma=0.002532 ! Fit after self-cal, rms=0.370965Jy sigma=0.002532 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0695557 Jy ! Component: 100 - total flux cleaned = 0.116429 Jy ! Total flux subtracted in 100 components = 0.116429 Jy ! Clean residual min=-0.017933 max=0.025713 Jy/beam ! Clean residual mean=0.000014 rms=0.004916 Jy/beam ! Combined flux in latest and established models = 1.00927 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 20 components and 1.00927 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.350529Jy sigma=0.002296 ! Fit after self-cal, rms=0.350482Jy sigma=0.002295 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.697 mas, bmaj=3.203 mas, bpa=27.08 degrees ! Estimated noise=402.451 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0505821 Jy ! Component: 100 - total flux cleaned = 0.0732739 Jy ! Component: 150 - total flux cleaned = 0.084887 Jy ! Component: 200 - total flux cleaned = 0.0910249 Jy ! Total flux subtracted in 200 components = 0.0910249 Jy ! Clean residual min=-0.017228 max=0.020385 Jy/beam ! Clean residual mean=0.000280 rms=0.004900 Jy/beam ! Combined flux in latest and established models = 1.1003 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 30 components and 1.1003 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.345674Jy sigma=0.002223 ! Fit after self-cal, rms=0.345619Jy sigma=0.002222 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 72 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.14 GC 1.08 GG 1.00* ! HN 0.94 KP 0.77 LA 1.16 MA 0.94 ! MK 1.02 NL 0.98 NY 1.17 ON 1.29 ! OV 1.04 PT 1.02 SC 0.87 TC 1.00* ! TS 1.08 WF 1.00* WZ 1.09 ! ! A total of 66 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 1.18 GC 1.03 GG 1.00* ! HN 0.98 KP 0.92 LA 1.19 MA 0.93 ! MK 0.93 NL 1.00 NY 1.15 ON 1.31 ! OV 1.37 PT 1.03 SC 0.90 TC 1.00* ! TS 1.10 WF 0.89 WZ 1.13 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.11 FD 1.00* GC 1.09 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 1.04 NL 1.00* NY 1.36 ON 1.00* ! OV 1.16 PT 1.00* SC 1.00* TC 1.00* ! TS 1.04 WF 1.00* WZ 1.00* ! ! ! Correcting IF 2. ! A total of 87 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 1.04 GC 0.93 GG 1.00* ! HN 1.03 KP 0.91 LA 1.08 MA 1.00* ! MK 1.05 NL 1.03 NY 0.89 ON 1.07 ! OV 1.06 PT 1.00 SC 0.90 TC 1.00* ! TS 1.07 WF 1.00* WZ 1.00* ! ! A total of 63 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.08 GC 0.93 GG 1.00* ! HN 1.03 KP 0.95 LA 1.13 MA 1.00* ! MK 1.00 NL 1.04 NY 0.87 ON 1.02 ! OV 1.12 PT 1.01 SC 0.92 TC 1.00* ! TS 1.10 WF 0.94 WZ 1.01 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.22 FD 1.00* GC 0.92 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 1.09 NL 1.00* NY 0.88 ON 1.00* ! OV 1.05 PT 1.00* SC 1.00* TC 1.00* ! TS 1.04 WF 1.00* WZ 1.00* ! ! ! Correcting IF 3. ! A total of 66 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.68 FD 0.99 GC 1.12 GG 1.00* ! HN 0.97 KP 0.92 LA 0.90 MA 0.93 ! MK 1.01 NL 1.05 NY 0.87 ON 1.25 ! OV 1.04 PT 1.02 SC 1.02 TC 1.00* ! TS 1.08 WF 1.00* WZ 1.24 ! ! A total of 78 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.16 FD 1.02 GC 1.14 GG 1.00* ! HN 1.02 KP 0.95 LA 0.91 MA 0.95 ! MK 0.97 NL 1.07 NY 0.85 ON 1.31 ! OV 1.14 PT 1.05 SC 1.06 TC 1.00* ! TS 1.11 WF 0.93 WZ 1.02 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.46 FD 1.00* GC 1.16 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 1.06 NL 1.00* NY 0.88 ON 1.00* ! OV 1.09 PT 1.00* SC 1.00* TC 1.00* ! TS 1.05 WF 1.00* WZ 1.00* ! ! ! Correcting IF 4. ! A total of 66 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.83 FD 0.88 GC 1.17 GG 1.00* ! HN 0.98 KP 0.92 LA 0.72 MA 0.94 ! MK 0.95 NL 0.99 NY 0.81 ON 1.36 ! OV 0.89 PT 0.98 SC 1.02 TC 1.00* ! TS 1.07 WF 1.00* WZ 1.11 ! ! A total of 87 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.80 FD 0.91 GC 1.17 GG 1.00* ! HN 1.00 KP 0.92 LA 0.75 MA 1.49 ! MK 0.94 NL 0.98 NY 0.79 ON 1.14 ! OV 1.10 PT 1.00 SC 1.10 TC 1.00* ! TS 1.13 WF 0.91 WZ 1.08 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 66 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.87 FD 1.00* GC 1.25 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 0.95 NL 1.00* NY 0.87 ON 1.00* ! OV 0.92 PT 1.00* SC 1.00* TC 1.00* ! TS 1.05 WF 1.00* WZ 1.00* ! ! ! Fit before self-cal, rms=0.345619Jy sigma=0.002222 ! Fit after self-cal, rms=0.303077Jy sigma=0.001742 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.734 mas, bmaj=3.221 mas, bpa=27.19 degrees ! Estimated noise=404.256 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00491687 Jy ! Component: 100 - total flux cleaned = -0.00560627 Jy ! Component: 150 - total flux cleaned = -0.00593735 Jy ! Component: 200 - total flux cleaned = -0.00612938 Jy ! Total flux subtracted in 200 components = -0.00612938 Jy ! Clean residual min=-0.010568 max=0.013651 Jy/beam ! Clean residual mean=0.000162 rms=0.003113 Jy/beam ! Combined flux in latest and established models = 1.09417 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 37 components and 1.09417 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.302922Jy sigma=0.001741 ! Fit after self-cal, rms=0.303212Jy sigma=0.001740 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.303212Jy sigma=0.001740 ! Fit after self-cal, rms=0.286622Jy sigma=0.001525 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.744 mas, bmaj=3.226 mas, bpa=26.97 degrees ! Estimated noise=404.3 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00208359 Jy ! Component: 100 - total flux cleaned = -0.0028301 Jy ! Component: 150 - total flux cleaned = -0.00283116 Jy ! Component: 200 - total flux cleaned = -0.00291814 Jy ! Total flux subtracted in 200 components = -0.00291814 Jy ! Clean residual min=-0.003288 max=0.007346 Jy/beam ! Clean residual mean=0.000307 rms=0.001100 Jy/beam ! Combined flux in latest and established models = 1.09125 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 40 components and 1.09125 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286451Jy sigma=0.001524 ! Fit after self-cal, rms=0.286645Jy sigma=0.001524 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (2.5, 3.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00766985 Jy ! Component: 100 - total flux cleaned = 0.0105732 Jy ! Component: 150 - total flux cleaned = 0.0111019 Jy ! Component: 200 - total flux cleaned = 0.0116461 Jy ! Total flux subtracted in 200 components = 0.0116461 Jy ! Clean residual min=-0.002983 max=0.004639 Jy/beam ! Clean residual mean=0.000247 rms=0.000953 Jy/beam ! Combined flux in latest and established models = 1.10289 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 56 components and 1.10289 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286431Jy sigma=0.001521 ! Fit after self-cal, rms=0.286379Jy sigma=0.001521 ! Inverting map ! Added new window around map position (11.5, 20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00447035 Jy ! Component: 100 - total flux cleaned = 0.00630773 Jy ! Component: 150 - total flux cleaned = 0.00681444 Jy ! Component: 200 - total flux cleaned = 0.00736278 Jy ! Total flux subtracted in 200 components = 0.00736278 Jy ! Clean residual min=-0.002950 max=0.004021 Jy/beam ! Clean residual mean=0.000214 rms=0.000902 Jy/beam ! Combined flux in latest and established models = 1.11026 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 73 components and 1.11026 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286276Jy sigma=0.001520 ! Fit after self-cal, rms=0.286269Jy sigma=0.001520 ! Inverting map ! Added new window around map position (-0.5, 11). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0032268 Jy ! Component: 100 - total flux cleaned = 0.00363808 Jy ! Component: 150 - total flux cleaned = 0.00407803 Jy ! Component: 200 - total flux cleaned = 0.00440104 Jy ! Total flux subtracted in 200 components = 0.00440104 Jy ! Clean residual min=-0.002866 max=0.003997 Jy/beam ! Clean residual mean=0.000192 rms=0.000867 Jy/beam ! Combined flux in latest and established models = 1.11466 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 91 components and 1.11466 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286181Jy sigma=0.001519 ! Fit after self-cal, rms=0.286175Jy sigma=0.001519 ! Inverting map ! Added new window around map position (30, 49.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00337397 Jy ! Component: 100 - total flux cleaned = 0.00417833 Jy ! Component: 150 - total flux cleaned = 0.00465471 Jy ! Component: 200 - total flux cleaned = 0.00509397 Jy ! Total flux subtracted in 200 components = 0.00509397 Jy ! Clean residual min=-0.002872 max=0.003604 Jy/beam ! Clean residual mean=0.000174 rms=0.000845 Jy/beam ! Combined flux in latest and established models = 1.11975 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 105 components and 1.11975 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286098Jy sigma=0.001518 ! Fit after self-cal, rms=0.286096Jy sigma=0.001518 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286096Jy sigma=0.001518 ! Fit after self-cal, rms=0.286202Jy sigma=0.001517 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.75 mas, bmaj=3.226 mas, bpa=26.48 degrees ! Estimated noise=404.158 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000666698 Jy ! Component: 100 - total flux cleaned = 4.622e-05 Jy ! Component: 150 - total flux cleaned = 0.000938759 Jy ! Component: 200 - total flux cleaned = 0.00153611 Jy ! Total flux subtracted in 200 components = 0.00153611 Jy ! Clean residual min=-0.002736 max=0.003402 Jy/beam ! Clean residual mean=0.000134 rms=0.000794 Jy/beam ! Combined flux in latest and established models = 1.12129 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 113 components and 1.12129 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286099Jy sigma=0.001516 ! Fit after self-cal, rms=0.286096Jy sigma=0.001516 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0012678 Jy ! Component: 100 - total flux cleaned = 0.00186168 Jy ! Component: 150 - total flux cleaned = 0.00206639 Jy ! Component: 200 - total flux cleaned = 0.00232325 Jy ! Total flux subtracted in 200 components = 0.00232325 Jy ! Clean residual min=-0.002703 max=0.003330 Jy/beam ! Clean residual mean=0.000123 rms=0.000779 Jy/beam ! Combined flux in latest and established models = 1.12361 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 120 components and 1.12361 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286057Jy sigma=0.001516 ! Fit after self-cal, rms=0.286055Jy sigma=0.001516 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 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.286055Jy sigma=0.001516 ! Fit after self-cal, rms=0.286374Jy sigma=0.001515 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.158 mas, bmaj=8.74 mas, bpa=-57.22 degrees ! Estimated noise=559.024 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00273418 Jy ! Component: 100 - total flux cleaned = 0.00314004 Jy ! Component: 150 - total flux cleaned = 0.00352836 Jy ! Component: 200 - total flux cleaned = 0.00395432 Jy ! Total flux subtracted in 200 components = 0.00395432 Jy ! Clean residual min=-0.003581 max=0.004218 Jy/beam ! Clean residual mean=0.000276 rms=0.001184 Jy/beam ! Combined flux in latest and established models = 1.12757 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 139 components and 1.12757 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286363Jy sigma=0.001515 ! Fit after self-cal, rms=0.286364Jy sigma=0.001515 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000730521 Jy ! Component: 100 - total flux cleaned = 0.000928892 Jy ! Component: 150 - total flux cleaned = 0.00110391 Jy ! Component: 200 - total flux cleaned = 0.00126397 Jy ! Total flux subtracted in 200 components = 0.00126397 Jy ! Clean residual min=-0.003631 max=0.003990 Jy/beam ! Clean residual mean=0.000258 rms=0.001157 Jy/beam ! Combined flux in latest and established models = 1.12883 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 150 components and 1.12883 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.286355Jy sigma=0.001515 ! Fit after self-cal, rms=0.286353Jy sigma=0.001515 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=2.755 mas, bmaj=3.226 mas, bpa=26.19 degrees ! Estimated noise=404.039 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.755 x 3.226 at 26.19 degrees (North through East) ! Clean map min=-0.0025852 max=0.98928 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.286353Jy sigma=0.001515 ! Fit after self-cal, rms=0.272379Jy sigma=0.001301 ! Inverting map and beam ! Estimated beam: bmin=2.756 mas, bmaj=3.224 mas, bpa=25.91 degrees ! Estimated noise=403.538 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00193677 Jy ! Component: 100 - total flux cleaned = -0.00214219 Jy ! Component: 150 - total flux cleaned = -0.00221916 Jy ! Component: 200 - total flux cleaned = -0.00215072 Jy ! Total flux subtracted in 200 components = -0.00215072 Jy ! Clean residual min=-0.002585 max=0.002313 Jy/beam ! Clean residual mean=0.000049 rms=0.000605 Jy/beam ! Combined flux in latest and established models = 1.12668 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 161 components and 1.12668 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272299Jy sigma=0.001301 ! Fit after self-cal, rms=0.272295Jy sigma=0.001301 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000175473 Jy ! Component: 100 - total flux cleaned = -1.20623e-05 Jy ! Component: 150 - total flux cleaned = -7.02713e-05 Jy ! Component: 200 - total flux cleaned = -0.000181419 Jy ! Total flux subtracted in 200 components = -0.000181419 Jy ! Clean residual min=-0.002587 max=0.002253 Jy/beam ! Clean residual mean=0.000049 rms=0.000595 Jy/beam ! Combined flux in latest and established models = 1.1265 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 171 components and 1.1265 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272264Jy sigma=0.001301 ! Fit after self-cal, rms=0.272262Jy sigma=0.001301 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272262Jy sigma=0.001301 ! Fit after self-cal, rms=0.272604Jy sigma=0.001300 print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.091 mas, bmaj=2.406 mas, bpa=13.1 degrees ! Estimated noise=840.311 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.745754 Jy ! Component: 100 - total flux cleaned = 0.908557 Jy ! Component: 150 - total flux cleaned = 0.961882 Jy ! Component: 200 - total flux cleaned = 0.996261 Jy ! Total flux subtracted in 200 components = 0.996261 Jy ! Clean residual min=-0.006384 max=0.018920 Jy/beam ! Clean residual mean=0.000060 rms=0.001789 Jy/beam ! Combined flux in latest and established models = 0.996261 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 26 model components to the UV plane model. ! The established model now contains 26 components and 0.996261 Jy ! Inverting map and beam ! Estimated beam: bmin=2.759 mas, bmaj=3.224 mas, bpa=25.63 degrees ! Estimated noise=403.31 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.0437299 Jy ! Component: 100 - total flux cleaned = 0.0671836 Jy ! Component: 150 - total flux cleaned = 0.0814587 Jy ! Component: 200 - total flux cleaned = 0.0911721 Jy ! Total flux subtracted in 200 components = 0.0911721 Jy ! Clean residual min=-0.002988 max=0.005607 Jy/beam ! Clean residual mean=0.000224 rms=0.000988 Jy/beam ! Combined flux in latest and established models = 1.08743 Jy ! Adding 63 model components to the UV plane model. ! The established model now contains 77 components and 1.08743 Jy ! Inverting map ! Added new window around map position (4.5, 13.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00725276 Jy ! Component: 100 - total flux cleaned = 0.0128958 Jy ! Component: 150 - total flux cleaned = 0.0173241 Jy ! Component: 200 - total flux cleaned = 0.0209599 Jy ! Total flux subtracted in 200 components = 0.0209599 Jy ! Clean residual min=-0.002693 max=0.002784 Jy/beam ! Clean residual mean=0.000120 rms=0.000673 Jy/beam ! Combined flux in latest and established models = 1.10839 Jy ! Adding 40 model components to the UV plane model. ! The established model now contains 104 components and 1.10839 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.147 mas, bmaj=8.73 mas, bpa=-56.5 degrees ! Estimated noise=557.973 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00597596 Jy ! Component: 100 - total flux cleaned = 0.00938533 Jy ! Component: 150 - total flux cleaned = 0.0118168 Jy ! Component: 200 - total flux cleaned = 0.0136191 Jy ! Total flux subtracted in 200 components = 0.0136191 Jy ! Clean residual min=-0.003009 max=0.003071 Jy/beam ! Clean residual mean=0.000170 rms=0.000817 Jy/beam ! Combined flux in latest and established models = 1.12201 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.00135827 Jy ! Component: 100 - total flux cleaned = 0.00238606 Jy ! Component: 150 - total flux cleaned = 0.00317814 Jy ! Component: 200 - total flux cleaned = 0.00366138 Jy ! Total flux subtracted in 200 components = 0.00366138 Jy ! Clean residual min=-0.002883 max=0.002932 Jy/beam ! Clean residual mean=0.000119 rms=0.000784 Jy/beam ! Combined flux in latest and established models = 1.12567 Jy ! Adding 86 model components to the UV plane model. ! The established model now contains 173 components and 1.12567 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.759 mas, bmaj=3.224 mas, bpa=25.63 degrees ! Estimated noise=403.31 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000267529 Jy ! Component: 100 - total flux cleaned = -0.000327576 Jy ! Component: 150 - total flux cleaned = -0.000212352 Jy ! Component: 200 - total flux cleaned = -4.99412e-05 Jy ! Total flux subtracted in 200 components = -4.99412e-05 Jy ! Clean residual min=-0.002695 max=0.002285 Jy/beam ! Clean residual mean=0.000038 rms=0.000575 Jy/beam ! Combined flux in latest and established models = 1.12562 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 191 components and 1.12562 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272552Jy sigma=0.001300 ! Fit after self-cal, rms=0.272539Jy sigma=0.001300 wmodel J1740+5211_S_2003_07_09_pus_map.mod ! Writing 191 model components to file: J1740+5211_S_2003_07_09_pus_map.mod wobs J1740+5211_S_2003_07_09_pus_uvs.fits ! Writing UV FITS file: J1740+5211_S_2003_07_09_pus_uvs.fits wwins J1740+5211_S_2003_07_09_pus_map.win ! wwins: Wrote 6 windows to J1740+5211_S_2003_07_09_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000672019 Jy ! Component: 100 - total flux cleaned = -0.000471107 Jy ! Component: 150 - total flux cleaned = -0.00039498 Jy ! Component: 200 - total flux cleaned = -0.000321441 Jy ! Component: 250 - total flux cleaned = -0.000214067 Jy ! Component: 300 - total flux cleaned = -0.000179053 Jy ! Component: 350 - total flux cleaned = -0.00017968 Jy ! Component: 400 - total flux cleaned = -0.000179675 Jy ! Component: 450 - total flux cleaned = -0.000113586 Jy ! Component: 500 - total flux cleaned = -1.63251e-05 Jy ! Component: 550 - total flux cleaned = -1.62256e-05 Jy ! Component: 600 - total flux cleaned = 0.000109849 Jy ! Component: 650 - total flux cleaned = 7.85437e-05 Jy ! Component: 700 - total flux cleaned = 7.87101e-05 Jy ! Component: 750 - total flux cleaned = 0.000139587 Jy ! Component: 800 - total flux cleaned = 0.000139515 Jy ! Component: 850 - total flux cleaned = 0.000169229 Jy ! Component: 900 - total flux cleaned = 0.000227892 Jy ! Component: 950 - total flux cleaned = 0.000198624 Jy ! Component: 1000 - total flux cleaned = 0.000169875 Jy ! Component: 1050 - total flux cleaned = 0.000227223 Jy ! Component: 1100 - total flux cleaned = 0.000198773 Jy ! Component: 1150 - total flux cleaned = 0.000227214 Jy ! Component: 1200 - total flux cleaned = 0.000255524 Jy ! Component: 1250 - total flux cleaned = 0.000367116 Jy ! Component: 1300 - total flux cleaned = 0.000284011 Jy ! Component: 1350 - total flux cleaned = 0.00039423 Jy ! Component: 1400 - total flux cleaned = 0.000366971 Jy ! Component: 1450 - total flux cleaned = 0.000421323 Jy ! Component: 1500 - total flux cleaned = 0.000448085 Jy ! Component: 1550 - total flux cleaned = 0.000475085 Jy ! Component: 1600 - total flux cleaned = 0.000501357 Jy ! Component: 1650 - total flux cleaned = 0.000474659 Jy ! Component: 1700 - total flux cleaned = 0.000579965 Jy ! Component: 1750 - total flux cleaned = 0.000448882 Jy ! Component: 1800 - total flux cleaned = 0.000527114 Jy ! Component: 1850 - total flux cleaned = 0.000501161 Jy ! Component: 1900 - total flux cleaned = 0.000501319 Jy ! Component: 1950 - total flux cleaned = 0.000527017 Jy ! Component: 2000 - total flux cleaned = 0.000501575 Jy ! Total flux subtracted in 2048 components = 0.000552391 Jy ! Clean residual min=-0.001475 max=0.001415 Jy/beam ! Clean residual mean=0.000035 rms=0.000483 Jy/beam ! Combined flux in latest and established models = 1.12618 Jy keep ! Adding 152 model components to the UV plane model. ! The established model now contains 343 components and 1.12618 Jy ![Exited script file: muppet] wmap "J1740+5211_S_2003_07_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.759 x 3.224 at 25.63 degrees (North through East) ! Clean map min=-0.002387 max=0.98728 Jy/beam ! Writing clean map to FITS file: J1740+5211_S_2003_07_09_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.987276 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000483318 device "J1740+5211_S_2003_07_09_pus\_map.eps/VPS" ! Attempting to open device: 'J1740+5211_S_2003_07_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.146864 0.146864 0.293728 0.587457 1.17491 2.34983 4.69965 9.39931 18.7986 37.5972 75.1945 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0014 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.987 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.759 x 3.224 at 25.63 degrees (North through East) ! Clean map min=-0.002387 max=0.98728 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J1740+5211_S_2003_07_09_pus_uvs.log closed on Mon Apr 8 11:53:55 2013