! Started logfile: /image_prod/ug002_uva//J1950-0436_S_dfm.log on Thu Mar 5 04:25:31 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//J1950-0436_S_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J1950-0436_S_uva.fits ! AN table 1: 90 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.990123 visibilities/baseline/integration-bin. ! Found source: J1950-0436 ! ! 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 9624 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 J1950-0436_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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=1.191778Jy sigma=7.269496 ! Fit after self-cal, rms=0.592579Jy sigma=4.016564 ! 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.913 mas, bmaj=12.84 mas, bpa=-18.42 degrees ! Estimated noise=1.67816 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.357969 Jy ! Component: 100 - total flux cleaned = 0.529333 Jy ! Total flux subtracted in 100 components = 0.529333 Jy ! Clean residual min=-0.046752 max=0.115145 Jy/beam ! Clean residual mean=0.000179 rms=0.015732 Jy/beam ! Combined flux in latest and established models = 0.529333 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 27 components and 0.529333 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.360727Jy sigma=1.934617 ! Fit after self-cal, rms=0.358442Jy sigma=1.923209 ! Inverting map ! Added new window around map position (-25.5, 51.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.127674 Jy ! Component: 100 - total flux cleaned = 0.215111 Jy ! Total flux subtracted in 100 components = 0.215111 Jy ! Clean residual min=-0.036025 max=0.062879 Jy/beam ! Clean residual mean=0.000135 rms=0.011098 Jy/beam ! Combined flux in latest and established models = 0.744445 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 50 components and 0.744445 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.264809Jy sigma=1.422478 ! Fit after self-cal, rms=0.262763Jy sigma=1.402335 ! 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.472 mas, bmaj=14.17 mas, bpa=-17.67 degrees ! Estimated noise=1.18015 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.112836 Jy ! Component: 100 - total flux cleaned = 0.189622 Jy ! Component: 150 - total flux cleaned = 0.247349 Jy ! Component: 200 - total flux cleaned = 0.29219 Jy ! Total flux subtracted in 200 components = 0.29219 Jy ! Clean residual min=-0.028263 max=0.076967 Jy/beam ! Clean residual mean=0.000119 rms=0.010402 Jy/beam ! Combined flux in latest and established models = 1.03664 Jy ! Performing phase self-cal ! Adding 92 model components to the UV plane model. ! The established model now contains 135 components and 1.03663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.174049Jy sigma=0.962081 ! Fit after self-cal, rms=0.169935Jy sigma=0.926074 ! Inverting map ! Added new window around map position (-16.5, 49). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0762907 Jy ! Component: 100 - total flux cleaned = 0.122422 Jy ! Component: 150 - total flux cleaned = 0.158059 Jy ! Component: 200 - total flux cleaned = 0.18752 Jy ! Total flux subtracted in 200 components = 0.18752 Jy ! Clean residual min=-0.016806 max=0.025450 Jy/beam ! Clean residual mean=0.000063 rms=0.005712 Jy/beam ! Combined flux in latest and established models = 1.22415 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 178 components and 1.22415 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.112333Jy sigma=0.631826 ! Fit after self-cal, rms=0.108227Jy sigma=0.604581 ! 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.472 x 14.17 at -17.67 degrees (North through East) ! Clean map min=-0.014779 max=0.46115 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 0.95 HN 1.04 KP 0.99 ! LA 0.97 MK 1.00 NL 0.98 OV 0.99 ! PT 0.94 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 0.93 HN 1.05 KP 1.00 ! LA 0.98 MK 0.98 NL 0.98 OV 0.99 ! PT 0.95 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 0.98 HN 1.03 KP 0.99 ! LA 0.98 MK 1.00 NL 0.95 OV 0.99 ! PT 0.95 ! ! ! Fit before self-cal, rms=0.108227Jy sigma=0.604581 ! Fit after self-cal, rms=0.095419Jy sigma=0.558669 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.493 mas, bmaj=14.39 mas, bpa=-17.73 degrees ! Estimated noise=1.15229 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.0252031 Jy ! Component: 100 - total flux cleaned = 0.0423412 Jy ! Component: 150 - total flux cleaned = 0.0537837 Jy ! Component: 200 - total flux cleaned = 0.0610901 Jy ! Total flux subtracted in 200 components = 0.0610901 Jy ! Clean residual min=-0.011158 max=0.015011 Jy/beam ! Clean residual mean=0.000006 rms=0.003229 Jy/beam ! Combined flux in latest and established models = 1.28525 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 216 components and 1.28525 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.085126Jy sigma=0.500554 ! Fit after self-cal, rms=0.084500Jy sigma=0.495911 ! 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.084500Jy sigma=0.495911 ! Fit after self-cal, rms=0.082057Jy sigma=0.478817 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.471 mas, bmaj=14.49 mas, bpa=-17.62 degrees ! Estimated noise=1.14759 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0120216 Jy ! Component: 100 - total flux cleaned = 0.0198253 Jy ! Component: 150 - total flux cleaned = 0.0248825 Jy ! Component: 200 - total flux cleaned = 0.0278527 Jy ! Total flux subtracted in 200 components = 0.0278527 Jy ! Clean residual min=-0.008223 max=0.012055 Jy/beam ! Clean residual mean=0.000009 rms=0.002441 Jy/beam ! Combined flux in latest and established models = 1.3131 Jy selfcal ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 259 components and 1.3131 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.078767Jy sigma=0.459299 ! Fit after self-cal, rms=0.078623Jy sigma=0.458138 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 (7, -128). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00971302 Jy ! Component: 100 - total flux cleaned = 0.0154946 Jy ! Component: 150 - total flux cleaned = 0.019295 Jy ! Component: 200 - total flux cleaned = 0.0208972 Jy ! Total flux subtracted in 200 components = 0.0208972 Jy ! Clean residual min=-0.006910 max=0.009373 Jy/beam ! Clean residual mean=0.000006 rms=0.002076 Jy/beam ! Combined flux in latest and established models = 1.33399 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 291 components and 1.33399 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.076254Jy sigma=0.444648 ! Fit after self-cal, rms=0.076151Jy sigma=0.443922 ! Inverting map ! Added new window around map position (-5.5, 46.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00997518 Jy ! Component: 100 - total flux cleaned = 0.0167196 Jy ! Component: 150 - total flux cleaned = 0.0220431 Jy ! Component: 200 - total flux cleaned = 0.0255832 Jy ! Total flux subtracted in 200 components = 0.0255832 Jy ! Clean residual min=-0.006754 max=0.007964 Jy/beam ! Clean residual mean=-0.000000 rms=0.001741 Jy/beam ! Combined flux in latest and established models = 1.35958 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 328 components and 1.35958 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.074051Jy sigma=0.432394 ! Fit after self-cal, rms=0.073950Jy sigma=0.432102 ! Inverting map ! Added new window around map position (-6.5, -6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00564986 Jy ! Component: 100 - total flux cleaned = 0.00821772 Jy ! Component: 150 - total flux cleaned = 0.0107168 Jy ! Component: 200 - total flux cleaned = 0.0124892 Jy ! Total flux subtracted in 200 components = 0.0124892 Jy ! Clean residual min=-0.006028 max=0.005671 Jy/beam ! Clean residual mean=-0.000003 rms=0.001554 Jy/beam ! Combined flux in latest and established models = 1.37207 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 363 components and 1.37207 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.072605Jy sigma=0.424706 ! Fit after self-cal, rms=0.072499Jy sigma=0.424424 ! 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.072499Jy sigma=0.424424 ! Fit after self-cal, rms=0.071667Jy sigma=0.421295 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.478 mas, bmaj=14.49 mas, bpa=-17.59 degrees ! Estimated noise=1.14775 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00432925 Jy ! Component: 100 - total flux cleaned = 0.00764945 Jy ! Component: 150 - total flux cleaned = 0.0102656 Jy ! Component: 200 - total flux cleaned = 0.0122269 Jy ! Total flux subtracted in 200 components = 0.0122269 Jy ! Clean residual min=-0.005197 max=0.005554 Jy/beam ! Clean residual mean=-0.000008 rms=0.001343 Jy/beam ! Combined flux in latest and established models = 1.38429 Jy selfcal ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 405 components and 1.38429 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.070736Jy sigma=0.416266 ! Fit after self-cal, rms=0.070680Jy sigma=0.416104 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.00196853 Jy ! Component: 100 - total flux cleaned = 0.00319807 Jy ! Component: 150 - total flux cleaned = 0.00410522 Jy ! Component: 200 - total flux cleaned = 0.00497062 Jy ! Total flux subtracted in 200 components = 0.00497062 Jy ! Clean residual min=-0.005019 max=0.005174 Jy/beam ! Clean residual mean=-0.000009 rms=0.001254 Jy/beam ! Combined flux in latest and established models = 1.38926 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 427 components and 1.38926 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.070166Jy sigma=0.413513 ! Fit after self-cal, rms=0.070131Jy sigma=0.413443 ! 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.070131Jy sigma=0.413443 ! Fit after self-cal, rms=0.069903Jy sigma=0.412905 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.397 mas, bmaj=15.19 mas, bpa=-11.65 degrees ! Estimated noise=1.45043 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.003497 Jy ! Component: 100 - total flux cleaned = 0.00574837 Jy ! Component: 150 - total flux cleaned = 0.00755049 Jy ! Component: 200 - total flux cleaned = 0.00904168 Jy ! Total flux subtracted in 200 components = 0.00904168 Jy ! Clean residual min=-0.005071 max=0.005763 Jy/beam ! Clean residual mean=-0.000016 rms=0.001477 Jy/beam ! Combined flux in latest and established models = 1.39831 Jy selfcal ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 471 components and 1.39831 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.069507Jy sigma=0.410662 ! Fit after self-cal, rms=0.069443Jy sigma=0.410462 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.00137199 Jy ! Component: 100 - total flux cleaned = 0.00217137 Jy ! Component: 150 - total flux cleaned = 0.00270181 Jy ! Component: 200 - total flux cleaned = 0.0034006 Jy ! Total flux subtracted in 200 components = 0.0034006 Jy ! Clean residual min=-0.004615 max=0.005463 Jy/beam ! Clean residual mean=-0.000016 rms=0.001364 Jy/beam ! Combined flux in latest and established models = 1.40171 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 495 components and 1.40171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.069209Jy sigma=0.409420 ! Fit after self-cal, rms=0.069176Jy sigma=0.409345 ! 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.478 mas, bmaj=14.48 mas, bpa=-17.59 degrees ! Estimated noise=1.1473 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.478 x 14.48 at -17.59 degrees (North through East) ! Clean map min=-0.0056893 max=0.44385 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.913 mas, bmaj=13.1 mas, bpa=-18.33 degrees ! Estimated noise=1.63853 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.347442 Jy ! Component: 100 - total flux cleaned = 0.52781 Jy ! Component: 150 - total flux cleaned = 0.649777 Jy ! Component: 200 - total flux cleaned = 0.740879 Jy ! Total flux subtracted in 200 components = 0.740879 Jy ! Clean residual min=-0.028176 max=0.054019 Jy/beam ! Clean residual mean=0.000113 rms=0.009908 Jy/beam ! Combined flux in latest and established models = 0.740878 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 50 model components to the UV plane model. ! The established model now contains 50 components and 0.740878 Jy ! Inverting map and beam ! Estimated beam: bmin=3.478 mas, bmaj=14.48 mas, bpa=-17.59 degrees ! Estimated noise=1.1473 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.113853 Jy ! Component: 100 - total flux cleaned = 0.198426 Jy ! Component: 150 - total flux cleaned = 0.264408 Jy ! Component: 200 - total flux cleaned = 0.31773 Jy ! Total flux subtracted in 200 components = 0.31773 Jy ! Clean residual min=-0.018768 max=0.032027 Jy/beam ! Clean residual mean=0.000103 rms=0.007233 Jy/beam ! Combined flux in latest and established models = 1.05861 Jy ! Adding 101 model components to the UV plane model. ! The established model now contains 134 components and 1.05861 Jy ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0436767 Jy ! Component: 100 - total flux cleaned = 0.0798397 Jy ! Component: 150 - total flux cleaned = 0.110177 Jy ! Component: 200 - total flux cleaned = 0.135988 Jy ! Total flux subtracted in 200 components = 0.135988 Jy ! Clean residual min=-0.011531 max=0.015922 Jy/beam ! Clean residual mean=0.000054 rms=0.004145 Jy/beam ! Combined flux in latest and established models = 1.1946 Jy ! Adding 78 model components to the UV plane model. ! The established model now contains 202 components and 1.1946 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.397 mas, bmaj=15.19 mas, bpa=-11.65 degrees ! Estimated noise=1.45043 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0360209 Jy ! Component: 100 - total flux cleaned = 0.0631432 Jy ! Component: 150 - total flux cleaned = 0.0845345 Jy ! Component: 200 - total flux cleaned = 0.101994 Jy ! Total flux subtracted in 200 components = 0.101994 Jy ! Clean residual min=-0.010705 max=0.012361 Jy/beam ! Clean residual mean=0.000038 rms=0.003602 Jy/beam ! Combined flux in latest and established models = 1.29659 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.0147231 Jy ! Component: 100 - total flux cleaned = 0.0273592 Jy ! Component: 150 - total flux cleaned = 0.0383434 Jy ! Component: 200 - total flux cleaned = 0.047928 Jy ! Total flux subtracted in 200 components = 0.047928 Jy ! Clean residual min=-0.007123 max=0.008337 Jy/beam ! Clean residual mean=0.000014 rms=0.002340 Jy/beam ! Combined flux in latest and established models = 1.34452 Jy ! Adding 195 model components to the UV plane model. ! The established model now contains 368 components and 1.34452 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.478 mas, bmaj=14.48 mas, bpa=-17.59 degrees ! Estimated noise=1.1473 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00714697 Jy ! Component: 100 - total flux cleaned = 0.0125802 Jy ! Component: 150 - total flux cleaned = 0.0174356 Jy ! Component: 200 - total flux cleaned = 0.0216162 Jy ! Total flux subtracted in 200 components = 0.0216162 Jy ! Clean residual min=-0.005110 max=0.005960 Jy/beam ! Clean residual mean=-0.000001 rms=0.001398 Jy/beam ! Combined flux in latest and established models = 1.36614 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 396 components and 1.36614 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.070811Jy sigma=0.418429 ! Fit after self-cal, rms=0.070585Jy sigma=0.416369 wmodel J1950-0436_S_map.mod ! Writing 396 model components to file: J1950-0436_S_map.mod wobs J1950-0436_S_uvs.fits ! Writing UV FITS file: J1950-0436_S_uvs.fits wwins J1950-0436_S_map.win ! wwins: Wrote 6 windows to J1950-0436_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.00181325 Jy ! Component: 100 - total flux cleaned = 0.00275868 Jy ! Component: 150 - total flux cleaned = 0.00324791 Jy ! Component: 200 - total flux cleaned = 0.00348303 Jy ! Component: 250 - total flux cleaned = 0.00370892 Jy ! Component: 300 - total flux cleaned = 0.00356217 Jy ! Component: 350 - total flux cleaned = 0.00349018 Jy ! Component: 400 - total flux cleaned = 0.00321574 Jy ! Component: 450 - total flux cleaned = 0.00315126 Jy ! Component: 500 - total flux cleaned = 0.00341259 Jy ! Component: 550 - total flux cleaned = 0.00322042 Jy ! Component: 600 - total flux cleaned = 0.00315892 Jy ! Component: 650 - total flux cleaned = 0.00322012 Jy ! Component: 700 - total flux cleaned = 0.00328001 Jy ! Component: 750 - total flux cleaned = 0.003338 Jy ! Component: 800 - total flux cleaned = 0.0032798 Jy ! Component: 850 - total flux cleaned = 0.00339297 Jy ! Component: 900 - total flux cleaned = 0.00344837 Jy ! Component: 950 - total flux cleaned = 0.00339332 Jy ! Component: 1000 - total flux cleaned = 0.00344793 Jy ! Component: 1050 - total flux cleaned = 0.0034485 Jy ! Component: 1100 - total flux cleaned = 0.00350118 Jy ! Component: 1150 - total flux cleaned = 0.00350118 Jy ! Component: 1200 - total flux cleaned = 0.00335014 Jy ! Component: 1250 - total flux cleaned = 0.00339899 Jy ! Component: 1300 - total flux cleaned = 0.00325256 Jy ! Component: 1350 - total flux cleaned = 0.00320659 Jy ! Component: 1400 - total flux cleaned = 0.00344185 Jy ! Component: 1450 - total flux cleaned = 0.00334689 Jy ! Component: 1500 - total flux cleaned = 0.00311655 Jy ! Component: 1550 - total flux cleaned = 0.00302564 Jy ! Component: 1600 - total flux cleaned = 0.00302552 Jy ! Component: 1650 - total flux cleaned = 0.00298159 Jy ! Component: 1700 - total flux cleaned = 0.00298205 Jy ! Component: 1750 - total flux cleaned = 0.00293909 Jy ! Component: 1800 - total flux cleaned = 0.00298178 Jy ! Component: 1850 - total flux cleaned = 0.00298186 Jy ! Component: 1900 - total flux cleaned = 0.00310445 Jy ! Component: 1950 - total flux cleaned = 0.00302215 Jy ! Component: 2000 - total flux cleaned = 0.00286213 Jy ! Component: 2050 - total flux cleaned = 0.00294121 Jy ! Component: 2100 - total flux cleaned = 0.00294105 Jy ! Component: 2150 - total flux cleaned = 0.0029022 Jy ! Component: 2200 - total flux cleaned = 0.00297779 Jy ! Component: 2250 - total flux cleaned = 0.00290208 Jy ! Component: 2300 - total flux cleaned = 0.00286548 Jy ! Component: 2350 - total flux cleaned = 0.00293848 Jy ! Component: 2400 - total flux cleaned = 0.00293838 Jy ! Component: 2450 - total flux cleaned = 0.0028671 Jy ! Component: 2500 - total flux cleaned = 0.00276176 Jy ! Component: 2550 - total flux cleaned = 0.00286669 Jy ! Component: 2600 - total flux cleaned = 0.00286613 Jy ! Component: 2650 - total flux cleaned = 0.00279865 Jy ! Component: 2700 - total flux cleaned = 0.00276495 Jy ! Component: 2750 - total flux cleaned = 0.00276558 Jy ! Component: 2800 - total flux cleaned = 0.00273279 Jy ! Component: 2850 - total flux cleaned = 0.00270128 Jy ! Component: 2900 - total flux cleaned = 0.00273395 Jy ! Component: 2950 - total flux cleaned = 0.00279757 Jy ! Component: 3000 - total flux cleaned = 0.00279771 Jy ! Component: 3050 - total flux cleaned = 0.00270406 Jy ! Component: 3100 - total flux cleaned = 0.00276574 Jy ! Component: 3150 - total flux cleaned = 0.00267474 Jy ! Component: 3200 - total flux cleaned = 0.00273524 Jy ! Component: 3250 - total flux cleaned = 0.00270512 Jy ! Component: 3300 - total flux cleaned = 0.00273423 Jy ! Component: 3350 - total flux cleaned = 0.00270479 Jy ! Component: 3400 - total flux cleaned = 0.00270485 Jy ! Component: 3450 - total flux cleaned = 0.00267611 Jy ! Component: 3500 - total flux cleaned = 0.00273309 Jy ! Component: 3550 - total flux cleaned = 0.00267649 Jy ! Component: 3600 - total flux cleaned = 0.00267648 Jy ! Component: 3650 - total flux cleaned = 0.00259435 Jy ! Component: 3700 - total flux cleaned = 0.0025674 Jy ! Component: 3750 - total flux cleaned = 0.00264847 Jy ! Component: 3800 - total flux cleaned = 0.0025411 Jy ! Component: 3850 - total flux cleaned = 0.00256766 Jy ! Component: 3900 - total flux cleaned = 0.00254145 Jy ! Component: 3950 - total flux cleaned = 0.00259373 Jy ! Component: 4000 - total flux cleaned = 0.00249056 Jy ! Component: 4050 - total flux cleaned = 0.00249049 Jy ! Total flux subtracted in 4096 components = 0.00249046 Jy ! Clean residual min=-0.001467 max=0.001532 Jy/beam ! Clean residual mean=0.000002 rms=0.000528 Jy/beam ! Combined flux in latest and established models = 1.36863 Jy keep ! Adding 342 model components to the UV plane model. ! The established model now contains 734 components and 1.36863 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J1950-0436_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.478 x 14.48 at -17.59 degrees (North through East) ! Clean map min=-0.0031851 max=0.44117 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J1950-0436_S_map.fits wmodel /image_prod/ug002_uva//J1950-0436_S_map.mod ! Writing 734 model components to file: /image_prod/ug002_uva//J1950-0436_S_map.mod wobs /image_prod/ug002_uva//J1950-0436_S_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J1950-0436_S_uvs.fits wwins /image_prod/ug002_uva//J1950-0436_S_map.win ! wwins: Wrote 7 windows to /image_prod/ug002_uva//J1950-0436_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.441173 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000522063 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.355006 0.355006 0.710012 1.42002 2.84005 5.6801 11.3602 22.7204 45.4408 90.8815 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0015 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.441 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 845.056 quit ! Quitting program ! Log file /image_prod/ug002_uva//J1950-0436_S_dfm.log closed on Thu Mar 5 04:25:36 2020