! Started logfile: J1505+0326_S_2010_10_13_pus_uvs.log on Mon Feb 28 00:02:43 2011 obs J1505+0326_S_2010_10_13_pus_uva.fits ! Reading UV FITS file: J1505+0326_S_2010_10_13_pus_uva.fits ! AN table 1: 15 integrations on 105 of 105 possible baselines. ! AN table 2: 113 integrations on 105 of 105 possible baselines. ! AN table 3: 34 integrations on 105 of 105 possible baselines. ! Apparent sampling: 0.1592 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1505+0326 ! ! 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.23696e+09 8e+06 1 8e+06 ! 02 2 2.26696e+09 8e+06 1 8e+06 ! 03 3 2.35696e+09 8e+06 1 8e+06 ! 04 4 2.37696e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 3545 lines of history. ! ! Reading 10832 visibilities. select I ! Polarization I is unavailable. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 512 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J1505+0326_S_2010_10_13_pus] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 512x512 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 104 telescope corrections were flagged in sub-array 1. ! A total of 1226 telescope corrections were flagged in sub-array 2. ! A total of 279 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 104 telescope corrections were flagged in sub-array 1. ! A total of 1226 telescope corrections were flagged in sub-array 2. ! A total of 279 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 104 telescope corrections were flagged in sub-array 1. ! A total of 1226 telescope corrections were flagged in sub-array 2. ! A total of 279 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 104 telescope corrections were flagged in sub-array 1. ! A total of 1226 telescope corrections were flagged in sub-array 2. ! A total of 279 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.508886Jy sigma=0.005036 ! Fit after self-cal, rms=0.476508Jy sigma=0.004774 ! 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.014 mas, bmaj=3.298 mas, bpa=19.1 degrees ! Estimated noise=1052.64 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.43358 Jy ! Component: 100 - total flux cleaned = 0.533792 Jy ! Total flux subtracted in 100 components = 0.533792 Jy ! Clean residual min=-0.013742 max=0.034319 Jy/beam ! Clean residual mean=-0.000006 rms=0.003707 Jy/beam ! Combined flux in latest and established models = 0.533792 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.533792 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180853Jy sigma=0.001596 ! Fit after self-cal, rms=0.180851Jy sigma=0.001596 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0350745 Jy ! Component: 100 - total flux cleaned = 0.0521327 Jy ! Total flux subtracted in 100 components = 0.0521327 Jy ! Clean residual min=-0.010982 max=0.015396 Jy/beam ! Clean residual mean=-0.000223 rms=0.002697 Jy/beam ! Combined flux in latest and established models = 0.585925 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 13 components and 0.585925 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.175153Jy sigma=0.001503 ! Fit after self-cal, rms=0.175141Jy sigma=0.001503 ! 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.489 mas, bmaj=4.16 mas, bpa=12.56 degrees ! Estimated noise=692.535 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 (-5.5, -4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0124592 Jy ! Component: 100 - total flux cleaned = 0.015022 Jy ! Component: 150 - total flux cleaned = 0.0152374 Jy ! Component: 200 - total flux cleaned = 0.0150376 Jy ! Total flux subtracted in 200 components = 0.0150376 Jy ! Clean residual min=-0.009107 max=0.009533 Jy/beam ! Clean residual mean=-0.002091 rms=0.001958 Jy/beam ! Combined flux in latest and established models = 0.600963 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 29 components and 0.600963 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174449Jy sigma=0.001487 ! Fit after self-cal, rms=0.174428Jy sigma=0.001486 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.01 HH 1.00* HN 1.04 ! KB 1.00* KK 1.00* KP 0.94 LA 1.03 ! MK 0.73 NL 1.00 NY 1.00* OV 0.98 ! PT 0.98 SC 0.98 TS 1.00* ! ! A total of 69 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 0.99 HH 1.02 HN 0.99 ! KB 0.88 KK 1.00* KP 0.96 LA 1.00* ! MK 1.00* NL 0.98 NY 1.02 OV 1.00* ! PT 0.96 SC 0.97 TS 0.96 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.97 FD 1.03 HH 1.00* HN 1.40 ! KB 1.42 KK 1.01 KP 0.96 LA 1.01 ! MK 0.97 NL 0.98 NY 1.00* OV 0.94 ! PT 0.99 SC 1.00* TS 0.92 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 0.98 HH 1.00* HN 1.04 ! KB 1.00* KK 1.00* KP 0.98 LA 0.99 ! MK 0.99 NL 1.01 NY 1.00* OV 1.03 ! PT 1.02 SC 0.99 TS 1.00* ! ! A total of 69 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 FD 0.98 HH 1.01 HN 1.00 ! KB 0.91 KK 1.00* KP 0.98 LA 1.00* ! MK 1.00* NL 1.01 NY 1.06 OV 1.00* ! PT 1.01 SC 0.97 TS 0.95 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.06 FD 0.99 HH 1.00* HN 1.47 ! KB 1.55 KK 0.97 KP 0.99 LA 1.00 ! MK 0.94 NL 0.99 NY 1.00* OV 1.02 ! PT 1.00 SC 1.00* TS 0.92 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.99 HH 1.00* HN 1.04 ! KB 1.00* KK 1.00* KP 1.02 LA 1.01 ! MK 1.01 NL 1.00 NY 1.00* OV 1.01 ! PT 0.99 SC 0.96 TS 1.00* ! ! A total of 69 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 0.98 HH 1.03 HN 1.00 ! KB 0.89 KK 1.00* KP 1.00 LA 1.00* ! MK 1.00* NL 1.00 NY 1.02 OV 1.00* ! PT 0.98 SC 0.97 TS 1.02 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.00 HH 1.00* HN 1.34 ! KB 1.45 KK 0.96 KP 1.02 LA 0.97 ! MK 0.95 NL 1.00 NY 1.00* OV 1.00 ! PT 0.99 SC 1.00* TS 0.92 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.79 FD 1.01 HH 1.00* HN 1.05 ! KB 1.00* KK 1.00* KP 0.97 LA 0.99 ! MK 1.01 NL 1.00 NY 1.00* OV 0.99 ! PT 1.01 SC 1.03 TS 1.00* ! ! A total of 69 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.89 FD 0.96 HH 0.99 HN 1.02 ! KB 0.99 KK 1.00* KP 0.97 LA 1.00* ! MK 1.00* NL 0.97 NY 1.04 OV 1.00* ! PT 0.99 SC 1.03 TS 1.14 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 0.97 HH 1.00* HN 1.43 ! KB 1.54 KK 0.97 KP 0.97 LA 1.02 ! MK 0.95 NL 0.97 NY 1.00* OV 0.98 ! PT 1.00 SC 1.00* TS 1.11 ! ! ! Fit before self-cal, rms=0.174428Jy sigma=0.001486 ! Fit after self-cal, rms=0.167681Jy sigma=0.001329 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.544 mas, bmaj=4.273 mas, bpa=11.55 degrees ! Estimated noise=693.934 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.00175797 Jy ! Component: 100 - total flux cleaned = -0.00174687 Jy ! Component: 150 - total flux cleaned = -0.0018368 Jy ! Component: 200 - total flux cleaned = -0.00202372 Jy ! Total flux subtracted in 200 components = -0.00202372 Jy ! Clean residual min=-0.004439 max=0.003492 Jy/beam ! Clean residual mean=-0.000338 rms=0.001124 Jy/beam ! Combined flux in latest and established models = 0.598939 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 42 components and 0.598939 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167472Jy sigma=0.001326 ! Fit after self-cal, rms=0.167577Jy sigma=0.001326 ! 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.167577Jy sigma=0.001326 ! Fit after self-cal, rms=0.167841Jy sigma=0.001318 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.542 mas, bmaj=4.286 mas, bpa=11.93 degrees ! Estimated noise=692.974 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000195597 Jy ! Component: 100 - total flux cleaned = 8.54991e-05 Jy ! Component: 150 - total flux cleaned = 0.000162907 Jy ! Component: 200 - total flux cleaned = 0.000235228 Jy ! Total flux subtracted in 200 components = 0.000235228 Jy ! Clean residual min=-0.003883 max=0.003053 Jy/beam ! Clean residual mean=-0.000281 rms=0.000899 Jy/beam ! Combined flux in latest and established models = 0.599174 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 51 components and 0.599174 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167717Jy sigma=0.001316 ! Fit after self-cal, rms=0.167815Jy sigma=0.001316 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.00022743 Jy ! Component: 100 - total flux cleaned = 0.0003469 Jy ! Component: 150 - total flux cleaned = 0.00034644 Jy ! Component: 200 - total flux cleaned = 0.000347047 Jy ! Total flux subtracted in 200 components = 0.000347047 Jy ! Clean residual min=-0.003871 max=0.003049 Jy/beam ! Clean residual mean=-0.000288 rms=0.000887 Jy/beam ! Combined flux in latest and established models = 0.599521 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 59 components and 0.599521 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167766Jy sigma=0.001316 ! Fit after self-cal, rms=0.167761Jy sigma=0.001316 ! 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.167761Jy sigma=0.001316 ! Fit after self-cal, rms=0.169167Jy sigma=0.001311 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.546 mas, bmaj=4.29 mas, bpa=11.88 degrees ! Estimated noise=693.01 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000288936 Jy ! Component: 100 - total flux cleaned = 0.000223131 Jy ! Component: 150 - total flux cleaned = 0.00027957 Jy ! Component: 200 - total flux cleaned = 0.000278756 Jy ! Total flux subtracted in 200 components = 0.000278756 Jy ! Clean residual min=-0.003198 max=0.003379 Jy/beam ! Clean residual mean=-0.000094 rms=0.000809 Jy/beam ! Combined flux in latest and established models = 0.5998 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 63 components and 0.5998 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169103Jy sigma=0.001311 ! Fit after self-cal, rms=0.169143Jy sigma=0.001311 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.000166426 Jy ! Component: 100 - total flux cleaned = 0.000211934 Jy ! Component: 150 - total flux cleaned = 0.000212789 Jy ! Component: 200 - total flux cleaned = 0.000212784 Jy ! Total flux subtracted in 200 components = 0.000212784 Jy ! Clean residual min=-0.003215 max=0.003379 Jy/beam ! Clean residual mean=-0.000100 rms=0.000804 Jy/beam ! Combined flux in latest and established models = 0.600013 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 65 components and 0.600013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169106Jy sigma=0.001310 ! Fit after self-cal, rms=0.169104Jy sigma=0.001310 ! 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.169104Jy sigma=0.001310 ! Fit after self-cal, rms=0.169427Jy sigma=0.001310 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.049 mas, bmaj=10.77 mas, bpa=-9.652 degrees ! Estimated noise=860.591 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00019071 Jy ! Component: 100 - total flux cleaned = 0.000241083 Jy ! Component: 150 - total flux cleaned = 0.000253241 Jy ! Component: 200 - total flux cleaned = 0.000265017 Jy ! Total flux subtracted in 200 components = 0.000265017 Jy ! Clean residual min=-0.003361 max=0.003297 Jy/beam ! Clean residual mean=-0.000118 rms=0.000957 Jy/beam ! Combined flux in latest and established models = 0.600278 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 72 components and 0.600278 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169420Jy sigma=0.001310 ! Fit after self-cal, rms=0.169424Jy sigma=0.001310 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 = 7.63516e-05 Jy ! Component: 100 - total flux cleaned = 8.80865e-05 Jy ! Component: 150 - total flux cleaned = 9.89179e-05 Jy ! Component: 200 - total flux cleaned = 0.000119304 Jy ! Total flux subtracted in 200 components = 0.000119304 Jy ! Clean residual min=-0.003348 max=0.003295 Jy/beam ! Clean residual mean=-0.000123 rms=0.000957 Jy/beam ! Combined flux in latest and established models = 0.600397 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 79 components and 0.600397 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169419Jy sigma=0.001310 ! Fit after self-cal, rms=0.169419Jy sigma=0.001310 ! 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.549 mas, bmaj=4.292 mas, bpa=11.84 degrees ! Estimated noise=692.957 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.549 x 4.292 at 11.84 degrees (North through East) ! Clean map min=-0.0029358 max=0.55606 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.169419Jy sigma=0.001310 ! Fit after self-cal, rms=0.171434Jy sigma=0.001177 ! Inverting map and beam ! Estimated beam: bmin=2.549 mas, bmaj=4.291 mas, bpa=11.85 degrees ! Estimated noise=692.555 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00051506 Jy ! Component: 100 - total flux cleaned = -0.000598081 Jy ! Component: 150 - total flux cleaned = -0.000598335 Jy ! Component: 200 - total flux cleaned = -0.000597805 Jy ! Total flux subtracted in 200 components = -0.000597805 Jy ! Clean residual min=-0.003082 max=0.003429 Jy/beam ! Clean residual mean=-0.000039 rms=0.000786 Jy/beam ! Combined flux in latest and established models = 0.599799 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 83 components and 0.599799 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171402Jy sigma=0.001177 ! Fit after self-cal, rms=0.171401Jy sigma=0.001177 ! Inverting map ! Added new window around map position (-59, -12). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00204918 Jy ! Component: 100 - total flux cleaned = 0.00203324 Jy ! Component: 150 - total flux cleaned = 0.00207972 Jy ! Component: 200 - total flux cleaned = 0.00225052 Jy ! Total flux subtracted in 200 components = 0.00225052 Jy ! Clean residual min=-0.003246 max=0.003236 Jy/beam ! Clean residual mean=-0.000080 rms=0.000768 Jy/beam ! Combined flux in latest and established models = 0.60205 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 96 components and 0.60205 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.171307Jy sigma=0.001176 ! Fit after self-cal, rms=0.171295Jy sigma=0.001176 ! 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.171295Jy sigma=0.001176 ! Fit after self-cal, rms=0.177814Jy sigma=0.001176 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.057 mas, bmaj=3.417 mas, bpa=18.41 degrees ! Estimated noise=1054.37 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.432211 Jy ! Component: 100 - total flux cleaned = 0.530992 Jy ! Component: 150 - total flux cleaned = 0.561679 Jy ! Component: 200 - total flux cleaned = 0.575371 Jy ! Total flux subtracted in 200 components = 0.575371 Jy ! Clean residual min=-0.005290 max=0.006881 Jy/beam ! Clean residual mean=0.000100 rms=0.001349 Jy/beam ! Combined flux in latest and established models = 0.575371 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 20 model components to the UV plane model. ! The established model now contains 20 components and 0.575371 Jy ! Inverting map and beam ! Estimated beam: bmin=2.552 mas, bmaj=4.294 mas, bpa=11.84 degrees ! Estimated noise=692.505 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.0109058 Jy ! Component: 100 - total flux cleaned = 0.0167949 Jy ! Component: 150 - total flux cleaned = 0.0207094 Jy ! Component: 200 - total flux cleaned = 0.0235627 Jy ! Total flux subtracted in 200 components = 0.0235627 Jy ! Clean residual min=-0.003251 max=0.003476 Jy/beam ! Clean residual mean=0.000009 rms=0.000803 Jy/beam ! Combined flux in latest and established models = 0.598934 Jy ! Adding 28 model components to the UV plane model. ! The established model now contains 42 components and 0.598934 Jy ! Inverting map ! Added new window around map position (-43.5, -34.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00176781 Jy ! Component: 100 - total flux cleaned = 0.00258054 Jy ! Component: 150 - total flux cleaned = 0.00322378 Jy ! Component: 200 - total flux cleaned = 0.00380956 Jy ! Total flux subtracted in 200 components = 0.00380956 Jy ! Clean residual min=-0.003152 max=0.002765 Jy/beam ! Clean residual mean=-0.000049 rms=0.000770 Jy/beam ! Combined flux in latest and established models = 0.602743 Jy ! Adding 26 model components to the UV plane model. ! The established model now contains 56 components and 0.602743 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.051 mas, bmaj=10.77 mas, bpa=-9.604 degrees ! Estimated noise=859.378 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00130426 Jy ! Component: 100 - total flux cleaned = 0.001861 Jy ! Component: 150 - total flux cleaned = 0.00222166 Jy ! Component: 200 - total flux cleaned = 0.00246869 Jy ! Total flux subtracted in 200 components = 0.00246869 Jy ! Clean residual min=-0.003348 max=0.002589 Jy/beam ! Clean residual mean=-0.000210 rms=0.000845 Jy/beam ! Combined flux in latest and established models = 0.605212 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.000214926 Jy ! Component: 100 - total flux cleaned = 0.000262916 Jy ! Component: 150 - total flux cleaned = 0.000308004 Jy ! Component: 200 - total flux cleaned = 0.000329904 Jy ! Total flux subtracted in 200 components = 0.000329904 Jy ! Clean residual min=-0.003317 max=0.002573 Jy/beam ! Clean residual mean=-0.000222 rms=0.000836 Jy/beam ! Combined flux in latest and established models = 0.605542 Jy ! Adding 44 model components to the UV plane model. ! The established model now contains 95 components and 0.605542 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.552 mas, bmaj=4.294 mas, bpa=11.84 degrees ! Estimated noise=692.505 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00114836 Jy ! Component: 100 - total flux cleaned = -0.00148185 Jy ! Component: 150 - total flux cleaned = -0.00147932 Jy ! Component: 200 - total flux cleaned = -0.00142241 Jy ! Total flux subtracted in 200 components = -0.00142241 Jy ! Clean residual min=-0.003058 max=0.002794 Jy/beam ! Clean residual mean=-0.000070 rms=0.000754 Jy/beam ! Combined flux in latest and established models = 0.604119 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 107 components and 0.60412 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177824Jy sigma=0.001176 ! Fit after self-cal, rms=0.177818Jy sigma=0.001176 wmodel J1505+0326_S_2010_10_13_pus_map.mod ! Writing 107 model components to file: J1505+0326_S_2010_10_13_pus_map.mod wobs J1505+0326_S_2010_10_13_pus_uvs.fits ! Writing UV FITS file: J1505+0326_S_2010_10_13_pus_uvs.fits wwins J1505+0326_S_2010_10_13_pus_map.win ! wwins: Wrote 4 windows to J1505+0326_S_2010_10_13_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.00112056 Jy ! Component: 100 - total flux cleaned = -0.00157868 Jy ! Component: 150 - total flux cleaned = -0.00177218 Jy ! Component: 200 - total flux cleaned = -0.00195912 Jy ! Component: 250 - total flux cleaned = -0.00218638 Jy ! Component: 300 - total flux cleaned = -0.00240859 Jy ! Component: 350 - total flux cleaned = -0.00266977 Jy ! Component: 400 - total flux cleaned = -0.0028411 Jy ! Component: 450 - total flux cleaned = -0.00305121 Jy ! Component: 500 - total flux cleaned = -0.00313407 Jy ! Component: 550 - total flux cleaned = -0.0032159 Jy ! Component: 600 - total flux cleaned = -0.00337695 Jy ! Component: 650 - total flux cleaned = -0.00349613 Jy ! Component: 700 - total flux cleaned = -0.00361347 Jy ! Component: 750 - total flux cleaned = -0.00365183 Jy ! Component: 800 - total flux cleaned = -0.00365181 Jy ! Component: 850 - total flux cleaned = -0.00372764 Jy ! Component: 900 - total flux cleaned = -0.00369028 Jy ! Component: 950 - total flux cleaned = -0.00372745 Jy ! Component: 1000 - total flux cleaned = -0.0037274 Jy ! Component: 1050 - total flux cleaned = -0.00372741 Jy ! Component: 1100 - total flux cleaned = -0.00369169 Jy ! Component: 1150 - total flux cleaned = -0.00379821 Jy ! Component: 1200 - total flux cleaned = -0.00376237 Jy ! Component: 1250 - total flux cleaned = -0.00372731 Jy ! Component: 1300 - total flux cleaned = -0.00372745 Jy ! Component: 1350 - total flux cleaned = -0.00372735 Jy ! Component: 1400 - total flux cleaned = -0.00369336 Jy ! Component: 1450 - total flux cleaned = -0.00376076 Jy ! Component: 1500 - total flux cleaned = -0.00372738 Jy ! Component: 1550 - total flux cleaned = -0.00376024 Jy ! Component: 1600 - total flux cleaned = -0.00376004 Jy ! Component: 1650 - total flux cleaned = -0.0037925 Jy ! Component: 1700 - total flux cleaned = -0.0037924 Jy ! Component: 1750 - total flux cleaned = -0.00385637 Jy ! Component: 1800 - total flux cleaned = -0.00385608 Jy ! Component: 1850 - total flux cleaned = -0.00385606 Jy ! Component: 1900 - total flux cleaned = -0.00385599 Jy ! Component: 1950 - total flux cleaned = -0.00385615 Jy ! Component: 2000 - total flux cleaned = -0.00391789 Jy ! Total flux subtracted in 2048 components = -0.00388714 Jy ! Clean residual min=-0.001881 max=0.001627 Jy/beam ! Clean residual mean=-0.000001 rms=0.000564 Jy/beam ! Combined flux in latest and established models = 0.600232 Jy keep ! Adding 138 model components to the UV plane model. ! The established model now contains 245 components and 0.600233 Jy ![Exited script file: muppet] wmap "J1505+0326_S_2010_10_13_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.552 x 4.294 at 11.84 degrees (North through East) ! Clean map min=-0.0027103 max=0.55559 Jy/beam ! Writing clean map to FITS file: J1505+0326_S_2010_10_13_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.555591 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000565137 device "J1505+0326_S_2010_10_13_pus\_map.eps/VPS" ! Attempting to open device: 'J1505+0326_S_2010_10_13_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.305155 0.305155 0.61031 1.22062 2.44124 4.88248 9.76496 19.5299 39.0598 78.1197 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0016 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.555 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.552 x 4.294 at 11.84 degrees (North through East) ! Clean map min=-0.0027103 max=0.55559 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1505+0326_S_2010_10_13_pus\_map.gif/vgif" ! Attempting to open device: 'J1505+0326_S_2010_10_13_pus_map.gif/vgif' mapplot cln