! Started logfile: J2225-0457_S_2003_07_09_pus_uvs.log on Mon Apr 8 12:58:52 2013 obs J2225-0457_S_2003_07_09_pus_uva.fits ! Reading UV FITS file: J2225-0457_S_2003_07_09_pus_uva.fits ! AN table 1: 295 integrations on 171 of 171 possible baselines. ! AN table 2: 136 integrations on 171 of 171 possible baselines. ! AN table 3: 383 integrations on 171 of 171 possible baselines. ! AN table 4: 47 integrations on 171 of 171 possible baselines. ! AN table 5: 151 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.05993 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2225-0457 ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 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 41484 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 = 2048 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J2225-0457_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 = 2048x2048 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 4582 telescope corrections were flagged in sub-array 1. ! A total of 2050 telescope corrections were flagged in sub-array 2. ! A total of 6098 telescope corrections were flagged in sub-array 3. ! A total of 802 telescope corrections were flagged in sub-array 4. ! A total of 2695 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 4662 telescope corrections were flagged in sub-array 1. ! A total of 2050 telescope corrections were flagged in sub-array 2. ! A total of 6188 telescope corrections were flagged in sub-array 3. ! A total of 802 telescope corrections were flagged in sub-array 4. ! A total of 2812 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 4582 telescope corrections were flagged in sub-array 1. ! A total of 2051 telescope corrections were flagged in sub-array 2. ! A total of 6126 telescope corrections were flagged in sub-array 3. ! A total of 802 telescope corrections were flagged in sub-array 4. ! A total of 2728 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 4582 telescope corrections were flagged in sub-array 1. ! A total of 2050 telescope corrections were flagged in sub-array 2. ! A total of 6126 telescope corrections were flagged in sub-array 3. ! A total of 802 telescope corrections were flagged in sub-array 4. ! A total of 2785 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=2.185769Jy sigma=0.015811 ! Fit after self-cal, rms=1.848311Jy sigma=0.013914 ! 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.217 mas, bmaj=4.549 mas, bpa=-0.4116 degrees ! Estimated noise=930.576 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 = 1.90874 Jy ! Component: 100 - total flux cleaned = 2.41894 Jy ! Total flux subtracted in 100 components = 2.41894 Jy ! Clean residual min=-0.068842 max=0.221085 Jy/beam ! Clean residual mean=0.000626 rms=0.016538 Jy/beam ! Combined flux in latest and established models = 2.41894 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 2.41894 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.952650Jy sigma=0.006032 ! Fit after self-cal, rms=0.951340Jy sigma=0.006028 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.2378 Jy ! Component: 100 - total flux cleaned = 0.368787 Jy ! Total flux subtracted in 100 components = 0.368787 Jy ! Clean residual min=-0.058670 max=0.147387 Jy/beam ! Clean residual mean=0.000547 rms=0.013838 Jy/beam ! Combined flux in latest and established models = 2.78772 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 22 components and 2.78772 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.861866Jy sigma=0.005095 ! Fit after self-cal, rms=0.861297Jy sigma=0.005091 ! Inverting map ! Added new window around map position (3, -1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.130712 Jy ! Component: 100 - total flux cleaned = 0.197082 Jy ! Total flux subtracted in 100 components = 0.197082 Jy ! Clean residual min=-0.055817 max=0.079983 Jy/beam ! Clean residual mean=0.000506 rms=0.012987 Jy/beam ! Combined flux in latest and established models = 2.9848 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 29 components and 2.9848 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.827743Jy sigma=0.004811 ! Fit after self-cal, rms=0.824496Jy sigma=0.004798 ! Inverting map ! Added new window around map position (6.5, 0.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0749921 Jy ! Component: 100 - total flux cleaned = 0.103613 Jy ! Total flux subtracted in 100 components = 0.103613 Jy ! Clean residual min=-0.052310 max=0.075011 Jy/beam ! Clean residual mean=0.000485 rms=0.012570 Jy/beam ! Combined flux in latest and established models = 3.08842 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 40 components and 3.08842 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.813119Jy sigma=0.004703 ! Fit after self-cal, rms=0.812215Jy sigma=0.004695 ! 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.623 mas, bmaj=5.46 mas, bpa=-3.598 degrees ! Estimated noise=536.81 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 (52.5, -20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.116208 Jy ! Component: 100 - total flux cleaned = 0.198893 Jy ! Component: 150 - total flux cleaned = 0.264581 Jy ! Component: 200 - total flux cleaned = 0.317782 Jy ! Total flux subtracted in 200 components = 0.317782 Jy ! Clean residual min=-0.053338 max=0.084926 Jy/beam ! Clean residual mean=0.006747 rms=0.015667 Jy/beam ! Combined flux in latest and established models = 3.4062 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 83 components and 3.4062 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.784867Jy sigma=0.004491 ! Fit after self-cal, rms=0.783690Jy sigma=0.004478 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 87 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.16 GC 1.00 GG 1.00* ! HN 0.90 KP 0.84 LA 1.18 MA 0.78 ! MK 1.01 NL 0.98 NY 1.00* ON 1.31 ! OV 0.89 PT 1.05 SC 0.91 TC 1.09 ! TS 1.07 WF 0.92 WZ 1.04 ! ! A total of 93 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 1.19 GC 0.99 GG 1.00* ! HN 1.00* KP 0.91 LA 1.19 MA 1.00* ! MK 1.00 NL 1.02 NY 1.00* ON 1.00* ! OV 0.95 PT 1.01 SC 1.00* TC 1.06 ! TS 1.06 WF 1.00* WZ 1.00* ! ! A total of 132 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 FD 1.16 GC 0.95 GG 1.00* ! HN 0.90 KP 0.87 LA 1.17 MA 0.77 ! MK 0.88 NL 0.99 NY 1.36 ON 1.27 ! OV 1.09 PT 1.03 SC 0.85 TC 1.13 ! TS 1.00* WF 0.94 WZ 1.08 ! ! A total of 39 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 0.98 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 0.94 NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 1.00* SC 1.00* TC 0.98 ! TS 1.35 WF 1.00* WZ 1.00* ! ! A total of 174 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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* ! ! ! Correcting IF 2. ! A total of 57 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.09 GC 0.89 GG 1.00* ! HN 1.00 KP 0.98 LA 1.12 MA 1.00* ! MK 1.08 NL 1.02 NY 1.00* ON 1.16 ! OV 0.98 PT 1.04 SC 0.93 TC 1.04 ! TS 1.04 WF 0.89 WZ 0.96 ! ! A total of 93 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 1.10 GC 0.89 GG 1.00* ! HN 1.00* KP 0.99 LA 1.12 MA 1.00* ! MK 1.08 NL 1.06 NY 1.00* ON 1.00* ! OV 1.04 PT 0.99 SC 1.00* TC 1.10 ! TS 1.03 WF 1.00* WZ 1.00* ! ! A total of 135 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.13 FD 1.08 GC 0.87 GG 1.00* ! HN 1.00 KP 0.97 LA 1.10 MA 1.00* ! MK 1.03 NL 1.02 NY 0.99 ON 1.00 ! OV 1.06 PT 1.02 SC 0.87 TC 1.09 ! TS 1.00* WF 0.89 WZ 1.00 ! ! A total of 39 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 0.86 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 0.98 NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 1.00* SC 1.00* TC 1.02 ! TS 1.34 WF 1.00* WZ 1.00* ! ! A total of 57 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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* ! ! ! Correcting IF 3. ! A total of 87 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.87 FD 1.01 GC 1.11 GG 1.00* ! HN 1.01 KP 0.96 LA 0.92 MA 0.85 ! MK 1.03 NL 1.05 NY 1.00* ON 1.45 ! OV 0.98 PT 1.08 SC 1.09 TC 0.96 ! TS 1.11 WF 0.90 WZ 1.06 ! ! A total of 93 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.19 FD 1.05 GC 1.06 GG 1.00* ! HN 1.00* KP 0.98 LA 0.90 MA 1.00* ! MK 1.04 NL 1.11 NY 1.00* ON 1.00* ! OV 1.03 PT 1.02 SC 1.00* TC 0.93 ! TS 1.04 WF 1.00* WZ 1.00* ! ! A total of 99 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 1.02 GC 1.25 GG 1.00* ! HN 1.00 KP 0.96 LA 0.91 MA 0.74 ! MK 1.17 NL 1.06 NY 0.97 ON 1.28 ! OV 1.11 PT 1.05 SC 0.99 TC 0.98 ! TS 1.00* WF 0.90 WZ 1.05 ! ! A total of 39 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 1.10 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 0.90 NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 1.00* SC 1.00* TC 0.87 ! TS 1.40 WF 1.00* WZ 1.00* ! ! A total of 141 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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* ! ! ! Correcting IF 4. ! A total of 87 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.22 FD 0.92 GC 1.14 GG 1.00* ! HN 0.84 KP 0.97 LA 0.74 MA 1.05 ! MK 0.97 NL 0.99 NY 1.00* ON 1.38 ! OV 0.81 PT 1.04 SC 1.13 TC 0.89 ! TS 1.02 WF 0.86 WZ 1.08 ! ! A total of 93 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.35 FD 0.96 GC 1.11 GG 1.00* ! HN 1.00* KP 0.98 LA 0.73 MA 1.00* ! MK 0.99 NL 1.03 NY 1.00* ON 1.00* ! OV 0.87 PT 0.97 SC 1.00* TC 0.96 ! TS 1.02 WF 1.00* WZ 1.00* ! ! A total of 99 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.12 FD 0.92 GC 1.20 GG 1.00* ! HN 0.81 KP 0.97 LA 0.74 MA 1.07 ! MK 1.03 NL 0.98 NY 0.92 ON 1.25 ! OV 0.92 PT 1.01 SC 0.99 TC 1.03 ! TS 1.00* WF 0.86 WZ 1.06 ! ! A total of 39 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00* FD 1.00* GC 1.16 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 0.87 NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 1.00* SC 1.00* TC 0.94 ! TS 1.36 WF 1.00* WZ 1.00* ! ! A total of 84 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! 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* ! ! ! Fit before self-cal, rms=0.783690Jy sigma=0.004478 ! Fit after self-cal, rms=0.630077Jy sigma=0.003044 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.667 mas, bmaj=5.76 mas, bpa=-3.881 degrees ! Estimated noise=531.826 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 (121, -34). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0614775 Jy ! Component: 100 - total flux cleaned = 0.107031 Jy ! Component: 150 - total flux cleaned = 0.144619 Jy ! Component: 200 - total flux cleaned = 0.17619 Jy ! Total flux subtracted in 200 components = 0.17619 Jy ! Clean residual min=-0.033512 max=0.049075 Jy/beam ! Clean residual mean=0.004275 rms=0.008738 Jy/beam ! Combined flux in latest and established models = 3.58239 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 111 components and 3.58239 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.619468Jy sigma=0.002952 ! Fit after self-cal, rms=0.619385Jy sigma=0.002946 ! Inverting map ! Added new window around map position (60, -21). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0464145 Jy ! Component: 100 - total flux cleaned = 0.0761527 Jy ! Component: 150 - total flux cleaned = 0.101611 Jy ! Component: 200 - total flux cleaned = 0.121959 Jy ! Total flux subtracted in 200 components = 0.121959 Jy ! Clean residual min=-0.031337 max=0.041844 Jy/beam ! Clean residual mean=0.003928 rms=0.008166 Jy/beam ! Combined flux in latest and established models = 3.70435 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 138 components and 3.70435 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.613472Jy sigma=0.002901 ! Fit after self-cal, rms=0.613066Jy sigma=0.002897 ! Inverting map ! Added new window around map position (128.5, -36). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0393313 Jy ! Component: 100 - total flux cleaned = 0.0636785 Jy ! Component: 150 - total flux cleaned = 0.0811328 Jy ! Component: 200 - total flux cleaned = 0.0929451 Jy ! Total flux subtracted in 200 components = 0.0929451 Jy ! Clean residual min=-0.029629 max=0.038827 Jy/beam ! Clean residual mean=0.003701 rms=0.007920 Jy/beam ! Combined flux in latest and established models = 3.79729 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 165 components and 3.79729 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.608693Jy sigma=0.002868 ! Fit after self-cal, rms=0.608531Jy sigma=0.002866 ! 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.608531Jy sigma=0.002866 ! Fit after self-cal, rms=0.571851Jy sigma=0.002504 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.671 mas, bmaj=5.752 mas, bpa=-3.944 degrees ! Estimated noise=535.564 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0280474 Jy ! Component: 100 - total flux cleaned = 0.0420982 Jy ! Component: 150 - total flux cleaned = 0.049974 Jy ! Component: 200 - total flux cleaned = 0.0557784 Jy ! Total flux subtracted in 200 components = 0.0557784 Jy ! Clean residual min=-0.025776 max=0.036658 Jy/beam ! Clean residual mean=0.001192 rms=0.006377 Jy/beam ! Combined flux in latest and established models = 3.85307 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 186 components and 3.85307 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.569339Jy sigma=0.002484 ! Fit after self-cal, rms=0.569055Jy sigma=0.002482 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 (44.5, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0376856 Jy ! Component: 100 - total flux cleaned = 0.0604985 Jy ! Component: 150 - total flux cleaned = 0.0746687 Jy ! Component: 200 - total flux cleaned = 0.0838892 Jy ! Total flux subtracted in 200 components = 0.0838892 Jy ! Clean residual min=-0.024105 max=0.029162 Jy/beam ! Clean residual mean=0.000891 rms=0.005879 Jy/beam ! Combined flux in latest and established models = 3.93696 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 233 components and 3.93696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.566227Jy sigma=0.002442 ! Fit after self-cal, rms=0.565857Jy sigma=0.002439 ! Inverting map ! Added new window around map position (14.5, -1.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0263189 Jy ! Component: 100 - total flux cleaned = 0.0365953 Jy ! Component: 150 - total flux cleaned = 0.0436739 Jy ! Component: 200 - total flux cleaned = 0.0482001 Jy ! Total flux subtracted in 200 components = 0.0482001 Jy ! Clean residual min=-0.023548 max=0.024973 Jy/beam ! Clean residual mean=0.000708 rms=0.005652 Jy/beam ! Combined flux in latest and established models = 3.98516 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 262 components and 3.98516 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.564170Jy sigma=0.002420 ! Fit after self-cal, rms=0.564164Jy sigma=0.002418 ! Inverting map ! Added new window around map position (41.5, -244.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0208639 Jy ! Component: 100 - total flux cleaned = 0.0300352 Jy ! Component: 150 - total flux cleaned = 0.035193 Jy ! Component: 200 - total flux cleaned = 0.0387424 Jy ! Total flux subtracted in 200 components = 0.0387424 Jy ! Clean residual min=-0.022736 max=0.023969 Jy/beam ! Clean residual mean=0.000645 rms=0.005498 Jy/beam ! Combined flux in latest and established models = 4.0239 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 285 components and 4.02391 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.562825Jy sigma=0.002406 ! Fit after self-cal, rms=0.562799Jy sigma=0.002406 ! Inverting map ! Added new window around map position (79.5, -22). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.02436 Jy ! Component: 100 - total flux cleaned = 0.0395127 Jy ! Component: 150 - total flux cleaned = 0.0479775 Jy ! Component: 200 - total flux cleaned = 0.0535629 Jy ! Total flux subtracted in 200 components = 0.0535629 Jy ! Clean residual min=-0.023090 max=0.021620 Jy/beam ! Clean residual mean=0.000450 rms=0.005292 Jy/beam ! Combined flux in latest and established models = 4.07747 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 321 components and 4.07747 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.561279Jy sigma=0.002391 ! Fit after self-cal, rms=0.561209Jy sigma=0.002390 ! 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.561209Jy sigma=0.002390 ! Fit after self-cal, rms=0.528612Jy sigma=0.001978 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.672 mas, bmaj=5.745 mas, bpa=-4.012 degrees ! Estimated noise=535.133 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0133058 Jy ! Component: 100 - total flux cleaned = 0.0228905 Jy ! Component: 150 - total flux cleaned = 0.0322224 Jy ! Component: 200 - total flux cleaned = 0.0387295 Jy ! Total flux subtracted in 200 components = 0.0387295 Jy ! Clean residual min=-0.015323 max=0.018393 Jy/beam ! Clean residual mean=0.000669 rms=0.004123 Jy/beam ! Combined flux in latest and established models = 4.1162 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 345 components and 4.1162 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.527717Jy sigma=0.001967 ! Fit after self-cal, rms=0.527689Jy sigma=0.001967 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 (144, -30). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0171572 Jy ! Component: 100 - total flux cleaned = 0.0280357 Jy ! Component: 150 - total flux cleaned = 0.0346303 Jy ! Component: 200 - total flux cleaned = 0.0393282 Jy ! Total flux subtracted in 200 components = 0.0393282 Jy ! Clean residual min=-0.015782 max=0.015919 Jy/beam ! Clean residual mean=0.000558 rms=0.003988 Jy/beam ! Combined flux in latest and established models = 4.15553 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 378 components and 4.15553 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.526636Jy sigma=0.001958 ! Fit after self-cal, rms=0.526606Jy sigma=0.001957 ! 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.526606Jy sigma=0.001957 ! Fit after self-cal, rms=0.526446Jy sigma=0.001951 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.464 mas, bmaj=11.66 mas, bpa=3.427 degrees ! Estimated noise=648.933 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0235227 Jy ! Component: 100 - total flux cleaned = 0.0411499 Jy ! Component: 150 - total flux cleaned = 0.0559019 Jy ! Component: 200 - total flux cleaned = 0.0686501 Jy ! Total flux subtracted in 200 components = 0.0686501 Jy ! Clean residual min=-0.024222 max=0.029418 Jy/beam ! Clean residual mean=0.001114 rms=0.007789 Jy/beam ! Combined flux in latest and established models = 4.22418 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 401 components and 4.22418 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.525983Jy sigma=0.001948 ! Fit after self-cal, rms=0.525958Jy sigma=0.001947 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.0112468 Jy ! Component: 100 - total flux cleaned = 0.0191201 Jy ! Component: 150 - total flux cleaned = 0.0263416 Jy ! Component: 200 - total flux cleaned = 0.0319246 Jy ! Total flux subtracted in 200 components = 0.0319246 Jy ! Clean residual min=-0.024166 max=0.029429 Jy/beam ! Clean residual mean=0.000872 rms=0.007778 Jy/beam ! Combined flux in latest and established models = 4.2561 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 413 components and 4.2561 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.525837Jy sigma=0.001947 ! Fit after self-cal, rms=0.525831Jy sigma=0.001947 ! 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.672 mas, bmaj=5.744 mas, bpa=-4.03 degrees ! Estimated noise=535.185 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.672 x 5.744 at -4.03 degrees (North through East) ! Clean map min=-0.014745 max=2.5081 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.525831Jy sigma=0.001947 ! Fit after self-cal, rms=0.522562Jy sigma=0.001432 ! Inverting map and beam ! Estimated beam: bmin=2.675 mas, bmaj=5.745 mas, bpa=-4.073 degrees ! Estimated noise=534.941 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00402879 Jy ! Component: 100 - total flux cleaned = -0.00401895 Jy ! Component: 150 - total flux cleaned = -0.0040198 Jy ! Component: 200 - total flux cleaned = -0.00285021 Jy ! Total flux subtracted in 200 components = -0.00285021 Jy ! Clean residual min=-0.013608 max=0.014946 Jy/beam ! Clean residual mean=0.000276 rms=0.003616 Jy/beam ! Combined flux in latest and established models = 4.25325 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 429 components and 4.25325 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.522163Jy sigma=0.001425 ! Fit after self-cal, rms=0.522167Jy sigma=0.001425 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00115851 Jy ! Component: 100 - total flux cleaned = 0.00223956 Jy ! Component: 150 - total flux cleaned = 0.00302158 Jy ! Component: 200 - total flux cleaned = 0.00426797 Jy ! Total flux subtracted in 200 components = 0.00426797 Jy ! Clean residual min=-0.013428 max=0.014699 Jy/beam ! Clean residual mean=0.000261 rms=0.003588 Jy/beam ! Combined flux in latest and established models = 4.25752 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 439 components and 4.25752 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.521896Jy sigma=0.001420 ! Fit after self-cal, rms=0.521898Jy sigma=0.001420 ! 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.521898Jy sigma=0.001420 ! Fit after self-cal, rms=0.563382Jy sigma=0.001419 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.249 mas, bmaj=4.662 mas, bpa=-0.2886 degrees ! Estimated noise=944.393 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.9298 Jy ! Component: 100 - total flux cleaned = 2.45384 Jy ! Component: 150 - total flux cleaned = 2.70778 Jy ! Component: 200 - total flux cleaned = 2.8613 Jy ! Total flux subtracted in 200 components = 2.8613 Jy ! Clean residual min=-0.028250 max=0.082444 Jy/beam ! Clean residual mean=0.000332 rms=0.007287 Jy/beam ! Combined flux in latest and established models = 2.8613 Jy uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 27 model components to the UV plane model. ! The established model now contains 27 components and 2.8613 Jy ! Inverting map and beam ! Estimated beam: bmin=2.678 mas, bmaj=5.75 mas, bpa=-4.103 degrees ! Estimated noise=534.597 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.155383 Jy ! Component: 100 - total flux cleaned = 0.268997 Jy ! Component: 150 - total flux cleaned = 0.358732 Jy ! Component: 200 - total flux cleaned = 0.432495 Jy ! Total flux subtracted in 200 components = 0.432495 Jy ! Clean residual min=-0.028143 max=0.045652 Jy/beam ! Clean residual mean=0.003393 rms=0.009334 Jy/beam ! Combined flux in latest and established models = 3.2938 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 100 components and 3.2938 Jy ! Inverting map ! Added new window around map position (0.5, -4.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0635736 Jy ! Component: 100 - total flux cleaned = 0.119712 Jy ! Component: 150 - total flux cleaned = 0.169995 Jy ! Component: 200 - total flux cleaned = 0.215582 Jy ! Total flux subtracted in 200 components = 0.215582 Jy ! Clean residual min=-0.021748 max=0.030251 Jy/beam ! Clean residual mean=0.002650 rms=0.007318 Jy/beam ! Combined flux in latest and established models = 3.50938 Jy ! Adding 105 model components to the UV plane model. ! The established model now contains 188 components and 3.50938 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=7.471 mas, bmaj=11.66 mas, bpa=3.398 degrees ! Estimated noise=646.338 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0880341 Jy ! Component: 100 - total flux cleaned = 0.161347 Jy ! Component: 150 - total flux cleaned = 0.224215 Jy ! Component: 200 - total flux cleaned = 0.279092 Jy ! Total flux subtracted in 200 components = 0.279092 Jy ! Clean residual min=-0.025774 max=0.042477 Jy/beam ! Clean residual mean=0.004126 rms=0.011897 Jy/beam ! Combined flux in latest and established models = 3.78847 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.0484619 Jy ! Component: 100 - total flux cleaned = 0.0915683 Jy ! Component: 150 - total flux cleaned = 0.130246 Jy ! Component: 200 - total flux cleaned = 0.165149 Jy ! Total flux subtracted in 200 components = 0.165149 Jy ! Clean residual min=-0.021527 max=0.032680 Jy/beam ! Clean residual mean=0.002836 rms=0.009717 Jy/beam ! Combined flux in latest and established models = 3.95362 Jy ! Adding 254 model components to the UV plane model. ! The established model now contains 428 components and 3.95362 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.678 mas, bmaj=5.75 mas, bpa=-4.103 degrees ! Estimated noise=534.597 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0223411 Jy ! Component: 100 - total flux cleaned = 0.0403365 Jy ! Component: 150 - total flux cleaned = 0.0567007 Jy ! Component: 200 - total flux cleaned = 0.0717593 Jy ! Total flux subtracted in 200 components = 0.0717593 Jy ! Clean residual min=-0.012886 max=0.016337 Jy/beam ! Clean residual mean=0.000959 rms=0.004121 Jy/beam ! Combined flux in latest and established models = 4.02538 Jy selfcal ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 464 components and 4.02538 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.566632Jy sigma=0.001436 ! Fit after self-cal, rms=0.565793Jy sigma=0.001428 wmodel J2225-0457_S_2003_07_09_pus_map.mod ! Writing 464 model components to file: J2225-0457_S_2003_07_09_pus_map.mod wobs J2225-0457_S_2003_07_09_pus_uvs.fits ! Writing UV FITS file: J2225-0457_S_2003_07_09_pus_uvs.fits wwins J2225-0457_S_2003_07_09_pus_map.win ! wwins: Wrote 13 windows to J2225-0457_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=8192 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.00722715 Jy ! Component: 100 - total flux cleaned = 0.013823 Jy ! Component: 150 - total flux cleaned = 0.0200876 Jy ! Component: 200 - total flux cleaned = 0.0251607 Jy ! Component: 250 - total flux cleaned = 0.0301015 Jy ! Component: 300 - total flux cleaned = 0.034927 Jy ! Component: 350 - total flux cleaned = 0.0398754 Jy ! Component: 400 - total flux cleaned = 0.0445058 Jy ! Component: 450 - total flux cleaned = 0.049052 Jy ! Component: 500 - total flux cleaned = 0.0535221 Jy ! Component: 550 - total flux cleaned = 0.058126 Jy ! Component: 600 - total flux cleaned = 0.061834 Jy ! Component: 650 - total flux cleaned = 0.0658964 Jy ! Component: 700 - total flux cleaned = 0.0700996 Jy ! Component: 750 - total flux cleaned = 0.073654 Jy ! Component: 800 - total flux cleaned = 0.0775534 Jy ! Component: 850 - total flux cleaned = 0.0813959 Jy ! Component: 900 - total flux cleaned = 0.0848123 Jy ! Component: 950 - total flux cleaned = 0.0881851 Jy ! Component: 1000 - total flux cleaned = 0.0915175 Jy ! Component: 1050 - total flux cleaned = 0.0942603 Jy ! Component: 1100 - total flux cleaned = 0.0978767 Jy ! Component: 1150 - total flux cleaned = 0.100377 Jy ! Component: 1200 - total flux cleaned = 0.103385 Jy ! Component: 1250 - total flux cleaned = 0.106532 Jy ! Component: 1300 - total flux cleaned = 0.109646 Jy ! Component: 1350 - total flux cleaned = 0.112214 Jy ! Component: 1400 - total flux cleaned = 0.115266 Jy ! Component: 1450 - total flux cleaned = 0.11812 Jy ! Component: 1500 - total flux cleaned = 0.120945 Jy ! Component: 1550 - total flux cleaned = 0.123741 Jy ! Component: 1600 - total flux cleaned = 0.12651 Jy ! Component: 1650 - total flux cleaned = 0.128931 Jy ! Component: 1700 - total flux cleaned = 0.131488 Jy ! Component: 1750 - total flux cleaned = 0.13402 Jy ! Component: 1800 - total flux cleaned = 0.136998 Jy ! Component: 1850 - total flux cleaned = 0.139172 Jy ! Component: 1900 - total flux cleaned = 0.141483 Jy ! Component: 1950 - total flux cleaned = 0.144229 Jy ! Component: 2000 - total flux cleaned = 0.146348 Jy ! Component: 2050 - total flux cleaned = 0.148746 Jy ! Component: 2100 - total flux cleaned = 0.150975 Jy ! Component: 2150 - total flux cleaned = 0.153185 Jy ! Component: 2200 - total flux cleaned = 0.154939 Jy ! Component: 2250 - total flux cleaned = 0.15726 Jy ! Component: 2300 - total flux cleaned = 0.159704 Jy ! Component: 2350 - total flux cleaned = 0.161416 Jy ! Component: 2400 - total flux cleaned = 0.163539 Jy ! Component: 2450 - total flux cleaned = 0.165505 Jy ! Component: 2500 - total flux cleaned = 0.167595 Jy ! Component: 2550 - total flux cleaned = 0.169256 Jy ! Component: 2600 - total flux cleaned = 0.171179 Jy ! Component: 2650 - total flux cleaned = 0.173088 Jy ! Component: 2700 - total flux cleaned = 0.174982 Jy ! Component: 2750 - total flux cleaned = 0.176862 Jy ! Component: 2800 - total flux cleaned = 0.178329 Jy ! Component: 2850 - total flux cleaned = 0.180448 Jy ! Component: 2900 - total flux cleaned = 0.181762 Jy ! Component: 2950 - total flux cleaned = 0.183329 Jy ! Component: 3000 - total flux cleaned = 0.185016 Jy ! Component: 3050 - total flux cleaned = 0.186818 Jy ! Component: 3100 - total flux cleaned = 0.188353 Jy ! Component: 3150 - total flux cleaned = 0.190005 Jy ! Component: 3200 - total flux cleaned = 0.19152 Jy ! Component: 3250 - total flux cleaned = 0.193149 Jy ! Component: 3300 - total flux cleaned = 0.194642 Jy ! Component: 3350 - total flux cleaned = 0.196125 Jy ! Component: 3400 - total flux cleaned = 0.197721 Jy ! Component: 3450 - total flux cleaned = 0.198941 Jy ! Component: 3500 - total flux cleaned = 0.200396 Jy ! Component: 3550 - total flux cleaned = 0.20172 Jy ! Component: 3600 - total flux cleaned = 0.203395 Jy ! Component: 3650 - total flux cleaned = 0.20494 Jy ! Component: 3700 - total flux cleaned = 0.20624 Jy ! Component: 3750 - total flux cleaned = 0.207765 Jy ! Component: 3800 - total flux cleaned = 0.209513 Jy ! Component: 3850 - total flux cleaned = 0.21044 Jy ! Component: 3900 - total flux cleaned = 0.211937 Jy ! Component: 3950 - total flux cleaned = 0.213196 Jy ! Component: 4000 - total flux cleaned = 0.214447 Jy ! Component: 4050 - total flux cleaned = 0.216144 Jy ! Component: 4100 - total flux cleaned = 0.216931 Jy ! Component: 4150 - total flux cleaned = 0.218272 Jy ! Component: 4200 - total flux cleaned = 0.219494 Jy ! Component: 4250 - total flux cleaned = 0.220819 Jy ! Component: 4300 - total flux cleaned = 0.221917 Jy ! Component: 4350 - total flux cleaned = 0.223228 Jy ! Component: 4400 - total flux cleaned = 0.224204 Jy ! Component: 4450 - total flux cleaned = 0.225391 Jy ! Component: 4500 - total flux cleaned = 0.226678 Jy ! Component: 4550 - total flux cleaned = 0.227638 Jy ! Component: 4600 - total flux cleaned = 0.228911 Jy ! Component: 4650 - total flux cleaned = 0.229966 Jy ! Component: 4700 - total flux cleaned = 0.231014 Jy ! Component: 4750 - total flux cleaned = 0.232473 Jy ! Component: 4800 - total flux cleaned = 0.233406 Jy ! Component: 4850 - total flux cleaned = 0.234644 Jy ! Component: 4900 - total flux cleaned = 0.235771 Jy ! Component: 4950 - total flux cleaned = 0.236891 Jy ! Component: 5000 - total flux cleaned = 0.237905 Jy ! Component: 5050 - total flux cleaned = 0.238711 Jy ! Component: 5100 - total flux cleaned = 0.240015 Jy ! Component: 5150 - total flux cleaned = 0.241012 Jy ! Component: 5200 - total flux cleaned = 0.241707 Jy ! Component: 5250 - total flux cleaned = 0.24299 Jy ! Component: 5300 - total flux cleaned = 0.243972 Jy ! Component: 5350 - total flux cleaned = 0.244852 Jy ! Component: 5400 - total flux cleaned = 0.246018 Jy ! Component: 5450 - total flux cleaned = 0.247082 Jy ! Component: 5500 - total flux cleaned = 0.248236 Jy ! Component: 5550 - total flux cleaned = 0.249193 Jy ! Component: 5600 - total flux cleaned = 0.249955 Jy ! Component: 5650 - total flux cleaned = 0.250903 Jy ! Component: 5700 - total flux cleaned = 0.252129 Jy ! Component: 5750 - total flux cleaned = 0.25288 Jy ! Component: 5800 - total flux cleaned = 0.253721 Jy ! Component: 5850 - total flux cleaned = 0.255024 Jy ! Component: 5900 - total flux cleaned = 0.256042 Jy ! Component: 5950 - total flux cleaned = 0.256688 Jy ! Component: 6000 - total flux cleaned = 0.257514 Jy ! Component: 6050 - total flux cleaned = 0.258792 Jy ! Component: 6100 - total flux cleaned = 0.259337 Jy ! Component: 6150 - total flux cleaned = 0.260514 Jy ! Component: 6200 - total flux cleaned = 0.261054 Jy ! Component: 6250 - total flux cleaned = 0.261952 Jy ! Component: 6300 - total flux cleaned = 0.262935 Jy ! Component: 6350 - total flux cleaned = 0.263468 Jy ! Component: 6400 - total flux cleaned = 0.26462 Jy ! Component: 6450 - total flux cleaned = 0.26559 Jy ! Component: 6500 - total flux cleaned = 0.266117 Jy ! Component: 6550 - total flux cleaned = 0.266904 Jy ! Component: 6600 - total flux cleaned = 0.268036 Jy ! Component: 6650 - total flux cleaned = 0.268816 Jy ! Component: 6700 - total flux cleaned = 0.269592 Jy ! Component: 6750 - total flux cleaned = 0.270366 Jy ! Component: 6800 - total flux cleaned = 0.271393 Jy ! Component: 6850 - total flux cleaned = 0.272075 Jy ! Component: 6900 - total flux cleaned = 0.272924 Jy ! Component: 6950 - total flux cleaned = 0.274024 Jy ! Component: 7000 - total flux cleaned = 0.274445 Jy ! Component: 7050 - total flux cleaned = 0.275116 Jy ! Component: 7100 - total flux cleaned = 0.276035 Jy ! Component: 7150 - total flux cleaned = 0.276701 Jy ! Component: 7200 - total flux cleaned = 0.277281 Jy ! Component: 7250 - total flux cleaned = 0.278107 Jy ! Component: 7300 - total flux cleaned = 0.278929 Jy ! Component: 7350 - total flux cleaned = 0.279748 Jy ! Component: 7400 - total flux cleaned = 0.280565 Jy ! Component: 7450 - total flux cleaned = 0.281296 Jy ! Component: 7500 - total flux cleaned = 0.282025 Jy ! Component: 7550 - total flux cleaned = 0.28259 Jy ! Component: 7600 - total flux cleaned = 0.283554 Jy ! Component: 7650 - total flux cleaned = 0.284355 Jy ! Component: 7700 - total flux cleaned = 0.284993 Jy ! Component: 7750 - total flux cleaned = 0.285788 Jy ! Component: 7800 - total flux cleaned = 0.286579 Jy ! Component: 7850 - total flux cleaned = 0.287368 Jy ! Component: 7900 - total flux cleaned = 0.287918 Jy ! Component: 7950 - total flux cleaned = 0.288622 Jy ! Component: 8000 - total flux cleaned = 0.289714 Jy ! Component: 8050 - total flux cleaned = 0.290259 Jy ! Component: 8100 - total flux cleaned = 0.291111 Jy ! Component: 8150 - total flux cleaned = 0.291959 Jy ! Total flux subtracted in 8192 components = 0.292344 Jy ! Clean residual min=-0.004619 max=0.003940 Jy/beam ! Clean residual mean=0.000311 rms=0.001439 Jy/beam ! Combined flux in latest and established models = 4.31772 Jy keep ! Adding 864 model components to the UV plane model. ! The established model now contains 1322 components and 4.31772 Jy ![Exited script file: muppet] wmap "J2225-0457_S_2003_07_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.678 x 5.75 at -4.103 degrees (North through East) ! Clean map min=-0.0068852 max=2.5344 Jy/beam ! Writing clean map to FITS file: J2225-0457_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 2.53444 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00143007 device "J2225-0457_S_2003_07_09_pus\_map.eps/VPS" ! Attempting to open device: 'J2225-0457_S_2003_07_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.169276 0.169276 0.338552 0.677104 1.35421 2.70842 5.41683 10.8337 21.6673 43.3347 86.6693 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0042 print "peak=",int(1000*peak_flux)/1000 ! peak= 2.534 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.678 x 5.75 at -4.103 degrees (North through East) ! Clean map min=-0.0068852 max=2.5344 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J2225-0457_S_2003_07_09_pus_uvs.log closed on Mon Apr 8 13:01:24 2013