! Started logfile: /image_prod/ug002_uva//J2230+6946_S_l_dfm.log on Thu Mar 5 05:44:48 2020 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.22000 obs /image_prod/ug002_uva//J2230+6946_S_l_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J2230+6946_S_l_uva.fits ! AN table 1: 12 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.657738 visibilities/baseline/integration-bin. ! Found source: J2230+6946 ! ! There are 3 IFs, and a total of 3 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 ! ! Polarization(s): RR ! ! Read 36 lines of history. ! ! Reading 663 visibilities. select rr, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 ! Restricting channel ranges to the available 3 channels. ! Selecting polarization: RR, channels: 1..3 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 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 ![@/opt64/bin/pima_mupet_01.dfm J2230+6946_S_l] 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 print "===============================================================" ! =============================================================== ! print "=== Difmap script: pima_mupet_01.dvm version of 2016.05.19 === " ! === Difmap script: pima_mupet_01.dvm version of 2016.05.19 === ! 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 22 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 22 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 22 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.160114Jy sigma=7.700175 ! Fit after self-cal, rms=0.418800Jy sigma=2.444684 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.78 mas, bmaj=6.65 mas, bpa=-70.31 degrees ! Estimated noise=5.01178 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.478417 Jy ! Component: 100 - total flux cleaned = 0.590066 Jy ! Total flux subtracted in 100 components = 0.590066 Jy ! Clean residual min=-0.026441 max=0.055344 Jy/beam ! Clean residual mean=-0.000011 rms=0.009953 Jy/beam ! Combined flux in latest and established models = 0.590066 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.590066 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.156962Jy sigma=0.892863 ! Fit after self-cal, rms=0.156960Jy sigma=0.892816 ! 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.351 mas, bmaj=6.737 mas, bpa=-70.76 degrees ! Estimated noise=4.39409 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.0557634 Jy ! Component: 100 - total flux cleaned = 0.0878199 Jy ! Component: 150 - total flux cleaned = 0.10875 Jy ! Component: 200 - total flux cleaned = 0.122453 Jy ! Total flux subtracted in 200 components = 0.122453 Jy ! Clean residual min=-0.028641 max=0.037593 Jy/beam ! Clean residual mean=-0.000006 rms=0.007842 Jy/beam ! Combined flux in latest and established models = 0.712518 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 19 components and 0.712518 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.137513Jy sigma=0.702162 ! Fit after self-cal, rms=0.137510Jy sigma=0.702074 ! 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.351 x 6.737 at -70.76 degrees (North through East) ! Clean map min=-0.027797 max=0.61966 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.01 HN 1.00 KP 1.02 ! MK 0.96 NL 0.96 OV 0.89 SC 1.23 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.00 HN 0.95 KP 1.04 ! MK 0.94 NL 1.13 OV 0.89 SC 3.63 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 FD 1.00 HN 0.96 KP 1.00 ! MK 0.93 NL 0.92 OV 0.89 SC 1.27 ! ! ! Fit before self-cal, rms=0.137510Jy sigma=0.702074 ! Fit after self-cal, rms=0.117378Jy sigma=0.417207 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.426 mas, bmaj=6.606 mas, bpa=-70.72 degrees ! Estimated noise=4.3837 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.0102444 Jy ! Component: 100 - total flux cleaned = -0.0140981 Jy ! Component: 150 - total flux cleaned = -0.0154428 Jy ! Component: 200 - total flux cleaned = -0.0165332 Jy ! Total flux subtracted in 200 components = -0.0165332 Jy ! Clean residual min=-0.011297 max=0.010585 Jy/beam ! Clean residual mean=-0.000007 rms=0.002382 Jy/beam ! Combined flux in latest and established models = 0.695985 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 29 components and 0.695985 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.116100Jy sigma=0.405902 ! Fit after self-cal, rms=0.116636Jy sigma=0.405628 ! 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. ! ! Fit before self-cal, rms=0.116636Jy sigma=0.405628 ! Fit after self-cal, rms=0.117784Jy sigma=0.397752 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.432 mas, bmaj=6.633 mas, bpa=-70.68 degrees ! Estimated noise=4.37915 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00540084 Jy ! Component: 100 - total flux cleaned = -0.00821673 Jy ! Component: 150 - total flux cleaned = -0.0086574 Jy ! Component: 200 - total flux cleaned = -0.00911819 Jy ! Total flux subtracted in 200 components = -0.00911819 Jy ! Clean residual min=-0.010116 max=0.008930 Jy/beam ! Clean residual mean=-0.000003 rms=0.002011 Jy/beam ! Combined flux in latest and established models = 0.686867 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 37 components and 0.686867 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117238Jy sigma=0.395152 ! Fit after self-cal, rms=0.117254Jy sigma=0.395140 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.000149048 Jy ! Component: 100 - total flux cleaned = -0.000360172 Jy ! Component: 150 - total flux cleaned = -0.000541083 Jy ! Component: 200 - total flux cleaned = -0.000645465 Jy ! Total flux subtracted in 200 components = -0.000645465 Jy ! Clean residual min=-0.010022 max=0.008823 Jy/beam ! Clean residual mean=-0.000003 rms=0.001993 Jy/beam ! Combined flux in latest and established models = 0.686221 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 45 components and 0.686222 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117129Jy sigma=0.394751 ! Fit after self-cal, rms=0.117123Jy sigma=0.394742 ! 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. ! ! Fit before self-cal, rms=0.117123Jy sigma=0.394742 ! Fit after self-cal, rms=0.118166Jy sigma=0.393298 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.442 mas, bmaj=6.636 mas, bpa=-70.64 degrees ! Estimated noise=4.37722 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00324353 Jy ! Component: 100 - total flux cleaned = -0.00486605 Jy ! Component: 150 - total flux cleaned = -0.00514566 Jy ! Component: 200 - total flux cleaned = -0.00544428 Jy ! Total flux subtracted in 200 components = -0.00544428 Jy ! Clean residual min=-0.009955 max=0.008762 Jy/beam ! Clean residual mean=-0.000003 rms=0.001881 Jy/beam ! Combined flux in latest and established models = 0.680777 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 48 components and 0.680777 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117926Jy sigma=0.392527 ! Fit after self-cal, rms=0.117938Jy sigma=0.392522 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.000261543 Jy ! Component: 100 - total flux cleaned = -0.000398725 Jy ! Component: 150 - total flux cleaned = -0.00044025 Jy ! Component: 200 - total flux cleaned = -0.000561254 Jy ! Total flux subtracted in 200 components = -0.000561254 Jy ! Clean residual min=-0.009853 max=0.008638 Jy/beam ! Clean residual mean=-0.000003 rms=0.001875 Jy/beam ! Combined flux in latest and established models = 0.680216 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 55 components and 0.680216 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.117866Jy sigma=0.392377 ! Fit after self-cal, rms=0.117858Jy sigma=0.392371 ! 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. ! ! Fit before self-cal, rms=0.117858Jy sigma=0.392371 ! Fit after self-cal, rms=0.118540Jy sigma=0.391832 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.25 mas, bmaj=8.766 mas, bpa=88.12 degrees ! Estimated noise=5.2496 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00306866 Jy ! Component: 100 - total flux cleaned = -0.00433984 Jy ! Component: 150 - total flux cleaned = -0.00466184 Jy ! Component: 200 - total flux cleaned = -0.00468822 Jy ! Total flux subtracted in 200 components = -0.00468822 Jy ! Clean residual min=-0.010212 max=0.009040 Jy/beam ! Clean residual mean=-0.000005 rms=0.002321 Jy/beam ! Combined flux in latest and established models = 0.675528 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 61 components and 0.675528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.118503Jy sigma=0.391840 ! Fit after self-cal, rms=0.118492Jy sigma=0.391739 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.000452645 Jy ! Component: 100 - total flux cleaned = -0.000453532 Jy ! Component: 150 - total flux cleaned = -0.000402918 Jy ! Component: 200 - total flux cleaned = -0.000358601 Jy ! Total flux subtracted in 200 components = -0.000358601 Jy ! Clean residual min=-0.010168 max=0.009160 Jy/beam ! Clean residual mean=-0.000004 rms=0.002320 Jy/beam ! Combined flux in latest and established models = 0.675169 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 67 components and 0.675169 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.118475Jy sigma=0.391757 ! Fit after self-cal, rms=0.118475Jy sigma=0.391727 ! 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.451 mas, bmaj=6.635 mas, bpa=-70.63 degrees ! Estimated noise=4.37592 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.451 x 6.635 at -70.63 degrees (North through East) ! Clean map min=-0.0091905 max=0.61782 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. ! ! Fit before self-cal, rms=0.118475Jy sigma=0.391727 ! Fit after self-cal, rms=0.116243Jy sigma=0.363679 ! Inverting map and beam ! Estimated beam: bmin=3.463 mas, bmaj=6.632 mas, bpa=-70.64 degrees ! Estimated noise=4.37635 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00104751 Jy ! Component: 100 - total flux cleaned = -0.00129322 Jy ! Component: 150 - total flux cleaned = -0.00151209 Jy ! Component: 200 - total flux cleaned = -0.00165966 Jy ! Total flux subtracted in 200 components = -0.00165966 Jy ! Clean residual min=-0.009156 max=0.008655 Jy/beam ! Clean residual mean=-0.000003 rms=0.001806 Jy/beam ! Combined flux in latest and established models = 0.67351 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 73 components and 0.67351 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.116081Jy sigma=0.363328 ! Fit after self-cal, rms=0.116095Jy sigma=0.363322 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000145589 Jy ! Component: 100 - total flux cleaned = -0.000271734 Jy ! Component: 150 - total flux cleaned = -0.000310779 Jy ! Component: 200 - total flux cleaned = -0.000347774 Jy ! Total flux subtracted in 200 components = -0.000347774 Jy ! Clean residual min=-0.009169 max=0.008553 Jy/beam ! Clean residual mean=-0.000002 rms=0.001803 Jy/beam ! Combined flux in latest and established models = 0.673162 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 74 components and 0.673162 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.116004Jy sigma=0.363195 ! Fit after self-cal, rms=0.116013Jy sigma=0.363194 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.116013Jy sigma=0.363194 ! Fit after self-cal, rms=0.116015Jy sigma=0.362996 print "============== Clearing model and starting over ==========" ! ============== Clearing model and starting over ========== clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.853 mas, bmaj=6.509 mas, bpa=-70.29 degrees ! Estimated noise=5.05116 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.482139 Jy ! Component: 100 - total flux cleaned = 0.589147 Jy ! Component: 150 - total flux cleaned = 0.625658 Jy ! Component: 200 - total flux cleaned = 0.642131 Jy ! Total flux subtracted in 200 components = 0.642131 Jy ! Clean residual min=-0.009029 max=0.009181 Jy/beam ! Clean residual mean=0.000000 rms=0.002419 Jy/beam ! Combined flux in latest and established models = 0.642131 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.642131 Jy ! Inverting map and beam ! Estimated beam: bmin=3.469 mas, bmaj=6.629 mas, bpa=-70.65 degrees ! Estimated noise=4.3763 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.0119604 Jy ! Component: 100 - total flux cleaned = 0.0186266 Jy ! Component: 150 - total flux cleaned = 0.0228154 Jy ! Component: 200 - total flux cleaned = 0.0256896 Jy ! Total flux subtracted in 200 components = 0.0256896 Jy ! Clean residual min=-0.008806 max=0.008265 Jy/beam ! Clean residual mean=-0.000001 rms=0.001874 Jy/beam ! Combined flux in latest and established models = 0.66782 Jy ! Adding 33 model components to the UV plane model. ! The established model now contains 41 components and 0.66782 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.264 mas, bmaj=8.769 mas, bpa=87.93 degrees ! Estimated noise=5.23027 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00255402 Jy ! Component: 100 - total flux cleaned = 0.00297256 Jy ! Component: 150 - total flux cleaned = 0.00315396 Jy ! Component: 200 - total flux cleaned = 0.00327012 Jy ! Total flux subtracted in 200 components = 0.00327012 Jy ! Clean residual min=-0.009503 max=0.008828 Jy/beam ! Clean residual mean=-0.000004 rms=0.002274 Jy/beam ! Combined flux in latest and established models = 0.67109 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 = 6.8988e-05 Jy ! Component: 100 - total flux cleaned = 6.81582e-05 Jy ! Component: 150 - total flux cleaned = 9.87477e-05 Jy ! Component: 200 - total flux cleaned = 9.72003e-05 Jy ! Total flux subtracted in 200 components = 9.72003e-05 Jy ! Clean residual min=-0.009500 max=0.008768 Jy/beam ! Clean residual mean=-0.000004 rms=0.002266 Jy/beam ! Combined flux in latest and established models = 0.671188 Jy ! Adding 19 model components to the UV plane model. ! The established model now contains 56 components and 0.671187 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.469 mas, bmaj=6.629 mas, bpa=-70.65 degrees ! Estimated noise=4.3763 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000994757 Jy ! Component: 100 - total flux cleaned = 0.00141071 Jy ! Component: 150 - total flux cleaned = 0.00171566 Jy ! Component: 200 - total flux cleaned = 0.00202041 Jy ! Total flux subtracted in 200 components = 0.00202041 Jy ! Clean residual min=-0.008673 max=0.008647 Jy/beam ! Clean residual mean=-0.000002 rms=0.001830 Jy/beam ! Combined flux in latest and established models = 0.673208 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 65 components and 0.673208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.115311Jy sigma=0.363955 ! Fit after self-cal, rms=0.115425Jy sigma=0.363920 wmodel J2230+6946_S_l_map.mod ! Writing 65 model components to file: J2230+6946_S_l_map.mod wobs J2230+6946_S_l_uvs.fits ! Writing UV FITS file: J2230+6946_S_l_uvs.fits wwins J2230+6946_S_l_map.win ! wwins: Wrote 1 windows to J2230+6946_S_l_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.000284684 Jy ! Component: 100 - total flux cleaned = 0.00141807 Jy ! Component: 150 - total flux cleaned = 0.00252178 Jy ! Component: 200 - total flux cleaned = 0.00322126 Jy ! Component: 250 - total flux cleaned = 0.00387979 Jy ! Component: 300 - total flux cleaned = 0.00414514 Jy ! Component: 350 - total flux cleaned = 0.00431177 Jy ! Component: 400 - total flux cleaned = 0.00463188 Jy ! Component: 450 - total flux cleaned = 0.00501603 Jy ! Component: 500 - total flux cleaned = 0.00524057 Jy ! Component: 550 - total flux cleaned = 0.00538477 Jy ! Component: 600 - total flux cleaned = 0.00559123 Jy ! Component: 650 - total flux cleaned = 0.00578986 Jy ! Component: 700 - total flux cleaned = 0.00572585 Jy ! Component: 750 - total flux cleaned = 0.00591425 Jy ! Component: 800 - total flux cleaned = 0.00603509 Jy ! Component: 850 - total flux cleaned = 0.00609299 Jy ! Component: 900 - total flux cleaned = 0.0060927 Jy ! Component: 950 - total flux cleaned = 0.00592395 Jy ! Component: 1000 - total flux cleaned = 0.00586846 Jy ! Component: 1050 - total flux cleaned = 0.00565306 Jy ! Component: 1100 - total flux cleaned = 0.00544199 Jy ! Component: 1150 - total flux cleaned = 0.0052869 Jy ! Component: 1200 - total flux cleaned = 0.00503428 Jy ! Component: 1250 - total flux cleaned = 0.00488595 Jy ! Component: 1300 - total flux cleaned = 0.00483769 Jy ! Component: 1350 - total flux cleaned = 0.00464704 Jy ! Component: 1400 - total flux cleaned = 0.004694 Jy ! Component: 1450 - total flux cleaned = 0.00450987 Jy ! Component: 1500 - total flux cleaned = 0.00428407 Jy ! Component: 1550 - total flux cleaned = 0.00415112 Jy ! Component: 1600 - total flux cleaned = 0.00410735 Jy ! Component: 1650 - total flux cleaned = 0.00389105 Jy ! Component: 1700 - total flux cleaned = 0.00389134 Jy ! Component: 1750 - total flux cleaned = 0.00384933 Jy ! Component: 1800 - total flux cleaned = 0.00368404 Jy ! Component: 1850 - total flux cleaned = 0.0037652 Jy ! Component: 1900 - total flux cleaned = 0.00356382 Jy ! Component: 1950 - total flux cleaned = 0.0035644 Jy ! Component: 2000 - total flux cleaned = 0.00356443 Jy ! Component: 2050 - total flux cleaned = 0.00352624 Jy ! Component: 2100 - total flux cleaned = 0.00352607 Jy ! Component: 2150 - total flux cleaned = 0.00348842 Jy ! Component: 2200 - total flux cleaned = 0.00348854 Jy ! Component: 2250 - total flux cleaned = 0.0035249 Jy ! Component: 2300 - total flux cleaned = 0.00338021 Jy ! Component: 2350 - total flux cleaned = 0.0033089 Jy ! Component: 2400 - total flux cleaned = 0.0032026 Jy ! Component: 2450 - total flux cleaned = 0.00302838 Jy ! Component: 2500 - total flux cleaned = 0.00295917 Jy ! Component: 2550 - total flux cleaned = 0.00285705 Jy ! Component: 2600 - total flux cleaned = 0.00275591 Jy ! Component: 2650 - total flux cleaned = 0.00265577 Jy ! Component: 2700 - total flux cleaned = 0.0026226 Jy ! Component: 2750 - total flux cleaned = 0.00255673 Jy ! Component: 2800 - total flux cleaned = 0.002427 Jy ! Component: 2850 - total flux cleaned = 0.00226604 Jy ! Component: 2900 - total flux cleaned = 0.00220217 Jy ! Component: 2950 - total flux cleaned = 0.00204461 Jy ! Component: 3000 - total flux cleaned = 0.00191992 Jy ! Component: 3050 - total flux cleaned = 0.00182755 Jy ! Component: 3100 - total flux cleaned = 0.00179667 Jy ! Component: 3150 - total flux cleaned = 0.0016445 Jy ! Component: 3200 - total flux cleaned = 0.00158441 Jy ! Component: 3250 - total flux cleaned = 0.00143517 Jy ! Component: 3300 - total flux cleaned = 0.00140538 Jy ! Component: 3350 - total flux cleaned = 0.00122919 Jy ! Component: 3400 - total flux cleaned = 0.00119997 Jy ! Component: 3450 - total flux cleaned = 0.00105548 Jy ! Component: 3500 - total flux cleaned = 0.00102714 Jy ! Component: 3550 - total flux cleaned = 0.000970039 Jy ! Component: 3600 - total flux cleaned = 0.00088525 Jy ! Component: 3650 - total flux cleaned = 0.000801187 Jy ! Component: 3700 - total flux cleaned = 0.000717953 Jy ! Component: 3750 - total flux cleaned = 0.000607579 Jy ! Component: 3800 - total flux cleaned = 0.000552933 Jy ! Component: 3850 - total flux cleaned = 0.000552771 Jy ! Component: 3900 - total flux cleaned = 0.000363708 Jy ! Component: 3950 - total flux cleaned = 0.000390556 Jy ! Component: 4000 - total flux cleaned = 0.000284145 Jy ! Component: 4050 - total flux cleaned = 0.00031103 Jy ! Total flux subtracted in 4096 components = 0.000337353 Jy ! Clean residual min=-0.001847 max=0.001611 Jy/beam ! Clean residual mean=-0.000001 rms=0.000463 Jy/beam ! Combined flux in latest and established models = 0.673545 Jy keep ! Adding 275 model components to the UV plane model. ! The established model now contains 340 components and 0.673545 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J2230+6946_S_l_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.469 x 6.629 at -70.65 degrees (North through East) ! Clean map min=-0.0045626 max=0.6219 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J2230+6946_S_l_map.fits wmodel /image_prod/ug002_uva//J2230+6946_S_l_map.mod ! Writing 340 model components to file: /image_prod/ug002_uva//J2230+6946_S_l_map.mod wobs /image_prod/ug002_uva//J2230+6946_S_l_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J2230+6946_S_l_uvs.fits wwins /image_prod/ug002_uva//J2230+6946_S_l_map.win ! wwins: Wrote 2 windows to /image_prod/ug002_uva//J2230+6946_S_l_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.621903 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000433035 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.208892 0.208892 0.417784 0.835567 1.67113 3.34227 6.68454 13.3691 26.7382 53.4763 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0012 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.621 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1436.15 quit ! Quitting program ! Log file /image_prod/ug002_uva//J2230+6946_S_l_dfm.log closed on Thu Mar 5 05:44:52 2020