! Started logfile: /image_prod/ug002_uva//J0122+2502_S_dfm.log on Wed Mar 4 19:02:07 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//J0122+2502_S_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J0122+2502_S_uva.fits ! AN table 1: 75 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.975407 visibilities/baseline/integration-bin. ! Found source: J0122+2502 ! ! 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 9876 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 J0122+2502_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 9 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.051271Jy sigma=8.258932 ! Fit after self-cal, rms=0.688983Jy sigma=5.422488 ! 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.802 mas, bmaj=6.468 mas, bpa=-5.294 degrees ! Estimated noise=1.47262 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.236286 Jy ! Component: 100 - total flux cleaned = 0.304906 Jy ! Total flux subtracted in 100 components = 0.304906 Jy ! Clean residual min=-0.014236 max=0.042065 Jy/beam ! Clean residual mean=0.000008 rms=0.004268 Jy/beam ! Combined flux in latest and established models = 0.304906 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 0.304906 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.108171Jy sigma=0.802087 ! Fit after self-cal, rms=0.108026Jy sigma=0.800422 ! Inverting map ! Added new window around map position (-6, -8.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.040641 Jy ! Component: 100 - total flux cleaned = 0.0651587 Jy ! Total flux subtracted in 100 components = 0.0651587 Jy ! Clean residual min=-0.009894 max=0.017541 Jy/beam ! Clean residual mean=0.000022 rms=0.002997 Jy/beam ! Combined flux in latest and established models = 0.370064 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 13 components and 0.370064 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.088486Jy sigma=0.641048 ! Fit after self-cal, rms=0.087256Jy sigma=0.634068 ! 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.302 mas, bmaj=7.248 mas, bpa=-4.64 degrees ! Estimated noise=0.929349 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.0220202 Jy ! Component: 100 - total flux cleaned = 0.0355551 Jy ! Component: 150 - total flux cleaned = 0.0441019 Jy ! Component: 200 - total flux cleaned = 0.0469705 Jy ! Total flux subtracted in 200 components = 0.0469705 Jy ! Clean residual min=-0.010251 max=0.015934 Jy/beam ! Clean residual mean=0.000024 rms=0.002984 Jy/beam ! Combined flux in latest and established models = 0.417035 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 39 components and 0.417035 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.080559Jy sigma=0.584188 ! Fit after self-cal, rms=0.080073Jy sigma=0.581001 ! 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.302 x 7.248 at -4.64 degrees (North through East) ! Clean map min=-0.010283 max=0.30766 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.10 HN 1.00 KP 1.02 ! LA 1.00 MK 0.98 NL 0.99 OV 0.87 ! PT 0.99 SC 1.00 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.08 HN 0.99 KP 1.06 ! LA 1.00 MK 1.00 NL 0.96 OV 0.88 ! PT 1.01 SC 0.93 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.09 HN 1.03 KP 1.07 ! LA 1.01 MK 1.02 NL 1.04 OV 0.86 ! PT 1.04 SC 0.99 ! ! ! Fit before self-cal, rms=0.080073Jy sigma=0.581001 ! Fit after self-cal, rms=0.072216Jy sigma=0.533269 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.3 mas, bmaj=7.193 mas, bpa=-4.154 degrees ! Estimated noise=0.916728 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 (-21, -32). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.014145 Jy ! Component: 100 - total flux cleaned = 0.0235241 Jy ! Component: 150 - total flux cleaned = 0.0306312 Jy ! Component: 200 - total flux cleaned = 0.0362779 Jy ! Total flux subtracted in 200 components = 0.0362779 Jy ! Clean residual min=-0.006760 max=0.008993 Jy/beam ! Clean residual mean=0.000018 rms=0.001761 Jy/beam ! Combined flux in latest and established models = 0.453313 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 91 components and 0.453313 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.068835Jy sigma=0.502212 ! Fit after self-cal, rms=0.068645Jy sigma=0.500742 ! Inverting map ! Added new window around map position (44, 114). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00810622 Jy ! Component: 100 - total flux cleaned = 0.0132451 Jy ! Component: 150 - total flux cleaned = 0.0172138 Jy ! Component: 200 - total flux cleaned = 0.020062 Jy ! Total flux subtracted in 200 components = 0.020062 Jy ! Clean residual min=-0.006042 max=0.007603 Jy/beam ! Clean residual mean=0.000009 rms=0.001578 Jy/beam ! Combined flux in latest and established models = 0.473375 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 127 components and 0.473375 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.067265Jy sigma=0.490646 ! Fit after self-cal, rms=0.067151Jy sigma=0.489742 ! 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.067151Jy sigma=0.489742 ! Fit after self-cal, rms=0.061085Jy sigma=0.431394 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.301 mas, bmaj=7.224 mas, bpa=-4.105 degrees ! Estimated noise=0.922651 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00378215 Jy ! Component: 100 - total flux cleaned = 0.00581808 Jy ! Component: 150 - total flux cleaned = 0.00739316 Jy ! Component: 200 - total flux cleaned = 0.00863595 Jy ! Total flux subtracted in 200 components = 0.00863595 Jy ! Clean residual min=-0.002483 max=0.003647 Jy/beam ! Clean residual mean=-0.000001 rms=0.000738 Jy/beam ! Combined flux in latest and established models = 0.482011 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 149 components and 0.482011 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060681Jy sigma=0.428054 ! Fit after self-cal, rms=0.060695Jy sigma=0.427897 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 (-12.5, -22.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00283026 Jy ! Component: 100 - total flux cleaned = 0.00355311 Jy ! Component: 150 - total flux cleaned = 0.00404899 Jy ! Component: 200 - total flux cleaned = 0.00439171 Jy ! Total flux subtracted in 200 components = 0.00439171 Jy ! Clean residual min=-0.002493 max=0.003606 Jy/beam ! Clean residual mean=-0.000003 rms=0.000686 Jy/beam ! Combined flux in latest and established models = 0.486402 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 170 components and 0.486402 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060488Jy sigma=0.426254 ! Fit after self-cal, rms=0.060478Jy sigma=0.426168 ! Inverting map ! Added new window around map position (-29.5, -39). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00326726 Jy ! Component: 100 - total flux cleaned = 0.00466463 Jy ! Component: 150 - total flux cleaned = 0.00557069 Jy ! Component: 200 - total flux cleaned = 0.00614828 Jy ! Total flux subtracted in 200 components = 0.00614828 Jy ! Clean residual min=-0.002178 max=0.002417 Jy/beam ! Clean residual mean=-0.000001 rms=0.000611 Jy/beam ! Combined flux in latest and established models = 0.492551 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 197 components and 0.492551 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060277Jy sigma=0.424154 ! Fit after self-cal, rms=0.060256Jy sigma=0.424034 ! 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.060256Jy sigma=0.424034 ! Fit after self-cal, rms=0.060070Jy sigma=0.422763 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.292 mas, bmaj=7.214 mas, bpa=-4.1 degrees ! Estimated noise=0.92127 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00116138 Jy ! Component: 100 - total flux cleaned = 0.0020548 Jy ! Component: 150 - total flux cleaned = 0.00275375 Jy ! Component: 200 - total flux cleaned = 0.00321858 Jy ! Total flux subtracted in 200 components = 0.00321858 Jy ! Clean residual min=-0.001965 max=0.001929 Jy/beam ! Clean residual mean=-0.000001 rms=0.000531 Jy/beam ! Combined flux in latest and established models = 0.495769 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 219 components and 0.495769 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059935Jy sigma=0.421595 ! Fit after self-cal, rms=0.059930Jy sigma=0.421553 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.000522095 Jy ! Component: 100 - total flux cleaned = 0.000708642 Jy ! Component: 150 - total flux cleaned = 0.000974294 Jy ! Component: 200 - total flux cleaned = 0.00118466 Jy ! Total flux subtracted in 200 components = 0.00118466 Jy ! Clean residual min=-0.001845 max=0.001907 Jy/beam ! Clean residual mean=-0.000000 rms=0.000509 Jy/beam ! Combined flux in latest and established models = 0.496954 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 235 components and 0.496954 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059867Jy sigma=0.421054 ! Fit after self-cal, rms=0.059864Jy sigma=0.421039 ! 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.059864Jy sigma=0.421039 ! Fit after self-cal, rms=0.059824Jy sigma=0.420846 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.199 mas, bmaj=10.53 mas, bpa=-5.767 degrees ! Estimated noise=1.15812 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00148456 Jy ! Component: 100 - total flux cleaned = 0.00187343 Jy ! Component: 150 - total flux cleaned = 0.00226533 Jy ! Component: 200 - total flux cleaned = 0.00260506 Jy ! Total flux subtracted in 200 components = 0.00260506 Jy ! Clean residual min=-0.002211 max=0.001983 Jy/beam ! Clean residual mean=0.000001 rms=0.000607 Jy/beam ! Combined flux in latest and established models = 0.499559 Jy selfcal ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 265 components and 0.499559 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059761Jy sigma=0.420345 ! Fit after self-cal, rms=0.059752Jy sigma=0.420289 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.000706392 Jy ! Component: 100 - total flux cleaned = 0.00101692 Jy ! Component: 150 - total flux cleaned = 0.00125499 Jy ! Component: 200 - total flux cleaned = 0.0015384 Jy ! Total flux subtracted in 200 components = 0.0015384 Jy ! Clean residual min=-0.002036 max=0.001856 Jy/beam ! Clean residual mean=0.000002 rms=0.000573 Jy/beam ! Combined flux in latest and established models = 0.501097 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 286 components and 0.501097 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059727Jy sigma=0.420109 ! Fit after self-cal, rms=0.059725Jy sigma=0.420085 ! 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.287 mas, bmaj=7.208 mas, bpa=-4.103 degrees ! Estimated noise=0.921199 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.287 x 7.208 at -4.103 degrees (North through East) ! Clean map min=-0.0016855 max=0.30401 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.795 mas, bmaj=6.427 mas, bpa=-4.788 degrees ! Estimated noise=1.46119 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.234089 Jy ! Component: 100 - total flux cleaned = 0.306066 Jy ! Component: 150 - total flux cleaned = 0.346774 Jy ! Component: 200 - total flux cleaned = 0.373029 Jy ! Total flux subtracted in 200 components = 0.373029 Jy ! Clean residual min=-0.006832 max=0.014260 Jy/beam ! Clean residual mean=-0.000009 rms=0.002044 Jy/beam ! Combined flux in latest and established models = 0.373029 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 20 model components to the UV plane model. ! The established model now contains 20 components and 0.373029 Jy ! Inverting map and beam ! Estimated beam: bmin=3.287 mas, bmaj=7.208 mas, bpa=-4.103 degrees ! Estimated noise=0.921199 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) ! Added new window around map position (-3, -4.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0225829 Jy ! Component: 100 - total flux cleaned = 0.0383418 Jy ! Component: 150 - total flux cleaned = 0.0507926 Jy ! Component: 200 - total flux cleaned = 0.0609622 Jy ! Total flux subtracted in 200 components = 0.0609622 Jy ! Clean residual min=-0.004407 max=0.006148 Jy/beam ! Clean residual mean=0.000002 rms=0.001411 Jy/beam ! Combined flux in latest and established models = 0.433991 Jy ! Adding 43 model components to the UV plane model. ! The established model now contains 57 components and 0.433991 Jy ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00845172 Jy ! Component: 100 - total flux cleaned = 0.01553 Jy ! Component: 150 - total flux cleaned = 0.0215265 Jy ! Component: 200 - total flux cleaned = 0.0267162 Jy ! Total flux subtracted in 200 components = 0.0267162 Jy ! Clean residual min=-0.003123 max=0.003476 Jy/beam ! Clean residual mean=0.000002 rms=0.000907 Jy/beam ! Combined flux in latest and established models = 0.460707 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 112 components and 0.460707 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.199 mas, bmaj=10.53 mas, bpa=-5.767 degrees ! Estimated noise=1.15812 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00803887 Jy ! Component: 100 - total flux cleaned = 0.0138858 Jy ! Component: 150 - total flux cleaned = 0.0184358 Jy ! Component: 200 - total flux cleaned = 0.0220654 Jy ! Total flux subtracted in 200 components = 0.0220654 Jy ! Clean residual min=-0.002924 max=0.002949 Jy/beam ! Clean residual mean=0.000005 rms=0.000841 Jy/beam ! Combined flux in latest and established models = 0.482773 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.00296889 Jy ! Component: 100 - total flux cleaned = 0.0054398 Jy ! Component: 150 - total flux cleaned = 0.00752588 Jy ! Component: 200 - total flux cleaned = 0.00931765 Jy ! Total flux subtracted in 200 components = 0.00931765 Jy ! Clean residual min=-0.002220 max=0.002252 Jy/beam ! Clean residual mean=0.000004 rms=0.000636 Jy/beam ! Combined flux in latest and established models = 0.49209 Jy ! Adding 168 model components to the UV plane model. ! The established model now contains 262 components and 0.49209 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.287 mas, bmaj=7.208 mas, bpa=-4.103 degrees ! Estimated noise=0.921199 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00144996 Jy ! Component: 100 - total flux cleaned = 0.00275772 Jy ! Component: 150 - total flux cleaned = 0.00371697 Jy ! Component: 200 - total flux cleaned = 0.0044028 Jy ! Total flux subtracted in 200 components = 0.0044028 Jy ! Clean residual min=-0.001893 max=0.002125 Jy/beam ! Clean residual mean=0.000002 rms=0.000499 Jy/beam ! Combined flux in latest and established models = 0.496493 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 292 components and 0.496493 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059924Jy sigma=0.421416 ! Fit after self-cal, rms=0.059893Jy sigma=0.421172 wmodel J0122+2502_S_map.mod ! Writing 292 model components to file: J0122+2502_S_map.mod wobs J0122+2502_S_uvs.fits ! Writing UV FITS file: J0122+2502_S_uvs.fits wwins J0122+2502_S_map.win ! wwins: Wrote 7 windows to J0122+2502_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.000234417 Jy ! Component: 100 - total flux cleaned = 0.000300965 Jy ! Component: 150 - total flux cleaned = 0.000301509 Jy ! Component: 200 - total flux cleaned = 0.000240686 Jy ! Component: 250 - total flux cleaned = 0.000300095 Jy ! Component: 300 - total flux cleaned = 0.000415229 Jy ! Component: 350 - total flux cleaned = 0.000471649 Jy ! Component: 400 - total flux cleaned = 0.000610205 Jy ! Component: 450 - total flux cleaned = 0.000691382 Jy ! Component: 500 - total flux cleaned = 0.000798166 Jy ! Component: 550 - total flux cleaned = 0.000929093 Jy ! Component: 600 - total flux cleaned = 0.00095482 Jy ! Component: 650 - total flux cleaned = 0.00110719 Jy ! Component: 700 - total flux cleaned = 0.00115723 Jy ! Component: 750 - total flux cleaned = 0.00125567 Jy ! Component: 800 - total flux cleaned = 0.00137738 Jy ! Component: 850 - total flux cleaned = 0.00149697 Jy ! Component: 900 - total flux cleaned = 0.00159124 Jy ! Component: 950 - total flux cleaned = 0.00166127 Jy ! Component: 1000 - total flux cleaned = 0.0016845 Jy ! Component: 1050 - total flux cleaned = 0.00184412 Jy ! Component: 1100 - total flux cleaned = 0.00188916 Jy ! Component: 1150 - total flux cleaned = 0.00193367 Jy ! Component: 1200 - total flux cleaned = 0.00199945 Jy ! Component: 1250 - total flux cleaned = 0.00195565 Jy ! Component: 1300 - total flux cleaned = 0.00199834 Jy ! Component: 1350 - total flux cleaned = 0.00197701 Jy ! Component: 1400 - total flux cleaned = 0.00197693 Jy ! Component: 1450 - total flux cleaned = 0.00201856 Jy ! Component: 1500 - total flux cleaned = 0.00210052 Jy ! Component: 1550 - total flux cleaned = 0.00203972 Jy ! Component: 1600 - total flux cleaned = 0.00205989 Jy ! Component: 1650 - total flux cleaned = 0.00207984 Jy ! Component: 1700 - total flux cleaned = 0.00211934 Jy ! Component: 1750 - total flux cleaned = 0.00217801 Jy ! Component: 1800 - total flux cleaned = 0.00227459 Jy ! Component: 1850 - total flux cleaned = 0.00223634 Jy ! Component: 1900 - total flux cleaned = 0.00227422 Jy ! Component: 1950 - total flux cleaned = 0.00229297 Jy ! Component: 2000 - total flux cleaned = 0.00223727 Jy ! Component: 2050 - total flux cleaned = 0.00231123 Jy ! Component: 2100 - total flux cleaned = 0.00227444 Jy ! Component: 2150 - total flux cleaned = 0.00234706 Jy ! Component: 2200 - total flux cleaned = 0.00222112 Jy ! Component: 2250 - total flux cleaned = 0.00231043 Jy ! Component: 2300 - total flux cleaned = 0.00227512 Jy ! Component: 2350 - total flux cleaned = 0.0022401 Jy ! Component: 2400 - total flux cleaned = 0.00225756 Jy ! Component: 2450 - total flux cleaned = 0.00230894 Jy ! Component: 2500 - total flux cleaned = 0.00227495 Jy ! Component: 2550 - total flux cleaned = 0.00232558 Jy ! Component: 2600 - total flux cleaned = 0.00230882 Jy ! Component: 2650 - total flux cleaned = 0.00230889 Jy ! Component: 2700 - total flux cleaned = 0.00234211 Jy ! Component: 2750 - total flux cleaned = 0.00224392 Jy ! Component: 2800 - total flux cleaned = 0.00229259 Jy ! Component: 2850 - total flux cleaned = 0.00227668 Jy ! Component: 2900 - total flux cleaned = 0.00230887 Jy ! Component: 2950 - total flux cleaned = 0.0023725 Jy ! Component: 3000 - total flux cleaned = 0.00237246 Jy ! Component: 3050 - total flux cleaned = 0.00240362 Jy ! Component: 3100 - total flux cleaned = 0.00240356 Jy ! Component: 3150 - total flux cleaned = 0.00241897 Jy ! Component: 3200 - total flux cleaned = 0.0024189 Jy ! Component: 3250 - total flux cleaned = 0.00243408 Jy ! Component: 3300 - total flux cleaned = 0.00237401 Jy ! Component: 3350 - total flux cleaned = 0.00243375 Jy ! Component: 3400 - total flux cleaned = 0.00240411 Jy ! Component: 3450 - total flux cleaned = 0.0024336 Jy ! Component: 3500 - total flux cleaned = 0.00246279 Jy ! Component: 3550 - total flux cleaned = 0.00247729 Jy ! Component: 3600 - total flux cleaned = 0.00249152 Jy ! Component: 3650 - total flux cleaned = 0.00247728 Jy ! Component: 3700 - total flux cleaned = 0.00254809 Jy ! Component: 3750 - total flux cleaned = 0.00253395 Jy ! Component: 3800 - total flux cleaned = 0.00249198 Jy ! Component: 3850 - total flux cleaned = 0.00247806 Jy ! Component: 3900 - total flux cleaned = 0.00251927 Jy ! Component: 3950 - total flux cleaned = 0.00246456 Jy ! Component: 4000 - total flux cleaned = 0.0024782 Jy ! Component: 4050 - total flux cleaned = 0.00245126 Jy ! Total flux subtracted in 4096 components = 0.00251812 Jy ! Clean residual min=-0.000865 max=0.000785 Jy/beam ! Clean residual mean=0.000000 rms=0.000254 Jy/beam ! Combined flux in latest and established models = 0.499011 Jy keep ! Adding 297 model components to the UV plane model. ! The established model now contains 589 components and 0.499011 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J0122+2502_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.287 x 7.208 at -4.103 degrees (North through East) ! Clean map min=-0.0013881 max=0.3029 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J0122+2502_S_map.fits wmodel /image_prod/ug002_uva//J0122+2502_S_map.mod ! Writing 589 model components to file: /image_prod/ug002_uva//J0122+2502_S_map.mod wobs /image_prod/ug002_uva//J0122+2502_S_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J0122+2502_S_uvs.fits wwins /image_prod/ug002_uva//J0122+2502_S_map.win ! wwins: Wrote 8 windows to /image_prod/ug002_uva//J0122+2502_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.302902 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000251668 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.249258 0.249258 0.498515 0.99703 1.99406 3.98812 7.97624 15.9525 31.905 63.8099 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0007 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.302 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1203.57 quit ! Quitting program ! Log file /image_prod/ug002_uva//J0122+2502_S_dfm.log closed on Wed Mar 4 19:02:12 2020