! Started logfile: /scr/pima/bp177a_uvs/J0921+6215_X_uvs.log on Fri Feb 7 04:26:05 2014 obs /scr/pima/bp177a_uvs/J0921+6215_X_uva.fits ! Reading UV FITS file: /scr/pima/bp177a_uvs/J0921+6215_X_uva.fits ! AN table 1: 5 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.9 visibilities/baseline/integration-bin. ! Found source: J0921+6215 ! ! There are 8 IFs, and a total of 8 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 7.392e+09 3.2e+07 1 3.2e+07 ! 02 2 7.424e+09 3.2e+07 1 3.2e+07 ! 03 3 7.456e+09 3.2e+07 1 3.2e+07 ! 04 4 7.552e+09 3.2e+07 1 3.2e+07 ! 05 5 7.744e+09 3.2e+07 1 3.2e+07 ! 06 6 7.776e+09 3.2e+07 1 3.2e+07 ! 07 7 7.808e+09 3.2e+07 1 3.2e+07 ! 08 8 7.84e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 1008 visibilities. select RR ! Selecting polarization: RR, channels: 1..8 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: 5..5 ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.2 print "taper_size =", taper_size ! taper_size = 70 ![@/home/lpetrov/bin/01_muppet /scr/pima/bp177a_uvs/J0921+6215_X] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 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 = 1024x1024 pixels with 0.200x0.200 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 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.315696Jy sigma=45.399994 ! Fit after self-cal, rms=0.346659Jy sigma=12.160398 ! 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=0.6987 mas, bmaj=3.056 mas, bpa=-49.08 degrees ! Estimated noise=1.01687 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.556735 Jy ! Component: 100 - total flux cleaned = 0.727324 Jy ! Total flux subtracted in 100 components = 0.727324 Jy ! Clean residual min=-0.046646 max=0.078190 Jy/beam ! Clean residual mean=-0.000047 rms=0.014471 Jy/beam ! Combined flux in latest and established models = 0.727325 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.727325 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.202204Jy sigma=7.078280 ! Fit after self-cal, rms=0.201245Jy sigma=7.057686 ! 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=0.9365 mas, bmaj=2.867 mas, bpa=-50.29 degrees ! Estimated noise=0.666703 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.120754 Jy ! Component: 100 - total flux cleaned = 0.192879 Jy ! Component: 150 - total flux cleaned = 0.24227 Jy ! Component: 200 - total flux cleaned = 0.27763 Jy ! Total flux subtracted in 200 components = 0.27763 Jy ! Clean residual min=-0.033789 max=0.035132 Jy/beam ! Clean residual mean=-0.000023 rms=0.009347 Jy/beam ! Combined flux in latest and established models = 1.00495 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 43 components and 1.00495 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.098600Jy sigma=3.280125 ! Fit after self-cal, rms=0.098000Jy sigma=3.253925 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9365 x 2.867 at -50.29 degrees (North through East) ! Clean map min=-0.029946 max=0.73901 Jy/beam if ( peak(flux,max)/imstat(rms) < lim_self ) selflim lim_val ! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2. ! Selfcal phase corrections will not be limited. gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.18 FD 1.00 HN 0.89 LA 0.94 ! MK 0.99 OV 0.92 PT 0.87 SC 0.99 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.08 HN 0.93 LA 0.92 ! MK 1.03 OV 0.97 PT 0.88 SC 1.03 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.07 HN 1.00 LA 0.96 ! MK 1.06 OV 0.96 PT 0.91 SC 1.07 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 0.99 HN 0.97 LA 0.96 ! MK 0.96 OV 0.95 PT 0.87 SC 1.00 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.06 HN 1.08 LA 1.03 ! MK 1.00 OV 0.99 PT 0.91 SC 1.06 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.12 FD 1.02 HN 0.99 LA 0.94 ! MK 0.99 OV 1.02 PT 0.88 SC 1.03 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* LA 1.00* ! MK 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! ! Fit before self-cal, rms=0.098000Jy sigma=3.253925 ! Fit after self-cal, rms=0.067937Jy sigma=2.210036 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9428 mas, bmaj=2.89 mas, bpa=-50.52 degrees ! Estimated noise=0.655477 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.00459976 Jy ! Component: 100 - total flux cleaned = 0.00515522 Jy ! Component: 150 - total flux cleaned = 0.00515765 Jy ! Component: 200 - total flux cleaned = 0.00561826 Jy ! Total flux subtracted in 200 components = 0.00561826 Jy ! Clean residual min=-0.014299 max=0.013979 Jy/beam ! Clean residual mean=-0.000025 rms=0.003656 Jy/beam ! Combined flux in latest and established models = 1.01057 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 55 components and 1.01057 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.065477Jy sigma=2.131590 ! Fit after self-cal, rms=0.065494Jy sigma=2.129676 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.065494Jy sigma=2.129676 ! Fit after self-cal, rms=0.065255Jy sigma=2.121257 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9403 mas, bmaj=2.889 mas, bpa=-50.51 degrees ! Estimated noise=0.655518 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00588129 Jy ! Component: 100 - total flux cleaned = 0.00693213 Jy ! Component: 150 - total flux cleaned = 0.00769315 Jy ! Component: 200 - total flux cleaned = 0.00872338 Jy ! Total flux subtracted in 200 components = 0.00872338 Jy ! Clean residual min=-0.013644 max=0.013297 Jy/beam ! Clean residual mean=-0.000026 rms=0.003506 Jy/beam ! Combined flux in latest and established models = 1.0193 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 62 components and 1.0193 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.063942Jy sigma=2.081403 ! Fit after self-cal, rms=0.063939Jy sigma=2.080753 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.000807802 Jy ! Component: 100 - total flux cleaned = 0.00151206 Jy ! Component: 150 - total flux cleaned = 0.00215941 Jy ! Component: 200 - total flux cleaned = 0.00274759 Jy ! Total flux subtracted in 200 components = 0.00274759 Jy ! Clean residual min=-0.013162 max=0.013008 Jy/beam ! Clean residual mean=-0.000024 rms=0.003421 Jy/beam ! Combined flux in latest and established models = 1.02204 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 69 components and 1.02204 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.063363Jy sigma=2.064622 ! Fit after self-cal, rms=0.063372Jy sigma=2.064347 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.063372Jy sigma=2.064347 ! Fit after self-cal, rms=0.063071Jy sigma=2.049626 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9382 mas, bmaj=2.887 mas, bpa=-50.47 degrees ! Estimated noise=0.65657 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00548499 Jy ! Component: 100 - total flux cleaned = 0.0088939 Jy ! Component: 150 - total flux cleaned = 0.00971509 Jy ! Component: 200 - total flux cleaned = 0.0103416 Jy ! Total flux subtracted in 200 components = 0.0103416 Jy ! Clean residual min=-0.012912 max=0.013020 Jy/beam ! Clean residual mean=-0.000026 rms=0.003443 Jy/beam ! Combined flux in latest and established models = 1.03239 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 70 components and 1.03239 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.062260Jy sigma=2.027289 ! Fit after self-cal, rms=0.062272Jy sigma=2.027077 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.000813675 Jy ! Component: 100 - total flux cleaned = 0.0013296 Jy ! Component: 150 - total flux cleaned = 0.00161796 Jy ! Component: 200 - total flux cleaned = 0.00188843 Jy ! Total flux subtracted in 200 components = 0.00188843 Jy ! Clean residual min=-0.012600 max=0.012763 Jy/beam ! Clean residual mean=-0.000024 rms=0.003393 Jy/beam ! Combined flux in latest and established models = 1.03427 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 72 components and 1.03427 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.061945Jy sigma=2.018043 ! Fit after self-cal, rms=0.061967Jy sigma=2.017878 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.061967Jy sigma=2.017878 ! Fit after self-cal, rms=0.061737Jy sigma=2.007170 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.083 mas, bmaj=3.091 mas, bpa=-51.04 degrees ! Estimated noise=0.81665 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00740329 Jy ! Component: 100 - total flux cleaned = 0.0113242 Jy ! Component: 150 - total flux cleaned = 0.0122012 Jy ! Component: 200 - total flux cleaned = 0.0127598 Jy ! Total flux subtracted in 200 components = 0.0127598 Jy ! Clean residual min=-0.017834 max=0.017089 Jy/beam ! Clean residual mean=-0.000050 rms=0.005371 Jy/beam ! Combined flux in latest and established models = 1.04703 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 81 components and 1.04703 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.061268Jy sigma=1.995720 ! Fit after self-cal, rms=0.061344Jy sigma=1.994671 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.000359211 Jy ! Component: 100 - total flux cleaned = 0.000518022 Jy ! Component: 150 - total flux cleaned = 0.000595143 Jy ! Component: 200 - total flux cleaned = 0.0007485 Jy ! Total flux subtracted in 200 components = 0.0007485 Jy ! Clean residual min=-0.017540 max=0.017186 Jy/beam ! Clean residual mean=-0.000049 rms=0.005342 Jy/beam ! Combined flux in latest and established models = 1.04778 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 87 components and 1.04778 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.061237Jy sigma=1.989444 ! Fit after self-cal, rms=0.061251Jy sigma=1.989121 ! 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=0.9368 mas, bmaj=2.886 mas, bpa=-50.44 degrees ! Estimated noise=0.657392 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9368 x 2.886 at -50.44 degrees (North through East) ! Clean map min=-0.019227 max=0.72533 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.061251Jy sigma=1.989121 ! Fit after self-cal, rms=0.044443Jy sigma=1.341946 ! Inverting map and beam ! Estimated beam: bmin=0.9344 mas, bmaj=2.893 mas, bpa=-50.51 degrees ! Estimated noise=0.659167 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00244143 Jy ! Component: 100 - total flux cleaned = -0.00408767 Jy ! Component: 150 - total flux cleaned = -0.00502306 Jy ! Component: 200 - total flux cleaned = -0.00571182 Jy ! Total flux subtracted in 200 components = -0.00571182 Jy ! Clean residual min=-0.009897 max=0.009828 Jy/beam ! Clean residual mean=0.000007 rms=0.002567 Jy/beam ! Combined flux in latest and established models = 1.04207 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 98 components and 1.04207 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.043772Jy sigma=1.323076 ! Fit after self-cal, rms=0.043784Jy sigma=1.322251 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000182858 Jy ! Component: 100 - total flux cleaned = -0.00034827 Jy ! Component: 150 - total flux cleaned = -0.000665732 Jy ! Component: 200 - total flux cleaned = -0.000956999 Jy ! Total flux subtracted in 200 components = -0.000956999 Jy ! Clean residual min=-0.009576 max=0.009731 Jy/beam ! Clean residual mean=0.000007 rms=0.002543 Jy/beam ! Combined flux in latest and established models = 1.04111 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 106 components and 1.04111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.043558Jy sigma=1.317889 ! Fit after self-cal, rms=0.043617Jy sigma=1.317260 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.043617Jy sigma=1.317260 ! Fit after self-cal, rms=0.043494Jy sigma=1.311325 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=0.6978 mas, bmaj=3.068 mas, bpa=-49.14 degrees ! Estimated noise=1.00446 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.549859 Jy ! Component: 100 - total flux cleaned = 0.719123 Jy ! Component: 150 - total flux cleaned = 0.800237 Jy ! Component: 200 - total flux cleaned = 0.848696 Jy ! Total flux subtracted in 200 components = 0.848696 Jy ! Clean residual min=-0.024966 max=0.032593 Jy/beam ! Clean residual mean=-0.000025 rms=0.006798 Jy/beam ! Combined flux in latest and established models = 0.848696 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 22 model components to the UV plane model. ! The established model now contains 22 components and 0.848696 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9348 mas, bmaj=2.896 mas, bpa=-50.5 degrees ! Estimated noise=0.660384 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.0589067 Jy ! Component: 100 - total flux cleaned = 0.0975847 Jy ! Component: 150 - total flux cleaned = 0.124733 Jy ! Component: 200 - total flux cleaned = 0.144061 Jy ! Total flux subtracted in 200 components = 0.144061 Jy ! Clean residual min=-0.015782 max=0.018858 Jy/beam ! Clean residual mean=-0.000019 rms=0.004494 Jy/beam ! Combined flux in latest and established models = 0.992758 Jy ! Adding 32 model components to the UV plane model. ! The established model now contains 48 components and 0.992758 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.084 mas, bmaj=3.112 mas, bpa=-51.44 degrees ! Estimated noise=0.820091 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0242657 Jy ! Component: 100 - total flux cleaned = 0.0359901 Jy ! Component: 150 - total flux cleaned = 0.0381144 Jy ! Component: 200 - total flux cleaned = 0.0383691 Jy ! Total flux subtracted in 200 components = 0.0383691 Jy ! Clean residual min=-0.014235 max=0.016310 Jy/beam ! Clean residual mean=0.000002 rms=0.004603 Jy/beam ! Combined flux in latest and established models = 1.03113 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.000450205 Jy ! Component: 100 - total flux cleaned = 0.000654817 Jy ! Component: 150 - total flux cleaned = 0.00101907 Jy ! Component: 200 - total flux cleaned = 0.00118697 Jy ! Total flux subtracted in 200 components = 0.00118697 Jy ! Clean residual min=-0.013884 max=0.015284 Jy/beam ! Clean residual mean=0.000006 rms=0.004303 Jy/beam ! Combined flux in latest and established models = 1.03231 Jy ! Adding 35 model components to the UV plane model. ! The established model now contains 76 components and 1.03231 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.9348 mas, bmaj=2.896 mas, bpa=-50.5 degrees ! Estimated noise=0.660384 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00303636 Jy ! Component: 100 - total flux cleaned = -0.00287517 Jy ! Component: 150 - total flux cleaned = -0.00255316 Jy ! Component: 200 - total flux cleaned = -0.00225966 Jy ! Total flux subtracted in 200 components = -0.00225966 Jy ! Clean residual min=-0.011927 max=0.012006 Jy/beam ! Clean residual mean=0.000002 rms=0.002753 Jy/beam ! Combined flux in latest and established models = 1.03005 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 83 components and 1.03005 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.045796Jy sigma=1.378321 ! Fit after self-cal, rms=0.045656Jy sigma=1.373936 wmodel /scr/pima/bp177a_uvs/J0921+6215_X_map.mod ! Writing 83 model components to file: /scr/pima/bp177a_uvs/J0921+6215_X_map.mod wobs /scr/pima/bp177a_uvs/J0921+6215_X_uvs.fits ! Writing UV FITS file: /scr/pima/bp177a_uvs/J0921+6215_X_uvs.fits wwins /scr/pima/bp177a_uvs/J0921+6215_X_map.win ! wwins: Wrote 1 windows to /scr/pima/bp177a_uvs/J0921+6215_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000302416 Jy ! Component: 100 - total flux cleaned = 0.00273012 Jy ! Component: 150 - total flux cleaned = 0.0049348 Jy ! Component: 200 - total flux cleaned = 0.00695398 Jy ! Component: 250 - total flux cleaned = 0.00883441 Jy ! Component: 300 - total flux cleaned = 0.0108431 Jy ! Component: 350 - total flux cleaned = 0.0124845 Jy ! Component: 400 - total flux cleaned = 0.0137014 Jy ! Component: 450 - total flux cleaned = 0.0146439 Jy ! Component: 500 - total flux cleaned = 0.0150466 Jy ! Component: 550 - total flux cleaned = 0.0157107 Jy ! Component: 600 - total flux cleaned = 0.0163411 Jy ! Component: 650 - total flux cleaned = 0.0163386 Jy ! Component: 700 - total flux cleaned = 0.0163383 Jy ! Component: 750 - total flux cleaned = 0.0164964 Jy ! Component: 800 - total flux cleaned = 0.0165701 Jy ! Component: 850 - total flux cleaned = 0.0162753 Jy ! Component: 900 - total flux cleaned = 0.0161359 Jy ! Component: 950 - total flux cleaned = 0.0160007 Jy ! Component: 1000 - total flux cleaned = 0.0157348 Jy ! Component: 1050 - total flux cleaned = 0.0154807 Jy ! Component: 1100 - total flux cleaned = 0.0152966 Jy ! Component: 1150 - total flux cleaned = 0.0149997 Jy ! Component: 1200 - total flux cleaned = 0.0148879 Jy ! Component: 1250 - total flux cleaned = 0.0147235 Jy ! Component: 1300 - total flux cleaned = 0.0146173 Jy ! Component: 1350 - total flux cleaned = 0.0145661 Jy ! Component: 1400 - total flux cleaned = 0.0145158 Jy ! Component: 1450 - total flux cleaned = 0.0146135 Jy ! Component: 1500 - total flux cleaned = 0.0148488 Jy ! Component: 1550 - total flux cleaned = 0.0148919 Jy ! Component: 1600 - total flux cleaned = 0.0149354 Jy ! Component: 1650 - total flux cleaned = 0.0149799 Jy ! Component: 1700 - total flux cleaned = 0.0151095 Jy ! Component: 1750 - total flux cleaned = 0.0152343 Jy ! Component: 1800 - total flux cleaned = 0.0152741 Jy ! Component: 1850 - total flux cleaned = 0.015115 Jy ! Component: 1900 - total flux cleaned = 0.015233 Jy ! Component: 1950 - total flux cleaned = 0.0154613 Jy ! Component: 2000 - total flux cleaned = 0.0155354 Jy ! Component: 2050 - total flux cleaned = 0.0156815 Jy ! Component: 2100 - total flux cleaned = 0.015824 Jy ! Component: 2150 - total flux cleaned = 0.0160352 Jy ! Component: 2200 - total flux cleaned = 0.016207 Jy ! Component: 2250 - total flux cleaned = 0.0162746 Jy ! Component: 2300 - total flux cleaned = 0.0163082 Jy ! Component: 2350 - total flux cleaned = 0.0164061 Jy ! Component: 2400 - total flux cleaned = 0.0165009 Jy ! Component: 2450 - total flux cleaned = 0.0164682 Jy ! Component: 2500 - total flux cleaned = 0.016499 Jy ! Component: 2550 - total flux cleaned = 0.0165901 Jy ! Component: 2600 - total flux cleaned = 0.0165601 Jy ! Component: 2650 - total flux cleaned = 0.0165895 Jy ! Component: 2700 - total flux cleaned = 0.0165898 Jy ! Component: 2750 - total flux cleaned = 0.0166475 Jy ! Component: 2800 - total flux cleaned = 0.0167046 Jy ! Component: 2850 - total flux cleaned = 0.0166765 Jy ! Component: 2900 - total flux cleaned = 0.016621 Jy ! Component: 2950 - total flux cleaned = 0.0166211 Jy ! Component: 3000 - total flux cleaned = 0.0165666 Jy ! Component: 3050 - total flux cleaned = 0.0165391 Jy ! Component: 3100 - total flux cleaned = 0.0165123 Jy ! Component: 3150 - total flux cleaned = 0.0165124 Jy ! Component: 3200 - total flux cleaned = 0.0164594 Jy ! Component: 3250 - total flux cleaned = 0.0164594 Jy ! Component: 3300 - total flux cleaned = 0.0164592 Jy ! Component: 3350 - total flux cleaned = 0.016433 Jy ! Component: 3400 - total flux cleaned = 0.016433 Jy ! Component: 3450 - total flux cleaned = 0.016433 Jy ! Component: 3500 - total flux cleaned = 0.0164588 Jy ! Component: 3550 - total flux cleaned = 0.0164849 Jy ! Component: 3600 - total flux cleaned = 0.016537 Jy ! Component: 3650 - total flux cleaned = 0.0165372 Jy ! Component: 3700 - total flux cleaned = 0.0165896 Jy ! Component: 3750 - total flux cleaned = 0.016616 Jy ! Component: 3800 - total flux cleaned = 0.0166161 Jy ! Component: 3850 - total flux cleaned = 0.0166425 Jy ! Component: 3900 - total flux cleaned = 0.0167221 Jy ! Component: 3950 - total flux cleaned = 0.0167223 Jy ! Component: 4000 - total flux cleaned = 0.0167764 Jy ! Component: 4050 - total flux cleaned = 0.0168033 Jy ! Total flux subtracted in 4096 components = 0.0168849 Jy ! Clean residual min=-0.002015 max=0.001574 Jy/beam ! Clean residual mean=0.000016 rms=0.000452 Jy/beam ! Combined flux in latest and established models = 1.04694 Jy keep ! Adding 295 model components to the UV plane model. ! The established model now contains 378 components and 1.04694 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/scr/pima/bp177a_uvs/J0921+6215_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9348 x 2.896 at -50.5 degrees (North through East) ! Clean map min=-0.010747 max=0.73222 Jy/beam ! Writing clean map to FITS file: /scr/pima/bp177a_uvs/J0921+6215_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.732218 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000442801 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.181422 0.181422 0.362844 0.725687 1.45137 2.90275 5.8055 11.611 23.222 46.444 92.888 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0013 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.732 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1653.6 quit ! Quitting program ! Log file /scr/pima/bp177a_uvs/J0921+6215_X_uvs.log closed on Fri Feb 7 04:26:09 2014