! Started logfile: /scr/pima/bp177a_uvs/J0533+4822_X_uvs.log on Fri Feb 7 04:24:17 2014 obs /scr/pima/bp177a_uvs/J0533+4822_X_uva.fits ! Reading UV FITS file: /scr/pima/bp177a_uvs/J0533+4822_X_uva.fits ! AN table 1: 3 integrations on 28 of 28 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J0533+4822 ! ! 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 672 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/J0533+4822_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. ! ! 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=1.243092Jy sigma=46.529366 ! Fit after self-cal, rms=0.271855Jy sigma=10.521632 ! 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.7589 mas, bmaj=2.739 mas, bpa=-49.24 degrees ! Estimated noise=1.08827 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.57256 Jy ! Component: 100 - total flux cleaned = 0.722788 Jy ! Total flux subtracted in 100 components = 0.722788 Jy ! Clean residual min=-0.046721 max=0.073452 Jy/beam ! Clean residual mean=0.000044 rms=0.012517 Jy/beam ! Combined flux in latest and established models = 0.722788 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.722788 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.177731Jy sigma=6.540949 ! Fit after self-cal, rms=0.177754Jy sigma=6.537968 ! 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.9111 mas, bmaj=2.765 mas, bpa=-49.32 degrees ! Estimated noise=0.766155 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.101146 Jy ! Component: 100 - total flux cleaned = 0.155895 Jy ! Component: 150 - total flux cleaned = 0.194247 Jy ! Component: 200 - total flux cleaned = 0.222171 Jy ! Total flux subtracted in 200 components = 0.222171 Jy ! Clean residual min=-0.033350 max=0.030049 Jy/beam ! Clean residual mean=0.000098 rms=0.008176 Jy/beam ! Combined flux in latest and established models = 0.944959 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 42 components and 0.944959 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.101914Jy sigma=3.674774 ! Fit after self-cal, rms=0.101251Jy sigma=3.666905 ! 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.9111 x 2.765 at -49.32 degrees (North through East) ! Clean map min=-0.0324 max=0.77631 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.16 FD 1.00 HN 0.92 LA 0.94 ! MK 1.02 OV 0.94 PT 0.87 SC 0.93 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.91 FD 1.10 HN 0.98 LA 0.91 ! MK 1.06 OV 0.98 PT 0.90 SC 0.98 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.09 HN 1.03 LA 0.94 ! MK 1.09 OV 0.98 PT 0.91 SC 1.01 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 0.99 HN 0.99 LA 0.95 ! MK 1.01 OV 0.96 PT 0.88 SC 0.96 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.07 HN 1.11 LA 1.00 ! MK 1.04 OV 1.02 PT 0.92 SC 1.00 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 1.03 HN 1.02 LA 0.93 ! MK 1.02 OV 1.03 PT 0.89 SC 0.99 ! ! ! 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.101251Jy sigma=3.666905 ! Fit after self-cal, rms=0.074681Jy sigma=2.617385 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9218 mas, bmaj=2.791 mas, bpa=-49.33 degrees ! Estimated noise=0.757318 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.00707255 Jy ! Component: 100 - total flux cleaned = 0.00925729 Jy ! Component: 150 - total flux cleaned = 0.00992952 Jy ! Component: 200 - total flux cleaned = 0.0106913 Jy ! Total flux subtracted in 200 components = 0.0106913 Jy ! Clean residual min=-0.014353 max=0.018490 Jy/beam ! Clean residual mean=0.000015 rms=0.003624 Jy/beam ! Combined flux in latest and established models = 0.95565 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 54 components and 0.95565 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.073590Jy sigma=2.582763 ! Fit after self-cal, rms=0.073380Jy sigma=2.581069 ! Inverting map ! Added new window around map position (4.4, -9.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00678129 Jy ! Component: 100 - total flux cleaned = 0.00451314 Jy ! Component: 150 - total flux cleaned = 0.00300979 Jy ! Component: 200 - total flux cleaned = 0.00172781 Jy ! Total flux subtracted in 200 components = 0.00172781 Jy ! Clean residual min=-0.010512 max=0.012100 Jy/beam ! Clean residual mean=0.000009 rms=0.002914 Jy/beam ! Combined flux in latest and established models = 0.957378 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 84 components and 0.957378 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.068751Jy sigma=2.405775 ! Fit after self-cal, rms=0.068462Jy sigma=2.400525 ! 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.068462Jy sigma=2.400525 ! Fit after self-cal, rms=0.058557Jy sigma=2.033387 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9212 mas, bmaj=2.79 mas, bpa=-49.34 degrees ! Estimated noise=0.757228 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00127672 Jy ! Component: 100 - total flux cleaned = 0.0012738 Jy ! Component: 150 - total flux cleaned = 0.000890877 Jy ! Component: 200 - total flux cleaned = 0.000897421 Jy ! Total flux subtracted in 200 components = 0.000897421 Jy ! Clean residual min=-0.009862 max=0.010305 Jy/beam ! Clean residual mean=-0.000014 rms=0.002499 Jy/beam ! Combined flux in latest and established models = 0.958275 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 96 components and 0.958276 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.057309Jy sigma=1.981803 ! Fit after self-cal, rms=0.057216Jy sigma=1.981038 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.000324983 Jy ! Component: 100 - total flux cleaned = 0.000165968 Jy ! Component: 150 - total flux cleaned = 0.000170048 Jy ! Component: 200 - total flux cleaned = 0.000306775 Jy ! Total flux subtracted in 200 components = 0.000306775 Jy ! Clean residual min=-0.009712 max=0.009313 Jy/beam ! Clean residual mean=-0.000017 rms=0.002345 Jy/beam ! Combined flux in latest and established models = 0.958582 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 105 components and 0.958582 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.056542Jy sigma=1.955340 ! Fit after self-cal, rms=0.056549Jy sigma=1.954888 ! 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.056549Jy sigma=1.954888 ! Fit after self-cal, rms=0.056208Jy sigma=1.947459 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9207 mas, bmaj=2.786 mas, bpa=-49.36 degrees ! Estimated noise=0.758102 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00320028 Jy ! Component: 100 - total flux cleaned = 0.00391958 Jy ! Component: 150 - total flux cleaned = 0.00483832 Jy ! Component: 200 - total flux cleaned = 0.00532182 Jy ! Total flux subtracted in 200 components = 0.00532182 Jy ! Clean residual min=-0.009661 max=0.009232 Jy/beam ! Clean residual mean=-0.000015 rms=0.002353 Jy/beam ! Combined flux in latest and established models = 0.963904 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 118 components and 0.963904 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.055561Jy sigma=1.919852 ! Fit after self-cal, rms=0.055556Jy sigma=1.919550 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.000350787 Jy ! Component: 100 - total flux cleaned = 0.000798928 Jy ! Component: 150 - total flux cleaned = 0.00111715 Jy ! Component: 200 - total flux cleaned = 0.00152397 Jy ! Total flux subtracted in 200 components = 0.00152397 Jy ! Clean residual min=-0.009471 max=0.008779 Jy/beam ! Clean residual mean=-0.000017 rms=0.002256 Jy/beam ! Combined flux in latest and established models = 0.965428 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 126 components and 0.965428 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.055194Jy sigma=1.903581 ! Fit after self-cal, rms=0.055197Jy sigma=1.903452 ! 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.055197Jy sigma=1.903452 ! Fit after self-cal, rms=0.043210Jy sigma=1.362691 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.08 mas, bmaj=2.96 mas, bpa=-52.85 degrees ! Estimated noise=0.92876 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00136972 Jy ! Component: 100 - total flux cleaned = -0.000807779 Jy ! Component: 150 - total flux cleaned = -0.00056018 Jy ! Component: 200 - total flux cleaned = -0.000676912 Jy ! Total flux subtracted in 200 components = -0.000676912 Jy ! Clean residual min=-0.009971 max=0.006932 Jy/beam ! Clean residual mean=0.000002 rms=0.002482 Jy/beam ! Combined flux in latest and established models = 0.964751 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 139 components and 0.964751 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.042815Jy sigma=1.343691 ! Fit after self-cal, rms=0.042852Jy sigma=1.342354 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 = -1.10969e-06 Jy ! Component: 100 - total flux cleaned = -0.000101066 Jy ! Component: 150 - total flux cleaned = -0.000195494 Jy ! Component: 200 - total flux cleaned = -0.000288603 Jy ! Total flux subtracted in 200 components = -0.000288603 Jy ! Clean residual min=-0.009381 max=0.006523 Jy/beam ! Clean residual mean=0.000003 rms=0.002340 Jy/beam ! Combined flux in latest and established models = 0.964462 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 148 components and 0.964462 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.042659Jy sigma=1.333747 ! Fit after self-cal, rms=0.042689Jy sigma=1.333068 ! 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.9214 mas, bmaj=2.791 mas, bpa=-49.38 degrees ! Estimated noise=0.759825 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9214 x 2.791 at -49.38 degrees (North through East) ! Clean map min=-0.0099368 max=0.76598 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.042689Jy sigma=1.333068 ! Fit after self-cal, rms=0.043294Jy sigma=1.180700 ! Inverting map and beam ! Estimated beam: bmin=0.9208 mas, bmaj=2.791 mas, bpa=-49.42 degrees ! Estimated noise=0.758645 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00161336 Jy ! Component: 100 - total flux cleaned = -0.00255172 Jy ! Component: 150 - total flux cleaned = -0.00287809 Jy ! Component: 200 - total flux cleaned = -0.0029483 Jy ! Total flux subtracted in 200 components = -0.0029483 Jy ! Clean residual min=-0.007339 max=0.005459 Jy/beam ! Clean residual mean=0.000001 rms=0.001511 Jy/beam ! Combined flux in latest and established models = 0.961514 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 167 components and 0.961514 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.042849Jy sigma=1.164460 ! Fit after self-cal, rms=0.042832Jy sigma=1.164303 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -7.23505e-05 Jy ! Component: 100 - total flux cleaned = -0.000136358 Jy ! Component: 150 - total flux cleaned = -0.000315046 Jy ! Component: 200 - total flux cleaned = -0.000426073 Jy ! Total flux subtracted in 200 components = -0.000426073 Jy ! Clean residual min=-0.007188 max=0.005418 Jy/beam ! Clean residual mean=0.000001 rms=0.001473 Jy/beam ! Combined flux in latest and established models = 0.961088 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 173 components and 0.961088 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.042630Jy sigma=1.156968 ! Fit after self-cal, rms=0.042627Jy sigma=1.156857 ! 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.042627Jy sigma=1.156857 ! Fit after self-cal, rms=0.042689Jy sigma=1.153634 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.7672 mas, bmaj=2.752 mas, bpa=-49.39 degrees ! Estimated noise=1.09019 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.573521 Jy ! Component: 100 - total flux cleaned = 0.71748 Jy ! Component: 150 - total flux cleaned = 0.785458 Jy ! Component: 200 - total flux cleaned = 0.822522 Jy ! Total flux subtracted in 200 components = 0.822522 Jy ! Clean residual min=-0.019243 max=0.025981 Jy/beam ! Clean residual mean=0.000023 rms=0.005224 Jy/beam ! Combined flux in latest and established models = 0.822522 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 16 model components to the UV plane model. ! The established model now contains 16 components and 0.822522 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9214 mas, bmaj=2.792 mas, bpa=-49.45 degrees ! Estimated noise=0.758908 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.0396309 Jy ! Component: 100 - total flux cleaned = 0.0650239 Jy ! Component: 150 - total flux cleaned = 0.0797289 Jy ! Component: 200 - total flux cleaned = 0.0920043 Jy ! Total flux subtracted in 200 components = 0.0920043 Jy ! Clean residual min=-0.011296 max=0.016574 Jy/beam ! Clean residual mean=0.000034 rms=0.003807 Jy/beam ! Combined flux in latest and established models = 0.914526 Jy ! Adding 38 model components to the UV plane model. ! The established model now contains 48 components and 0.914526 Jy ! Inverting map ! Added new window around map position (-0.2, 3). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0159994 Jy ! Component: 100 - total flux cleaned = 0.0261196 Jy ! Component: 150 - total flux cleaned = 0.0318748 Jy ! Component: 200 - total flux cleaned = 0.0370119 Jy ! Total flux subtracted in 200 components = 0.0370119 Jy ! Clean residual min=-0.007551 max=0.008283 Jy/beam ! Clean residual mean=0.000016 rms=0.002147 Jy/beam ! Combined flux in latest and established models = 0.951538 Jy ! Adding 32 model components to the UV plane model. ! The established model now contains 73 components and 0.951538 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.078 mas, bmaj=2.964 mas, bpa=-53.12 degrees ! Estimated noise=0.927775 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00918179 Jy ! Component: 100 - total flux cleaned = 0.0117373 Jy ! Component: 150 - total flux cleaned = 0.0135986 Jy ! Component: 200 - total flux cleaned = 0.0151409 Jy ! Total flux subtracted in 200 components = 0.0151409 Jy ! Clean residual min=-0.007748 max=0.008034 Jy/beam ! Clean residual mean=0.000015 rms=0.002213 Jy/beam ! Combined flux in latest and established models = 0.966679 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.000586299 Jy ! Component: 100 - total flux cleaned = 0.000709755 Jy ! Component: 150 - total flux cleaned = 0.00082591 Jy ! Component: 200 - total flux cleaned = 0.000931772 Jy ! Total flux subtracted in 200 components = 0.000931772 Jy ! Clean residual min=-0.006555 max=0.006993 Jy/beam ! Clean residual mean=0.000011 rms=0.001929 Jy/beam ! Combined flux in latest and established models = 0.967611 Jy ! Adding 57 model components to the UV plane model. ! The established model now contains 123 components and 0.967611 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.9214 mas, bmaj=2.792 mas, bpa=-49.45 degrees ! Estimated noise=0.758908 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00276184 Jy ! Component: 100 - total flux cleaned = -0.00434092 Jy ! Component: 150 - total flux cleaned = -0.00502066 Jy ! Component: 200 - total flux cleaned = -0.00526213 Jy ! Total flux subtracted in 200 components = -0.00526213 Jy ! Clean residual min=-0.005122 max=0.005029 Jy/beam ! Clean residual mean=0.000006 rms=0.001280 Jy/beam ! Combined flux in latest and established models = 0.962349 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 144 components and 0.962349 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.042024Jy sigma=1.129429 ! Fit after self-cal, rms=0.041965Jy sigma=1.118936 wmodel /scr/pima/bp177a_uvs/J0533+4822_X_map.mod ! Writing 144 model components to file: /scr/pima/bp177a_uvs/J0533+4822_X_map.mod wobs /scr/pima/bp177a_uvs/J0533+4822_X_uvs.fits ! Writing UV FITS file: /scr/pima/bp177a_uvs/J0533+4822_X_uvs.fits wwins /scr/pima/bp177a_uvs/J0533+4822_X_map.win ! wwins: Wrote 3 windows to /scr/pima/bp177a_uvs/J0533+4822_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.000718248 Jy ! Component: 100 - total flux cleaned = 0.000324493 Jy ! Component: 150 - total flux cleaned = 0.000251645 Jy ! Component: 200 - total flux cleaned = -2.01092e-05 Jy ! Component: 250 - total flux cleaned = -0.00014846 Jy ! Component: 300 - total flux cleaned = -0.000269804 Jy ! Component: 350 - total flux cleaned = -0.000327765 Jy ! Component: 400 - total flux cleaned = -0.000380707 Jy ! Component: 450 - total flux cleaned = -0.000537247 Jy ! Component: 500 - total flux cleaned = -0.000536211 Jy ! Component: 550 - total flux cleaned = -0.000536472 Jy ! Component: 600 - total flux cleaned = -0.000536849 Jy ! Component: 650 - total flux cleaned = -0.000537127 Jy ! Component: 700 - total flux cleaned = -0.000323132 Jy ! Component: 750 - total flux cleaned = -0.00040678 Jy ! Component: 800 - total flux cleaned = -0.000447791 Jy ! Component: 850 - total flux cleaned = -0.000487859 Jy ! Component: 900 - total flux cleaned = -0.00048707 Jy ! Component: 950 - total flux cleaned = -0.000339872 Jy ! Component: 1000 - total flux cleaned = -0.000375152 Jy ! Component: 1050 - total flux cleaned = -0.000374893 Jy ! Component: 1100 - total flux cleaned = -0.000306699 Jy ! Component: 1150 - total flux cleaned = -0.000373542 Jy ! Component: 1200 - total flux cleaned = -0.000340241 Jy ! Component: 1250 - total flux cleaned = -0.000372086 Jy ! Component: 1300 - total flux cleaned = -0.000434475 Jy ! Component: 1350 - total flux cleaned = -0.000403253 Jy ! Component: 1400 - total flux cleaned = -0.000403153 Jy ! Component: 1450 - total flux cleaned = -0.000343301 Jy ! Component: 1500 - total flux cleaned = -0.000461164 Jy ! Component: 1550 - total flux cleaned = -0.000547567 Jy ! Component: 1600 - total flux cleaned = -0.000518586 Jy ! Component: 1650 - total flux cleaned = -0.000518772 Jy ! Component: 1700 - total flux cleaned = -0.0005468 Jy ! Component: 1750 - total flux cleaned = -0.000601909 Jy ! Component: 1800 - total flux cleaned = -0.000601817 Jy ! Component: 1850 - total flux cleaned = -0.000601771 Jy ! Component: 1900 - total flux cleaned = -0.00060193 Jy ! Component: 1950 - total flux cleaned = -0.000654585 Jy ! Component: 2000 - total flux cleaned = -0.000759036 Jy ! Component: 2050 - total flux cleaned = -0.00081079 Jy ! Component: 2100 - total flux cleaned = -0.000862148 Jy ! Component: 2150 - total flux cleaned = -0.000963908 Jy ! Component: 2200 - total flux cleaned = -0.00098912 Jy ! Component: 2250 - total flux cleaned = -0.00101433 Jy ! Component: 2300 - total flux cleaned = -0.00111413 Jy ! Component: 2350 - total flux cleaned = -0.00116346 Jy ! Component: 2400 - total flux cleaned = -0.00116353 Jy ! Component: 2450 - total flux cleaned = -0.00113913 Jy ! Component: 2500 - total flux cleaned = -0.0012854 Jy ! Component: 2550 - total flux cleaned = -0.00128534 Jy ! Component: 2600 - total flux cleaned = -0.00138229 Jy ! Component: 2650 - total flux cleaned = -0.0014065 Jy ! Component: 2700 - total flux cleaned = -0.00152729 Jy ! Component: 2750 - total flux cleaned = -0.00162362 Jy ! Component: 2800 - total flux cleaned = -0.00171972 Jy ! Component: 2850 - total flux cleaned = -0.00179165 Jy ! Component: 2900 - total flux cleaned = -0.00193557 Jy ! Component: 2950 - total flux cleaned = -0.00198344 Jy ! Component: 3000 - total flux cleaned = -0.00205534 Jy ! Component: 3050 - total flux cleaned = -0.00212702 Jy ! Component: 3100 - total flux cleaned = -0.00219875 Jy ! Component: 3150 - total flux cleaned = -0.00222262 Jy ! Component: 3200 - total flux cleaned = -0.00227035 Jy ! Component: 3250 - total flux cleaned = -0.00231816 Jy ! Component: 3300 - total flux cleaned = -0.00241376 Jy ! Component: 3350 - total flux cleaned = -0.00248561 Jy ! Component: 3400 - total flux cleaned = -0.00253347 Jy ! Component: 3450 - total flux cleaned = -0.00255762 Jy ! Component: 3500 - total flux cleaned = -0.00260583 Jy ! Component: 3550 - total flux cleaned = -0.00260593 Jy ! Component: 3600 - total flux cleaned = -0.00265436 Jy ! Component: 3650 - total flux cleaned = -0.00265424 Jy ! Component: 3700 - total flux cleaned = -0.00267839 Jy ! Component: 3750 - total flux cleaned = -0.00270286 Jy ! Component: 3800 - total flux cleaned = -0.00272728 Jy ! Component: 3850 - total flux cleaned = -0.00275185 Jy ! Component: 3900 - total flux cleaned = -0.00277672 Jy ! Component: 3950 - total flux cleaned = -0.00277648 Jy ! Component: 4000 - total flux cleaned = -0.00280135 Jy ! Component: 4050 - total flux cleaned = -0.00285164 Jy ! Total flux subtracted in 4096 components = -0.00285156 Jy ! Clean residual min=-0.001743 max=0.001474 Jy/beam ! Clean residual mean=0.000006 rms=0.000397 Jy/beam ! Combined flux in latest and established models = 0.959497 Jy keep ! Adding 249 model components to the UV plane model. ! The established model now contains 393 components and 0.959498 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/scr/pima/bp177a_uvs/J0533+4822_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9214 x 2.792 at -49.45 degrees (North through East) ! Clean map min=-0.006301 max=0.76596 Jy/beam ! Writing clean map to FITS file: /scr/pima/bp177a_uvs/J0533+4822_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.765956 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000398366 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.156027 0.156027 0.312054 0.624108 1.24822 2.49643 4.99286 9.98572 19.9714 39.9429 79.8858 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0011 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.765 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1922.74 quit ! Quitting program ! Log file /scr/pima/bp177a_uvs/J0533+4822_X_uvs.log closed on Fri Feb 7 04:24:21 2014