! Started logfile: .//J1351+0830_C_dfm.log on Mon Jan 20 18:09:52 2020 float field_size field_size = 1024 float field_cell field_cell = 0.3000 float taper_size taper_size = 40.0 float freq freq = 4.12800 obs .//J1351+0830_C_uva.fits ! Reading UV FITS file: .//J1351+0830_C_uva.fits ! AN table 1: 15 integrations on 45 of 45 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J1351+0830 ! ! There are 8 IFs, and a total of 8 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 4.128e+09 3.2e+07 1 3.2e+07 ! 02 2 4.16e+09 3.2e+07 1 3.2e+07 ! 03 3 4.192e+09 3.2e+07 1 3.2e+07 ! 04 4 4.224e+09 3.2e+07 1 3.2e+07 ! 05 5 4.416e+09 3.2e+07 1 3.2e+07 ! 06 6 4.512e+09 3.2e+07 1 3.2e+07 ! 07 7 4.544e+09 3.2e+07 1 3.2e+07 ! 08 8 4.576e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 94 lines of history. ! ! Reading 5400 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 8 channels. ! Selecting polarization: RR, channels: 1..8 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: 5..5 ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 print "field_size =", field_size ! field_size = 1024 print "field_cell =", field_cell ! field_cell = 0.3 print "taper_size =", taper_size ! taper_size = 40 ![@/opt64/bin/pima_mupet_01.dfm J1351+0830_C] 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.300x0.300 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.044445Jy sigma=22.455231 ! Fit after self-cal, rms=0.727127Jy sigma=15.302222 ! 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=1.476 mas, bmaj=3.663 mas, bpa=-11.79 degrees ! Estimated noise=0.703213 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.199754 Jy ! Component: 100 - total flux cleaned = 0.267874 Jy ! Total flux subtracted in 100 components = 0.267874 Jy ! Clean residual min=-0.020115 max=0.042954 Jy/beam ! Clean residual mean=-0.000003 rms=0.006946 Jy/beam ! Combined flux in latest and established models = 0.267874 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.267874 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.121897Jy sigma=2.559900 ! Fit after self-cal, rms=0.121630Jy sigma=2.556765 ! Inverting map ! Added new window around map position (-2.4, -2.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0432863 Jy ! Component: 100 - total flux cleaned = 0.0715301 Jy ! Total flux subtracted in 100 components = 0.0715301 Jy ! Clean residual min=-0.019981 max=0.031144 Jy/beam ! Clean residual mean=0.000000 rms=0.005702 Jy/beam ! Combined flux in latest and established models = 0.339404 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 24 components and 0.339404 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.098194Jy sigma=1.967319 ! Fit after self-cal, rms=0.097445Jy sigma=1.945590 ! 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=1.822 mas, bmaj=4.107 mas, bpa=-10.59 degrees ! Estimated noise=0.457451 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.0324762 Jy ! Component: 100 - total flux cleaned = 0.051775 Jy ! Component: 150 - total flux cleaned = 0.0613161 Jy ! Component: 200 - total flux cleaned = 0.0656277 Jy ! Total flux subtracted in 200 components = 0.0656277 Jy ! Clean residual min=-0.018802 max=0.020571 Jy/beam ! Clean residual mean=0.000002 rms=0.004350 Jy/beam ! Combined flux in latest and established models = 0.405032 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 42 components and 0.405032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.085869Jy sigma=1.637955 ! Fit after self-cal, rms=0.085699Jy sigma=1.631384 ! Inverting map print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.822 x 4.107 at -10.59 degrees (North through East) ! Clean map min=-0.017843 max=0.27066 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.01 HN 0.84 KP 0.96 ! LA 0.99 MK 0.90 NL 0.94 OV 0.95 ! PT 1.01 SC 0.80 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.00 HN 0.84 KP 0.97 ! LA 0.98 MK 0.91 NL 0.92 OV 0.96 ! PT 1.01 SC 0.79 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.01 HN 0.84 KP 0.96 ! LA 0.98 MK 0.90 NL 0.92 OV 0.96 ! PT 1.02 SC 0.81 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.01 HN 0.84 KP 0.96 ! LA 0.97 MK 0.91 NL 0.91 OV 0.98 ! PT 0.99 SC 0.81 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.01 HN 1.25 KP 0.97 ! LA 0.97 MK 1.02 NL 0.93 OV 0.98 ! PT 1.00 SC 1.38 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.00 HN 1.97 KP 0.97 ! LA 0.97 MK 1.14 NL 1.01 OV 0.96 ! PT 1.01 SC 2.74 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 0.99 HN 2.58 KP 0.98 ! LA 0.95 MK 1.17 NL 1.06 OV 0.97 ! PT 1.00 SC 3.42 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 1.00 HN 3.67 KP 0.97 ! LA 0.96 MK 1.25 NL 1.10 OV 0.96 ! PT 0.99 SC 6.41 ! ! ! Fit before self-cal, rms=0.085699Jy sigma=1.631384 ! Fit after self-cal, rms=0.370226Jy sigma=1.122682 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.88 mas, bmaj=4.272 mas, bpa=-12.03 degrees ! Estimated noise=0.467768 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.00102556 Jy ! Component: 100 - total flux cleaned = -0.00262542 Jy ! Component: 150 - total flux cleaned = -0.00291312 Jy ! Component: 200 - total flux cleaned = -0.00304141 Jy ! Total flux subtracted in 200 components = -0.00304141 Jy ! Clean residual min=-0.011852 max=0.012765 Jy/beam ! Clean residual mean=-0.000009 rms=0.002539 Jy/beam ! Combined flux in latest and established models = 0.40199 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 55 components and 0.40199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.369558Jy sigma=1.106995 ! Fit after self-cal, rms=0.371781Jy sigma=1.104519 ! Inverting map ! Added new window around map position (-24, 13.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00932818 Jy ! Component: 100 - total flux cleaned = 0.0107674 Jy ! Component: 150 - total flux cleaned = 0.0104755 Jy ! Component: 200 - total flux cleaned = 0.0102203 Jy ! Total flux subtracted in 200 components = 0.0102203 Jy ! Clean residual min=-0.010510 max=0.010395 Jy/beam ! Clean residual mean=-0.000008 rms=0.002292 Jy/beam ! Combined flux in latest and established models = 0.412211 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 73 components and 0.412211 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.368675Jy sigma=1.057956 ! Fit after self-cal, rms=0.367938Jy sigma=1.055738 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.367938Jy sigma=1.055738 ! Fit after self-cal, rms=0.306605Jy sigma=1.053083 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.872 mas, bmaj=4.262 mas, bpa=-12.11 degrees ! Estimated noise=0.467047 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00324687 Jy ! Component: 100 - total flux cleaned = 0.00303091 Jy ! Component: 150 - total flux cleaned = 0.0030311 Jy ! Component: 200 - total flux cleaned = 0.00302898 Jy ! Total flux subtracted in 200 components = 0.00302898 Jy ! Clean residual min=-0.010193 max=0.010463 Jy/beam ! Clean residual mean=-0.000007 rms=0.002245 Jy/beam ! Combined flux in latest and established models = 0.41524 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 85 components and 0.41524 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.304853Jy sigma=1.043365 ! Fit after self-cal, rms=0.304374Jy sigma=1.042602 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 (-58.5, 8.1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0059523 Jy ! Component: 100 - total flux cleaned = 0.00813204 Jy ! Component: 150 - total flux cleaned = 0.008765 Jy ! Component: 200 - total flux cleaned = 0.00916401 Jy ! Total flux subtracted in 200 components = 0.00916401 Jy ! Clean residual min=-0.009241 max=0.009716 Jy/beam ! Clean residual mean=-0.000005 rms=0.002158 Jy/beam ! Combined flux in latest and established models = 0.424404 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 103 components and 0.424404 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.302141Jy sigma=1.021219 ! Fit after self-cal, rms=0.301746Jy sigma=1.018766 ! Inverting map ! Added new window around map position (4.2, 37.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00072197 Jy ! Component: 100 - total flux cleaned = 0.000192096 Jy ! Component: 150 - total flux cleaned = 0.000877589 Jy ! Component: 200 - total flux cleaned = 0.00151034 Jy ! Total flux subtracted in 200 components = 0.00151034 Jy ! Clean residual min=-0.007869 max=0.008581 Jy/beam ! Clean residual mean=-0.000005 rms=0.002032 Jy/beam ! Combined flux in latest and established models = 0.425914 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 113 components and 0.425914 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.299276Jy sigma=0.988822 ! Fit after self-cal, rms=0.298982Jy sigma=0.987712 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.298982Jy sigma=0.987712 ! Fit after self-cal, rms=0.215209Jy sigma=0.978028 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.864 mas, bmaj=4.24 mas, bpa=-11.91 degrees ! Estimated noise=0.464421 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00422468 Jy ! Component: 100 - total flux cleaned = 0.00485759 Jy ! Component: 150 - total flux cleaned = 0.00523241 Jy ! Component: 200 - total flux cleaned = 0.00505325 Jy ! Total flux subtracted in 200 components = 0.00505325 Jy ! Clean residual min=-0.007256 max=0.008355 Jy/beam ! Clean residual mean=-0.000003 rms=0.001965 Jy/beam ! Combined flux in latest and established models = 0.430967 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 129 components and 0.430967 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.213372Jy sigma=0.970203 ! Fit after self-cal, rms=0.212962Jy sigma=0.969628 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 (-9.3, 17.1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00289824 Jy ! Component: 100 - total flux cleaned = 0.00348675 Jy ! Component: 150 - total flux cleaned = 0.00477791 Jy ! Component: 200 - total flux cleaned = 0.00494769 Jy ! Total flux subtracted in 200 components = 0.00494769 Jy ! Clean residual min=-0.007168 max=0.007711 Jy/beam ! Clean residual mean=-0.000003 rms=0.001907 Jy/beam ! Combined flux in latest and established models = 0.435915 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 147 components and 0.435915 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.210799Jy sigma=0.955108 ! Fit after self-cal, rms=0.210610Jy sigma=0.953749 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 40 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.210610Jy sigma=0.953749 ! Fit after self-cal, rms=0.174081Jy sigma=0.951059 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.029 mas, bmaj=5.806 mas, bpa=4.898 degrees ! Estimated noise=0.534686 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00240206 Jy ! Component: 100 - total flux cleaned = -0.00339523 Jy ! Component: 150 - total flux cleaned = -0.00395747 Jy ! Component: 200 - total flux cleaned = -0.00406116 Jy ! Total flux subtracted in 200 components = -0.00406116 Jy ! Clean residual min=-0.006860 max=0.007272 Jy/beam ! Clean residual mean=-0.000001 rms=0.001969 Jy/beam ! Combined flux in latest and established models = 0.431854 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 168 components and 0.431854 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.173397Jy sigma=0.945177 ! Fit after self-cal, rms=0.173280Jy sigma=0.944537 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.00012027 Jy ! Component: 100 - total flux cleaned = -0.00030863 Jy ! Component: 150 - total flux cleaned = -0.000486059 Jy ! Component: 200 - total flux cleaned = -0.00074222 Jy ! Total flux subtracted in 200 components = -0.00074222 Jy ! Clean residual min=-0.006835 max=0.007108 Jy/beam ! Clean residual mean=-0.000000 rms=0.001945 Jy/beam ! Combined flux in latest and established models = 0.431111 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 178 components and 0.431111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.172660Jy sigma=0.943638 ! Fit after self-cal, rms=0.172632Jy sigma=0.943256 ! 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=1.864 mas, bmaj=4.228 mas, bpa=-11.73 degrees ! Estimated noise=0.462866 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.864 x 4.228 at -11.73 degrees (North through East) ! Clean map min=-0.013261 max=0.26921 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=1.492 mas, bmaj=3.726 mas, bpa=-12.68 degrees ! Estimated noise=0.713319 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.201232 Jy ! Component: 100 - total flux cleaned = 0.27016 Jy ! Component: 150 - total flux cleaned = 0.311298 Jy ! Component: 200 - total flux cleaned = 0.339516 Jy ! Total flux subtracted in 200 components = 0.339516 Jy ! Clean residual min=-0.011613 max=0.016956 Jy/beam ! Clean residual mean=-0.000004 rms=0.003280 Jy/beam ! Combined flux in latest and established models = 0.339516 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 24 model components to the UV plane model. ! The established model now contains 24 components and 0.339516 Jy ! Inverting map and beam ! Estimated beam: bmin=1.864 mas, bmaj=4.228 mas, bpa=-11.73 degrees ! Estimated noise=0.462866 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.0297669 Jy ! Component: 100 - total flux cleaned = 0.04964 Jy ! Component: 150 - total flux cleaned = 0.0650565 Jy ! Component: 200 - total flux cleaned = 0.0762226 Jy ! Total flux subtracted in 200 components = 0.0762226 Jy ! Clean residual min=-0.010099 max=0.013583 Jy/beam ! Clean residual mean=-0.000006 rms=0.002548 Jy/beam ! Combined flux in latest and established models = 0.415738 Jy ! Adding 32 model components to the UV plane model. ! The established model now contains 43 components and 0.415738 Jy ! Inverting map ! Added new window around map position (-46.5, -66). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0115866 Jy ! Component: 100 - total flux cleaned = 0.0150857 Jy ! Component: 150 - total flux cleaned = 0.0175563 Jy ! Component: 200 - total flux cleaned = 0.0198029 Jy ! Total flux subtracted in 200 components = 0.0198029 Jy ! Clean residual min=-0.008179 max=0.010608 Jy/beam ! Clean residual mean=-0.000001 rms=0.002162 Jy/beam ! Combined flux in latest and established models = 0.435541 Jy ! Adding 31 model components to the UV plane model. ! The established model now contains 64 components and 0.435541 Jy ! Inverting map ! Added new window around map position (-19.5, -1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00931058 Jy ! Component: 100 - total flux cleaned = 0.012912 Jy ! Component: 150 - total flux cleaned = 0.0148563 Jy ! Component: 200 - total flux cleaned = 0.0164207 Jy ! Total flux subtracted in 200 components = 0.0164207 Jy ! Clean residual min=-0.007173 max=0.008958 Jy/beam ! Clean residual mean=0.000002 rms=0.001948 Jy/beam ! Combined flux in latest and established models = 0.451962 Jy ! Adding 35 model components to the UV plane model. ! The established model now contains 82 components and 0.451962 Jy ! Inverting map ! Added new window around map position (32.4, 56.1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00539813 Jy ! Component: 100 - total flux cleaned = 0.00831714 Jy ! Component: 150 - total flux cleaned = 0.00991699 Jy ! Component: 200 - total flux cleaned = 0.0114356 Jy ! Total flux subtracted in 200 components = 0.0114356 Jy ! Clean residual min=-0.006518 max=0.008542 Jy/beam ! Clean residual mean=0.000002 rms=0.001823 Jy/beam ! Combined flux in latest and established models = 0.463397 Jy ! Adding 46 model components to the UV plane model. ! The established model now contains 105 components and 0.463397 Jy ! Inverting map ! Added new window around map position (43.5, 29.1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00632483 Jy ! Component: 100 - total flux cleaned = 0.00949285 Jy ! Component: 150 - total flux cleaned = 0.0120809 Jy ! Component: 200 - total flux cleaned = 0.0132329 Jy ! Total flux subtracted in 200 components = 0.0132329 Jy ! Clean residual min=-0.006281 max=0.007376 Jy/beam ! Clean residual mean=0.000002 rms=0.001709 Jy/beam ! Combined flux in latest and established models = 0.47663 Jy ! Adding 46 model components to the UV plane model. ! The established model now contains 120 components and 0.47663 Jy ! Inverting map ! Added new window around map position (-14.7, 43.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00459837 Jy ! Component: 100 - total flux cleaned = 0.00614188 Jy ! Component: 150 - total flux cleaned = 0.00760631 Jy ! Component: 200 - total flux cleaned = 0.00813384 Jy ! Total flux subtracted in 200 components = 0.00813384 Jy ! Clean residual min=-0.005732 max=0.006560 Jy/beam ! Clean residual mean=0.000003 rms=0.001637 Jy/beam ! Combined flux in latest and established models = 0.484764 Jy ! Adding 41 model components to the UV plane model. ! The established model now contains 140 components and 0.484764 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 40 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.029 mas, bmaj=5.806 mas, bpa=4.898 degrees ! Estimated noise=0.534686 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00137796 Jy ! Component: 100 - total flux cleaned = 0.00124675 Jy ! Component: 150 - total flux cleaned = 0.00163374 Jy ! Component: 200 - total flux cleaned = 0.00211997 Jy ! Total flux subtracted in 200 components = 0.00211997 Jy ! Clean residual min=-0.006167 max=0.005757 Jy/beam ! Clean residual mean=0.000006 rms=0.001599 Jy/beam ! Combined flux in latest and established models = 0.486884 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.000588201 Jy ! Component: 100 - total flux cleaned = 0.00115488 Jy ! Component: 150 - total flux cleaned = 0.001482 Jy ! Component: 200 - total flux cleaned = 0.00211597 Jy ! Total flux subtracted in 200 components = 0.00211597 Jy ! Clean residual min=-0.006017 max=0.005409 Jy/beam ! Clean residual mean=0.000006 rms=0.001539 Jy/beam ! Combined flux in latest and established models = 0.489 Jy ! Adding 43 model components to the UV plane model. ! The established model now contains 175 components and 0.489 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.864 mas, bmaj=4.228 mas, bpa=-11.73 degrees ! Estimated noise=0.462866 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00229135 Jy ! Component: 100 - total flux cleaned = 0.00326649 Jy ! Component: 150 - total flux cleaned = 0.00358248 Jy ! Component: 200 - total flux cleaned = 0.00402837 Jy ! Total flux subtracted in 200 components = 0.00402837 Jy ! Clean residual min=-0.005317 max=0.006593 Jy/beam ! Clean residual mean=0.000003 rms=0.001519 Jy/beam ! Combined flux in latest and established models = 0.493028 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 194 components and 0.493028 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.166893Jy sigma=0.887211 ! Fit after self-cal, rms=0.166593Jy sigma=0.882472 wmodel J1351+0830_C_map.mod ! Writing 194 model components to file: J1351+0830_C_map.mod wobs J1351+0830_C_uvs.fits ! Writing UV FITS file: J1351+0830_C_uvs.fits wwins J1351+0830_C_map.win ! wwins: Wrote 11 windows to J1351+0830_C_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.00287851 Jy ! Component: 100 - total flux cleaned = 0.00546935 Jy ! Component: 150 - total flux cleaned = 0.00781624 Jy ! Component: 200 - total flux cleaned = 0.00950442 Jy ! Component: 250 - total flux cleaned = 0.0110474 Jy ! Component: 300 - total flux cleaned = 0.0127217 Jy ! Component: 350 - total flux cleaned = 0.0142716 Jy ! Component: 400 - total flux cleaned = 0.0159597 Jy ! Component: 450 - total flux cleaned = 0.0173685 Jy ! Component: 500 - total flux cleaned = 0.0189169 Jy ! Component: 550 - total flux cleaned = 0.0201196 Jy ! Component: 600 - total flux cleaned = 0.0214627 Jy ! Component: 650 - total flux cleaned = 0.0227849 Jy ! Component: 700 - total flux cleaned = 0.0240872 Jy ! Component: 750 - total flux cleaned = 0.0252202 Jy ! Component: 800 - total flux cleaned = 0.0264122 Jy ! Component: 850 - total flux cleaned = 0.0275156 Jy ! Component: 900 - total flux cleaned = 0.0287494 Jy ! Component: 950 - total flux cleaned = 0.0296803 Jy ! Component: 1000 - total flux cleaned = 0.0307423 Jy ! Component: 1050 - total flux cleaned = 0.031792 Jy ! Component: 1100 - total flux cleaned = 0.0327601 Jy ! Component: 1150 - total flux cleaned = 0.0336483 Jy ! Component: 1200 - total flux cleaned = 0.0345282 Jy ! Component: 1250 - total flux cleaned = 0.0354643 Jy ! Component: 1300 - total flux cleaned = 0.0362587 Jy ! Component: 1350 - total flux cleaned = 0.037177 Jy ! Component: 1400 - total flux cleaned = 0.0380863 Jy ! Component: 1450 - total flux cleaned = 0.0388583 Jy ! Component: 1500 - total flux cleaned = 0.0396866 Jy ! Component: 1550 - total flux cleaned = 0.0405712 Jy ! Component: 1600 - total flux cleaned = 0.0413228 Jy ! Component: 1650 - total flux cleaned = 0.0421294 Jy ! Component: 1700 - total flux cleaned = 0.0428671 Jy ! Component: 1750 - total flux cleaned = 0.043599 Jy ! Component: 1800 - total flux cleaned = 0.0442642 Jy ! Component: 1850 - total flux cleaned = 0.0449245 Jy ! Component: 1900 - total flux cleaned = 0.0458763 Jy ! Component: 1950 - total flux cleaned = 0.0464078 Jy ! Component: 2000 - total flux cleaned = 0.0471111 Jy ! Component: 2050 - total flux cleaned = 0.0476927 Jy ! Component: 2100 - total flux cleaned = 0.0485585 Jy ! Component: 2150 - total flux cleaned = 0.0491309 Jy ! Component: 2200 - total flux cleaned = 0.0497564 Jy ! Component: 2250 - total flux cleaned = 0.0503215 Jy ! Component: 2300 - total flux cleaned = 0.0510507 Jy ! Component: 2350 - total flux cleaned = 0.0516081 Jy ! Component: 2400 - total flux cleaned = 0.0522176 Jy ! Component: 2450 - total flux cleaned = 0.0529872 Jy ! Component: 2500 - total flux cleaned = 0.053643 Jy ! Component: 2550 - total flux cleaned = 0.0542946 Jy ! Component: 2600 - total flux cleaned = 0.0550495 Jy ! Component: 2650 - total flux cleaned = 0.0557457 Jy ! Component: 2700 - total flux cleaned = 0.056438 Jy ! Component: 2750 - total flux cleaned = 0.0570197 Jy ! Component: 2800 - total flux cleaned = 0.0577039 Jy ! Component: 2850 - total flux cleaned = 0.0583834 Jy ! Component: 2900 - total flux cleaned = 0.0591632 Jy ! Component: 2950 - total flux cleaned = 0.0598341 Jy ! Component: 3000 - total flux cleaned = 0.0605528 Jy ! Component: 3050 - total flux cleaned = 0.0612165 Jy ! Component: 3100 - total flux cleaned = 0.0618758 Jy ! Component: 3150 - total flux cleaned = 0.062481 Jy ! Component: 3200 - total flux cleaned = 0.0631829 Jy ! Component: 3250 - total flux cleaned = 0.0638808 Jy ! Component: 3300 - total flux cleaned = 0.0645747 Jy ! Component: 3350 - total flux cleaned = 0.0650183 Jy ! Component: 3400 - total flux cleaned = 0.0657537 Jy ! Component: 3450 - total flux cleaned = 0.0662896 Jy ! Component: 3500 - total flux cleaned = 0.0669682 Jy ! Component: 3550 - total flux cleaned = 0.0675946 Jy ! Component: 3600 - total flux cleaned = 0.0681218 Jy ! Component: 3650 - total flux cleaned = 0.068694 Jy ! Component: 3700 - total flux cleaned = 0.0693104 Jy ! Component: 3750 - total flux cleaned = 0.0698291 Jy ! Component: 3800 - total flux cleaned = 0.0702983 Jy ! Component: 3850 - total flux cleaned = 0.0708581 Jy ! Component: 3900 - total flux cleaned = 0.0713686 Jy ! Component: 3950 - total flux cleaned = 0.0718762 Jy ! Component: 4000 - total flux cleaned = 0.0724731 Jy ! Component: 4050 - total flux cleaned = 0.0728386 Jy ! Total flux subtracted in 4096 components = 0.0733842 Jy ! Clean residual min=-0.002564 max=0.002454 Jy/beam ! Clean residual mean=0.000000 rms=0.000900 Jy/beam ! Combined flux in latest and established models = 0.566413 Jy keep ! Adding 268 model components to the UV plane model. ! The established model now contains 462 components and 0.566412 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap .//J1351+0830_C_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.864 x 4.228 at -11.73 degrees (North through East) ! Clean map min=-0.0076609 max=0.26685 Jy/beam ! Writing clean map to FITS file: .//J1351+0830_C_map.fits wmodel .//J1351+0830_C_map.mod ! Writing 462 model components to file: .//J1351+0830_C_map.mod wobs .//J1351+0830_C_uvs.fits ! Writing UV FITS file: .//J1351+0830_C_uvs.fits wwins .//J1351+0830_C_map.win ! wwins: Wrote 12 windows to .//J1351+0830_C_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.266851 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000904493 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.01685 1.01685 2.0337 4.0674 8.1348 16.2696 32.5392 65.0784 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0027 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.266 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 295.028 quit ! Quitting program ! Log file .//J1351+0830_C_dfm.log closed on Mon Jan 20 18:09:56 2020