! Started logfile: /scr/pima/bg219b_uvs/J0643-2451_S_uvs.log on Thu Apr 7 01:15:15 2016 obs /scr/pima/bg219b_uvs/J0643-2451_S_uva.fits ! Reading UV FITS file: /scr/pima/bg219b_uvs/J0643-2451_S_uva.fits ! AN table 1: 47 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.770095 visibilities/baseline/integration-bin. ! Found source: J0643-2451 ! ! 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.22e+09 3.2e+07 1 3.2e+07 ! 02 2 2.252e+09 3.2e+07 1 3.2e+07 ! 03 3 2.284e+09 3.2e+07 1 3.2e+07 ! 04 4 2.348e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 54 lines of history. ! ! Reading 5212 visibilities. select rr, 1,2,3,4,6,7,8 ! Restricting channel ranges to the available 4 channels. ! 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 = 1024 print "field_cell =", field_cell ! field_cell = 0.5 print "taper_size =", taper_size ! taper_size = 20 ![@/home/lpetrov/bin/01_muppet /scr/pima/bg219b_uvs/J0643-2451_S] 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.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 25 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 25 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 25 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 25 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.054506Jy sigma=18.092732 ! Fit after self-cal, rms=0.744023Jy sigma=12.837829 ! 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=3.086 mas, bmaj=12.83 mas, bpa=17.62 degrees ! Estimated noise=0.895412 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.187407 Jy ! Component: 100 - total flux cleaned = 0.247468 Jy ! Total flux subtracted in 100 components = 0.247468 Jy ! Clean residual min=-0.012910 max=0.027663 Jy/beam ! Clean residual mean=0.000049 rms=0.004487 Jy/beam ! Combined flux in latest and established models = 0.247468 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 15 components and 0.247468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146523Jy sigma=1.823138 ! Fit after self-cal, rms=0.146420Jy sigma=1.821543 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0313161 Jy ! Component: 100 - total flux cleaned = 0.0489584 Jy ! Total flux subtracted in 100 components = 0.0489584 Jy ! Clean residual min=-0.010227 max=0.017336 Jy/beam ! Clean residual mean=0.000084 rms=0.003277 Jy/beam ! Combined flux in latest and established models = 0.296427 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 34 components and 0.296427 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136188Jy sigma=1.639757 ! Fit after self-cal, rms=0.136118Jy sigma=1.638971 ! 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=3.881 mas, bmaj=14.07 mas, bpa=16.46 degrees ! Estimated noise=0.593768 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.012815 Jy ! Component: 100 - total flux cleaned = 0.0195103 Jy ! Component: 150 - total flux cleaned = 0.0204177 Jy ! Component: 200 - total flux cleaned = 0.0202603 Jy ! Total flux subtracted in 200 components = 0.0202603 Jy ! Clean residual min=-0.012284 max=0.022608 Jy/beam ! Clean residual mean=0.000165 rms=0.004351 Jy/beam ! Combined flux in latest and established models = 0.316687 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 53 components and 0.316687 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133889Jy sigma=1.621466 ! Fit after self-cal, rms=0.133857Jy sigma=1.621060 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.881 x 14.07 at 16.46 degrees (North through East) ! Clean map min=-0.012024 max=0.25896 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* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! ! Fit before self-cal, rms=0.133857Jy sigma=1.621060 ! Fit after self-cal, rms=0.133857Jy sigma=1.621060 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.881 mas, bmaj=14.07 mas, bpa=16.46 degrees ! Estimated noise=0.593768 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 (54, 20.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0190634 Jy ! Component: 100 - total flux cleaned = 0.0273786 Jy ! Component: 150 - total flux cleaned = 0.0309385 Jy ! Component: 200 - total flux cleaned = 0.0320005 Jy ! Total flux subtracted in 200 components = 0.0320005 Jy ! Clean residual min=-0.011535 max=0.017517 Jy/beam ! Clean residual mean=0.000126 rms=0.003654 Jy/beam ! Combined flux in latest and established models = 0.348688 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 99 components and 0.348688 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130129Jy sigma=1.567496 ! Fit after self-cal, rms=0.129847Jy sigma=1.560541 ! 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.129847Jy sigma=1.560541 ! Fit after self-cal, rms=0.128705Jy sigma=1.516944 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.861 mas, bmaj=14.08 mas, bpa=16.69 degrees ! Estimated noise=0.59487 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00193126 Jy ! Component: 100 - total flux cleaned = -0.00345788 Jy ! Component: 150 - total flux cleaned = -0.00404183 Jy ! Component: 200 - total flux cleaned = -0.00443384 Jy ! Total flux subtracted in 200 components = -0.00443384 Jy ! Clean residual min=-0.009599 max=0.009487 Jy/beam ! Clean residual mean=0.000055 rms=0.002570 Jy/beam ! Combined flux in latest and established models = 0.344254 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 124 components and 0.344254 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128268Jy sigma=1.513055 ! Fit after self-cal, rms=0.128270Jy sigma=1.512462 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.00117036 Jy ! Component: 100 - total flux cleaned = -0.00151114 Jy ! Component: 150 - total flux cleaned = -0.00182187 Jy ! Component: 200 - total flux cleaned = -0.00211104 Jy ! Total flux subtracted in 200 components = -0.00211104 Jy ! Clean residual min=-0.009363 max=0.009311 Jy/beam ! Clean residual mean=0.000060 rms=0.002561 Jy/beam ! Combined flux in latest and established models = 0.342143 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 136 components and 0.342143 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128080Jy sigma=1.511447 ! Fit after self-cal, rms=0.128080Jy sigma=1.511370 ! 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.128080Jy sigma=1.511370 ! Fit after self-cal, rms=0.127793Jy sigma=1.502558 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.861 mas, bmaj=14.03 mas, bpa=16.6 degrees ! Estimated noise=0.596043 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000485645 Jy ! Component: 100 - total flux cleaned = 0.000343599 Jy ! Component: 150 - total flux cleaned = 0.000213099 Jy ! Component: 200 - total flux cleaned = 2.82393e-05 Jy ! Total flux subtracted in 200 components = 2.82393e-05 Jy ! Clean residual min=-0.007924 max=0.008615 Jy/beam ! Clean residual mean=0.000055 rms=0.002341 Jy/beam ! Combined flux in latest and established models = 0.342171 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 157 components and 0.342171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127624Jy sigma=1.502058 ! Fit after self-cal, rms=0.127638Jy sigma=1.501947 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.000115068 Jy ! Component: 100 - total flux cleaned = -0.000346475 Jy ! Component: 150 - total flux cleaned = -0.000510515 Jy ! Component: 200 - total flux cleaned = -0.000719276 Jy ! Total flux subtracted in 200 components = -0.000719276 Jy ! Clean residual min=-0.007834 max=0.008575 Jy/beam ! Clean residual mean=0.000055 rms=0.002331 Jy/beam ! Combined flux in latest and established models = 0.341452 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 171 components and 0.341452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127541Jy sigma=1.501472 ! Fit after self-cal, rms=0.127541Jy sigma=1.501451 ! 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.127541Jy sigma=1.501451 ! Fit after self-cal, rms=0.127526Jy sigma=1.501324 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.693 mas, bmaj=19.96 mas, bpa=7.384 degrees ! Estimated noise=0.685403 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000756819 Jy ! Component: 100 - total flux cleaned = -0.00112871 Jy ! Component: 150 - total flux cleaned = -0.00136168 Jy ! Component: 200 - total flux cleaned = -0.00163805 Jy ! Total flux subtracted in 200 components = -0.00163805 Jy ! Clean residual min=-0.010827 max=0.010888 Jy/beam ! Clean residual mean=0.000091 rms=0.003246 Jy/beam ! Combined flux in latest and established models = 0.339814 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 180 components and 0.339814 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127383Jy sigma=1.502170 ! Fit after self-cal, rms=0.127372Jy sigma=1.501927 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.000303621 Jy ! Component: 100 - total flux cleaned = -0.000689148 Jy ! Component: 150 - total flux cleaned = -0.00110769 Jy ! Component: 200 - total flux cleaned = -0.00160996 Jy ! Total flux subtracted in 200 components = -0.00160996 Jy ! Clean residual min=-0.010705 max=0.010954 Jy/beam ! Clean residual mean=0.000091 rms=0.003221 Jy/beam ! Combined flux in latest and established models = 0.338204 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 185 components and 0.338204 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127311Jy sigma=1.502785 ! Fit after self-cal, rms=0.127308Jy sigma=1.502691 ! 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=3.859 mas, bmaj=14.01 mas, bpa=16.59 degrees ! Estimated noise=0.596083 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.859 x 14.01 at 16.59 degrees (North through East) ! Clean map min=-0.0072991 max=0.25721 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.078 mas, bmaj=12.81 mas, bpa=17.67 degrees ! Estimated noise=0.890305 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.186079 Jy ! Component: 100 - total flux cleaned = 0.246535 Jy ! Component: 150 - total flux cleaned = 0.277881 Jy ! Component: 200 - total flux cleaned = 0.296944 Jy ! Total flux subtracted in 200 components = 0.296944 Jy ! Clean residual min=-0.008838 max=0.011645 Jy/beam ! Clean residual mean=0.000025 rms=0.002542 Jy/beam ! Combined flux in latest and established models = 0.296944 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 32 model components to the UV plane model. ! The established model now contains 32 components and 0.296944 Jy ! Inverting map and beam ! Estimated beam: bmin=3.859 mas, bmaj=14.01 mas, bpa=16.59 degrees ! Estimated noise=0.596083 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.0182323 Jy ! Component: 100 - total flux cleaned = 0.0299274 Jy ! Component: 150 - total flux cleaned = 0.0378728 Jy ! Component: 200 - total flux cleaned = 0.0415365 Jy ! Total flux subtracted in 200 components = 0.0415365 Jy ! Clean residual min=-0.008718 max=0.010365 Jy/beam ! Clean residual mean=0.000048 rms=0.002526 Jy/beam ! Combined flux in latest and established models = 0.33848 Jy ! Adding 45 model components to the UV plane model. ! The established model now contains 72 components and 0.33848 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=7.693 mas, bmaj=19.96 mas, bpa=7.384 degrees ! Estimated noise=0.685403 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.002842 Jy ! Component: 100 - total flux cleaned = 0.00422646 Jy ! Component: 150 - total flux cleaned = 0.00409062 Jy ! Component: 200 - total flux cleaned = 0.00381395 Jy ! Total flux subtracted in 200 components = 0.00381395 Jy ! Clean residual min=-0.010349 max=0.012072 Jy/beam ! Clean residual mean=0.000093 rms=0.003366 Jy/beam ! Combined flux in latest and established models = 0.342294 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.000375647 Jy ! Component: 100 - total flux cleaned = -0.000818583 Jy ! Component: 150 - total flux cleaned = -0.00111751 Jy ! Component: 200 - total flux cleaned = -0.00130166 Jy ! Total flux subtracted in 200 components = -0.00130166 Jy ! Clean residual min=-0.010521 max=0.011759 Jy/beam ! Clean residual mean=0.000095 rms=0.003293 Jy/beam ! Combined flux in latest and established models = 0.340993 Jy ! Adding 71 model components to the UV plane model. ! The established model now contains 133 components and 0.340993 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.859 mas, bmaj=14.01 mas, bpa=16.59 degrees ! Estimated noise=0.596083 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000461031 Jy ! Component: 100 - total flux cleaned = -0.000206214 Jy ! Component: 150 - total flux cleaned = -0.000550874 Jy ! Component: 200 - total flux cleaned = -0.000648848 Jy ! Total flux subtracted in 200 components = -0.000648848 Jy ! Clean residual min=-0.007679 max=0.008966 Jy/beam ! Clean residual mean=0.000060 rms=0.002351 Jy/beam ! Combined flux in latest and established models = 0.340344 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 149 components and 0.340344 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127761Jy sigma=1.502662 ! Fit after self-cal, rms=0.127737Jy sigma=1.502163 wmodel /scr/pima/bg219b_uvs/J0643-2451_S_map.mod ! Writing 149 model components to file: /scr/pima/bg219b_uvs/J0643-2451_S_map.mod wobs /scr/pima/bg219b_uvs/J0643-2451_S_uvs.fits ! Writing UV FITS file: /scr/pima/bg219b_uvs/J0643-2451_S_uvs.fits wwins /scr/pima/bg219b_uvs/J0643-2451_S_map.win ! wwins: Wrote 2 windows to /scr/pima/bg219b_uvs/J0643-2451_S_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.00322855 Jy ! Component: 100 - total flux cleaned = 0.00282495 Jy ! Component: 150 - total flux cleaned = 0.00154911 Jy ! Component: 200 - total flux cleaned = 0.000215663 Jy ! Component: 250 - total flux cleaned = -0.00106052 Jy ! Component: 300 - total flux cleaned = -0.00217523 Jy ! Component: 350 - total flux cleaned = -0.00356593 Jy ! Component: 400 - total flux cleaned = -0.00449332 Jy ! Component: 450 - total flux cleaned = -0.00529131 Jy ! Component: 500 - total flux cleaned = -0.00606411 Jy ! Component: 550 - total flux cleaned = -0.00634488 Jy ! Component: 600 - total flux cleaned = -0.00680155 Jy ! Component: 650 - total flux cleaned = -0.00706857 Jy ! Component: 700 - total flux cleaned = -0.00732851 Jy ! Component: 750 - total flux cleaned = -0.00749917 Jy ! Component: 800 - total flux cleaned = -0.00766496 Jy ! Component: 850 - total flux cleaned = -0.00798904 Jy ! Component: 900 - total flux cleaned = -0.00822631 Jy ! Component: 950 - total flux cleaned = -0.00837993 Jy ! Component: 1000 - total flux cleaned = -0.00845596 Jy ! Component: 1050 - total flux cleaned = -0.00860603 Jy ! Component: 1100 - total flux cleaned = -0.00897115 Jy ! Component: 1150 - total flux cleaned = -0.00897208 Jy ! Component: 1200 - total flux cleaned = -0.00925539 Jy ! Component: 1250 - total flux cleaned = -0.00946349 Jy ! Component: 1300 - total flux cleaned = -0.00987372 Jy ! Component: 1350 - total flux cleaned = -0.00994031 Jy ! Component: 1400 - total flux cleaned = -0.0100731 Jy ! Component: 1450 - total flux cleaned = -0.0103338 Jy ! Component: 1500 - total flux cleaned = -0.0105898 Jy ! Component: 1550 - total flux cleaned = -0.010716 Jy ! Component: 1600 - total flux cleaned = -0.0108403 Jy ! Component: 1650 - total flux cleaned = -0.0109617 Jy ! Component: 1700 - total flux cleaned = -0.0110213 Jy ! Component: 1750 - total flux cleaned = -0.0110211 Jy ! Component: 1800 - total flux cleaned = -0.0110802 Jy ! Component: 1850 - total flux cleaned = -0.0111378 Jy ! Component: 1900 - total flux cleaned = -0.0113081 Jy ! Component: 1950 - total flux cleaned = -0.0112501 Jy ! Component: 2000 - total flux cleaned = -0.0110832 Jy ! Component: 2050 - total flux cleaned = -0.0111936 Jy ! Component: 2100 - total flux cleaned = -0.0111931 Jy ! Component: 2150 - total flux cleaned = -0.0111928 Jy ! Component: 2200 - total flux cleaned = -0.011034 Jy ! Component: 2250 - total flux cleaned = -0.0110866 Jy ! Component: 2300 - total flux cleaned = -0.0109322 Jy ! Component: 2350 - total flux cleaned = -0.0111373 Jy ! Component: 2400 - total flux cleaned = -0.0109863 Jy ! Component: 2450 - total flux cleaned = -0.0109868 Jy ! Component: 2500 - total flux cleaned = -0.0109375 Jy ! Component: 2550 - total flux cleaned = -0.0109861 Jy ! Component: 2600 - total flux cleaned = -0.0108883 Jy ! Component: 2650 - total flux cleaned = -0.0107922 Jy ! Component: 2700 - total flux cleaned = -0.0105562 Jy ! Component: 2750 - total flux cleaned = -0.0106037 Jy ! Component: 2800 - total flux cleaned = -0.0107881 Jy ! Component: 2850 - total flux cleaned = -0.0107417 Jy ! Component: 2900 - total flux cleaned = -0.0108772 Jy ! Component: 2950 - total flux cleaned = -0.0107427 Jy ! Component: 3000 - total flux cleaned = -0.0108764 Jy ! Component: 3050 - total flux cleaned = -0.0107444 Jy ! Component: 3100 - total flux cleaned = -0.0107008 Jy ! Component: 3150 - total flux cleaned = -0.0105706 Jy ! Component: 3200 - total flux cleaned = -0.0103134 Jy ! Component: 3250 - total flux cleaned = -0.010271 Jy ! Component: 3300 - total flux cleaned = -0.0101868 Jy ! Component: 3350 - total flux cleaned = -0.00985296 Jy ! Component: 3400 - total flux cleaned = -0.0097704 Jy ! Component: 3450 - total flux cleaned = -0.00956513 Jy ! Component: 3500 - total flux cleaned = -0.0094838 Jy ! Component: 3550 - total flux cleaned = -0.00924154 Jy ! Component: 3600 - total flux cleaned = -0.00924153 Jy ! Component: 3650 - total flux cleaned = -0.00900273 Jy ! Component: 3700 - total flux cleaned = -0.0089239 Jy ! Component: 3750 - total flux cleaned = -0.00884562 Jy ! Component: 3800 - total flux cleaned = -0.00872909 Jy ! Component: 3850 - total flux cleaned = -0.0085743 Jy ! Component: 3900 - total flux cleaned = -0.00849719 Jy ! Component: 3950 - total flux cleaned = -0.00853527 Jy ! Component: 4000 - total flux cleaned = -0.0084217 Jy ! Component: 4050 - total flux cleaned = -0.00838431 Jy ! Total flux subtracted in 4096 components = -0.00842156 Jy ! Clean residual min=-0.002394 max=0.002110 Jy/beam ! Clean residual mean=0.000023 rms=0.000786 Jy/beam ! Combined flux in latest and established models = 0.331922 Jy keep ! Adding 433 model components to the UV plane model. ! The established model now contains 582 components and 0.331922 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/scr/pima/bg219b_uvs/J0643-2451_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.859 x 14.01 at 16.59 degrees (North through East) ! Clean map min=-0.0060796 max=0.25583 Jy/beam ! Writing clean map to FITS file: /scr/pima/bg219b_uvs/J0643-2451_S_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.25583 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000787027 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.92291 0.92291 1.84582 3.69164 7.38328 14.7666 29.5331 59.0662 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0023 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.255 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 325.058 quit ! Quitting program ! Log file /scr/pima/bg219b_uvs/J0643-2451_S_uvs.log closed on Thu Apr 7 01:15:19 2016