! Started logfile: /vlbi/uf001/uva//J0341+3352_S_dfm.log on Sun Jan 13 21:53:13 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/uva//J0341+3352_S_uva.fits ! Reading UV FITS file: /vlbi/uf001/uva//J0341+3352_S_uva.fits ! AN table 1: 42 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.857143 visibilities/baseline/integration-bin. ! Found source: J0341+3352 ! ! 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 40 lines of history. ! ! Reading 3888 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 J0341+3352_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 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 78 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 78 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 78 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.007102Jy sigma=21.621574 ! Fit after self-cal, rms=0.924677Jy sigma=19.850258 ! 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.724 mas, bmaj=5.419 mas, bpa=-3.193 degrees ! Estimated noise=0.790039 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.0575708 Jy ! Component: 100 - total flux cleaned = 0.0752252 Jy ! Total flux subtracted in 100 components = 0.0752252 Jy ! Clean residual min=-0.004650 max=0.007664 Jy/beam ! Clean residual mean=0.000013 rms=0.001328 Jy/beam ! Combined flux in latest and established models = 0.0752252 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 12 components and 0.0752252 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.063003Jy sigma=1.250787 ! Fit after self-cal, rms=0.062579Jy sigma=1.242467 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00911076 Jy ! Component: 100 - total flux cleaned = 0.0142825 Jy ! Total flux subtracted in 100 components = 0.0142825 Jy ! Clean residual min=-0.003618 max=0.004465 Jy/beam ! Clean residual mean=0.000011 rms=0.001009 Jy/beam ! Combined flux in latest and established models = 0.0895076 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 29 components and 0.0895076 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.061179Jy sigma=1.208067 ! Fit after self-cal, rms=0.061153Jy sigma=1.207720 ! 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.194 mas, bmaj=6.134 mas, bpa=-2.629 degrees ! Estimated noise=0.582512 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.00378594 Jy ! Component: 100 - total flux cleaned = 0.0059127 Jy ! Component: 150 - total flux cleaned = 0.00721169 Jy ! Component: 200 - total flux cleaned = 0.00767485 Jy ! Total flux subtracted in 200 components = 0.00767485 Jy ! Clean residual min=-0.003763 max=0.003413 Jy/beam ! Clean residual mean=0.000015 rms=0.000820 Jy/beam ! Combined flux in latest and established models = 0.0971825 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 45 components and 0.0971825 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060931Jy sigma=1.201880 ! Fit after self-cal, rms=0.060924Jy sigma=1.201791 ! 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.194 x 6.134 at -2.629 degrees (North through East) ! Clean map min=-0.0036752 max=0.07525 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.77 FD 0.98 HN 1.05 LA 0.91 ! MK 1.06 NL 1.16 OV 0.95 PT 0.99 ! SC 1.06 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 0.97 HN 0.97 LA 0.96 ! MK 1.00 NL 1.07 OV 0.99 PT 1.04 ! SC 0.99 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.00 HN 1.03 LA 0.96 ! MK 0.95 NL 1.00 OV 0.97 PT 1.05 ! SC 1.04 ! ! ! Fit before self-cal, rms=0.060924Jy sigma=1.201791 ! Fit after self-cal, rms=0.059884Jy sigma=1.191913 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.209 mas, bmaj=6.153 mas, bpa=-2.611 degrees ! Estimated noise=0.580251 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.000619007 Jy ! Component: 100 - total flux cleaned = -0.000639171 Jy ! Component: 150 - total flux cleaned = -0.000639699 Jy ! Component: 200 - total flux cleaned = -0.00065663 Jy ! Total flux subtracted in 200 components = -0.00065663 Jy ! Clean residual min=-0.003519 max=0.002898 Jy/beam ! Clean residual mean=0.000018 rms=0.000735 Jy/beam ! Combined flux in latest and established models = 0.0965258 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 50 components and 0.0965258 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059868Jy sigma=1.191658 ! Fit after self-cal, rms=0.059900Jy sigma=1.191327 ! 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.059900Jy sigma=1.191327 ! Fit after self-cal, rms=0.060014Jy sigma=1.183106 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.202 mas, bmaj=6.14 mas, bpa=-2.88 degrees ! Estimated noise=0.57837 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000191058 Jy ! Component: 100 - total flux cleaned = 0.000172164 Jy ! Component: 150 - total flux cleaned = 0.000172099 Jy ! Component: 200 - total flux cleaned = 0.000172075 Jy ! Total flux subtracted in 200 components = 0.000172075 Jy ! Clean residual min=-0.002738 max=0.002647 Jy/beam ! Clean residual mean=0.000009 rms=0.000635 Jy/beam ! Combined flux in latest and established models = 0.0966979 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 54 components and 0.0966979 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060005Jy sigma=1.182920 ! Fit after self-cal, rms=0.060054Jy sigma=1.182709 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 = 2.70632e-05 Jy ! Component: 100 - total flux cleaned = 2.69138e-05 Jy ! Component: 150 - total flux cleaned = 3.82684e-05 Jy ! Component: 200 - total flux cleaned = 3.83758e-05 Jy ! Total flux subtracted in 200 components = 3.83758e-05 Jy ! Clean residual min=-0.002756 max=0.002623 Jy/beam ! Clean residual mean=0.000009 rms=0.000634 Jy/beam ! Combined flux in latest and established models = 0.0967363 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 56 components and 0.0967363 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060050Jy sigma=1.182616 ! Fit after self-cal, rms=0.060049Jy sigma=1.182612 ! 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.060049Jy sigma=1.182612 ! Fit after self-cal, rms=0.060192Jy sigma=1.182559 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.202 mas, bmaj=6.152 mas, bpa=-2.915 degrees ! Estimated noise=0.57849 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -2.31039e-05 Jy ! Component: 100 - total flux cleaned = -2.31367e-05 Jy ! Component: 150 - total flux cleaned = -2.29879e-05 Jy ! Component: 200 - total flux cleaned = -1.33363e-05 Jy ! Total flux subtracted in 200 components = -1.33363e-05 Jy ! Clean residual min=-0.002744 max=0.002650 Jy/beam ! Clean residual mean=0.000009 rms=0.000632 Jy/beam ! Combined flux in latest and established models = 0.0967229 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 56 components and 0.0967229 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060188Jy sigma=1.182474 ! Fit after self-cal, rms=0.060191Jy sigma=1.182469 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.52962e-06 Jy ! Component: 100 - total flux cleaned = 7.49701e-06 Jy ! Component: 150 - total flux cleaned = 7.47284e-06 Jy ! Component: 200 - total flux cleaned = 1.62295e-05 Jy ! Total flux subtracted in 200 components = 1.62295e-05 Jy ! Clean residual min=-0.002756 max=0.002657 Jy/beam ! Clean residual mean=0.000009 rms=0.000632 Jy/beam ! Combined flux in latest and established models = 0.0967392 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 57 components and 0.0967392 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060189Jy sigma=1.182409 ! Fit after self-cal, rms=0.060187Jy sigma=1.182404 ! 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.060187Jy sigma=1.182404 ! Fit after self-cal, rms=0.060299Jy sigma=1.182392 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.994 mas, bmaj=9.105 mas, bpa=-5.047 degrees ! Estimated noise=0.743055 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000583238 Jy ! Component: 100 - total flux cleaned = 0.000670164 Jy ! Component: 150 - total flux cleaned = 0.000721827 Jy ! Component: 200 - total flux cleaned = 0.000764994 Jy ! Total flux subtracted in 200 components = 0.000764994 Jy ! Clean residual min=-0.002865 max=0.003224 Jy/beam ! Clean residual mean=0.000020 rms=0.000760 Jy/beam ! Combined flux in latest and established models = 0.0975042 Jy selfcal ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 60 components and 0.0975042 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060310Jy sigma=1.182744 ! Fit after self-cal, rms=0.060299Jy sigma=1.182596 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 = 4.86185e-05 Jy ! Component: 100 - total flux cleaned = 4.74276e-05 Jy ! Component: 150 - total flux cleaned = 4.80813e-05 Jy ! Component: 200 - total flux cleaned = 4.75028e-05 Jy ! Total flux subtracted in 200 components = 4.75028e-05 Jy ! Clean residual min=-0.002860 max=0.003216 Jy/beam ! Clean residual mean=0.000020 rms=0.000757 Jy/beam ! Combined flux in latest and established models = 0.0975517 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 65 components and 0.0975517 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060309Jy sigma=1.182886 ! Fit after self-cal, rms=0.060304Jy sigma=1.182808 ! 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.203 mas, bmaj=6.162 mas, bpa=-2.917 degrees ! Estimated noise=0.57906 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.203 x 6.162 at -2.917 degrees (North through East) ! Clean map min=-0.002612 max=0.074724 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=2.737 mas, bmaj=5.454 mas, bpa=-3.246 degrees ! Estimated noise=0.783978 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0572606 Jy ! Component: 100 - total flux cleaned = 0.0753224 Jy ! Component: 150 - total flux cleaned = 0.0841423 Jy ! Component: 200 - total flux cleaned = 0.0892121 Jy ! Total flux subtracted in 200 components = 0.0892121 Jy ! Clean residual min=-0.003717 max=0.003936 Jy/beam ! Clean residual mean=0.000007 rms=0.000895 Jy/beam ! Combined flux in latest and established models = 0.0892121 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 25 model components to the UV plane model. ! The established model now contains 25 components and 0.0892121 Jy ! Inverting map and beam ! Estimated beam: bmin=3.203 mas, bmaj=6.162 mas, bpa=-2.917 degrees ! Estimated noise=0.57906 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.0039089 Jy ! Component: 100 - total flux cleaned = 0.00589099 Jy ! Component: 150 - total flux cleaned = 0.00632166 Jy ! Component: 200 - total flux cleaned = 0.00665251 Jy ! Total flux subtracted in 200 components = 0.00665251 Jy ! Clean residual min=-0.002785 max=0.002718 Jy/beam ! Clean residual mean=0.000009 rms=0.000647 Jy/beam ! Combined flux in latest and established models = 0.0958646 Jy ! Adding 21 model components to the UV plane model. ! The established model now contains 38 components and 0.0958646 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=6.994 mas, bmaj=9.105 mas, bpa=-5.047 degrees ! Estimated noise=0.743055 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00104648 Jy ! Component: 100 - total flux cleaned = 0.00110744 Jy ! Component: 150 - total flux cleaned = 0.00117443 Jy ! Component: 200 - total flux cleaned = 0.00119351 Jy ! Total flux subtracted in 200 components = 0.00119351 Jy ! Clean residual min=-0.002920 max=0.003155 Jy/beam ! Clean residual mean=0.000020 rms=0.000765 Jy/beam ! Combined flux in latest and established models = 0.0970581 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 = 4.82885e-05 Jy ! Component: 100 - total flux cleaned = 7.68531e-05 Jy ! Component: 150 - total flux cleaned = 0.000115545 Jy ! Component: 200 - total flux cleaned = 0.00013935 Jy ! Total flux subtracted in 200 components = 0.00013935 Jy ! Clean residual min=-0.002863 max=0.003168 Jy/beam ! Clean residual mean=0.000020 rms=0.000759 Jy/beam ! Combined flux in latest and established models = 0.0971975 Jy ! Adding 22 model components to the UV plane model. ! The established model now contains 52 components and 0.0971975 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.203 mas, bmaj=6.162 mas, bpa=-2.917 degrees ! Estimated noise=0.57906 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000451612 Jy ! Component: 100 - total flux cleaned = -0.000667586 Jy ! Component: 150 - total flux cleaned = -0.000669554 Jy ! Component: 200 - total flux cleaned = -0.000669595 Jy ! Total flux subtracted in 200 components = -0.000669595 Jy ! Clean residual min=-0.002795 max=0.002727 Jy/beam ! Clean residual mean=0.000008 rms=0.000635 Jy/beam ! Combined flux in latest and established models = 0.0965279 Jy selfcal ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 59 components and 0.0965279 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060300Jy sigma=1.182779 ! Fit after self-cal, rms=0.060299Jy sigma=1.182713 wmodel J0341+3352_S_map.mod ! Writing 59 model components to file: J0341+3352_S_map.mod wobs J0341+3352_S_uvs.fits ! Writing UV FITS file: J0341+3352_S_uvs.fits wwins J0341+3352_S_map.win ! wwins: Wrote 1 windows to J0341+3352_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.00010139 Jy ! Component: 100 - total flux cleaned = 0.000284708 Jy ! Component: 150 - total flux cleaned = 0.000329816 Jy ! Component: 200 - total flux cleaned = 0.00012599 Jy ! Component: 250 - total flux cleaned = 4.7779e-05 Jy ! Component: 300 - total flux cleaned = -0.000102975 Jy ! Component: 350 - total flux cleaned = -0.000211864 Jy ! Component: 400 - total flux cleaned = -0.000141714 Jy ! Component: 450 - total flux cleaned = -0.000107489 Jy ! Component: 500 - total flux cleaned = -4.08548e-05 Jy ! Component: 550 - total flux cleaned = -7.35514e-05 Jy ! Component: 600 - total flux cleaned = -0.000136951 Jy ! Component: 650 - total flux cleaned = -0.000230434 Jy ! Component: 700 - total flux cleaned = -4.85119e-05 Jy ! Component: 750 - total flux cleaned = -4.88562e-05 Jy ! Component: 800 - total flux cleaned = 9.15918e-06 Jy ! Component: 850 - total flux cleaned = 8.55459e-06 Jy ! Component: 900 - total flux cleaned = 3.69383e-05 Jy ! Component: 950 - total flux cleaned = -4.61317e-05 Jy ! Component: 1000 - total flux cleaned = -1.84287e-05 Jy ! Component: 1050 - total flux cleaned = -7.19129e-05 Jy ! Component: 1100 - total flux cleaned = -7.16968e-05 Jy ! Component: 1150 - total flux cleaned = -7.20198e-05 Jy ! Component: 1200 - total flux cleaned = -2.05023e-05 Jy ! Component: 1250 - total flux cleaned = -2.10596e-05 Jy ! Component: 1300 - total flux cleaned = 2.86994e-05 Jy ! Component: 1350 - total flux cleaned = 5.29511e-05 Jy ! Component: 1400 - total flux cleaned = 3.79348e-06 Jy ! Component: 1450 - total flux cleaned = 2.82562e-05 Jy ! Component: 1500 - total flux cleaned = 7.57116e-05 Jy ! Component: 1550 - total flux cleaned = 9.85462e-05 Jy ! Component: 1600 - total flux cleaned = 9.81081e-05 Jy ! Component: 1650 - total flux cleaned = 2.99074e-05 Jy ! Component: 1700 - total flux cleaned = 9.77141e-05 Jy ! Component: 1750 - total flux cleaned = 0.000141977 Jy ! Component: 1800 - total flux cleaned = 0.00014176 Jy ! Component: 1850 - total flux cleaned = 0.000142113 Jy ! Component: 1900 - total flux cleaned = 0.000206266 Jy ! Component: 1950 - total flux cleaned = 0.000164486 Jy ! Component: 2000 - total flux cleaned = 0.000206546 Jy ! Component: 2050 - total flux cleaned = 0.000227464 Jy ! Component: 2100 - total flux cleaned = 0.000145716 Jy ! Component: 2150 - total flux cleaned = 0.000226714 Jy ! Component: 2200 - total flux cleaned = 0.000226738 Jy ! Component: 2250 - total flux cleaned = 0.000206842 Jy ! Component: 2300 - total flux cleaned = 0.00020712 Jy ! Component: 2350 - total flux cleaned = 0.000226792 Jy ! Component: 2400 - total flux cleaned = 0.000245795 Jy ! Component: 2450 - total flux cleaned = 0.000245615 Jy ! Component: 2500 - total flux cleaned = 0.00030201 Jy ! Component: 2550 - total flux cleaned = 0.000301874 Jy ! Component: 2600 - total flux cleaned = 0.000265313 Jy ! Component: 2650 - total flux cleaned = 0.000283283 Jy ! Component: 2700 - total flux cleaned = 0.000301003 Jy ! Component: 2750 - total flux cleaned = 0.000301036 Jy ! Component: 2800 - total flux cleaned = 0.000301041 Jy ! Component: 2850 - total flux cleaned = 0.000231509 Jy ! Component: 2900 - total flux cleaned = 0.000197287 Jy ! Component: 2950 - total flux cleaned = 0.000231602 Jy ! Component: 3000 - total flux cleaned = 0.000163417 Jy ! Component: 3050 - total flux cleaned = 6.24844e-05 Jy ! Component: 3100 - total flux cleaned = 6.26876e-05 Jy ! Component: 3150 - total flux cleaned = 1.30466e-05 Jy ! Component: 3200 - total flux cleaned = -1.9916e-05 Jy ! Component: 3250 - total flux cleaned = -8.50692e-05 Jy ! Component: 3300 - total flux cleaned = -0.000181807 Jy ! Component: 3350 - total flux cleaned = -0.000181794 Jy ! Component: 3400 - total flux cleaned = -0.000245179 Jy ! Component: 3450 - total flux cleaned = -0.000292396 Jy ! Component: 3500 - total flux cleaned = -0.000261137 Jy ! Component: 3550 - total flux cleaned = -0.00024559 Jy ! Component: 3600 - total flux cleaned = -0.000261173 Jy ! Component: 3650 - total flux cleaned = -0.000306705 Jy ! Component: 3700 - total flux cleaned = -0.000321763 Jy ! Component: 3750 - total flux cleaned = -0.00036683 Jy ! Component: 3800 - total flux cleaned = -0.000426444 Jy ! Component: 3850 - total flux cleaned = -0.000456021 Jy ! Component: 3900 - total flux cleaned = -0.00048545 Jy ! Component: 3950 - total flux cleaned = -0.000558531 Jy ! Component: 4000 - total flux cleaned = -0.000573008 Jy ! Component: 4050 - total flux cleaned = -0.000587414 Jy ! Total flux subtracted in 4096 components = -0.000615923 Jy ! Clean residual min=-0.000959 max=0.000856 Jy/beam ! Clean residual mean=0.000001 rms=0.000251 Jy/beam ! Combined flux in latest and established models = 0.095912 Jy keep ! Adding 281 model components to the UV plane model. ! The established model now contains 340 components and 0.095912 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /vlbi/uf001/uva//J0341+3352_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.203 x 6.162 at -2.917 degrees (North through East) ! Clean map min=-0.0017019 max=0.074687 Jy/beam ! Writing clean map to FITS file: /vlbi/uf001/uva//J0341+3352_S_map.fits wmodel /vlbi/uf001/uva//J0341+3352_S_map.mod ! Writing 340 model components to file: /vlbi/uf001/uva//J0341+3352_S_map.mod wobs /vlbi/uf001/uva//J0341+3352_S_uvs.fits ! Writing UV FITS file: /vlbi/uf001/uva//J0341+3352_S_uvs.fits wwins /vlbi/uf001/uva//J0341+3352_S_map.win ! wwins: Wrote 2 windows to /vlbi/uf001/uva//J0341+3352_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.0746872 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000247602 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.994557 0.994557 1.98911 3.97823 7.95646 15.9129 31.8258 63.6517 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0007 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.074 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 301.641 quit ! Quitting program ! Log file /vlbi/uf001/uva//J0341+3352_S_dfm.log closed on Sun Jan 13 21:53:17 2019