! Started logfile: /vlbi/uf001/uvm//J1823-2726_S_a_dfm.log on Sun Jan 13 22:14:11 2019 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 /vlbi/uf001/uvm//J1823-2726_S_a_uvm.fits ! Reading UV FITS file: /vlbi/uf001/uvm//J1823-2726_S_a_uvm.fits ! AN table 1: 36 integrations on 21 of 21 possible baselines. ! AN table 2: 37 integrations on 21 of 21 possible baselines. ! AN table 3: 148 integrations on 28 of 28 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J1823-2726 ! ! 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 44 lines of history. ! ! Reading 17031 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 J1823-2726_S_a] 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 7 telescope corrections were flagged in sub-array 2. ! A total of 32 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 32 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 32 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.065698Jy sigma=15.093034 ! Fit after self-cal, rms=0.831964Jy sigma=11.829966 ! 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=4.979 mas, bmaj=24.36 mas, bpa=-0.5448 degrees ! Estimated noise=0.919408 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.15322 Jy ! Component: 100 - total flux cleaned = 0.237361 Jy ! Total flux subtracted in 100 components = 0.237361 Jy ! Clean residual min=-0.028567 max=0.063173 Jy/beam ! Clean residual mean=-0.001709 rms=0.015606 Jy/beam ! Combined flux in latest and established models = 0.237361 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 61 components and 0.237361 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.261106Jy sigma=3.233853 ! Fit after self-cal, rms=0.259973Jy sigma=3.221181 ! 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=6.073 mas, bmaj=27.06 mas, bpa=-0.486 degrees ! Estimated noise=0.394124 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.0825306 Jy ! Component: 100 - total flux cleaned = 0.138202 Jy ! Component: 150 - total flux cleaned = 0.178433 Jy ! Component: 200 - total flux cleaned = 0.2078 Jy ! Total flux subtracted in 200 components = 0.2078 Jy ! Clean residual min=-0.034960 max=0.079238 Jy/beam ! Clean residual mean=-0.001624 rms=0.017488 Jy/beam ! Combined flux in latest and established models = 0.445161 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 144 components and 0.445161 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.219390Jy sigma=2.753030 ! Fit after self-cal, rms=0.217866Jy sigma=2.739515 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 6.073 x 27.06 at -0.486 degrees (North through East) ! Clean map min=-0.034716 max=0.1945 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.49 FD 1.29 HN 1.72 KP 1.02 ! LA 0.75 NL 2.14 OV 1.42 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.51 FD 2.43 KP 0.78 LA 0.63 ! NL 2.74 OV 1.47 PT 0.54 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.28 FD 1.61 HN 1.36 KP 0.88 ! LA 0.61 NL 2.49 OV 1.91 PT 0.63 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.68 FD 1.36 HN 1.45 KP 1.01 ! LA 0.68 NL 2.45 OV 1.57 ! ! Telescope amplitude corrections in sub-array 2: ! BR 2.00 FD 2.05 KP 0.74 LA 0.66 ! NL 2.60 OV 1.43 PT 0.59 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.75 HN 1.23 KP 0.90 ! LA 0.62 NL 2.58 OV 1.67 PT 0.63 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.43 FD 1.22 HN 1.73 KP 1.13 ! LA 0.80 NL 1.68 OV 1.24 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.49 FD 1.85 KP 0.79 LA 0.60 ! NL 2.47 OV 1.43 PT 0.57 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.41 FD 1.74 HN 1.41 KP 0.89 ! LA 0.62 NL 2.39 OV 1.84 PT 0.65 ! ! ! Fit before self-cal, rms=0.217866Jy sigma=2.739515 ! Fit after self-cal, rms=0.270081Jy sigma=1.792672 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=6.868 mas, bmaj=29.52 mas, bpa=0.246 degrees ! Estimated noise=0.483472 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.00418384 Jy ! Component: 100 - total flux cleaned = 0.00212288 Jy ! Component: 150 - total flux cleaned = 0.000503287 Jy ! Component: 200 - total flux cleaned = -0.000944975 Jy ! Total flux subtracted in 200 components = -0.000944975 Jy ! Clean residual min=-0.026332 max=0.043073 Jy/beam ! Clean residual mean=-0.000279 rms=0.009730 Jy/beam ! Combined flux in latest and established models = 0.444216 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 157 components and 0.444216 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.266326Jy sigma=1.795525 ! Fit after self-cal, rms=0.271099Jy sigma=1.768935 ! 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.271099Jy sigma=1.768935 ! Fit after self-cal, rms=0.294220Jy sigma=1.762226 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.363 mas, bmaj=30.15 mas, bpa=-0.3345 degrees ! Estimated noise=0.484032 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00413863 Jy ! Component: 100 - total flux cleaned = 0.00314584 Jy ! Component: 150 - total flux cleaned = 0.00231921 Jy ! Component: 200 - total flux cleaned = 0.00182706 Jy ! Total flux subtracted in 200 components = 0.00182706 Jy ! Clean residual min=-0.024279 max=0.039086 Jy/beam ! Clean residual mean=-0.000251 rms=0.009276 Jy/beam ! Combined flux in latest and established models = 0.446043 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 172 components and 0.446043 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.291561Jy sigma=1.766222 ! Fit after self-cal, rms=0.291826Jy sigma=1.765382 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.00153667 Jy ! Component: 100 - total flux cleaned = -0.00240687 Jy ! Component: 150 - total flux cleaned = -0.00320259 Jy ! Component: 200 - total flux cleaned = -0.00376563 Jy ! Total flux subtracted in 200 components = -0.00376563 Jy ! Clean residual min=-0.024305 max=0.039565 Jy/beam ! Clean residual mean=-0.000272 rms=0.009291 Jy/beam ! Combined flux in latest and established models = 0.442277 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 181 components and 0.442277 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.290732Jy sigma=1.766774 ! Fit after self-cal, rms=0.290789Jy sigma=1.766497 ! 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.290789Jy sigma=1.766497 ! Fit after self-cal, rms=0.366914Jy sigma=1.678631 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.25 mas, bmaj=32.08 mas, bpa=-0.7475 degrees ! Estimated noise=0.507663 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00267393 Jy ! Component: 100 - total flux cleaned = 0.000324874 Jy ! Component: 150 - total flux cleaned = -0.00135032 Jy ! Component: 200 - total flux cleaned = -0.00253581 Jy ! Total flux subtracted in 200 components = -0.00253581 Jy ! Clean residual min=-0.025981 max=0.038131 Jy/beam ! Clean residual mean=-0.000239 rms=0.008743 Jy/beam ! Combined flux in latest and established models = 0.439741 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 188 components and 0.439741 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.363978Jy sigma=1.679154 ! Fit after self-cal, rms=0.364480Jy sigma=1.677620 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (18, -19). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0173281 Jy ! Component: 100 - total flux cleaned = 0.0129577 Jy ! Component: 150 - total flux cleaned = 0.0094141 Jy ! Component: 200 - total flux cleaned = 0.00644133 Jy ! Total flux subtracted in 200 components = 0.00644133 Jy ! Clean residual min=-0.026539 max=0.029689 Jy/beam ! Clean residual mean=-0.000167 rms=0.007946 Jy/beam ! Combined flux in latest and established models = 0.446183 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 235 components and 0.446183 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.354391Jy sigma=1.621097 ! Fit after self-cal, rms=0.351913Jy sigma=1.602656 ! 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.351913Jy sigma=1.602656 ! Fit after self-cal, rms=0.339340Jy sigma=1.575116 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.556 mas, bmaj=34.2 mas, bpa=1.546 degrees ! Estimated noise=0.506133 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.016817 Jy ! Component: 100 - total flux cleaned = 0.0142914 Jy ! Component: 150 - total flux cleaned = 0.0121714 Jy ! Component: 200 - total flux cleaned = 0.0107833 Jy ! Total flux subtracted in 200 components = 0.0107833 Jy ! Clean residual min=-0.030317 max=0.031859 Jy/beam ! Clean residual mean=-0.000194 rms=0.008838 Jy/beam ! Combined flux in latest and established models = 0.456966 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 279 components and 0.456966 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.334978Jy sigma=1.541957 ! Fit after self-cal, rms=0.334244Jy sigma=1.533842 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.00270773 Jy ! Component: 100 - total flux cleaned = 0.000600226 Jy ! Component: 150 - total flux cleaned = -0.000953964 Jy ! Component: 200 - total flux cleaned = -0.00241576 Jy ! Total flux subtracted in 200 components = -0.00241576 Jy ! Clean residual min=-0.030651 max=0.032243 Jy/beam ! Clean residual mean=-0.000194 rms=0.008706 Jy/beam ! Combined flux in latest and established models = 0.45455 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 318 components and 0.45455 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.332187Jy sigma=1.528354 ! Fit after self-cal, rms=0.331901Jy sigma=1.526720 ! 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=6.442 mas, bmaj=32.39 mas, bpa=-0.6097 degrees ! Estimated noise=0.512724 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 6.442 x 32.39 at -0.6097 degrees (North through East) ! Clean map min=-0.032745 max=0.25948 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=5.081 mas, bmaj=29.23 mas, bpa=-0.5575 degrees ! Estimated noise=1.24984 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.203457 Jy ! Component: 100 - total flux cleaned = 0.299216 Jy ! Component: 150 - total flux cleaned = 0.359058 Jy ! Component: 200 - total flux cleaned = 0.400713 Jy ! Total flux subtracted in 200 components = 0.400713 Jy ! Clean residual min=-0.025687 max=0.030798 Jy/beam ! Clean residual mean=-0.000557 rms=0.008375 Jy/beam ! Combined flux in latest and established models = 0.400713 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 79 model components to the UV plane model. ! The established model now contains 79 components and 0.400713 Jy ! Inverting map and beam ! Estimated beam: bmin=6.442 mas, bmaj=32.39 mas, bpa=-0.6097 degrees ! Estimated noise=0.512724 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.0274498 Jy ! Component: 100 - total flux cleaned = 0.033127 Jy ! Component: 150 - total flux cleaned = 0.0359511 Jy ! Component: 200 - total flux cleaned = 0.0383693 Jy ! Total flux subtracted in 200 components = 0.0383693 Jy ! Clean residual min=-0.026234 max=0.024574 Jy/beam ! Clean residual mean=-0.000292 rms=0.007693 Jy/beam ! Combined flux in latest and established models = 0.439082 Jy ! Adding 71 model components to the UV plane model. ! The established model now contains 142 components and 0.439082 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=9.556 mas, bmaj=34.2 mas, bpa=1.546 degrees ! Estimated noise=0.506133 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00555064 Jy ! Component: 100 - total flux cleaned = 0.00702135 Jy ! Component: 150 - total flux cleaned = 0.00760895 Jy ! Component: 200 - total flux cleaned = 0.00713146 Jy ! Total flux subtracted in 200 components = 0.00713146 Jy ! Clean residual min=-0.031706 max=0.031757 Jy/beam ! Clean residual mean=-0.000278 rms=0.008941 Jy/beam ! Combined flux in latest and established models = 0.446214 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.000413882 Jy ! Component: 100 - total flux cleaned = -0.000415867 Jy ! Component: 150 - total flux cleaned = -3.93931e-05 Jy ! Component: 200 - total flux cleaned = 0.000288065 Jy ! Total flux subtracted in 200 components = 0.000288065 Jy ! Clean residual min=-0.031718 max=0.032282 Jy/beam ! Clean residual mean=-0.000260 rms=0.008866 Jy/beam ! Combined flux in latest and established models = 0.446502 Jy ! Adding 79 model components to the UV plane model. ! The established model now contains 208 components and 0.446502 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=6.442 mas, bmaj=32.39 mas, bpa=-0.6097 degrees ! Estimated noise=0.512724 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000378303 Jy ! Component: 100 - total flux cleaned = -0.000208136 Jy ! Component: 150 - total flux cleaned = -0.000482557 Jy ! Component: 200 - total flux cleaned = -0.0014873 Jy ! Total flux subtracted in 200 components = -0.0014873 Jy ! Clean residual min=-0.026201 max=0.027771 Jy/beam ! Clean residual mean=-0.000199 rms=0.007324 Jy/beam ! Combined flux in latest and established models = 0.445015 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 237 components and 0.445014 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.329609Jy sigma=1.530563 ! Fit after self-cal, rms=0.328951Jy sigma=1.525897 wmodel J1823-2726_S_a_map.mod ! Writing 237 model components to file: J1823-2726_S_a_map.mod wobs J1823-2726_S_a_uvs.fits ! Writing UV FITS file: J1823-2726_S_a_uvs.fits wwins J1823-2726_S_a_map.win ! wwins: Wrote 2 windows to J1823-2726_S_a_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.00411694 Jy ! Component: 100 - total flux cleaned = 0.00927087 Jy ! Component: 150 - total flux cleaned = 0.0107303 Jy ! Component: 200 - total flux cleaned = 0.0100746 Jy ! Component: 250 - total flux cleaned = 0.00761832 Jy ! Component: 300 - total flux cleaned = 0.00616389 Jy ! Component: 350 - total flux cleaned = 0.00476818 Jy ! Component: 400 - total flux cleaned = 0.0039609 Jy ! Component: 450 - total flux cleaned = 0.00317774 Jy ! Component: 500 - total flux cleaned = 0.00292323 Jy ! Component: 550 - total flux cleaned = 0.00292539 Jy ! Component: 600 - total flux cleaned = 0.00364108 Jy ! Component: 650 - total flux cleaned = 0.00340968 Jy ! Component: 700 - total flux cleaned = 0.00386991 Jy ! Component: 750 - total flux cleaned = 0.00431154 Jy ! Component: 800 - total flux cleaned = 0.00474308 Jy ! Component: 850 - total flux cleaned = 0.00474269 Jy ! Component: 900 - total flux cleaned = 0.00516344 Jy ! Component: 950 - total flux cleaned = 0.0061866 Jy ! Component: 1000 - total flux cleaned = 0.00658533 Jy ! Component: 1050 - total flux cleaned = 0.00638633 Jy ! Component: 1100 - total flux cleaned = 0.00658269 Jy ! Component: 1150 - total flux cleaned = 0.00734339 Jy ! Component: 1200 - total flux cleaned = 0.00733971 Jy ! Component: 1250 - total flux cleaned = 0.00752138 Jy ! Component: 1300 - total flux cleaned = 0.00787932 Jy ! Component: 1350 - total flux cleaned = 0.00734469 Jy ! Component: 1400 - total flux cleaned = 0.00804945 Jy ! Component: 1450 - total flux cleaned = 0.00770084 Jy ! Component: 1500 - total flux cleaned = 0.00735821 Jy ! Component: 1550 - total flux cleaned = 0.0066854 Jy ! Component: 1600 - total flux cleaned = 0.00684722 Jy ! Component: 1650 - total flux cleaned = 0.00602472 Jy ! Component: 1700 - total flux cleaned = 0.005219 Jy ! Component: 1750 - total flux cleaned = 0.00505964 Jy ! Component: 1800 - total flux cleaned = 0.00521768 Jy ! Component: 1850 - total flux cleaned = 0.00521818 Jy ! Component: 1900 - total flux cleaned = 0.00506818 Jy ! Component: 1950 - total flux cleaned = 0.00537117 Jy ! Component: 2000 - total flux cleaned = 0.00552009 Jy ! Component: 2050 - total flux cleaned = 0.00596129 Jy ! Component: 2100 - total flux cleaned = 0.00610649 Jy ! Component: 2150 - total flux cleaned = 0.00639714 Jy ! Component: 2200 - total flux cleaned = 0.0066834 Jy ! Component: 2250 - total flux cleaned = 0.00724989 Jy ! Component: 2300 - total flux cleaned = 0.00711157 Jy ! Component: 2350 - total flux cleaned = 0.00752613 Jy ! Component: 2400 - total flux cleaned = 0.00738828 Jy ! Component: 2450 - total flux cleaned = 0.00738912 Jy ! Component: 2500 - total flux cleaned = 0.00712007 Jy ! Component: 2550 - total flux cleaned = 0.00765125 Jy ! Component: 2600 - total flux cleaned = 0.00712615 Jy ! Component: 2650 - total flux cleaned = 0.00712822 Jy ! Component: 2700 - total flux cleaned = 0.00725656 Jy ! Component: 2750 - total flux cleaned = 0.00751088 Jy ! Component: 2800 - total flux cleaned = 0.0076371 Jy ! Component: 2850 - total flux cleaned = 0.00738693 Jy ! Component: 2900 - total flux cleaned = 0.00726395 Jy ! Component: 2950 - total flux cleaned = 0.0080038 Jy ! Component: 3000 - total flux cleaned = 0.00824699 Jy ! Component: 3050 - total flux cleaned = 0.00824714 Jy ! Component: 3100 - total flux cleaned = 0.00848635 Jy ! Component: 3150 - total flux cleaned = 0.0090806 Jy ! Component: 3200 - total flux cleaned = 0.00896502 Jy ! Component: 3250 - total flux cleaned = 0.00978241 Jy ! Component: 3300 - total flux cleaned = 0.00978229 Jy ! Component: 3350 - total flux cleaned = 0.00978297 Jy ! Component: 3400 - total flux cleaned = 0.0103521 Jy ! Component: 3450 - total flux cleaned = 0.0108037 Jy ! Component: 3500 - total flux cleaned = 0.010914 Jy ! Component: 3550 - total flux cleaned = 0.0111355 Jy ! Component: 3600 - total flux cleaned = 0.0111364 Jy ! Component: 3650 - total flux cleaned = 0.0113529 Jy ! Component: 3700 - total flux cleaned = 0.0113524 Jy ! Component: 3750 - total flux cleaned = 0.0114597 Jy ! Component: 3800 - total flux cleaned = 0.0119947 Jy ! Component: 3850 - total flux cleaned = 0.0122056 Jy ! Component: 3900 - total flux cleaned = 0.012205 Jy ! Component: 3950 - total flux cleaned = 0.0125161 Jy ! Component: 4000 - total flux cleaned = 0.0124125 Jy ! Component: 4050 - total flux cleaned = 0.0125149 Jy ! Total flux subtracted in 4096 components = 0.0126172 Jy ! Clean residual min=-0.006062 max=0.005478 Jy/beam ! Clean residual mean=-0.000038 rms=0.002548 Jy/beam ! Combined flux in latest and established models = 0.457632 Jy keep ! Adding 537 model components to the UV plane model. ! The established model now contains 774 components and 0.457631 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /vlbi/uf001/uvm//J1823-2726_S_a_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 6.442 x 32.39 at -0.6097 degrees (North through East) ! Clean map min=-0.034825 max=0.26922 Jy/beam ! Writing clean map to FITS file: /vlbi/uf001/uvm//J1823-2726_S_a_map.fits wmodel /vlbi/uf001/uvm//J1823-2726_S_a_map.mod ! Writing 774 model components to file: /vlbi/uf001/uvm//J1823-2726_S_a_map.mod wobs /vlbi/uf001/uvm//J1823-2726_S_a_uvs.fits ! Writing UV FITS file: /vlbi/uf001/uvm//J1823-2726_S_a_uvs.fits wwins /vlbi/uf001/uvm//J1823-2726_S_a_map.win ! wwins: Wrote 3 windows to /vlbi/uf001/uvm//J1823-2726_S_a_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.26922 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00254376 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -2.83459 2.83459 5.66918 11.3384 22.6767 45.3535 90.7069 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0076 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.269 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 105.835 quit ! Quitting program ! Log file /vlbi/uf001/uvm//J1823-2726_S_a_dfm.log closed on Sun Jan 13 22:14:17 2019