! Started logfile: J0340+4732_S_dfm.log on Fri Feb 15 19:43:15 2019 float field_size field_size = 1024 float field_cell field_cell = 0.5000 float taper_size taper_size = 20.0 float freq freq = 2.22000 obs J0340+4732_S_uvm.fits ! Reading UV FITS file: J0340+4732_S_uvm.fits ! AN table 1: 66 integrations on 36 of 36 possible baselines. ! AN table 2: 48 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.880511 visibilities/baseline/integration-bin. ! Found source: J0340+4732 ! ! There are 4 IFs, and a total of 4 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 ! 04 4 2.348e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 54 lines of history. ! ! Reading 15976 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 4 channels. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 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 J0340+4732_S] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff print "===============================================================" ! =============================================================== ! print "=== Difmap script: pima_mupet_01.dvm version of 2016.05.19 === " ! === Difmap script: pima_mupet_01.dvm version of 2016.05.19 === ! mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 45 telescope corrections were flagged in sub-array 1. ! A total of 19 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 45 telescope corrections were flagged in sub-array 1. ! A total of 19 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 45 telescope corrections were flagged in sub-array 1. ! A total of 19 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 330 telescope corrections were flagged in sub-array 1. ! A total of 480 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.052272Jy sigma=18.275396 ! Fit after self-cal, rms=0.709769Jy sigma=12.678044 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. sele I,1,3 ! Polarization I is unavailable. ! Selecting polarization: RR, channels: 1..3 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: (none) 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.928 mas, bmaj=5.554 mas, bpa=1.197 degrees ! Estimated noise=0.61184 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.228447 Jy ! Component: 100 - total flux cleaned = 0.28429 Jy ! Total flux subtracted in 100 components = 0.28429 Jy ! Clean residual min=-0.011635 max=0.021709 Jy/beam ! Clean residual mean=0.000036 rms=0.002681 Jy/beam ! Combined flux in latest and established models = 0.28429 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 0.28429 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.111720Jy sigma=1.256118 ! Fit after self-cal, rms=0.111723Jy sigma=1.256016 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0174419 Jy ! Component: 100 - total flux cleaned = 0.0173573 Jy ! Total flux subtracted in 100 components = 0.0173573 Jy ! Clean residual min=-0.010558 max=0.010780 Jy/beam ! Clean residual mean=0.000037 rms=0.002218 Jy/beam ! Combined flux in latest and established models = 0.301647 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 12 components and 0.301647 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.106527Jy sigma=1.184385 ! Fit after self-cal, rms=0.106528Jy sigma=1.184329 ! 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.548 mas, bmaj=6.445 mas, bpa=2.472 degrees ! Estimated noise=0.380573 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.000560458 Jy ! Component: 100 - total flux cleaned = -0.000576289 Jy ! Component: 150 - total flux cleaned = -0.00058013 Jy ! Component: 200 - total flux cleaned = -0.000578837 Jy ! Total flux subtracted in 200 components = -0.000578837 Jy ! Clean residual min=-0.011769 max=0.009520 Jy/beam ! Clean residual mean=0.000077 rms=0.002331 Jy/beam ! Combined flux in latest and established models = 0.301068 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 22 components and 0.301068 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.105173Jy sigma=1.172462 ! Fit after self-cal, rms=0.105169Jy sigma=1.172431 ! 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.548 x 6.445 at 2.472 degrees (North through East) ! Clean map min=-0.01152 max=0.29245 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 0.92 HN 1.01 KP 0.95 ! LA 0.93 MK 0.98 NL 1.06 OV 0.93 ! PT 0.93 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 1.15 HN 1.09 KP 1.09 ! LA 1.01 MK 1.00 NL 1.17 OV 1.09 ! PT 1.10 SC 0.99 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 0.92 HN 1.04 KP 0.92 ! LA 0.93 MK 0.98 NL 0.81 OV 0.96 ! PT 0.93 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.76 FD 1.14 HN 1.10 KP 1.08 ! LA 1.02 MK 1.01 NL 1.26 OV 1.08 ! PT 1.11 SC 0.89 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 0.92 HN 1.02 KP 0.92 ! LA 0.88 MK 0.98 NL 0.97 OV 0.94 ! PT 0.93 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.13 HN 1.09 KP 1.11 ! LA 1.00 MK 1.00 NL 1.15 OV 1.07 ! PT 1.10 SC 0.95 ! ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.105169Jy sigma=1.172431 ! Fit after self-cal, rms=0.076771Jy sigma=0.860902 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.541 mas, bmaj=6.383 mas, bpa=1.257 degrees ! Estimated noise=0.382149 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.002306 Jy ! Component: 100 - total flux cleaned = 0.00230921 Jy ! Component: 150 - total flux cleaned = 0.00231117 Jy ! Component: 200 - total flux cleaned = 0.00239301 Jy ! Total flux subtracted in 200 components = 0.00239301 Jy ! Clean residual min=-0.003894 max=0.003124 Jy/beam ! Clean residual mean=0.000013 rms=0.000825 Jy/beam ! Combined flux in latest and established models = 0.303461 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 26 components and 0.303461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.076317Jy sigma=0.850438 ! Fit after self-cal, rms=0.076298Jy sigma=0.850260 ! 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. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.076298Jy sigma=0.850260 ! Fit after self-cal, rms=0.081805Jy sigma=0.839334 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.538 mas, bmaj=6.392 mas, bpa=0.5449 degrees ! Estimated noise=0.384634 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00251378 Jy ! Component: 100 - total flux cleaned = 0.00241457 Jy ! Component: 150 - total flux cleaned = 0.00241272 Jy ! Component: 200 - total flux cleaned = 0.0024079 Jy ! Total flux subtracted in 200 components = 0.0024079 Jy ! Clean residual min=-0.003278 max=0.002722 Jy/beam ! Clean residual mean=0.000011 rms=0.000744 Jy/beam ! Combined flux in latest and established models = 0.305869 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 35 components and 0.305869 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.081354Jy sigma=0.835824 ! Fit after self-cal, rms=0.081356Jy sigma=0.835810 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -1.55886e-07 Jy ! Component: 100 - total flux cleaned = 3.98865e-05 Jy ! Component: 150 - total flux cleaned = 7.33629e-05 Jy ! Component: 200 - total flux cleaned = 9.98672e-05 Jy ! Total flux subtracted in 200 components = 9.98672e-05 Jy ! Clean residual min=-0.003323 max=0.002725 Jy/beam ! Clean residual mean=0.000011 rms=0.000740 Jy/beam ! Combined flux in latest and established models = 0.305969 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 40 components and 0.305969 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.081262Jy sigma=0.835424 ! Fit after self-cal, rms=0.081262Jy sigma=0.835420 ! 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. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.081262Jy sigma=0.835420 ! Fit after self-cal, rms=0.047796Jy sigma=0.760998 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.547 mas, bmaj=6.381 mas, bpa=0.1363 degrees ! Estimated noise=0.384729 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0012236 Jy ! Component: 100 - total flux cleaned = 0.00123077 Jy ! Component: 150 - total flux cleaned = 0.00128242 Jy ! Component: 200 - total flux cleaned = 0.00128411 Jy ! Total flux subtracted in 200 components = 0.00128411 Jy ! Clean residual min=-0.001625 max=0.001800 Jy/beam ! Clean residual mean=0.000003 rms=0.000364 Jy/beam ! Combined flux in latest and established models = 0.307253 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 47 components and 0.307253 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047651Jy sigma=0.758639 ! Fit after self-cal, rms=0.047629Jy sigma=0.758484 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 (85.5, -5.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000991494 Jy ! Component: 100 - total flux cleaned = 0.000778636 Jy ! Component: 150 - total flux cleaned = 0.000622978 Jy ! Component: 200 - total flux cleaned = 0.000484714 Jy ! Total flux subtracted in 200 components = 0.000484714 Jy ! Clean residual min=-0.001535 max=0.001443 Jy/beam ! Clean residual mean=0.000003 rms=0.000350 Jy/beam ! Combined flux in latest and established models = 0.307738 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 63 components and 0.307738 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047561Jy sigma=0.757457 ! Fit after self-cal, rms=0.047554Jy sigma=0.757392 ! 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. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047554Jy sigma=0.757392 ! Fit after self-cal, rms=0.047420Jy sigma=0.756626 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.919 mas, bmaj=10.31 mas, bpa=-7.358 degrees ! Estimated noise=0.470199 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000595608 Jy ! Component: 100 - total flux cleaned = 0.000614993 Jy ! Component: 150 - total flux cleaned = 0.000615465 Jy ! Component: 200 - total flux cleaned = 0.000615841 Jy ! Total flux subtracted in 200 components = 0.000615841 Jy ! Clean residual min=-0.001365 max=0.001200 Jy/beam ! Clean residual mean=0.000006 rms=0.000350 Jy/beam ! Combined flux in latest and established models = 0.308353 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 70 components and 0.308353 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047397Jy sigma=0.756271 ! Fit after self-cal, rms=0.047393Jy sigma=0.756246 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 = 7.05222e-05 Jy ! Component: 100 - total flux cleaned = 3.97319e-06 Jy ! Component: 150 - total flux cleaned = -1.87508e-05 Jy ! Component: 200 - total flux cleaned = -2.96409e-05 Jy ! Total flux subtracted in 200 components = -2.96409e-05 Jy ! Clean residual min=-0.001359 max=0.001183 Jy/beam ! Clean residual mean=0.000006 rms=0.000349 Jy/beam ! Combined flux in latest and established models = 0.308324 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 75 components and 0.308324 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047388Jy sigma=0.756174 ! Fit after self-cal, rms=0.047387Jy sigma=0.756165 ! 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.549 mas, bmaj=6.375 mas, bpa=-0.1162 degrees ! Estimated noise=0.384729 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.549 x 6.375 at -0.1162 degrees (North through East) ! Clean map min=-0.0014105 max=0.29577 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.916 mas, bmaj=5.538 mas, bpa=-1.339 degrees ! Estimated noise=0.617782 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.230655 Jy ! Component: 100 - total flux cleaned = 0.283096 Jy ! Component: 150 - total flux cleaned = 0.300962 Jy ! Component: 200 - total flux cleaned = 0.308065 Jy ! Total flux subtracted in 200 components = 0.308065 Jy ! Clean residual min=-0.002488 max=0.002905 Jy/beam ! Clean residual mean=-0.000000 rms=0.000618 Jy/beam ! Combined flux in latest and established models = 0.308066 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 10 model components to the UV plane model. ! The established model now contains 10 components and 0.308066 Jy ! Inverting map and beam ! Estimated beam: bmin=3.549 mas, bmaj=6.375 mas, bpa=-0.1162 degrees ! Estimated noise=0.384729 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.00227173 Jy ! Component: 100 - total flux cleaned = 0.0024938 Jy ! Component: 150 - total flux cleaned = 0.00238873 Jy ! Component: 200 - total flux cleaned = 0.00234505 Jy ! Total flux subtracted in 200 components = 0.00234505 Jy ! Clean residual min=-0.001618 max=0.001339 Jy/beam ! Clean residual mean=0.000002 rms=0.000338 Jy/beam ! Combined flux in latest and established models = 0.310411 Jy ! Adding 16 model components to the UV plane model. ! The established model now contains 22 components and 0.310411 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.919 mas, bmaj=10.31 mas, bpa=-7.358 degrees ! Estimated noise=0.470199 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000763381 Jy ! Component: 100 - total flux cleaned = -0.00104399 Jy ! Component: 150 - total flux cleaned = -0.00116005 Jy ! Component: 200 - total flux cleaned = -0.00122041 Jy ! Total flux subtracted in 200 components = -0.00122041 Jy ! Clean residual min=-0.001391 max=0.001104 Jy/beam ! Clean residual mean=0.000005 rms=0.000362 Jy/beam ! Combined flux in latest and established models = 0.30919 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 = -5.6607e-05 Jy ! Component: 100 - total flux cleaned = -8.37208e-05 Jy ! Component: 150 - total flux cleaned = -0.000121512 Jy ! Component: 200 - total flux cleaned = -0.000156784 Jy ! Total flux subtracted in 200 components = -0.000156784 Jy ! Clean residual min=-0.001389 max=0.001119 Jy/beam ! Clean residual mean=0.000005 rms=0.000359 Jy/beam ! Combined flux in latest and established models = 0.309034 Jy ! Adding 21 model components to the UV plane model. ! The established model now contains 40 components and 0.309034 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.549 mas, bmaj=6.375 mas, bpa=-0.1162 degrees ! Estimated noise=0.384729 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000115209 Jy ! Component: 100 - total flux cleaned = 5.98393e-05 Jy ! Component: 150 - total flux cleaned = -1.51848e-05 Jy ! Component: 200 - total flux cleaned = -6.03179e-05 Jy ! Total flux subtracted in 200 components = -6.03179e-05 Jy ! Clean residual min=-0.001435 max=0.001301 Jy/beam ! Clean residual mean=0.000002 rms=0.000327 Jy/beam ! Combined flux in latest and established models = 0.308973 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 53 components and 0.308973 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Not correcting unselected IF 4. ! ! Fit before self-cal, rms=0.047344Jy sigma=0.755653 ! Fit after self-cal, rms=0.047341Jy sigma=0.755626 wmodel J0340+4732_S_map.mod ! Writing 53 model components to file: J0340+4732_S_map.mod wobs J0340+4732_S_uvs.fits ! Writing UV FITS file: J0340+4732_S_uvs.fits wwins J0340+4732_S_map.win ! wwins: Wrote 2 windows to J0340+4732_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.000187379 Jy ! Component: 100 - total flux cleaned = -0.000234294 Jy ! Component: 150 - total flux cleaned = -0.000302166 Jy ! Component: 200 - total flux cleaned = -0.000345733 Jy ! Component: 250 - total flux cleaned = -0.000366581 Jy ! Component: 300 - total flux cleaned = -0.000366886 Jy ! Component: 350 - total flux cleaned = -0.000387362 Jy ! Component: 400 - total flux cleaned = -0.000426602 Jy ! Component: 450 - total flux cleaned = -0.000387846 Jy ! Component: 500 - total flux cleaned = -0.000500858 Jy ! Component: 550 - total flux cleaned = -0.000518837 Jy ! Component: 600 - total flux cleaned = -0.000537269 Jy ! Component: 650 - total flux cleaned = -0.000537594 Jy ! Component: 700 - total flux cleaned = -0.000502831 Jy ! Component: 750 - total flux cleaned = -0.000554855 Jy ! Component: 800 - total flux cleaned = -0.000486499 Jy ! Component: 850 - total flux cleaned = -0.000452755 Jy ! Component: 900 - total flux cleaned = -0.000436035 Jy ! Component: 950 - total flux cleaned = -0.000386665 Jy ! Component: 1000 - total flux cleaned = -0.000402911 Jy ! Component: 1050 - total flux cleaned = -0.000386635 Jy ! Component: 1100 - total flux cleaned = -0.000275763 Jy ! Component: 1150 - total flux cleaned = -0.000244457 Jy ! Component: 1200 - total flux cleaned = -0.000197893 Jy ! Component: 1250 - total flux cleaned = -0.000121049 Jy ! Component: 1300 - total flux cleaned = 1.56708e-05 Jy ! Component: 1350 - total flux cleaned = 9.10382e-05 Jy ! Component: 1400 - total flux cleaned = 0.000165657 Jy ! Component: 1450 - total flux cleaned = 0.000224827 Jy ! Component: 1500 - total flux cleaned = 0.000341759 Jy ! Component: 1550 - total flux cleaned = 0.000442988 Jy ! Component: 1600 - total flux cleaned = 0.000471525 Jy ! Component: 1650 - total flux cleaned = 0.000570911 Jy ! Component: 1700 - total flux cleaned = 0.000627149 Jy ! Component: 1750 - total flux cleaned = 0.000655024 Jy ! Component: 1800 - total flux cleaned = 0.00071032 Jy ! Component: 1850 - total flux cleaned = 0.000833772 Jy ! Component: 1900 - total flux cleaned = 0.00080651 Jy ! Component: 1950 - total flux cleaned = 0.000847018 Jy ! Component: 2000 - total flux cleaned = 0.0009407 Jy ! Component: 2050 - total flux cleaned = 0.000927479 Jy ! Component: 2100 - total flux cleaned = 0.00092759 Jy ! Component: 2150 - total flux cleaned = 0.00104503 Jy ! Component: 2200 - total flux cleaned = 0.000993195 Jy ! Component: 2250 - total flux cleaned = 0.000993381 Jy ! Component: 2300 - total flux cleaned = 0.00113366 Jy ! Component: 2350 - total flux cleaned = 0.00103223 Jy ! Component: 2400 - total flux cleaned = 0.00108242 Jy ! Component: 2450 - total flux cleaned = 0.00113234 Jy ! Component: 2500 - total flux cleaned = 0.00116955 Jy ! Component: 2550 - total flux cleaned = 0.00119411 Jy ! Component: 2600 - total flux cleaned = 0.0012552 Jy ! Component: 2650 - total flux cleaned = 0.00130361 Jy ! Component: 2700 - total flux cleaned = 0.00135177 Jy ! Component: 2750 - total flux cleaned = 0.00138745 Jy ! Component: 2800 - total flux cleaned = 0.00138745 Jy ! Component: 2850 - total flux cleaned = 0.00137567 Jy ! Component: 2900 - total flux cleaned = 0.0014692 Jy ! Component: 2950 - total flux cleaned = 0.00153873 Jy ! Component: 3000 - total flux cleaned = 0.00155018 Jy ! Component: 3050 - total flux cleaned = 0.00155018 Jy ! Component: 3100 - total flux cleaned = 0.00156161 Jy ! Component: 3150 - total flux cleaned = 0.00160678 Jy ! Component: 3200 - total flux cleaned = 0.00165171 Jy ! Component: 3250 - total flux cleaned = 0.00169624 Jy ! Component: 3300 - total flux cleaned = 0.00168529 Jy ! Component: 3350 - total flux cleaned = 0.00172917 Jy ! Component: 3400 - total flux cleaned = 0.00179444 Jy ! Component: 3450 - total flux cleaned = 0.00182699 Jy ! Component: 3500 - total flux cleaned = 0.00185927 Jy ! Component: 3550 - total flux cleaned = 0.0019126 Jy ! Component: 3600 - total flux cleaned = 0.00195501 Jy ! Component: 3650 - total flux cleaned = 0.00191273 Jy ! Component: 3700 - total flux cleaned = 0.00195467 Jy ! Component: 3750 - total flux cleaned = 0.00204837 Jy ! Component: 3800 - total flux cleaned = 0.00201725 Jy ! Component: 3850 - total flux cleaned = 0.00207902 Jy ! Component: 3900 - total flux cleaned = 0.00207907 Jy ! Component: 3950 - total flux cleaned = 0.0021401 Jy ! Component: 4000 - total flux cleaned = 0.00216028 Jy ! Component: 4050 - total flux cleaned = 0.00215016 Jy ! Total flux subtracted in 4096 components = 0.00218018 Jy ! Clean residual min=-0.000595 max=0.000581 Jy/beam ! Clean residual mean=-0.000000 rms=0.000183 Jy/beam ! Combined flux in latest and established models = 0.311153 Jy keep ! Adding 297 model components to the UV plane model. ! The established model now contains 350 components and 0.311153 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap J0340+4732_S_map.fits ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.549 x 6.375 at -0.1162 degrees (North through East) ! Clean map min=-0.0010141 max=0.29585 Jy/beam ! Writing clean map to FITS file: J0340+4732_S_map.fits wmodel J0340+4732_S_map.mod ! Writing 350 model components to file: J0340+4732_S_map.mod wobs J0340+4732_S_uvs.fits ! Writing UV FITS file: J0340+4732_S_uvs.fits wwins J0340+4732_S_map.win ! wwins: Wrote 3 windows to J0340+4732_S_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.295845 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000182366 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.184927 0.184927 0.369854 0.739708 1.47942 2.95883 5.91767 11.8353 23.6707 47.3413 94.6827 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0005 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.295 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1622.26 quit ! Quitting program ! Log file J0340+4732_S_dfm.log closed on Fri Feb 15 19:43:20 2019