! Started logfile: J2202+4216_S_2003_07_09_pus_uvs.log on Mon Apr 8 12:56:32 2013 obs J2202+4216_S_2003_07_09_pus_uva.fits ! Reading UV FITS file: J2202+4216_S_2003_07_09_pus_uva.fits ! AN table 1: 339 integrations on 171 of 171 possible baselines. ! AN table 2: 129 integrations on 171 of 171 possible baselines. ! AN table 3: 102 integrations on 171 of 171 possible baselines. ! AN table 4: 12 integrations on 171 of 171 possible baselines. ! AN table 5: 23 integrations on 171 of 171 possible baselines. ! AN table 6: 29 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.120363 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2202+4216 ! 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. ! 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 52196 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 J2202+4216_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 4513 telescope corrections were flagged in sub-array 1. ! A total of 2152 telescope corrections were flagged in sub-array 2. ! A total of 1445 telescope corrections were flagged in sub-array 3. ! A total of 173 telescope corrections were flagged in sub-array 4. ! A total of 342 telescope corrections were flagged in sub-array 5. ! A total of 237 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 4691 telescope corrections were flagged in sub-array 1. ! A total of 2193 telescope corrections were flagged in sub-array 2. ! A total of 1445 telescope corrections were flagged in sub-array 3. ! A total of 173 telescope corrections were flagged in sub-array 4. ! A total of 352 telescope corrections were flagged in sub-array 5. ! A total of 237 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 4470 telescope corrections were flagged in sub-array 1. ! A total of 2297 telescope corrections were flagged in sub-array 2. ! A total of 1446 telescope corrections were flagged in sub-array 3. ! A total of 173 telescope corrections were flagged in sub-array 4. ! A total of 342 telescope corrections were flagged in sub-array 5. ! A total of 237 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 4484 telescope corrections were flagged in sub-array 1. ! A total of 2297 telescope corrections were flagged in sub-array 2. ! A total of 1445 telescope corrections were flagged in sub-array 3. ! A total of 173 telescope corrections were flagged in sub-array 4. ! A total of 342 telescope corrections were flagged in sub-array 5. ! A total of 237 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.003015Jy sigma=0.006619 ! Fit after self-cal, rms=0.645437Jy sigma=0.003888 ! 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.099 mas, bmaj=2.974 mas, bpa=-4.038 degrees ! Estimated noise=1012.08 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.548118 Jy ! Component: 100 - total flux cleaned = 0.790171 Jy ! Total flux subtracted in 100 components = 0.790171 Jy ! Clean residual min=-0.025278 max=0.113778 Jy/beam ! Clean residual mean=0.000542 rms=0.009221 Jy/beam ! Combined flux in latest and established models = 0.790172 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 0.790172 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.564706Jy sigma=0.003616 ! Fit after self-cal, rms=0.555235Jy sigma=0.003601 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.145023 Jy ! Component: 100 - total flux cleaned = 0.238424 Jy ! Total flux subtracted in 100 components = 0.238424 Jy ! Clean residual min=-0.021196 max=0.061596 Jy/beam ! Clean residual mean=0.000313 rms=0.006550 Jy/beam ! Combined flux in latest and established models = 1.0286 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 24 components and 1.0286 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.491762Jy sigma=0.002884 ! Fit after self-cal, rms=0.489914Jy sigma=0.002879 ! Inverting map ! Added new window around map position (0.5, -4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0737675 Jy ! Component: 100 - total flux cleaned = 0.12613 Jy ! Total flux subtracted in 100 components = 0.12613 Jy ! Clean residual min=-0.018580 max=0.037915 Jy/beam ! Clean residual mean=0.000191 rms=0.005608 Jy/beam ! Combined flux in latest and established models = 1.15473 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 32 components and 1.15473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.468939Jy sigma=0.002607 ! Fit after self-cal, rms=0.468074Jy sigma=0.002604 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0444235 Jy ! Component: 100 - total flux cleaned = 0.0778397 Jy ! Total flux subtracted in 100 components = 0.0778397 Jy ! Clean residual min=-0.016561 max=0.031377 Jy/beam ! Clean residual mean=0.000116 rms=0.005204 Jy/beam ! Combined flux in latest and established models = 1.23257 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 37 components and 1.23257 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.459038Jy sigma=0.002475 ! Fit after self-cal, rms=0.458640Jy sigma=0.002473 ! Inverting map ! Added new window around map position (3, -8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0346568 Jy ! Component: 100 - total flux cleaned = 0.0607521 Jy ! Total flux subtracted in 100 components = 0.0607521 Jy ! Clean residual min=-0.016925 max=0.022087 Jy/beam ! Clean residual mean=0.000064 rms=0.004951 Jy/beam ! Combined flux in latest and established models = 1.29332 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 55 components and 1.29332 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.452940Jy sigma=0.002390 ! Fit after self-cal, rms=0.452542Jy sigma=0.002388 ! 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.554 mas, bmaj=3.646 mas, bpa=2.198 degrees ! Estimated noise=513.893 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.0679137 Jy ! Component: 100 - total flux cleaned = 0.114835 Jy ! Component: 150 - total flux cleaned = 0.148697 Jy ! Component: 200 - total flux cleaned = 0.173389 Jy ! Total flux subtracted in 200 components = 0.173389 Jy ! Clean residual min=-0.016291 max=0.027868 Jy/beam ! Clean residual mean=0.003351 rms=0.006535 Jy/beam ! Combined flux in latest and established models = 1.46671 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 77 components and 1.46671 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.443259Jy sigma=0.002238 ! Fit after self-cal, rms=0.442179Jy sigma=0.002231 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 153 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.13 GC 1.03 GG 1.00* ! HN 0.96 KP 0.82 LA 1.07 MA 1.49 ! MK 1.02 NL 0.97 NY 1.14 ON 1.27 ! OV 1.02 PT 0.95 SC 0.83 TC 0.82 ! TS 1.00* WF 0.90 WZ 0.98 ! ! A total of 15 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 0.97 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.55 ! MK 1.00* NL 1.00* NY 1.15 ON 1.21 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 0.98 WF 1.00* WZ 1.03 ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.98 FD 1.11 GC 1.04 GG 1.00* ! HN 0.93 KP 0.92 LA 1.13 MA 1.00* ! MK 0.95 NL 0.90 NY 1.19 ON 1.00* ! OV 2.64 PT 0.98 SC 3.75 TC 1.00* ! TS 1.06 WF 0.84 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.04 FD 1.16 GC 0.99 GG 1.00* ! HN 1.00* KP 1.00* LA 1.09 MA 1.00* ! MK 1.00* NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 0.94 SC 1.00* TC 1.00* ! TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 24 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 0.91 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.17 ! MK 1.00* NL 1.00* NY 1.15 ON 1.14 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 0.97 WF 1.00* WZ 1.07 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.03 FD 1.08 GC 1.02 GG 1.00* ! HN 1.63 KP 1.03 LA 1.11 MA 1.00* ! MK 1.00* NL 0.98 NY 1.08 ON 1.50 ! OV 0.84 PT 0.96 SC 1.00* TC 1.00* ! TS 1.02 WF 0.75 WZ 1.00* ! ! ! Correcting IF 2. ! A total of 150 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 1.04 GC 0.92 GG 1.00* ! HN 0.96 KP 0.86 LA 1.01 MA 1.00* ! MK 1.08 NL 0.99 NY 0.89 ON 1.08 ! OV 1.04 PT 0.94 SC 0.87 TC 0.84 ! TS 1.00* WF 0.88 WZ 0.99 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 0.87 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 1.00* NL 1.00* NY 0.80 ON 0.98 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 1.03 WF 1.00* WZ 1.06 ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 FD 1.01 GC 0.93 GG 1.00* ! HN 0.97 KP 0.98 LA 1.07 MA 1.00* ! MK 1.03 NL 0.99 NY 0.91 ON 1.00* ! OV 1.27 PT 0.97 SC 0.97 TC 1.00* ! TS 1.05 WF 0.85 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.09 FD 1.05 GC 0.97 GG 1.00* ! HN 1.00* KP 1.00* LA 1.01 MA 1.00* ! MK 1.00* NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 1.01 SC 1.00* TC 1.00* ! TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 24 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 0.87 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.00* ! MK 1.00* NL 1.00* NY 0.86 ON 1.04 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 0.96 WF 1.00* WZ 0.97 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 1.09 FD 1.00 GC 0.91 GG 1.00* ! HN 1.15 KP 1.00 LA 1.06 MA 1.00* ! MK 1.00* NL 1.03 NY 0.83 ON 1.05 ! OV 0.78 PT 0.94 SC 1.00* TC 1.00* ! TS 1.02 WF 0.86 WZ 1.00* ! ! ! Correcting IF 3. ! A total of 138 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.00 GC 1.18 GG 1.00* ! HN 0.97 KP 0.89 LA 0.84 MA 1.01 ! MK 1.03 NL 1.02 NY 0.89 ON 1.39 ! OV 1.01 PT 0.98 SC 0.97 TC 0.87 ! TS 1.00* WF 0.90 WZ 1.00 ! ! A total of 21 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.11 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.01 ! MK 1.00* NL 1.00* NY 0.83 ON 1.25 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 1.02 WF 1.00* WZ 1.10 ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.12 FD 0.98 GC 1.23 GG 1.00* ! HN 0.97 KP 0.96 LA 0.89 MA 1.00* ! MK 0.98 NL 1.00 NY 0.89 ON 1.00* ! OV 1.60 PT 1.01 SC 1.00* TC 1.00* ! TS 1.08 WF 0.86 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.25 FD 1.04 GC 1.07 GG 1.00* ! HN 1.00* KP 1.00* LA 0.87 MA 1.00* ! MK 1.00* NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 0.99 SC 1.00* TC 1.00* ! TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 24 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.14 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 0.98 ! MK 1.00* NL 1.00* NY 0.86 ON 1.18 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 0.97 WF 1.00* WZ 1.02 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.95 FD 0.97 GC 1.13 GG 1.00* ! HN 1.22 KP 0.95 LA 0.86 MA 1.00* ! MK 1.00* NL 1.06 NY 0.83 ON 1.36 ! OV 0.68 PT 0.99 SC 1.00* TC 1.00* ! TS 1.02 WF 0.87 WZ 1.00* ! ! ! Correcting IF 4. ! A total of 141 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.19 FD 0.90 GC 1.23 GG 1.00* ! HN 1.06 KP 0.86 LA 0.68 MA 1.69 ! MK 0.96 NL 0.96 NY 0.83 ON 1.32 ! OV 0.89 PT 0.93 SC 0.97 TC 0.74 ! TS 1.00* WF 0.87 WZ 0.98 ! ! A total of 21 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.04 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.62 ! MK 1.00* NL 1.00* NY 0.78 ON 1.23 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 1.08 WF 1.00* WZ 0.97 ! ! A total of 9 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.51 FD 0.86 GC 1.17 GG 1.00* ! HN 1.01 KP 0.98 LA 0.71 MA 1.00* ! MK 0.91 NL 0.98 NY 0.80 ON 1.00* ! OV 3.81 PT 0.95 SC 1.80 TC 1.00* ! TS 1.07 WF 0.87 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.10 FD 0.92 GC 1.19 GG 1.00* ! HN 1.00* KP 1.00* LA 0.69 MA 1.00* ! MK 1.00* NL 1.00* NY 1.00* ON 1.00* ! OV 1.00* PT 0.95 SC 1.00* TC 1.00* ! TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 24 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.00* FD 1.00* GC 1.10 GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* MA 1.96 ! MK 1.00* NL 1.00* NY 0.84 ON 1.26 ! OV 1.00* PT 1.00* SC 1.00* TC 1.00* ! TS 0.95 WF 1.00* WZ 1.00 ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.85 FD 0.87 GC 1.16 GG 1.00* ! HN 2.32 KP 0.99 LA 0.71 MA 1.00* ! MK 1.00* NL 0.97 NY 0.73 ON 1.34 ! OV 0.56 PT 0.93 SC 1.00* TC 1.00* ! TS 1.01 WF 0.85 WZ 1.00* ! ! ! Fit before self-cal, rms=0.442179Jy sigma=0.002231 ! Fit after self-cal, rms=0.437370Jy sigma=0.001582 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.583 mas, bmaj=3.654 mas, bpa=2.167 degrees ! Estimated noise=512.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 (12, -20.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0201629 Jy ! Component: 100 - total flux cleaned = 0.0324964 Jy ! Component: 150 - total flux cleaned = 0.0402277 Jy ! Component: 200 - total flux cleaned = 0.0455382 Jy ! Total flux subtracted in 200 components = 0.0455382 Jy ! Clean residual min=-0.010588 max=0.013505 Jy/beam ! Clean residual mean=0.000665 rms=0.003051 Jy/beam ! Combined flux in latest and established models = 1.51224 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 108 components and 1.51224 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.436275Jy sigma=0.001566 ! Fit after self-cal, rms=0.437189Jy sigma=0.001563 ! 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.437189Jy sigma=0.001563 ! Fit after self-cal, rms=0.427202Jy sigma=0.001469 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.585 mas, bmaj=3.657 mas, bpa=2.202 degrees ! Estimated noise=511.809 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.010823 Jy ! Component: 100 - total flux cleaned = 0.0160972 Jy ! Component: 150 - total flux cleaned = 0.0196491 Jy ! Component: 200 - total flux cleaned = 0.0226344 Jy ! Total flux subtracted in 200 components = 0.0226344 Jy ! Clean residual min=-0.005701 max=0.012100 Jy/beam ! Clean residual mean=0.001041 rms=0.001710 Jy/beam ! Combined flux in latest and established models = 1.53488 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 119 components and 1.53488 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.426764Jy sigma=0.001463 ! Fit after self-cal, rms=0.426700Jy sigma=0.001462 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 (6.5, -7.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0122625 Jy ! Component: 100 - total flux cleaned = 0.01654 Jy ! Component: 150 - total flux cleaned = 0.0186281 Jy ! Component: 200 - total flux cleaned = 0.0203362 Jy ! Total flux subtracted in 200 components = 0.0203362 Jy ! Clean residual min=-0.005762 max=0.008946 Jy/beam ! Clean residual mean=0.000810 rms=0.001580 Jy/beam ! Combined flux in latest and established models = 1.55522 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 140 components and 1.55522 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.426281Jy sigma=0.001455 ! Fit after self-cal, rms=0.426180Jy sigma=0.001454 ! Inverting map ! Added new window around map position (4.5, -16). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00780984 Jy ! Component: 100 - total flux cleaned = 0.0101223 Jy ! Component: 150 - total flux cleaned = 0.011684 Jy ! Component: 200 - total flux cleaned = 0.0125761 Jy ! Total flux subtracted in 200 components = 0.0125761 Jy ! Clean residual min=-0.005839 max=0.007767 Jy/beam ! Clean residual mean=0.000668 rms=0.001508 Jy/beam ! Combined flux in latest and established models = 1.56779 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 153 components and 1.56779 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.425929Jy sigma=0.001451 ! Fit after self-cal, rms=0.425875Jy sigma=0.001451 ! Inverting map ! Added new window around map position (-12.5, -21). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00799602 Jy ! Component: 100 - total flux cleaned = 0.0112546 Jy ! Component: 150 - total flux cleaned = 0.0135762 Jy ! Component: 200 - total flux cleaned = 0.0150619 Jy ! Total flux subtracted in 200 components = 0.0150619 Jy ! Clean residual min=-0.005330 max=0.005814 Jy/beam ! Clean residual mean=0.000492 rms=0.001392 Jy/beam ! Combined flux in latest and established models = 1.58285 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 177 components and 1.58285 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.425623Jy sigma=0.001449 ! Fit after self-cal, rms=0.425575Jy sigma=0.001448 ! 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.425575Jy sigma=0.001448 ! Fit after self-cal, rms=0.410426Jy sigma=0.001432 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.589 mas, bmaj=3.656 mas, bpa=2.053 degrees ! Estimated noise=511.844 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00659071 Jy ! Component: 100 - total flux cleaned = 0.0109421 Jy ! Component: 150 - total flux cleaned = 0.0131138 Jy ! Component: 200 - total flux cleaned = 0.014733 Jy ! Total flux subtracted in 200 components = 0.014733 Jy ! Clean residual min=-0.004214 max=0.004610 Jy/beam ! Clean residual mean=0.000514 rms=0.001064 Jy/beam ! Combined flux in latest and established models = 1.59759 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 189 components and 1.59759 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.410239Jy sigma=0.001430 ! Fit after self-cal, rms=0.410181Jy sigma=0.001430 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 (12, -12.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00461955 Jy ! Component: 100 - total flux cleaned = 0.00587893 Jy ! Component: 150 - total flux cleaned = 0.00697797 Jy ! Component: 200 - total flux cleaned = 0.00786579 Jy ! Total flux subtracted in 200 components = 0.00786579 Jy ! Clean residual min=-0.004061 max=0.004458 Jy/beam ! Clean residual mean=0.000427 rms=0.001021 Jy/beam ! Combined flux in latest and established models = 1.60545 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 202 components and 1.60545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.410052Jy sigma=0.001429 ! Fit after self-cal, rms=0.410035Jy sigma=0.001429 ! Inverting map ! Added new window around map position (25, -42.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00312469 Jy ! Component: 100 - total flux cleaned = 0.00500936 Jy ! Component: 150 - total flux cleaned = 0.00617049 Jy ! Component: 200 - total flux cleaned = 0.00695365 Jy ! Total flux subtracted in 200 components = 0.00695365 Jy ! Clean residual min=-0.004219 max=0.004161 Jy/beam ! Clean residual mean=0.000351 rms=0.000982 Jy/beam ! Combined flux in latest and established models = 1.61241 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 217 components and 1.61241 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.409930Jy sigma=0.001428 ! Fit after self-cal, rms=0.409916Jy sigma=0.001428 ! 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.409916Jy sigma=0.001428 ! Fit after self-cal, rms=0.408716Jy sigma=0.001427 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.247 mas, bmaj=8.756 mas, bpa=-82.25 degrees ! Estimated noise=670.408 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0071068 Jy ! Component: 100 - total flux cleaned = 0.0109927 Jy ! Component: 150 - total flux cleaned = 0.0131348 Jy ! Component: 200 - total flux cleaned = 0.0150269 Jy ! Total flux subtracted in 200 components = 0.0150269 Jy ! Clean residual min=-0.003631 max=0.005382 Jy/beam ! Clean residual mean=0.000908 rms=0.001396 Jy/beam ! Combined flux in latest and established models = 1.62743 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 251 components and 1.62743 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.408680Jy sigma=0.001427 ! Fit after self-cal, rms=0.408654Jy sigma=0.001427 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.00167733 Jy ! Component: 100 - total flux cleaned = 0.00312429 Jy ! Component: 150 - total flux cleaned = 0.00430149 Jy ! Component: 200 - total flux cleaned = 0.0050678 Jy ! Total flux subtracted in 200 components = 0.0050678 Jy ! Clean residual min=-0.003422 max=0.004955 Jy/beam ! Clean residual mean=0.000740 rms=0.001320 Jy/beam ! Combined flux in latest and established models = 1.6325 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 269 components and 1.6325 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.408651Jy sigma=0.001427 ! Fit after self-cal, rms=0.408633Jy sigma=0.001426 ! 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.59 mas, bmaj=3.654 mas, bpa=2.035 degrees ! Estimated noise=512.101 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.59 x 3.654 at 2.035 degrees (North through East) ! Clean map min=-0.0032737 max=0.76513 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.408633Jy sigma=0.001426 ! Fit after self-cal, rms=0.461297Jy sigma=0.001300 ! Inverting map and beam ! Estimated beam: bmin=2.588 mas, bmaj=3.642 mas, bpa=1.844 degrees ! Estimated noise=511.557 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00272512 Jy ! Component: 100 - total flux cleaned = 0.0031472 Jy ! Component: 150 - total flux cleaned = 0.00333374 Jy ! Component: 200 - total flux cleaned = 0.00350228 Jy ! Total flux subtracted in 200 components = 0.00350228 Jy ! Clean residual min=-0.003528 max=0.004053 Jy/beam ! Clean residual mean=0.000304 rms=0.000821 Jy/beam ! Combined flux in latest and established models = 1.636 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 287 components and 1.636 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.461231Jy sigma=0.001300 ! Fit after self-cal, rms=0.461229Jy sigma=0.001300 ! Inverting map ! Added new window around map position (-5.5, -45). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0026969 Jy ! Component: 100 - total flux cleaned = 0.00333229 Jy ! Component: 150 - total flux cleaned = 0.00382645 Jy ! Component: 200 - total flux cleaned = 0.0040538 Jy ! Total flux subtracted in 200 components = 0.0040538 Jy ! Clean residual min=-0.003608 max=0.003173 Jy/beam ! Clean residual mean=0.000257 rms=0.000801 Jy/beam ! Combined flux in latest and established models = 1.64006 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 298 components and 1.64006 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.461146Jy sigma=0.001299 ! Fit after self-cal, rms=0.461144Jy sigma=0.001299 ! 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.461144Jy sigma=0.001299 ! Fit after self-cal, rms=0.458651Jy sigma=0.001299 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.155 mas, bmaj=2.937 mas, bpa=-3.193 degrees ! Estimated noise=1040.98 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.565796 Jy ! Component: 100 - total flux cleaned = 0.820974 Jy ! Component: 150 - total flux cleaned = 0.963964 Jy ! Component: 200 - total flux cleaned = 1.05833 Jy ! Total flux subtracted in 200 components = 1.05833 Jy ! Clean residual min=-0.014618 max=0.053347 Jy/beam ! Clean residual mean=0.000423 rms=0.004966 Jy/beam ! Combined flux in latest and established models = 1.05833 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 30 model components to the UV plane model. ! The established model now contains 30 components and 1.05833 Jy ! Inverting map and beam ! Estimated beam: bmin=2.589 mas, bmaj=3.64 mas, bpa=1.841 degrees ! Estimated noise=511.462 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.137017 Jy ! Component: 100 - total flux cleaned = 0.226042 Jy ! Component: 150 - total flux cleaned = 0.288377 Jy ! Component: 200 - total flux cleaned = 0.333961 Jy ! Total flux subtracted in 200 components = 0.333961 Jy ! Clean residual min=-0.005296 max=0.026931 Jy/beam ! Clean residual mean=0.003137 rms=0.004192 Jy/beam ! Combined flux in latest and established models = 1.39229 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 94 components and 1.39229 Jy ! Inverting map ! Added new window around map position (0.5, -10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0356128 Jy ! Component: 100 - total flux cleaned = 0.0642838 Jy ! Component: 150 - total flux cleaned = 0.0877838 Jy ! Component: 200 - total flux cleaned = 0.107402 Jy ! Total flux subtracted in 200 components = 0.107402 Jy ! Clean residual min=-0.003380 max=0.012352 Jy/beam ! Clean residual mean=0.001886 rms=0.002214 Jy/beam ! Combined flux in latest and established models = 1.4997 Jy ! Adding 94 model components to the UV plane model. ! The established model now contains 164 components and 1.4997 Jy ! Inverting map ! Added new window around map position (-2, -4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0166236 Jy ! Component: 100 - total flux cleaned = 0.0307403 Jy ! Component: 150 - total flux cleaned = 0.0428205 Jy ! Component: 200 - total flux cleaned = 0.0532222 Jy ! Total flux subtracted in 200 components = 0.0532222 Jy ! Clean residual min=-0.002900 max=0.007647 Jy/beam ! Clean residual mean=0.001277 rms=0.001411 Jy/beam ! Combined flux in latest and established models = 1.55292 Jy ! Adding 76 model components to the UV plane model. ! The established model now contains 209 components and 1.55292 Jy ! Inverting map ! Added new window around map position (5, -23). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00925025 Jy ! Component: 100 - total flux cleaned = 0.0172281 Jy ! Component: 150 - total flux cleaned = 0.0242934 Jy ! Component: 200 - total flux cleaned = 0.0305828 Jy ! Total flux subtracted in 200 components = 0.0305828 Jy ! Clean residual min=-0.002846 max=0.004751 Jy/beam ! Clean residual mean=0.000930 rms=0.001054 Jy/beam ! Combined flux in latest and established models = 1.5835 Jy ! Adding 83 model components to the UV plane model. ! The established model now contains 254 components and 1.5835 Jy ! Inverting map ! Added new window around map position (-2.5, -23). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00577372 Jy ! Component: 100 - total flux cleaned = 0.0109345 Jy ! Component: 150 - total flux cleaned = 0.0156549 Jy ! Component: 200 - total flux cleaned = 0.0199793 Jy ! Total flux subtracted in 200 components = 0.0199793 Jy ! Clean residual min=-0.002811 max=0.003723 Jy/beam ! Clean residual mean=0.000704 rms=0.000889 Jy/beam ! Combined flux in latest and established models = 1.60348 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 290 components and 1.60348 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.243 mas, bmaj=8.749 mas, bpa=-82.56 degrees ! Estimated noise=670.357 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00891617 Jy ! Component: 100 - total flux cleaned = 0.0156038 Jy ! Component: 150 - total flux cleaned = 0.0210212 Jy ! Component: 200 - total flux cleaned = 0.0254877 Jy ! Total flux subtracted in 200 components = 0.0254877 Jy ! Clean residual min=-0.002782 max=0.005331 Jy/beam ! Clean residual mean=0.001240 rms=0.001183 Jy/beam ! Combined flux in latest and established models = 1.62897 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) ! Added new window around map position (10, -44). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00534777 Jy ! Component: 100 - total flux cleaned = 0.00886069 Jy ! Component: 150 - total flux cleaned = 0.0117108 Jy ! Component: 200 - total flux cleaned = 0.0141076 Jy ! Total flux subtracted in 200 components = 0.0141076 Jy ! Clean residual min=-0.002565 max=0.003858 Jy/beam ! Clean residual mean=0.000780 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 1.64308 Jy ! Adding 175 model components to the UV plane model. ! The established model now contains 439 components and 1.64308 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.589 mas, bmaj=3.64 mas, bpa=1.841 degrees ! Estimated noise=511.462 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00229744 Jy ! Component: 100 - total flux cleaned = 0.002893 Jy ! Component: 150 - total flux cleaned = 0.0031319 Jy ! Component: 200 - total flux cleaned = 0.00351235 Jy ! Total flux subtracted in 200 components = 0.00351235 Jy ! Clean residual min=-0.003046 max=0.002810 Jy/beam ! Clean residual mean=0.000222 rms=0.000696 Jy/beam ! Combined flux in latest and established models = 1.64659 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 469 components and 1.64659 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.458512Jy sigma=0.001297 ! Fit after self-cal, rms=0.458401Jy sigma=0.001296 wmodel J2202+4216_S_2003_07_09_pus_map.mod ! Writing 469 model components to file: J2202+4216_S_2003_07_09_pus_map.mod wobs J2202+4216_S_2003_07_09_pus_uvs.fits ! Writing UV FITS file: J2202+4216_S_2003_07_09_pus_uvs.fits wwins J2202+4216_S_2003_07_09_pus_map.win ! wwins: Wrote 15 windows to J2202+4216_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.000602779 Jy ! Component: 100 - total flux cleaned = -0.000453327 Jy ! Component: 150 - total flux cleaned = -0.000312931 Jy ! Component: 200 - total flux cleaned = -4.01507e-05 Jy ! Component: 250 - total flux cleaned = 0.000313978 Jy ! Component: 300 - total flux cleaned = 0.000746423 Jy ! Component: 350 - total flux cleaned = 0.00100128 Jy ! Component: 400 - total flux cleaned = 0.00141821 Jy ! Component: 450 - total flux cleaned = 0.00170536 Jy ! Component: 500 - total flux cleaned = 0.00202871 Jy ! Component: 550 - total flux cleaned = 0.00230786 Jy ! Component: 600 - total flux cleaned = 0.00262302 Jy ! Component: 650 - total flux cleaned = 0.00297301 Jy ! Component: 700 - total flux cleaned = 0.0032804 Jy ! Component: 750 - total flux cleaned = 0.00362238 Jy ! Component: 800 - total flux cleaned = 0.00392285 Jy ! Component: 850 - total flux cleaned = 0.00422025 Jy ! Component: 900 - total flux cleaned = 0.0045512 Jy ! Component: 950 - total flux cleaned = 0.004806 Jy ! Component: 1000 - total flux cleaned = 0.00502301 Jy ! Component: 1050 - total flux cleaned = 0.00534541 Jy ! Component: 1100 - total flux cleaned = 0.00566454 Jy ! Component: 1150 - total flux cleaned = 0.00587548 Jy ! Component: 1200 - total flux cleaned = 0.00622401 Jy ! Component: 1250 - total flux cleaned = 0.00643131 Jy ! Component: 1300 - total flux cleaned = 0.00667124 Jy ! Component: 1350 - total flux cleaned = 0.00690941 Jy ! Component: 1400 - total flux cleaned = 0.007112 Jy ! Component: 1450 - total flux cleaned = 0.00727957 Jy ! Component: 1500 - total flux cleaned = 0.00747946 Jy ! Component: 1550 - total flux cleaned = 0.00764486 Jy ! Component: 1600 - total flux cleaned = 0.00777629 Jy ! Component: 1650 - total flux cleaned = 0.00803766 Jy ! Component: 1700 - total flux cleaned = 0.00819981 Jy ! Component: 1750 - total flux cleaned = 0.00839311 Jy ! Component: 1800 - total flux cleaned = 0.00845724 Jy ! Component: 1850 - total flux cleaned = 0.00861648 Jy ! Component: 1900 - total flux cleaned = 0.0087433 Jy ! Component: 1950 - total flux cleaned = 0.00899553 Jy ! Component: 2000 - total flux cleaned = 0.00915211 Jy ! Total flux subtracted in 2048 components = 0.00918342 Jy ! Clean residual min=-0.001873 max=0.002236 Jy/beam ! Clean residual mean=0.000134 rms=0.000585 Jy/beam ! Combined flux in latest and established models = 1.65577 Jy keep ! Adding 150 model components to the UV plane model. ! The established model now contains 618 components and 1.65577 Jy ![Exited script file: muppet] wmap "J2202+4216_S_2003_07_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.589 x 3.64 at 1.841 degrees (North through East) ! Clean map min=-0.0026652 max=0.76163 Jy/beam ! Writing clean map to FITS file: J2202+4216_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.761635 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00058431 device "J2202+4216_S_2003_07_09_pus\_map.eps/VPS" ! Attempting to open device: 'J2202+4216_S_2003_07_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.230153 0.230153 0.460307 0.920614 1.84123 3.68245 7.36491 14.7298 29.4596 58.9193 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0017 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.761 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.589 x 3.64 at 1.841 degrees (North through East) ! Clean map min=-0.0026652 max=0.76163 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J2202+4216_S_2003_07_09_pus_uvs.log closed on Mon Apr 8 12:56:59 2013