! Started logfile: J2225-0457_S_2001_04_09_pus_uvs.log on Mon Oct 8 10:45:53 2007 obs J2225-0457_S_2001_04_09_pus_uva.fits ! Reading UV FITS file: J2225-0457_S_2001_04_09_pus_uva.fits ! AN table 1: 76 integrations on 190 of 190 possible baselines. ! AN table 2: 171 integrations on 190 of 190 possible baselines. ! AN table 3: 33 integrations on 190 of 190 possible baselines. ! AN table 4: 13 integrations on 190 of 190 possible baselines. ! AN table 5: 167 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.0505721 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2225-0457 ! ! 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.22997e+09 8e+06 1 8e+06 ! 02 2 2.25997e+09 8e+06 1 8e+06 ! 03 3 2.34997e+09 8e+06 1 8e+06 ! 04 4 2.36997e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1949 lines of history. ! ! Reading 17680 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_2001_04_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 1124 telescope corrections were flagged in sub-array 1. ! A total of 2792 telescope corrections were flagged in sub-array 2. ! A total of 625 telescope corrections were flagged in sub-array 3. ! A total of 236 telescope corrections were flagged in sub-array 4. ! A total of 2854 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 1124 telescope corrections were flagged in sub-array 1. ! A total of 2792 telescope corrections were flagged in sub-array 2. ! A total of 625 telescope corrections were flagged in sub-array 3. ! A total of 236 telescope corrections were flagged in sub-array 4. ! A total of 2854 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 1124 telescope corrections were flagged in sub-array 1. ! A total of 2792 telescope corrections were flagged in sub-array 2. ! A total of 625 telescope corrections were flagged in sub-array 3. ! A total of 236 telescope corrections were flagged in sub-array 4. ! A total of 2854 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 1124 telescope corrections were flagged in sub-array 1. ! A total of 2792 telescope corrections were flagged in sub-array 2. ! A total of 625 telescope corrections were flagged in sub-array 3. ! A total of 236 telescope corrections were flagged in sub-array 4. ! A total of 2854 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=2.073743Jy sigma=0.017292 ! Fit after self-cal, rms=1.393120Jy sigma=0.011475 ! 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.322 mas, bmaj=4.937 mas, bpa=-1.089 degrees ! Estimated noise=997.566 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.71648 Jy ! Component: 100 - total flux cleaned = 2.13834 Jy ! Total flux subtracted in 100 components = 2.13834 Jy ! Clean residual min=-0.061444 max=0.179967 Jy/beam ! Clean residual mean=0.000159 rms=0.014844 Jy/beam ! Combined flux in latest and established models = 2.13834 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 2.13834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.571966Jy sigma=0.004221 ! Fit after self-cal, rms=0.569857Jy sigma=0.004210 ! Inverting map ! Added new window around map position (2.5, -1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.215214 Jy ! Component: 100 - total flux cleaned = 0.350313 Jy ! Total flux subtracted in 100 components = 0.350313 Jy ! Clean residual min=-0.046017 max=0.096551 Jy/beam ! Clean residual mean=0.000089 rms=0.010562 Jy/beam ! Combined flux in latest and established models = 2.48865 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 26 components and 2.48865 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.462072Jy sigma=0.003140 ! Fit after self-cal, rms=0.457676Jy sigma=0.003107 ! Inverting map ! Added new window around map position (6, 1.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.107979 Jy ! Component: 100 - total flux cleaned = 0.177548 Jy ! Total flux subtracted in 100 components = 0.177548 Jy ! Clean residual min=-0.040687 max=0.062663 Jy/beam ! Clean residual mean=0.000053 rms=0.009092 Jy/beam ! Combined flux in latest and established models = 2.6662 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 42 components and 2.6662 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.429927Jy sigma=0.002808 ! Fit after self-cal, rms=0.427818Jy sigma=0.002790 ! Inverting map ! Added new window around map position (51, -17.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0678725 Jy ! Component: 100 - total flux cleaned = 0.119177 Jy ! Total flux subtracted in 100 components = 0.119177 Jy ! Clean residual min=-0.037954 max=0.044909 Jy/beam ! Clean residual mean=0.000033 rms=0.008350 Jy/beam ! Combined flux in latest and established models = 2.78538 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 53 components and 2.78538 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.414322Jy sigma=0.002650 ! Fit after self-cal, rms=0.413601Jy sigma=0.002644 ! 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.705 mas, bmaj=5.939 mas, bpa=-0.9208 degrees ! Estimated noise=683.506 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 (57, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0707656 Jy ! Component: 100 - total flux cleaned = 0.120075 Jy ! Component: 150 - total flux cleaned = 0.15919 Jy ! Component: 200 - total flux cleaned = 0.19167 Jy ! Total flux subtracted in 200 components = 0.19167 Jy ! Clean residual min=-0.045010 max=0.047907 Jy/beam ! Clean residual mean=0.000251 rms=0.010005 Jy/beam ! Combined flux in latest and established models = 2.97705 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 92 components and 2.97705 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.395807Jy sigma=0.002480 ! Fit after self-cal, rms=0.394688Jy sigma=0.002471 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 42 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.00 GC 1.07 GG 1.00* ! HH 1.00* HN 1.00* KK 1.00* KP 1.01 ! LA 0.93 MA 1.00* MC 1.00* MK 1.01 ! NL 1.01 NY 1.00* OV 1.04 PT 0.92 ! SC 1.00* TS 0.97 WF 1.00* WZ 1.00* ! ! A total of 243 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.96 GC 1.04 GG 1.00* ! HH 1.00* HN 0.97 KK 1.02 KP 0.94 ! LA 0.96 MA 1.00* MC 1.00* MK 0.95 ! NL 0.98 NY 1.00* OV 1.07 PT 0.95 ! SC 0.99 TS 1.21 WF 1.02 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HH 0.97 HN 1.00* KK 1.00* KP 1.00* ! LA 1.00* MA 1.10 MC 1.33 MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.00* ! SC 1.00* TS 1.00* WF 1.00* WZ 0.93 ! ! A total of 24 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MA 0.00* MC 0.00* MK 0.00* ! NL 0.00* NY 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 93 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.99 FD 1.04 GC 1.00* GG 1.00* ! HH 1.08 HN 0.97 KK 1.00* KP 0.81 ! LA 1.00 MA 2.77 MC 0.82 MK 1.00* ! NL 0.99 NY 1.00* OV 1.04 PT 0.98 ! SC 1.00 TS 1.00* WF 1.02 WZ 0.85 ! ! ! Correcting IF 2. ! A total of 42 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.99 GC 1.02 GG 1.00* ! HH 1.00* HN 1.00* KK 1.00* KP 0.87 ! LA 0.96 MA 1.00* MC 1.00* MK 1.01 ! NL 1.01 NY 1.00* OV 1.01 PT 0.95 ! SC 1.00* TS 1.01 WF 1.00* WZ 1.00* ! ! A total of 243 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.95 GC 1.01 GG 1.00* ! HH 1.00* HN 0.97 KK 1.01 KP 1.00 ! LA 0.98 MA 1.00* MC 1.00* MK 0.97 ! NL 0.98 NY 1.00* OV 1.02 PT 0.98 ! SC 0.98 TS 1.22 WF 1.02 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HH 0.95 HN 1.00* KK 1.00* KP 1.00* ! LA 1.00* MA 0.77 MC 1.01 MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.00* ! SC 1.00* TS 1.00* WF 1.00* WZ 0.90 ! ! A total of 24 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MA 0.00* MC 0.00* MK 0.00* ! NL 0.00* NY 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 93 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 1.02 GC 1.00* GG 1.00* ! HH 1.01 HN 0.97 KK 1.00* KP 0.92 ! LA 1.03 MA 2.78 MC 0.91 MK 1.00* ! NL 0.98 NY 1.00* OV 1.01 PT 1.01 ! SC 1.01 TS 1.00* WF 1.01 WZ 0.80 ! ! ! Correcting IF 3. ! A total of 42 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 0.99 GC 0.99 GG 1.00* ! HH 1.00* HN 1.00* KK 1.00* KP 0.98 ! LA 0.92 MA 1.00* MC 1.00* MK 0.98 ! NL 0.98 NY 1.00* OV 0.99 PT 0.95 ! SC 1.00* TS 1.01 WF 1.00* WZ 1.00* ! ! A total of 243 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 0.95 GC 0.98 GG 1.00* ! HH 1.00* HN 0.98 KK 1.06 KP 1.04 ! LA 0.92 MA 1.00* MC 1.00* MK 0.93 ! NL 0.98 NY 1.00* OV 1.01 PT 0.98 ! SC 0.97 TS 1.15 WF 1.00 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HH 0.86 HN 1.00* KK 1.00* KP 1.00* ! LA 1.00* MA 1.06 MC 0.93 MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.00* ! SC 1.00* TS 1.00* WF 1.00* WZ 0.93 ! ! A total of 24 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MA 0.00* MC 0.00* MK 0.00* ! NL 0.00* NY 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 93 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.02 GC 1.00* GG 1.00* ! HH 0.95 HN 0.99 KK 1.00* KP 0.99 ! LA 0.98 MA 0.93 MC 0.83 MK 1.00* ! NL 0.96 NY 1.00* OV 0.99 PT 0.99 ! SC 0.95 TS 1.00* WF 1.00 WZ 0.84 ! ! ! Correcting IF 4. ! A total of 42 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 0.99 GC 1.02 GG 1.00* ! HH 1.00* HN 1.00* KK 1.00* KP 0.97 ! LA 0.95 MA 1.00* MC 1.00* MK 0.96 ! NL 0.97 NY 1.00* OV 0.98 PT 0.95 ! SC 1.00* TS 0.96 WF 1.00* WZ 1.00* ! ! A total of 243 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 0.94 GC 1.01 GG 1.00* ! HH 1.00* HN 0.99 KK 1.02 KP 1.00 ! LA 0.95 MA 1.00* MC 1.00* MK 0.92 ! NL 0.97 NY 1.00* OV 1.01 PT 1.11 ! SC 0.99 TS 0.99 WF 1.01 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HH 0.94 HN 1.00* KK 1.00* KP 1.00* ! LA 1.00* MA 1.02 MC 2.43 MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.00* ! SC 1.00* TS 1.00* WF 1.00* WZ 0.99 ! ! A total of 24 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HH 0.00* HN 0.00* KK 0.00* KP 0.00* ! LA 0.00* MA 0.00* MC 0.00* MK 0.00* ! NL 0.00* NY 0.00* OV 0.00* PT 0.00* ! SC 0.00* TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 93 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.95 FD 1.02 GC 1.00* GG 1.00* ! HH 1.08 HN 0.99 KK 1.00* KP 0.96 ! LA 1.02 MA 1.61 MC 1.42 MK 1.00* ! NL 0.96 NY 1.00* OV 1.00 PT 1.03 ! SC 0.98 TS 1.00* WF 1.00 WZ 0.75 ! ! ! Fit before self-cal, rms=0.394688Jy sigma=0.002471 ! Fit after self-cal, rms=0.348618Jy sigma=0.002094 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.674 mas, bmaj=6.553 mas, bpa=-1.393 degrees ! Estimated noise=688.484 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 (173.5, 168). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0373394 Jy ! Component: 100 - total flux cleaned = 0.0582971 Jy ! Component: 150 - total flux cleaned = 0.0716023 Jy ! Component: 200 - total flux cleaned = 0.0801595 Jy ! Total flux subtracted in 200 components = 0.0801595 Jy ! Clean residual min=-0.031350 max=0.031368 Jy/beam ! Clean residual mean=0.000170 rms=0.007232 Jy/beam ! Combined flux in latest and established models = 3.05721 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 126 components and 3.05721 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.341906Jy sigma=0.002032 ! Fit after self-cal, rms=0.343826Jy sigma=0.002026 ! 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.343826Jy sigma=0.002026 ! Fit after self-cal, rms=0.330138Jy sigma=0.001765 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=6.546 mas, bpa=-1.268 degrees ! Estimated noise=692.312 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0236089 Jy ! Component: 100 - total flux cleaned = 0.036309 Jy ! Component: 150 - total flux cleaned = 0.0455741 Jy ! Component: 200 - total flux cleaned = 0.0507647 Jy ! Total flux subtracted in 200 components = 0.0507647 Jy ! Clean residual min=-0.020120 max=0.029926 Jy/beam ! Clean residual mean=0.000143 rms=0.005630 Jy/beam ! Combined flux in latest and established models = 3.10797 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 151 components and 3.10797 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.327364Jy sigma=0.001735 ! Fit after self-cal, rms=0.327583Jy sigma=0.001734 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 (129, -35). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.031181 Jy ! Component: 100 - total flux cleaned = 0.0514548 Jy ! Component: 150 - total flux cleaned = 0.0667857 Jy ! Component: 200 - total flux cleaned = 0.0796851 Jy ! Total flux subtracted in 200 components = 0.0796851 Jy ! Clean residual min=-0.017710 max=0.023443 Jy/beam ! Clean residual mean=0.000126 rms=0.005108 Jy/beam ! Combined flux in latest and established models = 3.18766 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 196 components and 3.18766 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.322971Jy sigma=0.001691 ! Fit after self-cal, rms=0.322658Jy sigma=0.001689 ! Inverting map ! Added new window around map position (46.5, -16). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0264165 Jy ! Component: 100 - total flux cleaned = 0.0423656 Jy ! Component: 150 - total flux cleaned = 0.0550124 Jy ! Component: 200 - total flux cleaned = 0.0661696 Jy ! Total flux subtracted in 200 components = 0.0661696 Jy ! Clean residual min=-0.016506 max=0.018556 Jy/beam ! Clean residual mean=0.000075 rms=0.004724 Jy/beam ! Combined flux in latest and established models = 3.25383 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 232 components and 3.25383 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.319698Jy sigma=0.001656 ! Fit after self-cal, rms=0.319577Jy sigma=0.001654 ! 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.319577Jy sigma=0.001654 ! Fit after self-cal, rms=0.318718Jy sigma=0.001647 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=6.545 mas, bpa=-1.269 degrees ! Estimated noise=691.932 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0131227 Jy ! Component: 100 - total flux cleaned = 0.0240402 Jy ! Component: 150 - total flux cleaned = 0.0318287 Jy ! Component: 200 - total flux cleaned = 0.0393243 Jy ! Total flux subtracted in 200 components = 0.0393243 Jy ! Clean residual min=-0.016678 max=0.017985 Jy/beam ! Clean residual mean=0.000037 rms=0.004494 Jy/beam ! Combined flux in latest and established models = 3.29315 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 251 components and 3.29315 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.317460Jy sigma=0.001636 ! Fit after self-cal, rms=0.317397Jy sigma=0.001635 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.00764226 Jy ! Component: 100 - total flux cleaned = 0.0139936 Jy ! Component: 150 - total flux cleaned = 0.0190386 Jy ! Component: 200 - total flux cleaned = 0.0237821 Jy ! Total flux subtracted in 200 components = 0.0237821 Jy ! Clean residual min=-0.016668 max=0.017386 Jy/beam ! Clean residual mean=0.000022 rms=0.004421 Jy/beam ! Combined flux in latest and established models = 3.31694 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 263 components and 3.31693 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316696Jy sigma=0.001630 ! Fit after self-cal, rms=0.316670Jy sigma=0.001630 ! 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.316670Jy sigma=0.001630 ! Fit after self-cal, rms=0.316551Jy sigma=0.001628 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.219 mas, bmaj=10.7 mas, bpa=0.9095 degrees ! Estimated noise=932.188 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0210122 Jy ! Component: 100 - total flux cleaned = 0.0367123 Jy ! Component: 150 - total flux cleaned = 0.0496904 Jy ! Component: 200 - total flux cleaned = 0.0605767 Jy ! Total flux subtracted in 200 components = 0.0605767 Jy ! Clean residual min=-0.028201 max=0.032066 Jy/beam ! Clean residual mean=-0.000036 rms=0.009610 Jy/beam ! Combined flux in latest and established models = 3.37751 Jy selfcal ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 288 components and 3.37751 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.316034Jy sigma=0.001629 ! Fit after self-cal, rms=0.315838Jy sigma=0.001627 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.0104857 Jy ! Component: 100 - total flux cleaned = 0.018302 Jy ! Component: 150 - total flux cleaned = 0.0236421 Jy ! Component: 200 - total flux cleaned = 0.0284534 Jy ! Total flux subtracted in 200 components = 0.0284534 Jy ! Clean residual min=-0.028812 max=0.032249 Jy/beam ! Clean residual mean=-0.000084 rms=0.009557 Jy/beam ! Combined flux in latest and established models = 3.40597 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 295 components and 3.40597 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.315789Jy sigma=0.001631 ! Fit after self-cal, rms=0.315751Jy sigma=0.001631 ! 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.684 mas, bmaj=6.545 mas, bpa=-1.271 degrees ! Estimated noise=692.231 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.684 x 6.545 at -1.271 degrees (North through East) ! Clean map min=-0.017692 max=2.2138 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.315751Jy sigma=0.001631 ! Fit after self-cal, rms=0.347921Jy sigma=0.001314 ! Inverting map and beam ! Estimated beam: bmin=2.682 mas, bmaj=6.551 mas, bpa=-1.267 degrees ! Estimated noise=691.683 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00244968 Jy ! Component: 100 - total flux cleaned = -0.00297467 Jy ! Component: 150 - total flux cleaned = -0.00250067 Jy ! Component: 200 - total flux cleaned = -0.00228052 Jy ! Total flux subtracted in 200 components = -0.00228052 Jy ! Clean residual min=-0.015737 max=0.015671 Jy/beam ! Clean residual mean=-0.000035 rms=0.004110 Jy/beam ! Combined flux in latest and established models = 3.40369 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 318 components and 3.40369 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.347580Jy sigma=0.001308 ! Fit after self-cal, rms=0.347563Jy sigma=0.001308 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.0911e-06 Jy ! Component: 100 - total flux cleaned = 0.000199571 Jy ! Component: 150 - total flux cleaned = 0.000568655 Jy ! Component: 200 - total flux cleaned = 0.00110464 Jy ! Total flux subtracted in 200 components = 0.00110464 Jy ! Clean residual min=-0.015279 max=0.015609 Jy/beam ! Clean residual mean=-0.000035 rms=0.004087 Jy/beam ! Combined flux in latest and established models = 3.40479 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 326 components and 3.40479 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.347335Jy sigma=0.001304 ! Fit after self-cal, rms=0.347313Jy sigma=0.001304 ! 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.347313Jy sigma=0.001304 ! Fit after self-cal, rms=0.376918Jy sigma=0.001304 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.394 mas, bmaj=5.291 mas, bpa=0.1708 degrees ! Estimated noise=991.867 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.70447 Jy ! Component: 100 - total flux cleaned = 2.14483 Jy ! Component: 150 - total flux cleaned = 2.36018 Jy ! Component: 200 - total flux cleaned = 2.49262 Jy ! Total flux subtracted in 200 components = 2.49262 Jy ! Clean residual min=-0.032630 max=0.070625 Jy/beam ! Clean residual mean=0.000105 rms=0.007088 Jy/beam ! Combined flux in latest and established models = 2.49262 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 19 model components to the UV plane model. ! The established model now contains 19 components and 2.49262 Jy ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=6.551 mas, bpa=-1.268 degrees ! Estimated noise=691.406 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.114908 Jy ! Component: 100 - total flux cleaned = 0.198553 Jy ! Component: 150 - total flux cleaned = 0.265765 Jy ! Component: 200 - total flux cleaned = 0.321422 Jy ! Total flux subtracted in 200 components = 0.321422 Jy ! Clean residual min=-0.026153 max=0.035405 Jy/beam ! Clean residual mean=0.000410 rms=0.007805 Jy/beam ! Combined flux in latest and established models = 2.81404 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 79 components and 2.81404 Jy ! Inverting map ! Added new window around map position (51, -11.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0471655 Jy ! Component: 100 - total flux cleaned = 0.0879989 Jy ! Component: 150 - total flux cleaned = 0.12444 Jy ! Component: 200 - total flux cleaned = 0.157256 Jy ! Total flux subtracted in 200 components = 0.157256 Jy ! Clean residual min=-0.020555 max=0.024421 Jy/beam ! Clean residual mean=0.000272 rms=0.006319 Jy/beam ! Combined flux in latest and established models = 2.9713 Jy ! Adding 73 model components to the UV plane model. ! The established model now contains 137 components and 2.9713 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.224 mas, bmaj=10.7 mas, bpa=0.896 degrees ! Estimated noise=929.998 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0652989 Jy ! Component: 100 - total flux cleaned = 0.118464 Jy ! Component: 150 - total flux cleaned = 0.163182 Jy ! Component: 200 - total flux cleaned = 0.201398 Jy ! Total flux subtracted in 200 components = 0.201398 Jy ! Clean residual min=-0.031662 max=0.031989 Jy/beam ! Clean residual mean=0.000328 rms=0.011585 Jy/beam ! Combined flux in latest and established models = 3.17269 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.0330639 Jy ! Component: 100 - total flux cleaned = 0.061928 Jy ! Component: 150 - total flux cleaned = 0.0872299 Jy ! Component: 200 - total flux cleaned = 0.109507 Jy ! Total flux subtracted in 200 components = 0.109507 Jy ! Clean residual min=-0.027261 max=0.030759 Jy/beam ! Clean residual mean=0.000123 rms=0.010158 Jy/beam ! Combined flux in latest and established models = 3.2822 Jy ! Adding 202 model components to the UV plane model. ! The established model now contains 320 components and 3.2822 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=6.551 mas, bpa=-1.268 degrees ! Estimated noise=691.406 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0164902 Jy ! Component: 100 - total flux cleaned = 0.0296326 Jy ! Component: 150 - total flux cleaned = 0.0411203 Jy ! Component: 200 - total flux cleaned = 0.0489749 Jy ! Total flux subtracted in 200 components = 0.0489749 Jy ! Clean residual min=-0.014570 max=0.017294 Jy/beam ! Clean residual mean=0.000012 rms=0.004295 Jy/beam ! Combined flux in latest and established models = 3.33118 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 350 components and 3.33118 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.378919Jy sigma=0.001325 ! Fit after self-cal, rms=0.378726Jy sigma=0.001318 wmodel J2225-0457_S_2001_04_09_pus_map.mod ! Writing 350 model components to file: J2225-0457_S_2001_04_09_pus_map.mod wobs J2225-0457_S_2001_04_09_pus_uvs.fits ! Writing UV FITS file: J2225-0457_S_2001_04_09_pus_uvs.fits wwins J2225-0457_S_2001_04_09_pus_map.win ! wwins: Wrote 9 windows to J2225-0457_S_2001_04_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.00759566 Jy ! Component: 100 - total flux cleaned = 0.0137427 Jy ! Component: 150 - total flux cleaned = 0.0177296 Jy ! Component: 200 - total flux cleaned = 0.0210622 Jy ! Component: 250 - total flux cleaned = 0.0235392 Jy ! Component: 300 - total flux cleaned = 0.0259462 Jy ! Component: 350 - total flux cleaned = 0.0282918 Jy ! Component: 400 - total flux cleaned = 0.0301298 Jy ! Component: 450 - total flux cleaned = 0.0321605 Jy ! Component: 500 - total flux cleaned = 0.033708 Jy ! Component: 550 - total flux cleaned = 0.0350158 Jy ! Component: 600 - total flux cleaned = 0.0365145 Jy ! Component: 650 - total flux cleaned = 0.0384032 Jy ! Component: 700 - total flux cleaned = 0.0398445 Jy ! Component: 750 - total flux cleaned = 0.0414664 Jy ! Component: 800 - total flux cleaned = 0.0432582 Jy ! Component: 850 - total flux cleaned = 0.0444338 Jy ! Component: 900 - total flux cleaned = 0.0457831 Jy ! Component: 950 - total flux cleaned = 0.0474945 Jy ! Component: 1000 - total flux cleaned = 0.048434 Jy ! Component: 1050 - total flux cleaned = 0.0502854 Jy ! Component: 1100 - total flux cleaned = 0.0517455 Jy ! Component: 1150 - total flux cleaned = 0.0531859 Jy ! Component: 1200 - total flux cleaned = 0.0542544 Jy ! Component: 1250 - total flux cleaned = 0.0553093 Jy ! Component: 1300 - total flux cleaned = 0.0565268 Jy ! Component: 1350 - total flux cleaned = 0.0580712 Jy ! Component: 1400 - total flux cleaned = 0.0589172 Jy ! Component: 1450 - total flux cleaned = 0.0600862 Jy ! Component: 1500 - total flux cleaned = 0.0607459 Jy ! Component: 1550 - total flux cleaned = 0.0618918 Jy ! Component: 1600 - total flux cleaned = 0.0627018 Jy ! Component: 1650 - total flux cleaned = 0.0636632 Jy ! Component: 1700 - total flux cleaned = 0.0647732 Jy ! Component: 1750 - total flux cleaned = 0.0661858 Jy ! Component: 1800 - total flux cleaned = 0.0672724 Jy ! Component: 1850 - total flux cleaned = 0.0685001 Jy ! Component: 1900 - total flux cleaned = 0.0698671 Jy ! Component: 1950 - total flux cleaned = 0.0704672 Jy ! Component: 2000 - total flux cleaned = 0.0713622 Jy ! Component: 2050 - total flux cleaned = 0.0725428 Jy ! Component: 2100 - total flux cleaned = 0.0734189 Jy ! Component: 2150 - total flux cleaned = 0.0739981 Jy ! Component: 2200 - total flux cleaned = 0.075003 Jy ! Component: 2250 - total flux cleaned = 0.0757139 Jy ! Component: 2300 - total flux cleaned = 0.0769815 Jy ! Component: 2350 - total flux cleaned = 0.0775391 Jy ! Component: 2400 - total flux cleaned = 0.0780928 Jy ! Component: 2450 - total flux cleaned = 0.0793256 Jy ! Component: 2500 - total flux cleaned = 0.080278 Jy ! Component: 2550 - total flux cleaned = 0.081085 Jy ! Component: 2600 - total flux cleaned = 0.081486 Jy ! Component: 2650 - total flux cleaned = 0.0822809 Jy ! Component: 2700 - total flux cleaned = 0.0826754 Jy ! Component: 2750 - total flux cleaned = 0.0835877 Jy ! Component: 2800 - total flux cleaned = 0.0843622 Jy ! Component: 2850 - total flux cleaned = 0.0848748 Jy ! Component: 2900 - total flux cleaned = 0.0848763 Jy ! Component: 2950 - total flux cleaned = 0.0863892 Jy ! Component: 3000 - total flux cleaned = 0.0868882 Jy ! Component: 3050 - total flux cleaned = 0.0871372 Jy ! Component: 3100 - total flux cleaned = 0.0875051 Jy ! Component: 3150 - total flux cleaned = 0.0881161 Jy ! Component: 3200 - total flux cleaned = 0.0888429 Jy ! Component: 3250 - total flux cleaned = 0.0894427 Jy ! Component: 3300 - total flux cleaned = 0.0900382 Jy ! Component: 3350 - total flux cleaned = 0.0905117 Jy ! Component: 3400 - total flux cleaned = 0.0910986 Jy ! Component: 3450 - total flux cleaned = 0.0919139 Jy ! Component: 3500 - total flux cleaned = 0.0923767 Jy ! Component: 3550 - total flux cleaned = 0.093066 Jy ! Component: 3600 - total flux cleaned = 0.093407 Jy ! Component: 3650 - total flux cleaned = 0.0938591 Jy ! Component: 3700 - total flux cleaned = 0.0944209 Jy ! Component: 3750 - total flux cleaned = 0.094979 Jy ! Component: 3800 - total flux cleaned = 0.0956426 Jy ! Component: 3850 - total flux cleaned = 0.0963009 Jy ! Component: 3900 - total flux cleaned = 0.096846 Jy ! Component: 3950 - total flux cleaned = 0.0971714 Jy ! Component: 4000 - total flux cleaned = 0.0974937 Jy ! Component: 4050 - total flux cleaned = 0.0980274 Jy ! Component: 4100 - total flux cleaned = 0.098664 Jy ! Component: 4150 - total flux cleaned = 0.0988747 Jy ! Component: 4200 - total flux cleaned = 0.0993971 Jy ! Component: 4250 - total flux cleaned = 0.100124 Jy ! Component: 4300 - total flux cleaned = 0.100536 Jy ! Component: 4350 - total flux cleaned = 0.100638 Jy ! Component: 4400 - total flux cleaned = 0.101248 Jy ! Component: 4450 - total flux cleaned = 0.102156 Jy ! Component: 4500 - total flux cleaned = 0.102557 Jy ! Component: 4550 - total flux cleaned = 0.102855 Jy ! Component: 4600 - total flux cleaned = 0.103349 Jy ! Component: 4650 - total flux cleaned = 0.103742 Jy ! Component: 4700 - total flux cleaned = 0.10423 Jy ! Component: 4750 - total flux cleaned = 0.104908 Jy ! Component: 4800 - total flux cleaned = 0.105485 Jy ! Component: 4850 - total flux cleaned = 0.105579 Jy ! Component: 4900 - total flux cleaned = 0.106244 Jy ! Component: 4950 - total flux cleaned = 0.106716 Jy ! Component: 5000 - total flux cleaned = 0.106527 Jy ! Component: 5050 - total flux cleaned = 0.107926 Jy ! Component: 5100 - total flux cleaned = 0.108203 Jy ! Component: 5150 - total flux cleaned = 0.107833 Jy ! Component: 5200 - total flux cleaned = 0.108566 Jy ! Component: 5250 - total flux cleaned = 0.10902 Jy ! Component: 5300 - total flux cleaned = 0.109563 Jy ! Component: 5350 - total flux cleaned = 0.109833 Jy ! Component: 5400 - total flux cleaned = 0.11037 Jy ! Component: 5450 - total flux cleaned = 0.110814 Jy ! Component: 5500 - total flux cleaned = 0.110991 Jy ! Component: 5550 - total flux cleaned = 0.111254 Jy ! Component: 5600 - total flux cleaned = 0.111864 Jy ! Component: 5650 - total flux cleaned = 0.112212 Jy ! Component: 5700 - total flux cleaned = 0.112902 Jy ! Component: 5750 - total flux cleaned = 0.113502 Jy ! Component: 5800 - total flux cleaned = 0.113927 Jy ! Component: 5850 - total flux cleaned = 0.114097 Jy ! Component: 5900 - total flux cleaned = 0.114433 Jy ! Component: 5950 - total flux cleaned = 0.114685 Jy ! Component: 6000 - total flux cleaned = 0.115183 Jy ! Component: 6050 - total flux cleaned = 0.115432 Jy ! Component: 6100 - total flux cleaned = 0.115843 Jy ! Component: 6150 - total flux cleaned = 0.116579 Jy ! Component: 6200 - total flux cleaned = 0.116661 Jy ! Component: 6250 - total flux cleaned = 0.117065 Jy ! Component: 6300 - total flux cleaned = 0.117145 Jy ! Component: 6350 - total flux cleaned = 0.117386 Jy ! Component: 6400 - total flux cleaned = 0.117704 Jy ! Component: 6450 - total flux cleaned = 0.117942 Jy ! Component: 6500 - total flux cleaned = 0.118179 Jy ! Component: 6550 - total flux cleaned = 0.118492 Jy ! Component: 6600 - total flux cleaned = 0.118881 Jy ! Component: 6650 - total flux cleaned = 0.119268 Jy ! Component: 6700 - total flux cleaned = 0.119345 Jy ! Component: 6750 - total flux cleaned = 0.119575 Jy ! Component: 6800 - total flux cleaned = 0.119498 Jy ! Component: 6850 - total flux cleaned = 0.119727 Jy ! Component: 6900 - total flux cleaned = 0.119878 Jy ! Component: 6950 - total flux cleaned = 0.120029 Jy ! Component: 7000 - total flux cleaned = 0.120328 Jy ! Component: 7050 - total flux cleaned = 0.120253 Jy ! Component: 7100 - total flux cleaned = 0.120401 Jy ! Component: 7150 - total flux cleaned = 0.120696 Jy ! Component: 7200 - total flux cleaned = 0.120916 Jy ! Component: 7250 - total flux cleaned = 0.121135 Jy ! Component: 7300 - total flux cleaned = 0.121425 Jy ! Component: 7350 - total flux cleaned = 0.121353 Jy ! Component: 7400 - total flux cleaned = 0.121569 Jy ! Component: 7450 - total flux cleaned = 0.121926 Jy ! Component: 7500 - total flux cleaned = 0.122139 Jy ! Component: 7550 - total flux cleaned = 0.12221 Jy ! Component: 7600 - total flux cleaned = 0.122279 Jy ! Component: 7650 - total flux cleaned = 0.12284 Jy ! Component: 7700 - total flux cleaned = 0.122491 Jy ! Component: 7750 - total flux cleaned = 0.122629 Jy ! Component: 7800 - total flux cleaned = 0.122836 Jy ! Component: 7850 - total flux cleaned = 0.123043 Jy ! Component: 7900 - total flux cleaned = 0.123316 Jy ! Component: 7950 - total flux cleaned = 0.123385 Jy ! Component: 8000 - total flux cleaned = 0.123453 Jy ! Component: 8050 - total flux cleaned = 0.123385 Jy ! Component: 8100 - total flux cleaned = 0.123855 Jy ! Component: 8150 - total flux cleaned = 0.12399 Jy ! Total flux subtracted in 8192 components = 0.123989 Jy ! Clean residual min=-0.003629 max=0.003834 Jy/beam ! Clean residual mean=0.000015 rms=0.001279 Jy/beam ! Combined flux in latest and established models = 3.45516 Jy keep ! Adding 868 model components to the UV plane model. ! The established model now contains 1218 components and 3.45516 Jy ![Exited script file: muppet] wmap "J2225-0457_S_2001_04_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.684 x 6.551 at -1.268 degrees (North through East) ! Clean map min=-0.0098232 max=2.2174 Jy/beam ! Writing clean map to FITS file: J2225-0457_S_2001_04_09_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 2.21736 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00126547 device "J2225-0457_S_2001_04_09_pus\_map.eps/VPS" ! Attempting to open device: 'J2225-0457_S_2001_04_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.171213 0.171213 0.342426 0.684851 1.3697 2.7394 5.47881 10.9576 21.9152 43.8305 87.6609 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0037 print "peak=",int(1000*peak_flux)/1000 ! peak= 2.217 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.684 x 6.551 at -1.268 degrees (North through East) ! Clean map min=-0.0098232 max=2.2174 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J2225-0457_S_2001_04_09_pus\_map.gif/vgif" ! Attempting to open device: 'J2225-0457_S_2001_04_09_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J2225-0457_S_2001_04_09_pus_uvs.log closed on Mon Oct 8 10:48:45 2007