! Started logfile: /image_prod/ug002_uva//J1619+2247_S_a_dfm.log on Thu Mar 5 02:38:15 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//J1619+2247_S_a_uva.fits ! Reading UV FITS file: /image_prod/ug002_uva//J1619+2247_S_a_uva.fits ! AN table 1: 18 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.722222 visibilities/baseline/integration-bin. ! Found source: J1619+2247 ! ! 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 1755 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 J1619+2247_S_a] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff print "===============================================================" ! =============================================================== ! print "=== Difmap script: pima_mupet_01.dvm version of 2016.05.19 === " ! === Difmap script: pima_mupet_01.dvm version of 2016.05.19 === ! mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 28 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 28 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 28 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.199813Jy sigma=7.814694 ! Fit after self-cal, rms=0.384938Jy sigma=2.594769 ! 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.884 mas, bmaj=5.984 mas, bpa=-7.169 degrees ! Estimated noise=3.06089 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.499493 Jy ! Component: 100 - total flux cleaned = 0.608415 Jy ! Total flux subtracted in 100 components = 0.608415 Jy ! Clean residual min=-0.022259 max=0.030377 Jy/beam ! Clean residual mean=-0.000002 rms=0.005927 Jy/beam ! Combined flux in latest and established models = 0.608415 Jy ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 0.608415 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.162761Jy sigma=0.819026 ! Fit after self-cal, rms=0.162758Jy sigma=0.819014 ! 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.359 mas, bmaj=6.644 mas, bpa=-6.468 degrees ! Estimated noise=2.61993 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.0234348 Jy ! Component: 100 - total flux cleaned = 0.0296166 Jy ! Component: 150 - total flux cleaned = 0.0296279 Jy ! Component: 200 - total flux cleaned = 0.029053 Jy ! Total flux subtracted in 200 components = 0.029053 Jy ! Clean residual min=-0.029077 max=0.025922 Jy/beam ! Clean residual mean=0.000033 rms=0.006462 Jy/beam ! Combined flux in latest and established models = 0.637468 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 27 components and 0.637468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.157279Jy sigma=0.799701 ! Fit after self-cal, rms=0.157253Jy sigma=0.799548 ! 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.359 x 6.644 at -6.468 degrees (North through East) ! Clean map min=-0.028488 max=0.63678 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.96 HN 0.98 KP 0.98 ! LA 1.05 MK 0.98 NL 1.10 OV 1.03 ! PT 1.00 SC 1.01 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 0.96 HN 0.99 KP 0.98 ! LA 1.03 MK 1.01 NL 1.05 OV 0.95 ! PT 1.01 SC 0.98 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.91 FD 0.99 HN 1.00 KP 0.99 ! LA 1.03 MK 1.00 NL 1.14 OV 0.98 ! PT 0.98 SC 1.03 ! ! ! Fit before self-cal, rms=0.157253Jy sigma=0.799548 ! Fit after self-cal, rms=0.160769Jy sigma=0.763903 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.357 mas, bmaj=6.627 mas, bpa=-6.288 degrees ! Estimated noise=2.63353 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.00920373 Jy ! Component: 100 - total flux cleaned = 0.0106175 Jy ! Component: 150 - total flux cleaned = 0.010405 Jy ! Component: 200 - total flux cleaned = 0.0104077 Jy ! Total flux subtracted in 200 components = 0.0104077 Jy ! Clean residual min=-0.033009 max=0.021984 Jy/beam ! Clean residual mean=0.000011 rms=0.006162 Jy/beam ! Combined flux in latest and established models = 0.647876 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 36 components and 0.647875 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.159529Jy sigma=0.760086 ! Fit after self-cal, rms=0.159522Jy sigma=0.760005 ! 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.159522Jy sigma=0.760005 ! Fit after self-cal, rms=0.064855Jy sigma=0.359270 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.359 mas, bmaj=6.644 mas, bpa=-6.167 degrees ! Estimated noise=2.63705 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00310598 Jy ! Component: 100 - total flux cleaned = 0.00297949 Jy ! Component: 150 - total flux cleaned = 0.00262843 Jy ! Component: 200 - total flux cleaned = 0.00221032 Jy ! Total flux subtracted in 200 components = 0.00221032 Jy ! Clean residual min=-0.005701 max=0.006646 Jy/beam ! Clean residual mean=-0.000005 rms=0.001367 Jy/beam ! Combined flux in latest and established models = 0.650086 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 48 components and 0.650086 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.064288Jy sigma=0.355498 ! Fit after self-cal, rms=0.064447Jy sigma=0.355030 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 (-83, -68). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00550141 Jy ! Component: 100 - total flux cleaned = 0.00770708 Jy ! Component: 150 - total flux cleaned = 0.00830994 Jy ! Component: 200 - total flux cleaned = 0.00848932 Jy ! Total flux subtracted in 200 components = 0.00848932 Jy ! Clean residual min=-0.005267 max=0.006278 Jy/beam ! Clean residual mean=-0.000012 rms=0.001263 Jy/beam ! Combined flux in latest and established models = 0.658575 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 67 components and 0.658575 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.063883Jy sigma=0.350804 ! Fit after self-cal, rms=0.063863Jy sigma=0.350581 ! Inverting map ! Added new window around map position (44, 54.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00334883 Jy ! Component: 100 - total flux cleaned = 0.0035671 Jy ! Component: 150 - total flux cleaned = 0.00385544 Jy ! Component: 200 - total flux cleaned = 0.00368101 Jy ! Total flux subtracted in 200 components = 0.00368101 Jy ! Clean residual min=-0.005230 max=0.005839 Jy/beam ! Clean residual mean=-0.000015 rms=0.001187 Jy/beam ! Combined flux in latest and established models = 0.662256 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 83 components and 0.662256 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.063289Jy sigma=0.347017 ! Fit after self-cal, rms=0.063266Jy sigma=0.346750 ! Inverting map ! Added new window around map position (19.5, -34.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00395695 Jy ! Component: 100 - total flux cleaned = 0.00508352 Jy ! Component: 150 - total flux cleaned = 0.00547614 Jy ! Component: 200 - total flux cleaned = 0.00593216 Jy ! Total flux subtracted in 200 components = 0.00593216 Jy ! Clean residual min=-0.005006 max=0.004298 Jy/beam ! Clean residual mean=-0.000017 rms=0.001101 Jy/beam ! Combined flux in latest and established models = 0.668188 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 102 components and 0.668188 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062718Jy sigma=0.343370 ! Fit after self-cal, rms=0.062717Jy sigma=0.343266 ! 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.062717Jy sigma=0.343266 ! Fit after self-cal, rms=0.062974Jy sigma=0.341354 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.353 mas, bmaj=6.646 mas, bpa=-6.158 degrees ! Estimated noise=2.639 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00126491 Jy ! Component: 100 - total flux cleaned = 0.00170966 Jy ! Component: 150 - total flux cleaned = 0.00203617 Jy ! Component: 200 - total flux cleaned = 0.00257415 Jy ! Total flux subtracted in 200 components = 0.00257415 Jy ! Clean residual min=-0.004200 max=0.004125 Jy/beam ! Clean residual mean=-0.000017 rms=0.001016 Jy/beam ! Combined flux in latest and established models = 0.670763 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 112 components and 0.670763 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062733Jy sigma=0.339791 ! Fit after self-cal, rms=0.062744Jy sigma=0.339741 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.000708225 Jy ! Component: 100 - total flux cleaned = 0.000991776 Jy ! Component: 150 - total flux cleaned = 0.0011953 Jy ! Component: 200 - total flux cleaned = 0.0013931 Jy ! Total flux subtracted in 200 components = 0.0013931 Jy ! Clean residual min=-0.004135 max=0.003968 Jy/beam ! Clean residual mean=-0.000018 rms=0.000995 Jy/beam ! Combined flux in latest and established models = 0.672156 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 131 components and 0.672156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062590Jy sigma=0.338829 ! Fit after self-cal, rms=0.062593Jy sigma=0.338811 ! 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.062593Jy sigma=0.338811 ! Fit after self-cal, rms=0.062653Jy sigma=0.338713 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.738 mas, bmaj=10.04 mas, bpa=-8.459 degrees ! Estimated noise=3.36085 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.731e-06 Jy ! Component: 100 - total flux cleaned = 0.000262969 Jy ! Component: 150 - total flux cleaned = 0.000768635 Jy ! Component: 200 - total flux cleaned = 0.00111193 Jy ! Total flux subtracted in 200 components = 0.00111193 Jy ! Clean residual min=-0.004867 max=0.005186 Jy/beam ! Clean residual mean=-0.000038 rms=0.001325 Jy/beam ! Combined flux in latest and established models = 0.673268 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 147 components and 0.673267 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062610Jy sigma=0.338651 ! Fit after self-cal, rms=0.062604Jy sigma=0.338603 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.000451698 Jy ! Component: 100 - total flux cleaned = 0.000540621 Jy ! Component: 150 - total flux cleaned = 0.000661916 Jy ! Component: 200 - total flux cleaned = 0.000775456 Jy ! Total flux subtracted in 200 components = 0.000775456 Jy ! Clean residual min=-0.004760 max=0.005201 Jy/beam ! Clean residual mean=-0.000039 rms=0.001303 Jy/beam ! Combined flux in latest and established models = 0.674043 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 155 components and 0.674043 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062582Jy sigma=0.338624 ! Fit after self-cal, rms=0.062582Jy sigma=0.338609 ! 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.351 mas, bmaj=6.647 mas, bpa=-6.16 degrees ! Estimated noise=2.64088 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.351 x 6.647 at -6.16 degrees (North through East) ! Clean map min=-0.0038924 max=0.64163 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 ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.062582Jy sigma=0.338609 ! Fit after self-cal, rms=0.059238Jy sigma=0.315341 ! Inverting map and beam ! Estimated beam: bmin=3.352 mas, bmaj=6.65 mas, bpa=-6.175 degrees ! Estimated noise=2.64291 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000853946 Jy ! Component: 100 - total flux cleaned = 0.00134946 Jy ! Component: 150 - total flux cleaned = 0.00160744 Jy ! Component: 200 - total flux cleaned = 0.00178904 Jy ! Total flux subtracted in 200 components = 0.00178904 Jy ! Clean residual min=-0.003992 max=0.003859 Jy/beam ! Clean residual mean=-0.000018 rms=0.000954 Jy/beam ! Combined flux in latest and established models = 0.675832 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 169 components and 0.675832 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.059087Jy sigma=0.314326 ! Fit after self-cal, rms=0.059090Jy sigma=0.314307 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000425396 Jy ! Component: 100 - total flux cleaned = 0.000539827 Jy ! Component: 150 - total flux cleaned = 0.00054058 Jy ! Component: 200 - total flux cleaned = 0.000595328 Jy ! Total flux subtracted in 200 components = 0.000595328 Jy ! Clean residual min=-0.003935 max=0.003790 Jy/beam ! Clean residual mean=-0.000019 rms=0.000942 Jy/beam ! Combined flux in latest and established models = 0.676427 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 174 components and 0.676427 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.058985Jy sigma=0.313672 ! Fit after self-cal, rms=0.058986Jy sigma=0.313665 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.058986Jy sigma=0.313665 ! Fit after self-cal, rms=0.059027Jy sigma=0.313579 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.885 mas, bmaj=5.982 mas, bpa=-7.123 degrees ! Estimated noise=3.05977 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.499259 Jy ! Component: 100 - total flux cleaned = 0.608248 Jy ! Component: 150 - total flux cleaned = 0.635706 Jy ! Component: 200 - total flux cleaned = 0.647809 Jy ! Total flux subtracted in 200 components = 0.647809 Jy ! Clean residual min=-0.007416 max=0.006756 Jy/beam ! Clean residual mean=-0.000010 rms=0.001614 Jy/beam ! Combined flux in latest and established models = 0.647809 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 18 model components to the UV plane model. ! The established model now contains 18 components and 0.647809 Jy ! Inverting map and beam ! Estimated beam: bmin=3.351 mas, bmaj=6.651 mas, bpa=-6.168 degrees ! Estimated noise=2.6449 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.00804422 Jy ! Component: 100 - total flux cleaned = 0.0132972 Jy ! Component: 150 - total flux cleaned = 0.0156806 Jy ! Component: 200 - total flux cleaned = 0.0171097 Jy ! Total flux subtracted in 200 components = 0.0171097 Jy ! Clean residual min=-0.006038 max=0.005378 Jy/beam ! Clean residual mean=-0.000007 rms=0.001232 Jy/beam ! Combined flux in latest and established models = 0.664919 Jy ! Adding 32 model components to the UV plane model. ! The established model now contains 41 components and 0.664919 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.738 mas, bmaj=10.05 mas, bpa=-8.439 degrees ! Estimated noise=3.36598 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.003701 Jy ! Component: 100 - total flux cleaned = 0.00587236 Jy ! Component: 150 - total flux cleaned = 0.00748538 Jy ! Component: 200 - total flux cleaned = 0.00822749 Jy ! Total flux subtracted in 200 components = 0.00822749 Jy ! Clean residual min=-0.004957 max=0.005638 Jy/beam ! Clean residual mean=-0.000028 rms=0.001366 Jy/beam ! Combined flux in latest and established models = 0.673146 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.000503995 Jy ! Component: 100 - total flux cleaned = 0.000887628 Jy ! Component: 150 - total flux cleaned = 0.00124265 Jy ! Component: 200 - total flux cleaned = 0.00163501 Jy ! Total flux subtracted in 200 components = 0.00163501 Jy ! Clean residual min=-0.004854 max=0.005513 Jy/beam ! Clean residual mean=-0.000031 rms=0.001324 Jy/beam ! Combined flux in latest and established models = 0.674781 Jy ! Adding 40 model components to the UV plane model. ! The established model now contains 79 components and 0.674781 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.351 mas, bmaj=6.651 mas, bpa=-6.168 degrees ! Estimated noise=2.6449 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000121326 Jy ! Component: 100 - total flux cleaned = -0.000116835 Jy ! Component: 150 - total flux cleaned = 0.00012939 Jy ! Component: 200 - total flux cleaned = 0.000246469 Jy ! Total flux subtracted in 200 components = 0.000246469 Jy ! Clean residual min=-0.004595 max=0.004592 Jy/beam ! Clean residual mean=-0.000015 rms=0.001087 Jy/beam ! Combined flux in latest and established models = 0.675027 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 98 components and 0.675027 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Fit before self-cal, rms=0.060009Jy sigma=0.321136 ! Fit after self-cal, rms=0.059942Jy sigma=0.320711 wmodel J1619+2247_S_a_map.mod ! Writing 98 model components to file: J1619+2247_S_a_map.mod wobs J1619+2247_S_a_uvs.fits ! Writing UV FITS file: J1619+2247_S_a_uvs.fits wwins J1619+2247_S_a_map.win ! wwins: Wrote 4 windows to J1619+2247_S_a_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000935542 Jy ! Component: 100 - total flux cleaned = -0.00109158 Jy ! Component: 150 - total flux cleaned = -0.000576319 Jy ! Component: 200 - total flux cleaned = -0.000152614 Jy ! Component: 250 - total flux cleaned = 0.000458658 Jy ! Component: 300 - total flux cleaned = 0.000522392 Jy ! Component: 350 - total flux cleaned = 0.000522648 Jy ! Component: 400 - total flux cleaned = 0.000523365 Jy ! Component: 450 - total flux cleaned = 0.000464944 Jy ! Component: 500 - total flux cleaned = 0.000584282 Jy ! Component: 550 - total flux cleaned = 0.000469027 Jy ! Component: 600 - total flux cleaned = 0.000413425 Jy ! Component: 650 - total flux cleaned = 0.000523968 Jy ! Component: 700 - total flux cleaned = 0.000522277 Jy ! Component: 750 - total flux cleaned = 0.000416072 Jy ! Component: 800 - total flux cleaned = 0.000467949 Jy ! Component: 850 - total flux cleaned = 0.000569328 Jy ! Component: 900 - total flux cleaned = 0.000519425 Jy ! Component: 950 - total flux cleaned = 0.000471835 Jy ! Component: 1000 - total flux cleaned = 0.000664147 Jy ! Component: 1050 - total flux cleaned = 0.000663825 Jy ! Component: 1100 - total flux cleaned = 0.00066361 Jy ! Component: 1150 - total flux cleaned = 0.000617903 Jy ! Component: 1200 - total flux cleaned = 0.000751836 Jy ! Component: 1250 - total flux cleaned = 0.000618365 Jy ! Component: 1300 - total flux cleaned = 0.000574751 Jy ! Component: 1350 - total flux cleaned = 0.000618569 Jy ! Component: 1400 - total flux cleaned = 0.000660702 Jy ! Component: 1450 - total flux cleaned = 0.000577364 Jy ! Component: 1500 - total flux cleaned = 0.000699975 Jy ! Component: 1550 - total flux cleaned = 0.000659752 Jy ! Component: 1600 - total flux cleaned = 0.000580512 Jy ! Component: 1650 - total flux cleaned = 0.000502109 Jy ! Component: 1700 - total flux cleaned = 0.00050283 Jy ! Component: 1750 - total flux cleaned = 0.000312823 Jy ! Component: 1800 - total flux cleaned = 0.0001639 Jy ! Component: 1850 - total flux cleaned = 9.02669e-05 Jy ! Component: 1900 - total flux cleaned = 1.74313e-05 Jy ! Component: 1950 - total flux cleaned = 1.70448e-05 Jy ! Component: 2000 - total flux cleaned = -0.000231708 Jy ! Component: 2050 - total flux cleaned = -0.000266541 Jy ! Component: 2100 - total flux cleaned = -0.00019728 Jy ! Component: 2150 - total flux cleaned = -0.000265655 Jy ! Component: 2200 - total flux cleaned = -9.68139e-05 Jy ! Component: 2250 - total flux cleaned = -0.000196562 Jy ! Component: 2300 - total flux cleaned = -3.16237e-05 Jy ! Component: 2350 - total flux cleaned = -0.000293322 Jy ! Component: 2400 - total flux cleaned = -0.000163149 Jy ! Component: 2450 - total flux cleaned = -0.00022644 Jy ! Component: 2500 - total flux cleaned = -0.000195216 Jy ! Component: 2550 - total flux cleaned = -0.00019487 Jy ! Component: 2600 - total flux cleaned = -0.000318489 Jy ! Component: 2650 - total flux cleaned = -0.000288385 Jy ! Component: 2700 - total flux cleaned = -0.000349006 Jy ! Component: 2750 - total flux cleaned = -0.000408833 Jy ! Component: 2800 - total flux cleaned = -0.000379562 Jy ! Component: 2850 - total flux cleaned = -0.000380276 Jy ! Component: 2900 - total flux cleaned = -0.000409272 Jy ! Component: 2950 - total flux cleaned = -0.000351343 Jy ! Component: 3000 - total flux cleaned = -0.000408196 Jy ! Component: 3050 - total flux cleaned = -0.000379914 Jy ! Component: 3100 - total flux cleaned = -0.000323771 Jy ! Component: 3150 - total flux cleaned = -0.000406824 Jy ! Component: 3200 - total flux cleaned = -0.000296624 Jy ! Component: 3250 - total flux cleaned = -0.000215076 Jy ! Component: 3300 - total flux cleaned = -0.000242274 Jy ! Component: 3350 - total flux cleaned = -0.000188456 Jy ! Component: 3400 - total flux cleaned = -0.000161656 Jy ! Component: 3450 - total flux cleaned = -0.000109352 Jy ! Component: 3500 - total flux cleaned = -0.000161547 Jy ! Component: 3550 - total flux cleaned = -0.000213419 Jy ! Component: 3600 - total flux cleaned = -0.000162132 Jy ! Component: 3650 - total flux cleaned = -0.000289366 Jy ! Component: 3700 - total flux cleaned = -0.000289748 Jy ! Component: 3750 - total flux cleaned = -0.000389837 Jy ! Component: 3800 - total flux cleaned = -0.000414652 Jy ! Component: 3850 - total flux cleaned = -0.000439509 Jy ! Component: 3900 - total flux cleaned = -0.000463875 Jy ! Component: 3950 - total flux cleaned = -0.000609458 Jy ! Component: 4000 - total flux cleaned = -0.000657642 Jy ! Component: 4050 - total flux cleaned = -0.000729023 Jy ! Total flux subtracted in 4096 components = -0.000847255 Jy ! Clean residual min=-0.001471 max=0.001404 Jy/beam ! Clean residual mean=-0.000003 rms=0.000428 Jy/beam ! Combined flux in latest and established models = 0.67418 Jy keep ! Adding 296 model components to the UV plane model. ! The established model now contains 394 components and 0.67418 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /image_prod/ug002_uva//J1619+2247_S_a_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.351 x 6.651 at -6.168 degrees (North through East) ! Clean map min=-0.002939 max=0.64082 Jy/beam ! Writing clean map to FITS file: /image_prod/ug002_uva//J1619+2247_S_a_map.fits wmodel /image_prod/ug002_uva//J1619+2247_S_a_map.mod ! Writing 394 model components to file: /image_prod/ug002_uva//J1619+2247_S_a_map.mod wobs /image_prod/ug002_uva//J1619+2247_S_a_uvs.fits ! Writing UV FITS file: /image_prod/ug002_uva//J1619+2247_S_a_uvs.fits wwins /image_prod/ug002_uva//J1619+2247_S_a_map.win ! wwins: Wrote 5 windows to /image_prod/ug002_uva//J1619+2247_S_a_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.640825 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000426303 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.199572 0.199572 0.399144 0.798289 1.59658 3.19316 6.38631 12.7726 25.5452 51.0905 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0012 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.64 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1503.21 quit ! Quitting program ! Log file /image_prod/ug002_uva//J1619+2247_S_a_dfm.log closed on Thu Mar 5 02:38:19 2020