! Started logfile: J1120+1420_S_2007_12_05_pus_uvs.log on Sun Sep 12 11:00:06 2010 obs J1120+1420_S_2007_12_05_pus_uva.fits ! Reading UV FITS file: J1120+1420_S_2007_12_05_pus_uva.fits ! AN table 1: 292 integrations on 136 of 136 possible baselines. ! Apparent sampling: 0.123892 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1120+1420 ! ! 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 1686 lines of history. ! ! Reading 19680 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 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 3140 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 3140 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 3140 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 3140 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=0.785808Jy sigma=0.008576 ! Fit after self-cal, rms=0.757472Jy sigma=0.008274 ! 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=7.203 mas, bmaj=9.707 mas, bpa=29.42 degrees ! Estimated noise=941.952 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.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.25576 Jy ! Component: 100 - total flux cleaned = 0.371815 Jy ! Total flux subtracted in 100 components = 0.371815 Jy ! Clean residual min=-0.093757 max=0.128686 Jy/beam ! Clean residual mean=-0.000036 rms=0.020884 Jy/beam ! Combined flux in latest and established models = 0.371815 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 19 components and 0.371815 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.447486Jy sigma=0.004552 ! Fit after self-cal, rms=0.444565Jy sigma=0.004526 ! Inverting map ! Added new window around map position (-65, 37). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.116582 Jy ! Component: 100 - total flux cleaned = 0.195615 Jy ! Total flux subtracted in 100 components = 0.195615 Jy ! Clean residual min=-0.075508 max=0.074126 Jy/beam ! Clean residual mean=-0.000029 rms=0.016411 Jy/beam ! Combined flux in latest and established models = 0.56743 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 39 components and 0.56743 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.388145Jy sigma=0.003868 ! Fit after self-cal, rms=0.381631Jy sigma=0.003817 ! 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=8.125 mas, bmaj=11.43 mas, bpa=27.77 degrees ! Estimated noise=483.874 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.148583 Jy ! Component: 100 - total flux cleaned = 0.256325 Jy ! Component: 150 - total flux cleaned = 0.34526 Jy ! Component: 200 - total flux cleaned = 0.420475 Jy ! Total flux subtracted in 200 components = 0.420475 Jy ! Clean residual min=-0.068318 max=0.087033 Jy/beam ! Clean residual mean=-0.000102 rms=0.022814 Jy/beam ! Combined flux in latest and established models = 0.987905 Jy ! Performing phase self-cal ! Adding 71 model components to the UV plane model. ! The established model now contains 106 components and 0.987905 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.285683Jy sigma=0.002752 ! Fit after self-cal, rms=0.275387Jy sigma=0.002667 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.75 FD 1.41 HH 1.00* HN 1.00* ! KB 1.00* KP 0.88 LA 0.71 MK 1.00* ! NL 1.08 NY 1.00* ON 1.00* OV 1.24 ! PT 0.87 SV 1.00* TS 1.00* WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.84 FD 1.29 HH 1.00* HN 1.00* ! KB 1.00* KP 0.90 LA 0.78 MK 1.00* ! NL 1.13 NY 1.00* ON 1.00* OV 1.21 ! PT 0.87 SV 1.00* TS 1.00* WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.18 HH 1.00* HN 1.00* ! KB 1.00* KP 0.85 LA 0.76 MK 1.00* ! NL 1.31 NY 1.00* ON 1.00* OV 1.27 ! PT 0.90 SV 1.00* TS 1.00* WZ 1.00* ! ZC 1.00* ! ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.87 FD 1.20 HH 1.00* HN 1.00* ! KB 1.00* KP 0.85 LA 0.77 MK 1.00* ! NL 1.23 NY 1.00* ON 1.00* OV 1.30 ! PT 0.78 SV 1.00* TS 1.00* WZ 1.00* ! ZC 1.00* ! ! ! Fit before self-cal, rms=0.275387Jy sigma=0.002667 ! Fit after self-cal, rms=0.220976Jy sigma=0.002038 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=8.122 mas, bmaj=12.23 mas, bpa=33.78 degrees ! Estimated noise=443.421 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 (-75.5, 41.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0128632 Jy ! Component: 100 - total flux cleaned = 0.015607 Jy ! Component: 150 - total flux cleaned = 0.0211232 Jy ! Component: 200 - total flux cleaned = 0.0236465 Jy ! Total flux subtracted in 200 components = 0.0236465 Jy ! Clean residual min=-0.028694 max=0.028262 Jy/beam ! Clean residual mean=-0.000000 rms=0.007446 Jy/beam ! Combined flux in latest and established models = 1.01155 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 151 components and 1.01155 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199067Jy sigma=0.001814 ! Fit after self-cal, rms=0.196157Jy sigma=0.001778 ! 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.196157Jy sigma=0.001778 ! Fit after self-cal, rms=0.188179Jy sigma=0.001727 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.014 mas, bmaj=12.3 mas, bpa=35.08 degrees ! Estimated noise=426.051 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00466602 Jy ! Component: 100 - total flux cleaned = 0.00391416 Jy ! Component: 150 - total flux cleaned = 0.00451637 Jy ! Component: 200 - total flux cleaned = 0.00260249 Jy ! Total flux subtracted in 200 components = 0.00260249 Jy ! Clean residual min=-0.022985 max=0.023866 Jy/beam ! Clean residual mean=-0.000001 rms=0.005446 Jy/beam ! Combined flux in latest and established models = 1.01415 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 178 components and 1.01415 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181270Jy sigma=0.001674 ! Fit after self-cal, rms=0.180199Jy sigma=0.001665 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 (-15, 7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0106948 Jy ! Component: 100 - total flux cleaned = 0.0126117 Jy ! Component: 150 - total flux cleaned = 0.0126138 Jy ! Component: 200 - total flux cleaned = 0.0126198 Jy ! Total flux subtracted in 200 components = 0.0126198 Jy ! Clean residual min=-0.017929 max=0.018419 Jy/beam ! Clean residual mean=0.000002 rms=0.004579 Jy/beam ! Combined flux in latest and established models = 1.02677 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 211 components and 1.02677 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174143Jy sigma=0.001613 ! Fit after self-cal, rms=0.173253Jy sigma=0.001606 ! 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.173253Jy sigma=0.001606 ! Fit after self-cal, rms=0.172491Jy sigma=0.001590 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.029 mas, bmaj=12.28 mas, bpa=35.05 degrees ! Estimated noise=420.218 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00597534 Jy ! Component: 100 - total flux cleaned = 0.0103311 Jy ! Component: 150 - total flux cleaned = 0.0111043 Jy ! Component: 200 - total flux cleaned = 0.00977047 Jy ! Total flux subtracted in 200 components = 0.00977047 Jy ! Clean residual min=-0.013953 max=0.014651 Jy/beam ! Clean residual mean=0.000001 rms=0.003676 Jy/beam ! Combined flux in latest and established models = 1.03654 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 247 components and 1.03654 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169731Jy sigma=0.001569 ! Fit after self-cal, rms=0.169412Jy sigma=0.001566 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.0032223 Jy ! Component: 100 - total flux cleaned = 0.00160662 Jy ! Component: 150 - total flux cleaned = 0.000108053 Jy ! Component: 200 - total flux cleaned = -0.00101717 Jy ! Total flux subtracted in 200 components = -0.00101717 Jy ! Clean residual min=-0.013137 max=0.013891 Jy/beam ! Clean residual mean=0.000002 rms=0.003398 Jy/beam ! Combined flux in latest and established models = 1.03553 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 268 components and 1.03553 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167860Jy sigma=0.001554 ! Fit after self-cal, rms=0.167594Jy sigma=0.001553 ! 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.167594Jy sigma=0.001553 ! Fit after self-cal, rms=0.163950Jy sigma=0.001550 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.368 mas, bmaj=13.11 mas, bpa=35.06 degrees ! Estimated noise=431.061 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00700093 Jy ! Component: 100 - total flux cleaned = 0.00785206 Jy ! Component: 150 - total flux cleaned = 0.00785835 Jy ! Component: 200 - total flux cleaned = 0.00786753 Jy ! Total flux subtracted in 200 components = 0.00786753 Jy ! Clean residual min=-0.011755 max=0.013237 Jy/beam ! Clean residual mean=0.000001 rms=0.003320 Jy/beam ! Combined flux in latest and established models = 1.04339 Jy selfcal ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 296 components and 1.04339 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162750Jy sigma=0.001541 ! Fit after self-cal, rms=0.162427Jy sigma=0.001539 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.00220359 Jy ! Component: 100 - total flux cleaned = 0.00172004 Jy ! Component: 150 - total flux cleaned = 0.00106637 Jy ! Component: 200 - total flux cleaned = 0.000657625 Jy ! Total flux subtracted in 200 components = 0.000657625 Jy ! Clean residual min=-0.011256 max=0.012662 Jy/beam ! Clean residual mean=0.000001 rms=0.003120 Jy/beam ! Combined flux in latest and established models = 1.04405 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 320 components and 1.04405 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161674Jy sigma=0.001533 ! Fit after self-cal, rms=0.161451Jy sigma=0.001532 ! 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=8.063 mas, bmaj=12.24 mas, bpa=34.77 degrees ! Estimated noise=418.083 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 8.063 x 12.24 at 34.77 degrees (North through East) ! Clean map min=-0.042767 max=0.29666 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=7.005 mas, bmaj=10.37 mas, bpa=36.76 degrees ! Estimated noise=869.972 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.228323 Jy ! Component: 100 - total flux cleaned = 0.375785 Jy ! Component: 150 - total flux cleaned = 0.478551 Jy ! Component: 200 - total flux cleaned = 0.547241 Jy ! Total flux subtracted in 200 components = 0.547241 Jy ! Clean residual min=-0.049071 max=0.057135 Jy/beam ! Clean residual mean=-0.000021 rms=0.012482 Jy/beam ! Combined flux in latest and established models = 0.547241 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 47 model components to the UV plane model. ! The established model now contains 47 components and 0.547241 Jy ! Inverting map and beam ! Estimated beam: bmin=8.063 mas, bmaj=12.24 mas, bpa=34.77 degrees ! Estimated noise=418.083 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.113998 Jy ! Component: 100 - total flux cleaned = 0.195034 Jy ! Component: 150 - total flux cleaned = 0.25826 Jy ! Component: 200 - total flux cleaned = 0.308972 Jy ! Total flux subtracted in 200 components = 0.308972 Jy ! Clean residual min=-0.033845 max=0.036248 Jy/beam ! Clean residual mean=-0.000040 rms=0.010194 Jy/beam ! Combined flux in latest and established models = 0.856213 Jy ! Adding 108 model components to the UV plane model. ! The established model now contains 151 components and 0.856213 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=9.368 mas, bmaj=13.11 mas, bpa=35.06 degrees ! Estimated noise=431.061 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0495437 Jy ! Component: 100 - total flux cleaned = 0.0802518 Jy ! Component: 150 - total flux cleaned = 0.102333 Jy ! Component: 200 - total flux cleaned = 0.121312 Jy ! Total flux subtracted in 200 components = 0.121312 Jy ! Clean residual min=-0.021394 max=0.020934 Jy/beam ! Clean residual mean=-0.000017 rms=0.005847 Jy/beam ! Combined flux in latest and established models = 0.977525 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.0150788 Jy ! Component: 100 - total flux cleaned = 0.0285707 Jy ! Component: 150 - total flux cleaned = 0.0394259 Jy ! Component: 200 - total flux cleaned = 0.0471905 Jy ! Total flux subtracted in 200 components = 0.0471905 Jy ! Clean residual min=-0.014560 max=0.012983 Jy/beam ! Clean residual mean=-0.000004 rms=0.003915 Jy/beam ! Combined flux in latest and established models = 1.02472 Jy ! Adding 142 model components to the UV plane model. ! The established model now contains 275 components and 1.02471 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=8.063 mas, bmaj=12.24 mas, bpa=34.77 degrees ! Estimated noise=418.083 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0053004 Jy ! Component: 100 - total flux cleaned = 0.00945126 Jy ! Component: 150 - total flux cleaned = 0.0128301 Jy ! Component: 200 - total flux cleaned = 0.0159178 Jy ! Total flux subtracted in 200 components = 0.0159178 Jy ! Clean residual min=-0.014116 max=0.011257 Jy/beam ! Clean residual mean=0.000000 rms=0.003216 Jy/beam ! Combined flux in latest and established models = 1.04063 Jy selfcal ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 313 components and 1.04063 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162695Jy sigma=0.001547 ! Fit after self-cal, rms=0.161605Jy sigma=0.001538 wmodel J1120+1420_S_2007_12_05_pus_map.mod ! Writing 313 model components to file: J1120+1420_S_2007_12_05_pus_map.mod wobs J1120+1420_S_2007_12_05_pus_uvs.fits ! Writing UV FITS file: J1120+1420_S_2007_12_05_pus_uvs.fits wwins J1120+1420_S_2007_12_05_pus_map.win ! wwins: Wrote 4 windows to J1120+1420_S_2007_12_05_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.00356287 Jy ! Component: 100 - total flux cleaned = -0.00475349 Jy ! Component: 150 - total flux cleaned = -0.00493715 Jy ! Component: 200 - total flux cleaned = -0.00546961 Jy ! Component: 250 - total flux cleaned = -0.0053058 Jy ! Component: 300 - total flux cleaned = -0.00498795 Jy ! Component: 350 - total flux cleaned = -0.00498789 Jy ! Component: 400 - total flux cleaned = -0.00469275 Jy ! Component: 450 - total flux cleaned = -0.00512482 Jy ! Component: 500 - total flux cleaned = -0.00540317 Jy ! Component: 550 - total flux cleaned = -0.00540672 Jy ! Component: 600 - total flux cleaned = -0.00617867 Jy ! Component: 650 - total flux cleaned = -0.00655333 Jy ! Component: 700 - total flux cleaned = -0.00715914 Jy ! Component: 750 - total flux cleaned = -0.00751411 Jy ! Component: 800 - total flux cleaned = -0.00808438 Jy ! Component: 850 - total flux cleaned = -0.00863501 Jy ! Component: 900 - total flux cleaned = -0.00852841 Jy ! Component: 950 - total flux cleaned = -0.00915843 Jy ! Component: 1000 - total flux cleaned = -0.00915447 Jy ! Component: 1050 - total flux cleaned = -0.00935004 Jy ! Component: 1100 - total flux cleaned = -0.00915474 Jy ! Component: 1150 - total flux cleaned = -0.00915582 Jy ! Component: 1200 - total flux cleaned = -0.00915689 Jy ! Component: 1250 - total flux cleaned = -0.00897692 Jy ! Component: 1300 - total flux cleaned = -0.00906551 Jy ! Component: 1350 - total flux cleaned = -0.00906363 Jy ! Component: 1400 - total flux cleaned = -0.00863969 Jy ! Component: 1450 - total flux cleaned = -0.00872297 Jy ! Component: 1500 - total flux cleaned = -0.00872118 Jy ! Component: 1550 - total flux cleaned = -0.00856289 Jy ! Component: 1600 - total flux cleaned = -0.00848543 Jy ! Component: 1650 - total flux cleaned = -0.00848675 Jy ! Component: 1700 - total flux cleaned = -0.00871353 Jy ! Component: 1750 - total flux cleaned = -0.00849176 Jy ! Component: 1800 - total flux cleaned = -0.00870797 Jy ! Component: 1850 - total flux cleaned = -0.008707 Jy ! Component: 1900 - total flux cleaned = -0.00877601 Jy ! Component: 1950 - total flux cleaned = -0.00863858 Jy ! Component: 2000 - total flux cleaned = -0.00884232 Jy ! Component: 2050 - total flux cleaned = -0.00884139 Jy ! Component: 2100 - total flux cleaned = -0.00864599 Jy ! Component: 2150 - total flux cleaned = -0.00877478 Jy ! Component: 2200 - total flux cleaned = -0.00877412 Jy ! Component: 2250 - total flux cleaned = -0.00877351 Jy ! Component: 2300 - total flux cleaned = -0.0088953 Jy ! Component: 2350 - total flux cleaned = -0.00877483 Jy ! Component: 2400 - total flux cleaned = -0.00865641 Jy ! Component: 2450 - total flux cleaned = -0.00871505 Jy ! Component: 2500 - total flux cleaned = -0.00848358 Jy ! Component: 2550 - total flux cleaned = -0.00853948 Jy ! Component: 2600 - total flux cleaned = -0.00831503 Jy ! Component: 2650 - total flux cleaned = -0.0081497 Jy ! Component: 2700 - total flux cleaned = -0.0080951 Jy ! Component: 2750 - total flux cleaned = -0.00793314 Jy ! Component: 2800 - total flux cleaned = -0.00777444 Jy ! Component: 2850 - total flux cleaned = -0.00756517 Jy ! Component: 2900 - total flux cleaned = -0.00756527 Jy ! Component: 2950 - total flux cleaned = -0.00731002 Jy ! Component: 3000 - total flux cleaned = -0.00741097 Jy ! Component: 3050 - total flux cleaned = -0.00721113 Jy ! Component: 3100 - total flux cleaned = -0.00725973 Jy ! Component: 3150 - total flux cleaned = -0.0069201 Jy ! Component: 3200 - total flux cleaned = -0.00687217 Jy ! Component: 3250 - total flux cleaned = -0.00682447 Jy ! Component: 3300 - total flux cleaned = -0.00649687 Jy ! Component: 3350 - total flux cleaned = -0.00631138 Jy ! Component: 3400 - total flux cleaned = -0.00617381 Jy ! Component: 3450 - total flux cleaned = -0.00599283 Jy ! Component: 3500 - total flux cleaned = -0.00576881 Jy ! Component: 3550 - total flux cleaned = -0.0057244 Jy ! Component: 3600 - total flux cleaned = -0.0054182 Jy ! Component: 3650 - total flux cleaned = -0.00528863 Jy ! Component: 3700 - total flux cleaned = -0.00507474 Jy ! Component: 3750 - total flux cleaned = -0.00503218 Jy ! Component: 3800 - total flux cleaned = -0.00486463 Jy ! Component: 3850 - total flux cleaned = -0.00469877 Jy ! Component: 3900 - total flux cleaned = -0.00465766 Jy ! Component: 3950 - total flux cleaned = -0.00445435 Jy ! Component: 4000 - total flux cleaned = -0.00429402 Jy ! Component: 4050 - total flux cleaned = -0.00417455 Jy ! Component: 4100 - total flux cleaned = -0.00413519 Jy ! Component: 4150 - total flux cleaned = -0.00405687 Jy ! Component: 4200 - total flux cleaned = -0.00374772 Jy ! Component: 4250 - total flux cleaned = -0.00386257 Jy ! Component: 4300 - total flux cleaned = -0.00363475 Jy ! Component: 4350 - total flux cleaned = -0.00355923 Jy ! Component: 4400 - total flux cleaned = -0.0034479 Jy ! Component: 4450 - total flux cleaned = -0.00341106 Jy ! Component: 4500 - total flux cleaned = -0.00337476 Jy ! Component: 4550 - total flux cleaned = -0.00323034 Jy ! Component: 4600 - total flux cleaned = -0.00312336 Jy ! Component: 4650 - total flux cleaned = -0.00312339 Jy ! Component: 4700 - total flux cleaned = -0.00298248 Jy ! Component: 4750 - total flux cleaned = -0.00305208 Jy ! Component: 4800 - total flux cleaned = -0.00291394 Jy ! Component: 4850 - total flux cleaned = -0.00284524 Jy ! Component: 4900 - total flux cleaned = -0.00277734 Jy ! Component: 4950 - total flux cleaned = -0.00271 Jy ! Component: 5000 - total flux cleaned = -0.00260969 Jy ! Component: 5050 - total flux cleaned = -0.00254363 Jy ! Component: 5100 - total flux cleaned = -0.00244526 Jy ! Component: 5150 - total flux cleaned = -0.00231464 Jy ! Component: 5200 - total flux cleaned = -0.00224999 Jy ! Component: 5250 - total flux cleaned = -0.0021212 Jy ! Component: 5300 - total flux cleaned = -0.00208935 Jy ! Component: 5350 - total flux cleaned = -0.0020578 Jy ! Component: 5400 - total flux cleaned = -0.00196331 Jy ! Component: 5450 - total flux cleaned = -0.00190101 Jy ! Component: 5500 - total flux cleaned = -0.00171531 Jy ! Component: 5550 - total flux cleaned = -0.00180772 Jy ! Component: 5600 - total flux cleaned = -0.00171655 Jy ! Component: 5650 - total flux cleaned = -0.0017469 Jy ! Component: 5700 - total flux cleaned = -0.00174691 Jy ! Component: 5750 - total flux cleaned = -0.00165735 Jy ! Component: 5800 - total flux cleaned = -0.00159819 Jy ! Component: 5850 - total flux cleaned = -0.0015981 Jy ! Component: 5900 - total flux cleaned = -0.00148113 Jy ! Component: 5950 - total flux cleaned = -0.00153924 Jy ! Component: 6000 - total flux cleaned = -0.00151021 Jy ! Component: 6050 - total flux cleaned = -0.00153883 Jy ! Component: 6100 - total flux cleaned = -0.001425 Jy ! Component: 6150 - total flux cleaned = -0.00145342 Jy ! Component: 6200 - total flux cleaned = -0.00145335 Jy ! Component: 6250 - total flux cleaned = -0.00145324 Jy ! Component: 6300 - total flux cleaned = -0.00134183 Jy ! Component: 6350 - total flux cleaned = -0.00134206 Jy ! Component: 6400 - total flux cleaned = -0.00139706 Jy ! Component: 6450 - total flux cleaned = -0.00136985 Jy ! Component: 6500 - total flux cleaned = -0.00145175 Jy ! Component: 6550 - total flux cleaned = -0.0014247 Jy ! Component: 6600 - total flux cleaned = -0.00147877 Jy ! Component: 6650 - total flux cleaned = -0.00158649 Jy ! Component: 6700 - total flux cleaned = -0.00161353 Jy ! Component: 6750 - total flux cleaned = -0.00169356 Jy ! Component: 6800 - total flux cleaned = -0.00169369 Jy ! Component: 6850 - total flux cleaned = -0.00179949 Jy ! Component: 6900 - total flux cleaned = -0.0018259 Jy ! Component: 6950 - total flux cleaned = -0.0019834 Jy ! Component: 7000 - total flux cleaned = -0.00198329 Jy ! Component: 7050 - total flux cleaned = -0.00200925 Jy ! Component: 7100 - total flux cleaned = -0.00219081 Jy ! Component: 7150 - total flux cleaned = -0.00221683 Jy ! Component: 7200 - total flux cleaned = -0.00224249 Jy ! Component: 7250 - total flux cleaned = -0.00237073 Jy ! Component: 7300 - total flux cleaned = -0.00239628 Jy ! Component: 7350 - total flux cleaned = -0.00244714 Jy ! Component: 7400 - total flux cleaned = -0.00249778 Jy ! Component: 7450 - total flux cleaned = -0.0025987 Jy ! Component: 7500 - total flux cleaned = -0.00257357 Jy ! Component: 7550 - total flux cleaned = -0.00267408 Jy ! Component: 7600 - total flux cleaned = -0.00272408 Jy ! Component: 7650 - total flux cleaned = -0.002749 Jy ! Component: 7700 - total flux cleaned = -0.00284847 Jy ! Component: 7750 - total flux cleaned = -0.00282366 Jy ! Component: 7800 - total flux cleaned = -0.00284855 Jy ! Component: 7850 - total flux cleaned = -0.00299652 Jy ! Component: 7900 - total flux cleaned = -0.00299644 Jy ! Component: 7950 - total flux cleaned = -0.00302084 Jy ! Component: 8000 - total flux cleaned = -0.00302095 Jy ! Component: 8050 - total flux cleaned = -0.00299653 Jy ! Component: 8100 - total flux cleaned = -0.00302083 Jy ! Component: 8150 - total flux cleaned = -0.00302093 Jy ! Total flux subtracted in 8192 components = -0.00304511 Jy ! Clean residual min=-0.001294 max=0.001352 Jy/beam ! Clean residual mean=-0.000001 rms=0.000432 Jy/beam ! Combined flux in latest and established models = 1.03759 Jy keep ! Adding 678 model components to the UV plane model. ! The established model now contains 990 components and 1.03759 Jy wmap "J1120+1420_S_2007_12_05_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 8.063 x 12.24 at 34.77 degrees (North through East) ! Clean map min=-0.032347 max=0.28588 Jy/beam ! Writing clean map to FITS file: J1120+1420_S_2007_12_05_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.285882 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000437052 device "J1120+1420_S_2007_12_05_pus\_map.eps/VPS" ! Attempting to open device: 'J1120+1420_S_2007_12_05_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.458635 0.458635 0.917269 1.83454 3.66908 7.33815 14.6763 29.3526 58.7052 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0013 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.285 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 8.063 x 12.24 at 34.77 degrees (North through East) ! Clean map min=-0.032347 max=0.28588 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1120+1420_S_2007_12_05_pus\_map.gif/vgif" ! Attempting to open device: 'J1120+1420_S_2007_12_05_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J1120+1420_S_2007_12_05_pus_uvs.log closed on Sun Sep 12 11:01:44 2010