! Started logfile: J2359+3850_S_2008_07_09_pus_uvs.log on Wed Oct 13 13:19:01 2010 obs J2359+3850_S_2008_07_09_pus_uva.fits ! Reading UV FITS file: J2359+3850_S_2008_07_09_pus_uva.fits ! AN table 1: 654 integrations on 105 of 105 possible baselines. ! AN table 2: 108 integrations on 105 of 105 possible baselines. ! AN table 3: 130 integrations on 105 of 105 possible baselines. ! Apparent sampling: 0.260709 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2359+3850 ! ! 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 1054 lines of history. ! ! Reading 97672 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 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 5127 telescope corrections were flagged in sub-array 1. ! A total of 982 telescope corrections were flagged in sub-array 2. ! A total of 1152 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 5052 telescope corrections were flagged in sub-array 1. ! A total of 982 telescope corrections were flagged in sub-array 2. ! A total of 1152 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 5052 telescope corrections were flagged in sub-array 1. ! A total of 982 telescope corrections were flagged in sub-array 2. ! A total of 1152 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 5052 telescope corrections were flagged in sub-array 1. ! A total of 982 telescope corrections were flagged in sub-array 2. ! A total of 1152 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.807645Jy sigma=0.005274 ! Fit after self-cal, rms=0.731376Jy sigma=0.004914 ! 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=1.888 mas, bmaj=3.076 mas, bpa=-15.17 degrees ! Estimated noise=818.737 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.238687 Jy ! Component: 100 - total flux cleaned = 0.306191 Jy ! Total flux subtracted in 100 components = 0.306191 Jy ! Clean residual min=-0.008171 max=0.030041 Jy/beam ! Clean residual mean=0.000062 rms=0.002407 Jy/beam ! Combined flux in latest and established models = 0.306191 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.306191 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.290639Jy sigma=0.001267 ! Fit after self-cal, rms=0.290380Jy sigma=0.001267 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0333436 Jy ! Component: 100 - total flux cleaned = 0.0527922 Jy ! Total flux subtracted in 100 components = 0.0527922 Jy ! Clean residual min=-0.005688 max=0.010073 Jy/beam ! Clean residual mean=0.000001 rms=0.001485 Jy/beam ! Combined flux in latest and established models = 0.358983 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 15 components and 0.358983 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.284808Jy sigma=0.001208 ! Fit after self-cal, rms=0.284478Jy sigma=0.001207 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0134667 Jy ! Component: 100 - total flux cleaned = 0.0222613 Jy ! Total flux subtracted in 100 components = 0.0222613 Jy ! Clean residual min=-0.004863 max=0.004962 Jy/beam ! Clean residual mean=-0.000023 rms=0.001300 Jy/beam ! Combined flux in latest and established models = 0.381244 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 19 components and 0.381244 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.283338Jy sigma=0.001200 ! Fit after self-cal, rms=0.283193Jy sigma=0.001200 ! 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.282 mas, bmaj=3.665 mas, bpa=-12.83 degrees ! Estimated noise=352.589 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.00657495 Jy ! Component: 100 - total flux cleaned = 0.0100163 Jy ! Component: 150 - total flux cleaned = 0.0108039 Jy ! Component: 200 - total flux cleaned = 0.0113952 Jy ! Total flux subtracted in 200 components = 0.0113952 Jy ! Clean residual min=-0.004388 max=0.003733 Jy/beam ! Clean residual mean=-0.000192 rms=0.001138 Jy/beam ! Combined flux in latest and established models = 0.39264 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 26 components and 0.392639 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.282817Jy sigma=0.001199 ! Fit after self-cal, rms=0.282726Jy sigma=0.001199 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 174 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.04 HH 1.00* HN 0.99 ! KK 1.02 KP 1.00 LA 1.01 MK 0.99 ! NL 1.00 NY 0.98 ON 0.86 PT 1.02 ! SC 1.04 WZ 0.93 ZC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.00* HH 1.00* HN 1.00* ! KK 1.05 KP 1.00* LA 1.00* MK 1.00* ! NL 1.00* NY 0.95 ON 0.81 PT 1.00* ! SC 1.00* WZ 0.92 ZC 0.98 ! ! A total of 12 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HH 1.00* HN 0.99 ! KK 1.00* KP 1.00* LA 1.00* MK 1.00* ! NL 1.13 NY 0.95 ON 0.87 PT 1.00* ! SC 1.01 WZ 0.86 ZC 1.01 ! ! ! Correcting IF 2. ! A total of 249 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 1.03 HH 1.00* HN 1.00 ! KK 0.92 KP 1.03 LA 1.02 MK 0.99 ! NL 0.96 NY 1.01 ON 0.91 PT 1.03 ! SC 1.03 WZ 0.99 ZC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.00* HH 1.00* HN 1.00* ! KK 0.99 KP 1.00* LA 1.00* MK 1.00* ! NL 1.00* NY 0.92 ON 0.84 PT 1.00* ! SC 1.00* WZ 0.97 ZC 0.99 ! ! A total of 12 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HH 1.00* HN 0.98 ! KK 1.00* KP 1.00* LA 1.00* MK 1.00* ! NL 1.03 NY 0.94 ON 0.90 PT 1.00* ! SC 1.00 WZ 1.00 ZC 0.99 ! ! ! Correcting IF 3. ! A total of 249 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.03 HH 1.00* HN 1.00 ! KK 1.02 KP 1.04 LA 1.02 MK 0.96 ! NL 1.03 NY 0.99 ON 0.90 PT 1.06 ! SC 1.05 WZ 0.99 ZC 1.04 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.06 FD 1.00* HH 1.00* HN 1.00* ! KK 1.06 KP 1.00* LA 1.00* MK 1.00* ! NL 1.00* NY 1.00 ON 0.83 PT 1.00* ! SC 1.00* WZ 0.99 ZC 1.03 ! ! A total of 12 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HH 1.00* HN 1.05 ! KK 1.00* KP 1.00* LA 1.00* MK 1.00* ! NL 1.09 NY 0.94 ON 0.84 PT 1.00* ! SC 1.04 WZ 0.89 ZC 1.06 ! ! ! Correcting IF 4. ! A total of 249 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.03 HH 1.00* HN 1.00 ! KK 1.01 KP 1.01 LA 1.02 MK 0.99 ! NL 0.99 NY 0.95 ON 0.91 PT 1.03 ! SC 1.08 WZ 0.97 ZC 0.81 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.13 FD 1.00* HH 1.00* HN 1.00* ! KK 1.01 KP 1.00* LA 1.00* MK 1.00* ! NL 1.00* NY 0.90 ON 0.96 PT 1.00* ! SC 1.00* WZ 1.01 ZC 0.79 ! ! A total of 12 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* HH 1.00* HN 1.02 ! KK 1.00* KP 1.00* LA 1.00* MK 1.00* ! NL 1.11 NY 0.87 ON 0.91 PT 1.00* ! SC 1.07 WZ 0.96 ZC 0.78 ! ! ! Fit before self-cal, rms=0.282726Jy sigma=0.001199 ! Fit after self-cal, rms=0.256984Jy sigma=0.001190 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.26 mas, bmaj=3.643 mas, bpa=-12.97 degrees ! Estimated noise=350.72 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.00219568 Jy ! Component: 100 - total flux cleaned = 0.00280337 Jy ! Component: 150 - total flux cleaned = 0.00321693 Jy ! Component: 200 - total flux cleaned = 0.00350248 Jy ! Total flux subtracted in 200 components = 0.00350248 Jy ! Clean residual min=-0.002989 max=0.002791 Jy/beam ! Clean residual mean=-0.000106 rms=0.000730 Jy/beam ! Combined flux in latest and established models = 0.396142 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 31 components and 0.396142 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.256923Jy sigma=0.001190 ! Fit after self-cal, rms=0.256753Jy sigma=0.001190 ! 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.256753Jy sigma=0.001190 ! Fit after self-cal, rms=0.255464Jy sigma=0.001186 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.255 mas, bmaj=3.664 mas, bpa=-12.97 degrees ! Estimated noise=351.275 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00190687 Jy ! Component: 100 - total flux cleaned = 0.00206716 Jy ! Component: 150 - total flux cleaned = 0.00215219 Jy ! Component: 200 - total flux cleaned = 0.00218922 Jy ! Total flux subtracted in 200 components = 0.00218922 Jy ! Clean residual min=-0.001655 max=0.001784 Jy/beam ! Clean residual mean=0.000044 rms=0.000452 Jy/beam ! Combined flux in latest and established models = 0.398331 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 35 components and 0.398331 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255437Jy sigma=0.001186 ! Fit after self-cal, rms=0.255424Jy sigma=0.001186 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 = 2.7869e-05 Jy ! Component: 100 - total flux cleaned = -8.89765e-06 Jy ! Component: 150 - total flux cleaned = -8.95745e-06 Jy ! Component: 200 - total flux cleaned = 1.63487e-05 Jy ! Total flux subtracted in 200 components = 1.63487e-05 Jy ! Clean residual min=-0.001599 max=0.001797 Jy/beam ! Clean residual mean=0.000042 rms=0.000449 Jy/beam ! Combined flux in latest and established models = 0.398348 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 38 components and 0.398347 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255413Jy sigma=0.001186 ! Fit after self-cal, rms=0.255410Jy sigma=0.001186 ! 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.255410Jy sigma=0.001186 ! Fit after self-cal, rms=0.255188Jy sigma=0.001186 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.25 mas, bmaj=3.664 mas, bpa=-12.95 degrees ! Estimated noise=351.22 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00105977 Jy ! Component: 100 - total flux cleaned = 0.00119294 Jy ! Component: 150 - total flux cleaned = 0.00122055 Jy ! Component: 200 - total flux cleaned = 0.00124491 Jy ! Total flux subtracted in 200 components = 0.00124491 Jy ! Clean residual min=-0.001530 max=0.001683 Jy/beam ! Clean residual mean=0.000046 rms=0.000432 Jy/beam ! Combined flux in latest and established models = 0.399592 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 38 components and 0.399592 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255178Jy sigma=0.001186 ! Fit after self-cal, rms=0.255181Jy sigma=0.001186 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 = 6.00415e-05 Jy ! Component: 100 - total flux cleaned = -7.31032e-06 Jy ! Component: 150 - total flux cleaned = -6.79572e-06 Jy ! Component: 200 - total flux cleaned = -6.56463e-06 Jy ! Total flux subtracted in 200 components = -6.56463e-06 Jy ! Clean residual min=-0.001495 max=0.001685 Jy/beam ! Clean residual mean=0.000046 rms=0.000430 Jy/beam ! Combined flux in latest and established models = 0.399586 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 39 components and 0.399586 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255178Jy sigma=0.001186 ! Fit after self-cal, rms=0.255176Jy sigma=0.001186 ! 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.255176Jy sigma=0.001186 ! Fit after self-cal, rms=0.255166Jy sigma=0.001186 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.377 mas, bmaj=8.045 mas, bpa=-14.92 degrees ! Estimated noise=533.316 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00127652 Jy ! Component: 100 - total flux cleaned = 0.00162432 Jy ! Component: 150 - total flux cleaned = 0.0016511 Jy ! Component: 200 - total flux cleaned = 0.00166593 Jy ! Total flux subtracted in 200 components = 0.00166593 Jy ! Clean residual min=-0.001966 max=0.003093 Jy/beam ! Clean residual mean=0.000161 rms=0.000783 Jy/beam ! Combined flux in latest and established models = 0.401252 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 42 components and 0.401252 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255169Jy sigma=0.001186 ! Fit after self-cal, rms=0.255163Jy sigma=0.001186 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 = 4.56886e-05 Jy ! Component: 100 - total flux cleaned = 4.48769e-05 Jy ! Component: 150 - total flux cleaned = 6.31705e-05 Jy ! Component: 200 - total flux cleaned = 8.57489e-05 Jy ! Total flux subtracted in 200 components = 8.57489e-05 Jy ! Clean residual min=-0.001964 max=0.003068 Jy/beam ! Clean residual mean=0.000160 rms=0.000782 Jy/beam ! Combined flux in latest and established models = 0.401338 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 44 components and 0.401338 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255165Jy sigma=0.001186 ! Fit after self-cal, rms=0.255162Jy sigma=0.001186 ! 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.247 mas, bmaj=3.664 mas, bpa=-12.93 degrees ! Estimated noise=351.213 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.247 x 3.664 at -12.93 degrees (North through East) ! Clean map min=-0.0013559 max=0.31981 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 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=1.864 mas, bmaj=3.093 mas, bpa=-15.34 degrees ! Estimated noise=806.266 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.23428 Jy ! Component: 100 - total flux cleaned = 0.303389 Jy ! Component: 150 - total flux cleaned = 0.338151 Jy ! Component: 200 - total flux cleaned = 0.359531 Jy ! Total flux subtracted in 200 components = 0.359531 Jy ! Clean residual min=-0.004834 max=0.011114 Jy/beam ! Clean residual mean=0.000040 rms=0.001328 Jy/beam ! Combined flux in latest and established models = 0.359531 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 13 model components to the UV plane model. ! The established model now contains 13 components and 0.359531 Jy ! Inverting map and beam ! Estimated beam: bmin=2.247 mas, bmaj=3.664 mas, bpa=-12.93 degrees ! Estimated noise=351.213 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.0187684 Jy ! Component: 100 - total flux cleaned = 0.0285753 Jy ! Component: 150 - total flux cleaned = 0.0341145 Jy ! Component: 200 - total flux cleaned = 0.0374241 Jy ! Total flux subtracted in 200 components = 0.0374241 Jy ! Clean residual min=-0.001850 max=0.002175 Jy/beam ! Clean residual mean=0.000068 rms=0.000464 Jy/beam ! Combined flux in latest and established models = 0.396955 Jy ! Adding 23 model components to the UV plane model. ! The established model now contains 29 components and 0.396955 Jy ! Inverting map ! Added new window around map position (-2, 0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00162244 Jy ! Component: 100 - total flux cleaned = 0.00240505 Jy ! Component: 150 - total flux cleaned = 0.0030083 Jy ! Component: 200 - total flux cleaned = 0.00353876 Jy ! Total flux subtracted in 200 components = 0.00353876 Jy ! Clean residual min=-0.001479 max=0.001631 Jy/beam ! Clean residual mean=0.000046 rms=0.000425 Jy/beam ! Combined flux in latest and established models = 0.400494 Jy ! Adding 13 model components to the UV plane model. ! The established model now contains 40 components and 0.400494 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.377 mas, bmaj=8.045 mas, bpa=-14.92 degrees ! Estimated noise=533.316 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00116781 Jy ! Component: 100 - total flux cleaned = 0.00170695 Jy ! Component: 150 - total flux cleaned = 0.00187492 Jy ! Component: 200 - total flux cleaned = 0.00192206 Jy ! Total flux subtracted in 200 components = 0.00192206 Jy ! Clean residual min=-0.002167 max=0.003057 Jy/beam ! Clean residual mean=0.000135 rms=0.000771 Jy/beam ! Combined flux in latest and established models = 0.402416 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 = 3.50974e-05 Jy ! Component: 100 - total flux cleaned = 5.73708e-05 Jy ! Component: 150 - total flux cleaned = 6.47816e-05 Jy ! Component: 200 - total flux cleaned = 6.83009e-05 Jy ! Total flux subtracted in 200 components = 6.83009e-05 Jy ! Clean residual min=-0.002180 max=0.003061 Jy/beam ! Clean residual mean=0.000133 rms=0.000771 Jy/beam ! Combined flux in latest and established models = 0.402485 Jy ! Adding 25 model components to the UV plane model. ! The established model now contains 62 components and 0.402485 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.247 mas, bmaj=3.664 mas, bpa=-12.93 degrees ! Estimated noise=351.213 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000694231 Jy ! Component: 100 - total flux cleaned = -0.00066451 Jy ! Component: 150 - total flux cleaned = -0.000608326 Jy ! Component: 200 - total flux cleaned = -0.00055646 Jy ! Total flux subtracted in 200 components = -0.00055646 Jy ! Clean residual min=-0.001492 max=0.001660 Jy/beam ! Clean residual mean=0.000038 rms=0.000421 Jy/beam ! Combined flux in latest and established models = 0.401928 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 71 components and 0.401928 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.255147Jy sigma=0.001186 ! Fit after self-cal, rms=0.255160Jy sigma=0.001186 wmodel J2359+3850_S_2008_07_09_pus_map.mod ! Writing 71 model components to file: J2359+3850_S_2008_07_09_pus_map.mod wobs J2359+3850_S_2008_07_09_pus_uvs.fits ! Writing UV FITS file: J2359+3850_S_2008_07_09_pus_uvs.fits wwins J2359+3850_S_2008_07_09_pus_map.win ! wwins: Wrote 2 windows to J2359+3850_S_2008_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.000706342 Jy ! Component: 100 - total flux cleaned = 0.00113217 Jy ! Component: 150 - total flux cleaned = 0.00143243 Jy ! Component: 200 - total flux cleaned = 0.00167006 Jy ! Component: 250 - total flux cleaned = 0.0019526 Jy ! Component: 300 - total flux cleaned = 0.00215364 Jy ! Component: 350 - total flux cleaned = 0.00235087 Jy ! Component: 400 - total flux cleaned = 0.002544 Jy ! Component: 450 - total flux cleaned = 0.00273338 Jy ! Component: 500 - total flux cleaned = 0.00285004 Jy ! Component: 550 - total flux cleaned = 0.0029879 Jy ! Component: 600 - total flux cleaned = 0.00316882 Jy ! Component: 650 - total flux cleaned = 0.00332499 Jy ! Component: 700 - total flux cleaned = 0.00345733 Jy ! Component: 750 - total flux cleaned = 0.00358841 Jy ! Component: 800 - total flux cleaned = 0.00369665 Jy ! Component: 850 - total flux cleaned = 0.00378224 Jy ! Component: 900 - total flux cleaned = 0.00388836 Jy ! Component: 950 - total flux cleaned = 0.00397243 Jy ! Component: 1000 - total flux cleaned = 0.00409737 Jy ! Component: 1050 - total flux cleaned = 0.00415943 Jy ! Component: 1100 - total flux cleaned = 0.00428241 Jy ! Component: 1150 - total flux cleaned = 0.00438402 Jy ! Component: 1200 - total flux cleaned = 0.00448471 Jy ! Component: 1250 - total flux cleaned = 0.00460468 Jy ! Component: 1300 - total flux cleaned = 0.00470385 Jy ! Component: 1350 - total flux cleaned = 0.00476287 Jy ! Component: 1400 - total flux cleaned = 0.00486054 Jy ! Component: 1450 - total flux cleaned = 0.00491858 Jy ! Component: 1500 - total flux cleaned = 0.00499535 Jy ! Component: 1550 - total flux cleaned = 0.0049954 Jy ! Component: 1600 - total flux cleaned = 0.0050524 Jy ! Component: 1650 - total flux cleaned = 0.00505265 Jy ! Component: 1700 - total flux cleaned = 0.00522145 Jy ! Component: 1750 - total flux cleaned = 0.00520277 Jy ! Component: 1800 - total flux cleaned = 0.00529557 Jy ! Component: 1850 - total flux cleaned = 0.00535076 Jy ! Component: 1900 - total flux cleaned = 0.00533236 Jy ! Component: 1950 - total flux cleaned = 0.00536881 Jy ! Component: 2000 - total flux cleaned = 0.00542309 Jy ! Total flux subtracted in 2048 components = 0.00544102 Jy ! Clean residual min=-0.001127 max=0.001063 Jy/beam ! Clean residual mean=0.000012 rms=0.000342 Jy/beam ! Combined flux in latest and established models = 0.407369 Jy keep ! Adding 163 model components to the UV plane model. ! The established model now contains 234 components and 0.407369 Jy wmap "J2359+3850_S_2008_07_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.247 x 3.664 at -12.93 degrees (North through East) ! Clean map min=-0.0013195 max=0.3199 Jy/beam ! Writing clean map to FITS file: J2359+3850_S_2008_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.319904 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000342012 device "J2359+3850_S_2008_07_09_pus\_map.eps/VPS" ! Attempting to open device: 'J2359+3850_S_2008_07_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.320732 0.320732 0.641465 1.28293 2.56586 5.13172 10.2634 20.5269 41.0537 82.1075 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.001 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.319 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.247 x 3.664 at -12.93 degrees (North through East) ! Clean map min=-0.0013195 max=0.3199 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J2359+3850_S_2008_07_09_pus\_map.gif/vgif" ! Attempting to open device: 'J2359+3850_S_2008_07_09_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J2359+3850_S_2008_07_09_pus_uvs.log closed on Wed Oct 13 13:19:15 2010