! Started logfile: J1153+8058_S_2001_03_12_pus_uvs.log on Wed Jan 2 18:52:48 2013 obs J1153+8058_S_2001_03_12_pus_uva.fits ! Reading UV FITS file: J1153+8058_S_2001_03_12_pus_uva.fits ! AN table 1: 446 integrations on 136 of 136 possible baselines. ! AN table 2: 134 integrations on 136 of 136 possible baselines. ! AN table 3: 35 integrations on 136 of 136 possible baselines. ! AN table 4: 39 integrations on 136 of 136 possible baselines. ! Apparent sampling: 0.277422 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1153+8058 ! ! 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.22497e+09 8e+06 1 8e+06 ! 02 2 2.24497e+09 8e+06 1 8e+06 ! 03 3 2.33497e+09 8e+06 1 8e+06 ! 04 4 2.36497e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1899 lines of history. ! ! Reading 98700 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 J1153+8058_S_2001_03_12_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 4871 telescope corrections were flagged in sub-array 1. ! A total of 1173 telescope corrections were flagged in sub-array 2. ! A total of 306 telescope corrections were flagged in sub-array 3. ! A total of 402 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 5284 telescope corrections were flagged in sub-array 1. ! A total of 1381 telescope corrections were flagged in sub-array 2. ! A total of 342 telescope corrections were flagged in sub-array 3. ! A total of 447 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 4822 telescope corrections were flagged in sub-array 1. ! A total of 1102 telescope corrections were flagged in sub-array 2. ! A total of 282 telescope corrections were flagged in sub-array 3. ! A total of 388 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 5251 telescope corrections were flagged in sub-array 1. ! A total of 1372 telescope corrections were flagged in sub-array 2. ! A total of 373 telescope corrections were flagged in sub-array 3. ! A total of 469 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.061622Jy sigma=0.008428 ! Fit after self-cal, rms=0.665191Jy sigma=0.005094 ! 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=1.931 mas, bpa=43.4 degrees ! Estimated noise=663.598 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.655556 Jy ! Component: 100 - total flux cleaned = 0.951977 Jy ! Total flux subtracted in 100 components = 0.951977 Jy ! Clean residual min=-0.024466 max=0.139629 Jy/beam ! Clean residual mean=0.000036 rms=0.006976 Jy/beam ! Combined flux in latest and established models = 0.951977 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.951977 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.579820Jy sigma=0.004473 ! Fit after self-cal, rms=0.570051Jy sigma=0.004409 ! Inverting map ! Added new window around map position (-0.5, 2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.176003 Jy ! Component: 100 - total flux cleaned = 0.271384 Jy ! Total flux subtracted in 100 components = 0.271384 Jy ! Clean residual min=-0.015246 max=0.056083 Jy/beam ! Clean residual mean=0.000031 rms=0.004229 Jy/beam ! Combined flux in latest and established models = 1.22336 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 11 components and 1.22336 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.466122Jy sigma=0.003458 ! Fit after self-cal, rms=0.464641Jy sigma=0.003447 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0642245 Jy ! Component: 100 - total flux cleaned = 0.106567 Jy ! Total flux subtracted in 100 components = 0.106567 Jy ! Clean residual min=-0.014187 max=0.045210 Jy/beam ! Clean residual mean=0.000030 rms=0.003686 Jy/beam ! Combined flux in latest and established models = 1.32993 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 19 components and 1.32993 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.440237Jy sigma=0.003218 ! Fit after self-cal, rms=0.440079Jy sigma=0.003216 ! Inverting map ! Added new window around map position (2, -3). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0422909 Jy ! Component: 100 - total flux cleaned = 0.0702469 Jy ! Total flux subtracted in 100 components = 0.0702469 Jy ! Clean residual min=-0.013444 max=0.018430 Jy/beam ! Clean residual mean=0.000028 rms=0.003437 Jy/beam ! Combined flux in latest and established models = 1.40017 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 31 components and 1.40017 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.428278Jy sigma=0.003106 ! Fit after self-cal, rms=0.428074Jy sigma=0.003104 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0247544 Jy ! Component: 100 - total flux cleaned = 0.0437983 Jy ! Total flux subtracted in 100 components = 0.0437983 Jy ! Clean residual min=-0.013483 max=0.016497 Jy/beam ! Clean residual mean=0.000027 rms=0.003330 Jy/beam ! Combined flux in latest and established models = 1.44397 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 41 components and 1.44397 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.422456Jy sigma=0.003053 ! Fit after self-cal, rms=0.422369Jy sigma=0.003052 ! 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.383 mas, bmaj=2.489 mas, bpa=34.38 degrees ! Estimated noise=328.681 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.056661 Jy ! Component: 100 - total flux cleaned = 0.0929836 Jy ! Component: 150 - total flux cleaned = 0.117328 Jy ! Component: 200 - total flux cleaned = 0.133865 Jy ! Total flux subtracted in 200 components = 0.133865 Jy ! Clean residual min=-0.014441 max=0.017332 Jy/beam ! Clean residual mean=0.000068 rms=0.004049 Jy/beam ! Combined flux in latest and established models = 1.57784 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 68 components and 1.57784 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.413156Jy sigma=0.002971 ! Fit after self-cal, rms=0.413045Jy sigma=0.002970 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 75 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.13 HH 1.00* HN 0.93 ! KK 1.02 KP 1.00* LA 0.93 MA 1.26 ! MC 1.04 MK 0.94 NL 0.92 ON 0.88 ! OV 1.01 PT 0.97 SC 1.08 TS 1.04 ! WZ 0.85 ! ! A total of 75 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.15 HH 1.00* HN 0.92 ! KK 0.96 KP 1.00* LA 0.96 MA 1.36 ! MC 1.09 MK 0.89 NL 0.95 ON 0.90 ! OV 1.08 PT 0.99 SC 1.12 TS 1.00 ! WZ 0.86 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.90 FD 1.06 HH 1.00* HN 0.85 ! KK 0.99 KP 1.00* LA 0.92 MA 1.05 ! MC 1.12 MK 0.92 NL 0.89 ON 0.86 ! OV 1.00* PT 0.96 SC 0.97 TS 1.02 ! WZ 0.86 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.95 FD 1.16 HH 1.00* HN 0.93 ! KK 0.95 KP 1.00* LA 0.96 MA 1.35 ! MC 1.05 MK 0.95 NL 0.94 ON 0.87 ! OV 1.07 PT 1.00 SC 1.10 TS 0.96 ! WZ 0.88 ! ! ! Correcting IF 2. ! A total of 237 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.04 HH 1.00* HN 0.99 ! KK 1.00* KP 1.04 LA 0.98 MA 0.60 ! MC 1.00* MK 0.99 NL 1.00 ON 1.00* ! OV 0.94 PT 1.03 SC 0.91 TS 0.98 ! WZ 0.82 ! ! A total of 21 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.92 FD 1.07 HH 1.00* HN 0.99 ! KK 1.00* KP 1.09 LA 1.01 MA 0.89 ! MC 1.00* MK 0.95 NL 1.02 ON 1.00* ! OV 0.98 PT 1.07 SC 0.92 TS 0.96 ! WZ 0.80 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.88 FD 0.98 HH 1.00* HN 0.91 ! KK 1.00* KP 1.02 LA 0.94 MA 0.94 ! MC 1.00* MK 0.94 NL 0.98 ON 1.00* ! OV 1.00* PT 1.02 SC 0.80 TS 1.00 ! WZ 0.78 ! ! A total of 9 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.93 FD 1.06 HH 1.00* HN 1.02 ! KK 1.00* KP 1.07 LA 0.98 MA 0.51 ! MC 1.00* MK 1.01 NL 1.03 ON 1.00* ! OV 0.99 PT 1.05 SC 0.94 TS 0.97 ! WZ 0.78 ! ! ! Correcting IF 3. ! A total of 72 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.98 HH 1.00* HN 1.28 ! KK 0.98 KP 1.02 LA 0.93 MA 0.87 ! MC 1.14 MK 1.05 NL 0.87 ON 0.85 ! OV 1.12 PT 0.85 SC 1.19 TS 1.09 ! WZ 0.91 ! ! A total of 78 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 0.69 HH 1.00* HN 1.17 ! KK 0.97 KP 1.16 LA 0.95 MA 0.74 ! MC 1.28 MK 1.02 NL 1.26 ON 0.85 ! OV 1.08 PT 1.29 SC 1.24 TS 0.97 ! WZ 0.88 ! ! A total of 6 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.07 HH 1.00* HN 1.24 ! KK 0.94 KP 1.02 LA 1.12 MA 0.89 ! MC 1.11 MK 1.03 NL 1.04 ON 0.83 ! OV 1.00* PT 0.97 SC 1.10 TS 1.01 ! WZ 0.84 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 2.03 FD 1.00 HH 1.00* HN 1.24 ! KK 1.08 KP 1.12 LA 0.72 MA 0.88 ! MC 1.25 MK 1.05 NL 0.94 ON 0.82 ! OV 1.08 PT 0.89 SC 1.29 TS 0.94 ! WZ 0.87 ! ! ! Correcting IF 4. ! A total of 210 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 0.90 HH 1.00* HN 1.20 ! KK 1.00* KP 1.09 LA 0.86 MA 1.00* ! MC 1.11 MK 0.99 NL 1.04 ON 1.00* ! OV 0.95 PT 1.05 SC 1.26 TS 1.04 ! WZ 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.93 HH 1.00* HN 1.19 ! KK 1.00* KP 1.14 LA 0.89 MA 1.00* ! MC 1.19 MK 0.94 NL 1.06 ON 1.00* ! OV 1.00 PT 0.93 SC 1.30 TS 1.04 ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.98 FD 0.86 HH 1.00* HN 1.10 ! KK 1.00* KP 1.07 LA 0.84 MA 1.00* ! MC 1.24 MK 0.96 NL 1.01 ON 1.00* ! OV 1.00* PT 1.04 SC 1.14 TS 1.05 ! WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 0.94 HH 1.00* HN 1.22 ! KK 1.00* KP 1.12 LA 0.88 MA 1.00* ! MC 0.88 MK 0.97 NL 1.07 ON 1.00* ! OV 0.99 PT 0.91 SC 1.31 TS 0.99 ! WZ 1.00* ! ! ! Fit before self-cal, rms=0.413045Jy sigma=0.002970 ! Fit after self-cal, rms=0.370100Jy sigma=0.002661 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.397 mas, bmaj=2.494 mas, bpa=38.87 degrees ! Estimated noise=329.628 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.00129372 Jy ! Component: 100 - total flux cleaned = -0.00152544 Jy ! Component: 150 - total flux cleaned = -0.00153545 Jy ! Component: 200 - total flux cleaned = -0.00170358 Jy ! Total flux subtracted in 200 components = -0.00170358 Jy ! Clean residual min=-0.012091 max=0.009145 Jy/beam ! Clean residual mean=-0.000166 rms=0.002247 Jy/beam ! Combined flux in latest and established models = 1.57613 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 74 components and 1.57613 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.369774Jy sigma=0.002659 ! Fit after self-cal, rms=0.370727Jy sigma=0.002656 ! 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.370727Jy sigma=0.002656 ! Fit after self-cal, rms=0.357761Jy sigma=0.002528 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.39 mas, bmaj=2.493 mas, bpa=40.17 degrees ! Estimated noise=330.264 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00280585 Jy ! Component: 100 - total flux cleaned = -0.00339919 Jy ! Component: 150 - total flux cleaned = -0.0041509 Jy ! Component: 200 - total flux cleaned = -0.00439362 Jy ! Total flux subtracted in 200 components = -0.00439362 Jy ! Clean residual min=-0.005865 max=0.007381 Jy/beam ! Clean residual mean=-0.000036 rms=0.001322 Jy/beam ! Combined flux in latest and established models = 1.57174 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 76 components and 1.57174 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.357566Jy sigma=0.002527 ! Fit after self-cal, rms=0.358072Jy sigma=0.002526 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 (5.5, -5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00761555 Jy ! Component: 100 - total flux cleaned = 0.00969342 Jy ! Component: 150 - total flux cleaned = 0.00983134 Jy ! Component: 200 - total flux cleaned = 0.0100649 Jy ! Total flux subtracted in 200 components = 0.0100649 Jy ! Clean residual min=-0.005554 max=0.004505 Jy/beam ! Clean residual mean=-0.000039 rms=0.001218 Jy/beam ! Combined flux in latest and established models = 1.5818 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 94 components and 1.5818 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.357869Jy sigma=0.002524 ! Fit after self-cal, rms=0.357857Jy sigma=0.002524 ! 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.357857Jy sigma=0.002524 ! Fit after self-cal, rms=0.358917Jy sigma=0.002523 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.391 mas, bmaj=2.491 mas, bpa=40.87 degrees ! Estimated noise=331.097 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00242923 Jy ! Component: 100 - total flux cleaned = 0.00301873 Jy ! Component: 150 - total flux cleaned = 0.00326912 Jy ! Component: 200 - total flux cleaned = 0.00360593 Jy ! Total flux subtracted in 200 components = 0.00360593 Jy ! Clean residual min=-0.005108 max=0.004450 Jy/beam ! Clean residual mean=-0.000035 rms=0.001173 Jy/beam ! Combined flux in latest and established models = 1.58541 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 99 components and 1.58541 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.358798Jy sigma=0.002522 ! Fit after self-cal, rms=0.358818Jy sigma=0.002522 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.000700539 Jy ! Component: 100 - total flux cleaned = 0.000988104 Jy ! Component: 150 - total flux cleaned = 0.00124587 Jy ! Component: 200 - total flux cleaned = 0.00147966 Jy ! Total flux subtracted in 200 components = 0.00147966 Jy ! Clean residual min=-0.005027 max=0.004459 Jy/beam ! Clean residual mean=-0.000035 rms=0.001163 Jy/beam ! Combined flux in latest and established models = 1.58689 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 102 components and 1.58689 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.358767Jy sigma=0.002522 ! Fit after self-cal, rms=0.358764Jy sigma=0.002522 ! 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.358764Jy sigma=0.002522 ! Fit after self-cal, rms=0.359502Jy sigma=0.002521 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.829 mas, bmaj=9.304 mas, bpa=45.81 degrees ! Estimated noise=501.775 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00406578 Jy ! Component: 100 - total flux cleaned = 0.00433868 Jy ! Component: 150 - total flux cleaned = 0.00460735 Jy ! Component: 200 - total flux cleaned = 0.00475871 Jy ! Total flux subtracted in 200 components = 0.00475871 Jy ! Clean residual min=-0.008132 max=0.007049 Jy/beam ! Clean residual mean=-0.000127 rms=0.002385 Jy/beam ! Combined flux in latest and established models = 1.59165 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 107 components and 1.59165 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.359487Jy sigma=0.002521 ! Fit after self-cal, rms=0.359483Jy sigma=0.002521 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.000271517 Jy ! Component: 100 - total flux cleaned = 0.00051833 Jy ! Component: 150 - total flux cleaned = 0.000663661 Jy ! Component: 200 - total flux cleaned = 0.000663116 Jy ! Total flux subtracted in 200 components = 0.000663116 Jy ! Clean residual min=-0.008025 max=0.007028 Jy/beam ! Clean residual mean=-0.000126 rms=0.002395 Jy/beam ! Combined flux in latest and established models = 1.59231 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 113 components and 1.59231 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.359483Jy sigma=0.002521 ! Fit after self-cal, rms=0.359479Jy sigma=0.002521 ! 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.392 mas, bmaj=2.491 mas, bpa=41.34 degrees ! Estimated noise=331.891 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.392 x 2.491 at 41.34 degrees (North through East) ! Clean map min=-0.0038483 max=0.9235 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.359479Jy sigma=0.002521 ! Fit after self-cal, rms=0.398034Jy sigma=0.002124 ! Inverting map and beam ! Estimated beam: bmin=2.393 mas, bmaj=2.495 mas, bpa=40.72 degrees ! Estimated noise=332.386 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00264697 Jy ! Component: 100 - total flux cleaned = -0.0028776 Jy ! Component: 150 - total flux cleaned = -0.00316339 Jy ! Component: 200 - total flux cleaned = -0.00332131 Jy ! Total flux subtracted in 200 components = -0.00332131 Jy ! Clean residual min=-0.003571 max=0.003739 Jy/beam ! Clean residual mean=-0.000013 rms=0.000914 Jy/beam ! Combined flux in latest and established models = 1.58899 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 118 components and 1.58899 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.397961Jy sigma=0.002124 ! Fit after self-cal, rms=0.397962Jy sigma=0.002124 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -7.62015e-05 Jy ! Component: 100 - total flux cleaned = -7.80851e-05 Jy ! Component: 150 - total flux cleaned = -1.91798e-05 Jy ! Component: 200 - total flux cleaned = -2.01292e-05 Jy ! Total flux subtracted in 200 components = -2.01292e-05 Jy ! Clean residual min=-0.003551 max=0.003706 Jy/beam ! Clean residual mean=-0.000013 rms=0.000907 Jy/beam ! Combined flux in latest and established models = 1.58897 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 122 components and 1.58897 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.397933Jy sigma=0.002123 ! Fit after self-cal, rms=0.397929Jy sigma=0.002123 ! 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.397929Jy sigma=0.002123 ! Fit after self-cal, rms=0.400609Jy sigma=0.002123 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.879 mas, bmaj=1.923 mas, bpa=52.11 degrees ! Estimated noise=668.706 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.67067 Jy ! Component: 100 - total flux cleaned = 0.984588 Jy ! Component: 150 - total flux cleaned = 1.15012 Jy ! Component: 200 - total flux cleaned = 1.25011 Jy ! Total flux subtracted in 200 components = 1.25011 Jy ! Clean residual min=-0.009155 max=0.053015 Jy/beam ! Clean residual mean=0.000014 rms=0.002985 Jy/beam ! Combined flux in latest and established models = 1.25011 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 12 model components to the UV plane model. ! The established model now contains 12 components and 1.25011 Jy ! Inverting map and beam ! Estimated beam: bmin=2.395 mas, bmaj=2.495 mas, bpa=41.11 degrees ! Estimated noise=332.633 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.131136 Jy ! Component: 100 - total flux cleaned = 0.202485 Jy ! Component: 150 - total flux cleaned = 0.243871 Jy ! Component: 200 - total flux cleaned = 0.270138 Jy ! Total flux subtracted in 200 components = 0.270138 Jy ! Clean residual min=-0.005097 max=0.014968 Jy/beam ! Clean residual mean=-0.000005 rms=0.001852 Jy/beam ! Combined flux in latest and established models = 1.52024 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 54 components and 1.52024 Jy ! Inverting map ! Added new window around map position (-0.5, -2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0177781 Jy ! Component: 100 - total flux cleaned = 0.0303712 Jy ! Component: 150 - total flux cleaned = 0.0395776 Jy ! Component: 200 - total flux cleaned = 0.0464632 Jy ! Total flux subtracted in 200 components = 0.0464632 Jy ! Clean residual min=-0.004021 max=0.004689 Jy/beam ! Clean residual mean=-0.000011 rms=0.001036 Jy/beam ! Combined flux in latest and established models = 1.56671 Jy ! Adding 47 model components to the UV plane model. ! The established model now contains 87 components and 1.56671 Jy ! Inverting map ! Added new window around map position (9.5, -11). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00556756 Jy ! Component: 100 - total flux cleaned = 0.00990783 Jy ! Component: 150 - total flux cleaned = 0.0133746 Jy ! Component: 200 - total flux cleaned = 0.0161585 Jy ! Total flux subtracted in 200 components = 0.0161585 Jy ! Clean residual min=-0.003656 max=0.003964 Jy/beam ! Clean residual mean=-0.000014 rms=0.000897 Jy/beam ! Combined flux in latest and established models = 1.58287 Jy ! Adding 50 model components to the UV plane model. ! The established model now contains 117 components and 1.58287 Jy ! Inverting map ! Added new window around map position (23.5, -15.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00339043 Jy ! Component: 100 - total flux cleaned = 0.00528516 Jy ! Component: 150 - total flux cleaned = 0.00641383 Jy ! Component: 200 - total flux cleaned = 0.00709813 Jy ! Total flux subtracted in 200 components = 0.00709813 Jy ! Clean residual min=-0.003529 max=0.003231 Jy/beam ! Clean residual mean=-0.000018 rms=0.000855 Jy/beam ! Combined flux in latest and established models = 1.58996 Jy ! Adding 46 model components to the UV plane model. ! The established model now contains 137 components and 1.58996 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.835 mas, bmaj=9.311 mas, bpa=44 degrees ! Estimated noise=500.377 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00281444 Jy ! Component: 100 - total flux cleaned = 0.00438971 Jy ! Component: 150 - total flux cleaned = 0.00501081 Jy ! Component: 200 - total flux cleaned = 0.00532593 Jy ! Total flux subtracted in 200 components = 0.00532593 Jy ! Clean residual min=-0.005581 max=0.004988 Jy/beam ! Clean residual mean=-0.000072 rms=0.001685 Jy/beam ! Combined flux in latest and established models = 1.59529 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.000206652 Jy ! Component: 100 - total flux cleaned = 0.000257262 Jy ! Component: 150 - total flux cleaned = 0.000304067 Jy ! Component: 200 - total flux cleaned = 0.000326508 Jy ! Total flux subtracted in 200 components = 0.000326508 Jy ! Clean residual min=-0.005595 max=0.005049 Jy/beam ! Clean residual mean=-0.000072 rms=0.001681 Jy/beam ! Combined flux in latest and established models = 1.59562 Jy ! Adding 41 model components to the UV plane model. ! The established model now contains 170 components and 1.59562 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.395 mas, bmaj=2.495 mas, bpa=41.11 degrees ! Estimated noise=332.633 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000459898 Jy ! Component: 100 - total flux cleaned = -0.00035235 Jy ! Component: 150 - total flux cleaned = -0.000304998 Jy ! Component: 200 - total flux cleaned = -0.000260316 Jy ! Total flux subtracted in 200 components = -0.000260316 Jy ! Clean residual min=-0.003341 max=0.003314 Jy/beam ! Clean residual mean=-0.000019 rms=0.000836 Jy/beam ! Combined flux in latest and established models = 1.59536 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 178 components and 1.59536 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.400479Jy sigma=0.002122 ! Fit after self-cal, rms=0.400450Jy sigma=0.002122 wmodel J1153+8058_S_2001_03_12_pus_map.mod ! Writing 178 model components to file: J1153+8058_S_2001_03_12_pus_map.mod wobs J1153+8058_S_2001_03_12_pus_uvs.fits ! Writing UV FITS file: J1153+8058_S_2001_03_12_pus_uvs.fits wwins J1153+8058_S_2001_03_12_pus_map.win ! wwins: Wrote 7 windows to J1153+8058_S_2001_03_12_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.000540419 Jy ! Component: 100 - total flux cleaned = -0.00120383 Jy ! Component: 150 - total flux cleaned = -0.00152108 Jy ! Component: 200 - total flux cleaned = -0.00198185 Jy ! Component: 250 - total flux cleaned = -0.00238143 Jy ! Component: 300 - total flux cleaned = -0.00272368 Jy ! Component: 350 - total flux cleaned = -0.00301169 Jy ! Component: 400 - total flux cleaned = -0.00324687 Jy ! Component: 450 - total flux cleaned = -0.00343134 Jy ! Component: 500 - total flux cleaned = -0.00352224 Jy ! Component: 550 - total flux cleaned = -0.00361232 Jy ! Component: 600 - total flux cleaned = -0.00370153 Jy ! Component: 650 - total flux cleaned = -0.00387652 Jy ! Component: 700 - total flux cleaned = -0.00387684 Jy ! Component: 750 - total flux cleaned = -0.00400628 Jy ! Component: 800 - total flux cleaned = -0.00421844 Jy ! Component: 850 - total flux cleaned = -0.00421789 Jy ! Component: 900 - total flux cleaned = -0.00434295 Jy ! Component: 950 - total flux cleaned = -0.00434303 Jy ! Component: 1000 - total flux cleaned = -0.00438412 Jy ! Component: 1050 - total flux cleaned = -0.00442535 Jy ! Component: 1100 - total flux cleaned = -0.0046668 Jy ! Component: 1150 - total flux cleaned = -0.00478594 Jy ! Component: 1200 - total flux cleaned = -0.00470624 Jy ! Component: 1250 - total flux cleaned = -0.00478453 Jy ! Component: 1300 - total flux cleaned = -0.00490189 Jy ! Component: 1350 - total flux cleaned = -0.00490173 Jy ! Component: 1400 - total flux cleaned = -0.00482491 Jy ! Component: 1450 - total flux cleaned = -0.0049783 Jy ! Component: 1500 - total flux cleaned = -0.00505407 Jy ! Component: 1550 - total flux cleaned = -0.00497842 Jy ! Component: 1600 - total flux cleaned = -0.00497881 Jy ! Component: 1650 - total flux cleaned = -0.00524013 Jy ! Component: 1700 - total flux cleaned = -0.00512896 Jy ! Component: 1750 - total flux cleaned = -0.00512876 Jy ! Component: 1800 - total flux cleaned = -0.00509219 Jy ! Component: 1850 - total flux cleaned = -0.00520165 Jy ! Component: 1900 - total flux cleaned = -0.00512908 Jy ! Component: 1950 - total flux cleaned = -0.00509301 Jy ! Component: 2000 - total flux cleaned = -0.00512878 Jy ! Total flux subtracted in 2048 components = -0.005093 Jy ! Clean residual min=-0.001790 max=0.002533 Jy/beam ! Clean residual mean=-0.000019 rms=0.000687 Jy/beam ! Combined flux in latest and established models = 1.59026 Jy keep ! Adding 178 model components to the UV plane model. ! The established model now contains 356 components and 1.59026 Jy ![Exited script file: muppet] wmap "J1153+8058_S_2001_03_12_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.395 x 2.495 at 41.11 degrees (North through East) ! Clean map min=-0.0025918 max=0.9285 Jy/beam ! Writing clean map to FITS file: J1153+8058_S_2001_03_12_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.928496 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000686292 device "J1153+8058_S_2001_03_12_pus\_map.eps/VPS" ! Attempting to open device: 'J1153+8058_S_2001_03_12_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.221743 0.221743 0.443486 0.886973 1.77395 3.54789 7.09578 14.1916 28.3831 56.7663 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.002 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.928 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.395 x 2.495 at 41.11 degrees (North through East) ! Clean map min=-0.0025918 max=0.9285 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J1153+8058_S_2001_03_12_pus_uvs.log closed on Wed Jan 2 18:53:13 2013