! Started logfile: /scr/pima/rv102_uvs/J1642+6856_X_dfm.log on Sat Nov 9 10:20:34 2019 float field_size field_size = 1024 float field_cell field_cell = 0.2000 float taper_size taper_size = 70.0 float freq freq = 8.40599 obs /scr/pima/rv102_uvs/J1642+6856_X_uva.fits ! Reading UV FITS file: /scr/pima/rv102_uvs/J1642+6856_X_uva.fits ! AN table 1: 72 integrations on 28 of 28 possible baselines. ! AN table 2: 11 integrations on 10 of 10 possible baselines. ! AN table 3: 11 integrations on 55 of 55 possible baselines. ! AN table 4: 74 integrations on 66 of 66 possible baselines. ! AN table 5: 74 integrations on 66 of 66 possible baselines. ! Apparent sampling: 0.962877 visibilities/baseline/integration-bin. ! Found source: J1642+6856 ! ! 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 8.40599e+09 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 79 lines of history. ! ! Reading 48140 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.2 print "taper_size =", taper_size ! taper_size = 70 ![@/opt64/bin/pima_mupet_01.dfm J1642+6856_X] 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.200x0.200 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 16 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 2. ! A total of 25 telescope corrections were flagged in sub-array 3. ! A total of 31 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 2. ! A total of 25 telescope corrections were flagged in sub-array 3. ! A total of 31 telescope corrections were flagged in sub-array 4. ! A total of 38 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 27 telescope corrections were flagged in sub-array 2. ! A total of 35 telescope corrections were flagged in sub-array 3. ! A total of 107 telescope corrections were flagged in sub-array 4. ! A total of 111 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 27 telescope corrections were flagged in sub-array 2. ! A total of 35 telescope corrections were flagged in sub-array 3. ! A total of 107 telescope corrections were flagged in sub-array 4. ! A total of 111 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.707156Jy sigma=20.053328 ! Fit after self-cal, rms=0.675547Jy sigma=7.926877 ! 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 ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.4872 mas, bmaj=0.5836 mas, bpa=-30.62 degrees ! Estimated noise=0.517721 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.723011 Jy ! Component: 100 - total flux cleaned = 1.04422 Jy ! Total flux subtracted in 100 components = 1.04422 Jy ! Clean residual min=-0.105698 max=0.194502 Jy/beam ! Clean residual mean=0.000357 rms=0.010611 Jy/beam ! Combined flux in latest and established models = 1.04422 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.04422 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.602920Jy sigma=6.720881 ! Fit after self-cal, rms=0.593925Jy sigma=6.625201 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.234972 Jy ! Component: 100 - total flux cleaned = 0.391811 Jy ! Total flux subtracted in 100 components = 0.391811 Jy ! Clean residual min=-0.070925 max=0.120329 Jy/beam ! Clean residual mean=0.000221 rms=0.006445 Jy/beam ! Combined flux in latest and established models = 1.43603 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 9 components and 1.43603 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.371504Jy sigma=4.153882 ! Fit after self-cal, rms=0.365189Jy sigma=4.083897 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-0.6, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.137117 Jy ! Component: 100 - total flux cleaned = 0.228001 Jy ! Total flux subtracted in 100 components = 0.228001 Jy ! Clean residual min=-0.034928 max=0.054154 Jy/beam ! Clean residual mean=0.000142 rms=0.004090 Jy/beam ! Combined flux in latest and established models = 1.66404 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 12 components and 1.66404 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.256428Jy sigma=2.859357 ! Fit after self-cal, rms=0.252811Jy sigma=2.819067 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0722605 Jy ! Component: 100 - total flux cleaned = 0.12155 Jy ! Total flux subtracted in 100 components = 0.12155 Jy ! Clean residual min=-0.026657 max=0.050468 Jy/beam ! Clean residual mean=0.000100 rms=0.003077 Jy/beam ! Combined flux in latest and established models = 1.78559 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 13 components and 1.78559 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210731Jy sigma=2.334641 ! Fit after self-cal, rms=0.209522Jy sigma=2.321830 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-0.2, -1). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0516505 Jy ! Component: 100 - total flux cleaned = 0.0814291 Jy ! Total flux subtracted in 100 components = 0.0814291 Jy ! Clean residual min=-0.022201 max=0.027222 Jy/beam ! Clean residual mean=0.000072 rms=0.002431 Jy/beam ! Combined flux in latest and established models = 1.86701 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 17 components and 1.86701 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.189340Jy sigma=2.071599 ! Fit after self-cal, rms=0.187415Jy sigma=2.052438 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0293987 Jy ! Component: 100 - total flux cleaned = 0.0429044 Jy ! Total flux subtracted in 100 components = 0.0429044 Jy ! Clean residual min=-0.017964 max=0.019096 Jy/beam ! Clean residual mean=0.000057 rms=0.002072 Jy/beam ! Combined flux in latest and established models = 1.90992 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 18 components and 1.90992 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179842Jy sigma=1.952235 ! Fit after self-cal, rms=0.179188Jy sigma=1.946087 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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 ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5797 mas, bmaj=0.7172 mas, bpa=-17.87 degrees ! Estimated noise=0.293714 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.0261548 Jy ! Component: 100 - total flux cleaned = 0.037491 Jy ! Component: 150 - total flux cleaned = 0.041712 Jy ! Component: 200 - total flux cleaned = 0.0440324 Jy ! Total flux subtracted in 200 components = 0.0440324 Jy ! Clean residual min=-0.012890 max=0.015085 Jy/beam ! Clean residual mean=0.000363 rms=0.002331 Jy/beam ! Combined flux in latest and established models = 1.95395 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 22 components and 1.95395 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.175369Jy sigma=1.887924 ! Fit after self-cal, rms=0.174828Jy sigma=1.882586 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-2.2, -6.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.011029 Jy ! Component: 100 - total flux cleaned = 0.0164899 Jy ! Component: 150 - total flux cleaned = 0.0194195 Jy ! Component: 200 - total flux cleaned = 0.0204147 Jy ! Total flux subtracted in 200 components = 0.0204147 Jy ! Clean residual min=-0.012109 max=0.012829 Jy/beam ! Clean residual mean=0.000310 rms=0.002141 Jy/beam ! Combined flux in latest and established models = 1.97437 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 38 components and 1.97437 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.173315Jy sigma=1.865637 ! Fit after self-cal, rms=0.173139Jy sigma=1.863264 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec print "============== Finished natural weighting clean ==============" ! ============== Finished natural weighting clean ============== ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5797 x 0.7172 at -17.87 degrees (North through East) ! Clean map min=-0.010272 max=1.0949 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! FD 0.99 HN 0.94 KP 1.02 LA 0.99 ! MK 0.97 NL 1.04 PT 1.03 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! MA 0.98 NY 1.01 ON 0.97 WF 0.92 ! WZ 0.96 ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.02 HN 1.03 KP 1.06 LA 1.01 ! MA 1.01 NL 1.01 NY 1.04 ON 0.99 ! PT 1.02 WF 0.90 WZ 0.92 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.01 HN 0.95 KP 1.00 LA 1.01 ! MA 1.02 MK 1.05 NY 1.03 ON 1.02 ! PT 1.00 SC 1.00 WF 1.06 WZ 1.01 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.99 HN 0.97 LA 1.00 MA 1.06 ! MK 0.98 NL 1.00 NY 0.99 ON 1.02 ! PT 0.94 SC 0.99 WF 1.03 WZ 0.97 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 0.99 HN 0.93 KP 1.01 LA 0.99 ! MK 0.96 NL 1.03 PT 1.04 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! MA 0.96 NY 1.01 ON 0.98 WF 0.86 ! WZ 0.91 ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.03 HN 1.03 KP 1.05 LA 1.00 ! MA 0.95 NL 1.01 NY 1.04 ON 1.01 ! PT 1.05 WF 0.87 WZ 0.89 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.01 HN 0.94 KP 1.01 LA 1.00 ! MA 0.98 MK 1.04 NY 1.03 ON 1.02 ! PT 1.01 SC 0.99 WF 1.08 WZ 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! FD 0.99 HN 0.95 LA 1.01 MA 1.02 ! MK 0.97 NL 0.98 NY 0.98 ON 1.03 ! PT 0.94 SC 0.99 WF 1.03 WZ 0.96 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 0.99 HN 0.97 KP 1.01 LA 1.00 ! MK 0.95 NL 1.04 PT 1.05 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! MA 1.06 NY 1.00* ON 1.01 WF 0.95 ! WZ 0.93 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.03 HN 1.01 KP 1.05 LA 1.00 ! MA 1.01 NL 1.03 NY 1.00* ON 1.03 ! PT 1.04 WF 0.87 WZ 0.91 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.01 HN 0.94 KP 1.02 LA 1.01 ! MA 1.04 MK 1.03 NY 1.00* ON 1.01 ! PT 1.01 SC 0.98 WF 1.08 WZ 0.98 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.98 HN 0.97 LA 1.01 MA 1.11 ! MK 0.98 NL 0.99 NY 1.00* ON 1.01 ! PT 0.94 SC 1.00 WF 1.02 WZ 0.95 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00 HN 0.97 KP 1.02 LA 1.00 ! MK 0.98 NL 1.05 PT 1.05 SC 0.99 ! ! Telescope amplitude corrections in sub-array 2: ! MA 1.15 NY 1.00* ON 1.00 WF 0.90 ! WZ 0.87 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.00 HN 0.99 KP 1.04 LA 1.01 ! MA 1.08 NL 1.02 NY 1.00* ON 1.01 ! PT 1.04 WF 0.87 WZ 0.87 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.01 HN 0.95 KP 1.02 LA 1.00 ! MA 1.11 MK 1.04 NY 1.00* ON 1.03 ! PT 1.02 SC 0.98 WF 1.07 WZ 0.95 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.99 HN 0.97 LA 1.01 MA 1.18 ! MK 1.00 NL 1.00 NY 1.00* ON 1.05 ! PT 0.94 SC 0.99 WF 1.03 WZ 0.92 ! ! ! Fit before self-cal, rms=0.173139Jy sigma=1.863264 ! Fit after self-cal, rms=0.163900Jy sigma=1.718601 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5788 mas, bmaj=0.7187 mas, bpa=-18.65 degrees ! Estimated noise=0.292867 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 (-1.6, -0.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0104379 Jy ! Component: 100 - total flux cleaned = 0.0160227 Jy ! Component: 150 - total flux cleaned = 0.0192335 Jy ! Component: 200 - total flux cleaned = 0.0219046 Jy ! Total flux subtracted in 200 components = 0.0219046 Jy ! Clean residual min=-0.006321 max=0.008084 Jy/beam ! Clean residual mean=0.000235 rms=0.001085 Jy/beam ! Combined flux in latest and established models = 1.99627 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 52 components and 1.99627 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162897Jy sigma=1.707835 ! Fit after self-cal, rms=0.162918Jy sigma=1.706923 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-1.6, -12.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00661924 Jy ! Component: 100 - total flux cleaned = 0.00954778 Jy ! Component: 150 - total flux cleaned = 0.0125685 Jy ! Component: 200 - total flux cleaned = 0.0150732 Jy ! Total flux subtracted in 200 components = 0.0150732 Jy ! Clean residual min=-0.005738 max=0.006490 Jy/beam ! Clean residual mean=0.000198 rms=0.000963 Jy/beam ! Combined flux in latest and established models = 2.01134 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 72 components and 2.01134 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162340Jy sigma=1.701180 ! Fit after self-cal, rms=0.162285Jy sigma=1.700668 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-2.2, -9). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0066931 Jy ! Component: 100 - total flux cleaned = 0.00946916 Jy ! Component: 150 - total flux cleaned = 0.011908 Jy ! Component: 200 - total flux cleaned = 0.0141684 Jy ! Total flux subtracted in 200 components = 0.0141684 Jy ! Clean residual min=-0.005005 max=0.004526 Jy/beam ! Clean residual mean=0.000160 rms=0.000835 Jy/beam ! Combined flux in latest and established models = 2.02551 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 89 components and 2.02551 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161827Jy sigma=1.696157 ! Fit after self-cal, rms=0.161788Jy sigma=1.695744 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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. ! ! Fit before self-cal, rms=0.161788Jy sigma=1.695744 ! Fit after self-cal, rms=0.161730Jy sigma=1.694664 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5787 mas, bmaj=0.7186 mas, bpa=-18.69 degrees ! Estimated noise=0.292824 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00282301 Jy ! Component: 100 - total flux cleaned = 0.00516203 Jy ! Component: 150 - total flux cleaned = 0.00721818 Jy ! Component: 200 - total flux cleaned = 0.00905403 Jy ! Total flux subtracted in 200 components = 0.00905403 Jy ! Clean residual min=-0.003895 max=0.004623 Jy/beam ! Clean residual mean=0.000133 rms=0.000763 Jy/beam ! Combined flux in latest and established models = 2.03457 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 94 components and 2.03457 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161497Jy sigma=1.692297 ! Fit after self-cal, rms=0.161484Jy sigma=1.692043 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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 (2.2, -1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00297196 Jy ! Component: 100 - total flux cleaned = 0.0033647 Jy ! Component: 150 - total flux cleaned = 0.00439241 Jy ! Component: 200 - total flux cleaned = 0.00575143 Jy ! Total flux subtracted in 200 components = 0.00575143 Jy ! Clean residual min=-0.003325 max=0.003797 Jy/beam ! Clean residual mean=0.000117 rms=0.000725 Jy/beam ! Combined flux in latest and established models = 2.04032 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 105 components and 2.04032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161268Jy sigma=1.689094 ! Fit after self-cal, rms=0.161255Jy sigma=1.688970 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (20.2, 15.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00260245 Jy ! Component: 100 - total flux cleaned = 0.0041264 Jy ! Component: 150 - total flux cleaned = 0.00514165 Jy ! Component: 200 - total flux cleaned = 0.0061859 Jy ! Total flux subtracted in 200 components = 0.0061859 Jy ! Clean residual min=-0.002966 max=0.003599 Jy/beam ! Clean residual mean=0.000103 rms=0.000687 Jy/beam ! Combined flux in latest and established models = 2.0465 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 118 components and 2.0465 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161110Jy sigma=1.687401 ! Fit after self-cal, rms=0.161105Jy sigma=1.687315 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (3.6, -9.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00291338 Jy ! Component: 100 - total flux cleaned = 0.0041199 Jy ! Component: 150 - total flux cleaned = 0.00512991 Jy ! Component: 200 - total flux cleaned = 0.00592482 Jy ! Total flux subtracted in 200 components = 0.00592482 Jy ! Clean residual min=-0.002959 max=0.003363 Jy/beam ! Clean residual mean=0.000088 rms=0.000659 Jy/beam ! Combined flux in latest and established models = 2.05243 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 134 components and 2.05243 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160981Jy sigma=1.685934 ! Fit after self-cal, rms=0.160971Jy sigma=1.685811 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (15.2, 26). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00238248 Jy ! Component: 100 - total flux cleaned = 0.0032204 Jy ! Component: 150 - total flux cleaned = 0.00419427 Jy ! Component: 200 - total flux cleaned = 0.00497657 Jy ! Total flux subtracted in 200 components = 0.00497657 Jy ! Clean residual min=-0.002777 max=0.003011 Jy/beam ! Clean residual mean=0.000078 rms=0.000631 Jy/beam ! Combined flux in latest and established models = 2.0574 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 144 components and 2.0574 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160871Jy sigma=1.684755 ! Fit after self-cal, rms=0.160866Jy sigma=1.684684 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (8.2, 23.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00173828 Jy ! Component: 100 - total flux cleaned = 0.00219147 Jy ! Component: 150 - total flux cleaned = 0.00314993 Jy ! Component: 200 - total flux cleaned = 0.00366443 Jy ! Total flux subtracted in 200 components = 0.00366443 Jy ! Clean residual min=-0.002701 max=0.002698 Jy/beam ! Clean residual mean=0.000069 rms=0.000613 Jy/beam ! Combined flux in latest and established models = 2.06107 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 152 components and 2.06107 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160789Jy sigma=1.683743 ! Fit after self-cal, rms=0.160787Jy sigma=1.683695 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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. ! ! Fit before self-cal, rms=0.160787Jy sigma=1.683695 ! Fit after self-cal, rms=0.160865Jy sigma=1.683287 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5787 mas, bmaj=0.7187 mas, bpa=-18.65 degrees ! Estimated noise=0.292892 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00131007 Jy ! Component: 100 - total flux cleaned = 0.0023576 Jy ! Component: 150 - total flux cleaned = 0.00332713 Jy ! Component: 200 - total flux cleaned = 0.00400054 Jy ! Total flux subtracted in 200 components = 0.00400054 Jy ! Clean residual min=-0.002625 max=0.002549 Jy/beam ! Clean residual mean=0.000062 rms=0.000596 Jy/beam ! Combined flux in latest and established models = 2.06507 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 159 components and 2.06507 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160803Jy sigma=1.682632 ! Fit after self-cal, rms=0.160802Jy sigma=1.682600 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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.000378629 Jy ! Component: 100 - total flux cleaned = 0.00114055 Jy ! Component: 150 - total flux cleaned = 0.00153213 Jy ! Component: 200 - total flux cleaned = 0.00191188 Jy ! Total flux subtracted in 200 components = 0.00191188 Jy ! Clean residual min=-0.002567 max=0.002482 Jy/beam ! Clean residual mean=0.000057 rms=0.000585 Jy/beam ! Combined flux in latest and established models = 2.06698 Jy ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 166 components and 2.06698 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160758Jy sigma=1.682121 ! Fit after self-cal, rms=0.160757Jy sigma=1.682095 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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. ! ! Fit before self-cal, rms=0.160757Jy sigma=1.682095 ! Fit after self-cal, rms=0.160792Jy sigma=1.682037 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=2.032 mas, bmaj=2.629 mas, bpa=-7.096 degrees ! Estimated noise=0.617787 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00464591 Jy ! Component: 100 - total flux cleaned = 0.00754656 Jy ! Component: 150 - total flux cleaned = 0.00903374 Jy ! Component: 200 - total flux cleaned = 0.0101073 Jy ! Total flux subtracted in 200 components = 0.0101073 Jy ! Clean residual min=-0.003668 max=0.003882 Jy/beam ! Clean residual mean=0.000158 rms=0.001128 Jy/beam ! Combined flux in latest and established models = 2.07709 Jy selfcal ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 175 components and 2.07709 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160761Jy sigma=1.682218 ! Fit after self-cal, rms=0.160754Jy sigma=1.682104 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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.00117792 Jy ! Component: 100 - total flux cleaned = 0.00142323 Jy ! Component: 150 - total flux cleaned = 0.00164738 Jy ! Component: 200 - total flux cleaned = 0.00180281 Jy ! Total flux subtracted in 200 components = 0.00180281 Jy ! Clean residual min=-0.003622 max=0.003694 Jy/beam ! Clean residual mean=0.000134 rms=0.001080 Jy/beam ! Combined flux in latest and established models = 2.07889 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 187 components and 2.07889 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160747Jy sigma=1.682239 ! Fit after self-cal, rms=0.160747Jy sigma=1.682209 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec 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 ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5788 mas, bmaj=0.7187 mas, bpa=-18.66 degrees ! Estimated noise=0.292945 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5788 x 0.7187 at -18.66 degrees (North through East) ! Clean map min=-0.0047949 max=1.0903 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. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160747Jy sigma=1.682209 ! Fit after self-cal, rms=0.154554Jy sigma=1.583704 ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5787 mas, bmaj=0.7186 mas, bpa=-18.72 degrees ! Estimated noise=0.293114 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000197001 Jy ! Component: 100 - total flux cleaned = -6.85712e-05 Jy ! Component: 150 - total flux cleaned = -0.000127151 Jy ! Component: 200 - total flux cleaned = -0.00035212 Jy ! Total flux subtracted in 200 components = -0.00035212 Jy ! Clean residual min=-0.002648 max=0.002272 Jy/beam ! Clean residual mean=0.000027 rms=0.000550 Jy/beam ! Combined flux in latest and established models = 2.07854 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 199 components and 2.07854 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154507Jy sigma=1.582998 ! Fit after self-cal, rms=0.154504Jy sigma=1.582974 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -4.75385e-06 Jy ! Component: 100 - total flux cleaned = -5.04028e-06 Jy ! Component: 150 - total flux cleaned = -0.000158289 Jy ! Component: 200 - total flux cleaned = -0.000303529 Jy ! Total flux subtracted in 200 components = -0.000303529 Jy ! Clean residual min=-0.002634 max=0.002267 Jy/beam ! Clean residual mean=0.000028 rms=0.000544 Jy/beam ! Combined flux in latest and established models = 2.07824 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 203 components and 2.07824 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154471Jy sigma=1.582519 ! Fit after self-cal, rms=0.154468Jy sigma=1.582500 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154468Jy sigma=1.582500 ! Fit after self-cal, rms=0.154504Jy sigma=1.582472 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 ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.4865 mas, bmaj=0.5864 mas, bpa=-31.53 degrees ! Estimated noise=0.516136 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.734372 Jy ! Component: 100 - total flux cleaned = 1.08534 Jy ! Component: 150 - total flux cleaned = 1.31408 Jy ! Component: 200 - total flux cleaned = 1.47479 Jy ! Total flux subtracted in 200 components = 1.47479 Jy ! Clean residual min=-0.044362 max=0.089524 Jy/beam ! Clean residual mean=0.000207 rms=0.005457 Jy/beam ! Combined flux in latest and established models = 1.47479 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 9 model components to the UV plane model. ! The established model now contains 9 components and 1.47479 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5788 mas, bmaj=0.7186 mas, bpa=-18.72 degrees ! Estimated noise=0.293177 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.191906 Jy ! Component: 100 - total flux cleaned = 0.307656 Jy ! Component: 150 - total flux cleaned = 0.381701 Jy ! Component: 200 - total flux cleaned = 0.429233 Jy ! Total flux subtracted in 200 components = 0.429233 Jy ! Clean residual min=-0.009273 max=0.025089 Jy/beam ! Clean residual mean=0.000488 rms=0.002290 Jy/beam ! Combined flux in latest and established models = 1.90402 Jy ! Adding 17 model components to the UV plane model. ! The established model now contains 17 components and 1.90402 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.030585 Jy ! Component: 100 - total flux cleaned = 0.0513205 Jy ! Component: 150 - total flux cleaned = 0.0664793 Jy ! Component: 200 - total flux cleaned = 0.0759741 Jy ! Total flux subtracted in 200 components = 0.0759741 Jy ! Clean residual min=-0.007707 max=0.010333 Jy/beam ! Clean residual mean=0.000279 rms=0.001257 Jy/beam ! Combined flux in latest and established models = 1.98 Jy ! Adding 27 model components to the UV plane model. ! The established model now contains 36 components and 1.98 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-0.2, 0.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00882691 Jy ! Component: 100 - total flux cleaned = 0.0157488 Jy ! Component: 150 - total flux cleaned = 0.0218614 Jy ! Component: 200 - total flux cleaned = 0.0270138 Jy ! Total flux subtracted in 200 components = 0.0270138 Jy ! Clean residual min=-0.004683 max=0.006670 Jy/beam ! Clean residual mean=0.000207 rms=0.000978 Jy/beam ! Combined flux in latest and established models = 2.00701 Jy ! Adding 35 model components to the UV plane model. ! The established model now contains 54 components and 2.00701 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-1, -2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00561917 Jy ! Component: 100 - total flux cleaned = 0.0098813 Jy ! Component: 150 - total flux cleaned = 0.0137715 Jy ! Component: 200 - total flux cleaned = 0.0171429 Jy ! Total flux subtracted in 200 components = 0.0171429 Jy ! Clean residual min=-0.003953 max=0.004873 Jy/beam ! Clean residual mean=0.000163 rms=0.000827 Jy/beam ! Combined flux in latest and established models = 2.02415 Jy ! Adding 44 model components to the UV plane model. ! The established model now contains 70 components and 2.02415 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-1.8, -5.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00465732 Jy ! Component: 100 - total flux cleaned = 0.00778477 Jy ! Component: 150 - total flux cleaned = 0.0107057 Jy ! Component: 200 - total flux cleaned = 0.0136046 Jy ! Total flux subtracted in 200 components = 0.0136046 Jy ! Clean residual min=-0.003661 max=0.004056 Jy/beam ! Clean residual mean=0.000129 rms=0.000728 Jy/beam ! Combined flux in latest and established models = 2.03776 Jy ! Adding 49 model components to the UV plane model. ! The established model now contains 86 components and 2.03776 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (-2.4, -10.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0034699 Jy ! Component: 100 - total flux cleaned = 0.00594168 Jy ! Component: 150 - total flux cleaned = 0.00838425 Jy ! Component: 200 - total flux cleaned = 0.010421 Jy ! Total flux subtracted in 200 components = 0.010421 Jy ! Clean residual min=-0.003302 max=0.003326 Jy/beam ! Clean residual mean=0.000102 rms=0.000656 Jy/beam ! Combined flux in latest and established models = 2.04818 Jy ! Adding 54 model components to the UV plane model. ! The established model now contains 101 components and 2.04818 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (1.2, -11.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00273393 Jy ! Component: 100 - total flux cleaned = 0.00438109 Jy ! Component: 150 - total flux cleaned = 0.00593672 Jy ! Component: 200 - total flux cleaned = 0.00720804 Jy ! Total flux subtracted in 200 components = 0.00720804 Jy ! Clean residual min=-0.002938 max=0.002978 Jy/beam ! Clean residual mean=0.000084 rms=0.000611 Jy/beam ! Combined flux in latest and established models = 2.05539 Jy ! Adding 56 model components to the UV plane model. ! The established model now contains 117 components and 2.05539 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Added new window around map position (0.8, -24.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00239995 Jy ! Component: 100 - total flux cleaned = 0.0035785 Jy ! Component: 150 - total flux cleaned = 0.00471099 Jy ! Component: 200 - total flux cleaned = 0.00579964 Jy ! Total flux subtracted in 200 components = 0.00579964 Jy ! Clean residual min=-0.002659 max=0.002565 Jy/beam ! Clean residual mean=0.000070 rms=0.000579 Jy/beam ! Combined flux in latest and established models = 2.06119 Jy ! Adding 57 model components to the UV plane model. ! The established model now contains 130 components and 2.06119 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=2.031 mas, bmaj=2.629 mas, bpa=-7.139 degrees ! Estimated noise=0.618574 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00632965 Jy ! Component: 100 - total flux cleaned = 0.0100452 Jy ! Component: 150 - total flux cleaned = 0.0128166 Jy ! Component: 200 - total flux cleaned = 0.0150228 Jy ! Total flux subtracted in 200 components = 0.0150228 Jy ! Clean residual min=-0.003241 max=0.003379 Jy/beam ! Clean residual mean=0.000165 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 2.07621 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.00180664 Jy ! Component: 100 - total flux cleaned = 0.00330388 Jy ! Component: 150 - total flux cleaned = 0.00402902 Jy ! Component: 200 - total flux cleaned = 0.0044944 Jy ! Total flux subtracted in 200 components = 0.0044944 Jy ! Clean residual min=-0.003063 max=0.003175 Jy/beam ! Clean residual mean=0.000109 rms=0.000887 Jy/beam ! Combined flux in latest and established models = 2.0807 Jy ! Adding 60 model components to the UV plane model. ! The established model now contains 183 components and 2.0807 Jy ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! Estimated beam: bmin=0.5788 mas, bmaj=0.7186 mas, bpa=-18.72 degrees ! Estimated noise=0.293177 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00121228 Jy ! Component: 100 - total flux cleaned = -0.00147564 Jy ! Component: 150 - total flux cleaned = -0.0017126 Jy ! Component: 200 - total flux cleaned = -0.00193633 Jy ! Total flux subtracted in 200 components = -0.00193633 Jy ! Clean residual min=-0.002387 max=0.002264 Jy/beam ! Clean residual mean=0.000026 rms=0.000527 Jy/beam ! Combined flux in latest and established models = 2.07877 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 199 components and 2.07877 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.154467Jy sigma=1.583984 ! Fit after self-cal, rms=0.154428Jy sigma=1.583216 wmodel J1642+6856_X_map.mod ! Writing 199 model components to file: J1642+6856_X_map.mod wobs J1642+6856_X_uvs.fits ! Writing UV FITS file: J1642+6856_X_uvs.fits wwins J1642+6856_X_map.win ! wwins: Wrote 18 windows to J1642+6856_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000455387 Jy ! Component: 100 - total flux cleaned = -0.000663531 Jy ! Component: 150 - total flux cleaned = -0.000584222 Jy ! Component: 200 - total flux cleaned = -0.000427825 Jy ! Component: 250 - total flux cleaned = -0.000312719 Jy ! Component: 300 - total flux cleaned = -0.000237152 Jy ! Component: 350 - total flux cleaned = -0.000163273 Jy ! Component: 400 - total flux cleaned = -1.7322e-05 Jy ! Component: 450 - total flux cleaned = 0.000126321 Jy ! Component: 500 - total flux cleaned = 0.00019683 Jy ! Component: 550 - total flux cleaned = 0.000337044 Jy ! Component: 600 - total flux cleaned = 0.000371603 Jy ! Component: 650 - total flux cleaned = 0.000440664 Jy ! Component: 700 - total flux cleaned = 0.000441214 Jy ! Component: 750 - total flux cleaned = 0.000610071 Jy ! Component: 800 - total flux cleaned = 0.000609891 Jy ! Component: 850 - total flux cleaned = 0.000775332 Jy ! Component: 900 - total flux cleaned = 0.000808091 Jy ! Component: 950 - total flux cleaned = 0.000775339 Jy ! Component: 1000 - total flux cleaned = 0.000937174 Jy ! Component: 1050 - total flux cleaned = 0.000969041 Jy ! Component: 1100 - total flux cleaned = 0.00100089 Jy ! Component: 1150 - total flux cleaned = 0.00122223 Jy ! Component: 1200 - total flux cleaned = 0.00122238 Jy ! Component: 1250 - total flux cleaned = 0.00140908 Jy ! Component: 1300 - total flux cleaned = 0.00140891 Jy ! Component: 1350 - total flux cleaned = 0.00156306 Jy ! Component: 1400 - total flux cleaned = 0.00174645 Jy ! Component: 1450 - total flux cleaned = 0.00192855 Jy ! Component: 1500 - total flux cleaned = 0.00207965 Jy ! Component: 1550 - total flux cleaned = 0.00210974 Jy ! Component: 1600 - total flux cleaned = 0.00219954 Jy ! Component: 1650 - total flux cleaned = 0.00225889 Jy ! Component: 1700 - total flux cleaned = 0.00228866 Jy ! Component: 1750 - total flux cleaned = 0.00234775 Jy ! Component: 1800 - total flux cleaned = 0.00246502 Jy ! Component: 1850 - total flux cleaned = 0.00266902 Jy ! Component: 1900 - total flux cleaned = 0.00275587 Jy ! Component: 1950 - total flux cleaned = 0.00281354 Jy ! Component: 2000 - total flux cleaned = 0.00284224 Jy ! Component: 2050 - total flux cleaned = 0.00289927 Jy ! Component: 2100 - total flux cleaned = 0.00304147 Jy ! Component: 2150 - total flux cleaned = 0.00312644 Jy ! Component: 2200 - total flux cleaned = 0.00312624 Jy ! Component: 2250 - total flux cleaned = 0.0031541 Jy ! Component: 2300 - total flux cleaned = 0.00332169 Jy ! Component: 2350 - total flux cleaned = 0.00334964 Jy ! Component: 2400 - total flux cleaned = 0.00334945 Jy ! Component: 2450 - total flux cleaned = 0.00354248 Jy ! Component: 2500 - total flux cleaned = 0.00362492 Jy ! Component: 2550 - total flux cleaned = 0.00357003 Jy ! Component: 2600 - total flux cleaned = 0.00378811 Jy ! Component: 2650 - total flux cleaned = 0.00378818 Jy ! Component: 2700 - total flux cleaned = 0.00376134 Jy ! Component: 2750 - total flux cleaned = 0.00394972 Jy ! Component: 2800 - total flux cleaned = 0.00411054 Jy ! Component: 2850 - total flux cleaned = 0.00411055 Jy ! Component: 2900 - total flux cleaned = 0.00416361 Jy ! Component: 2950 - total flux cleaned = 0.00429607 Jy ! Component: 3000 - total flux cleaned = 0.00442806 Jy ! Component: 3050 - total flux cleaned = 0.00442809 Jy ! Component: 3100 - total flux cleaned = 0.00458517 Jy ! Component: 3150 - total flux cleaned = 0.00455895 Jy ! Component: 3200 - total flux cleaned = 0.00463696 Jy ! Component: 3250 - total flux cleaned = 0.00468887 Jy ! Component: 3300 - total flux cleaned = 0.00486975 Jy ! Component: 3350 - total flux cleaned = 0.00494694 Jy ! Component: 3400 - total flux cleaned = 0.00510077 Jy ! Component: 3450 - total flux cleaned = 0.00504976 Jy ! Component: 3500 - total flux cleaned = 0.00522792 Jy ! Component: 3550 - total flux cleaned = 0.00522789 Jy ! Component: 3600 - total flux cleaned = 0.00535429 Jy ! Component: 3650 - total flux cleaned = 0.00555576 Jy ! Component: 3700 - total flux cleaned = 0.00560584 Jy ! Component: 3750 - total flux cleaned = 0.00563084 Jy ! Component: 3800 - total flux cleaned = 0.00560587 Jy ! Component: 3850 - total flux cleaned = 0.00570523 Jy ! Component: 3900 - total flux cleaned = 0.0059036 Jy ! Component: 3950 - total flux cleaned = 0.00597762 Jy ! Component: 4000 - total flux cleaned = 0.00600229 Jy ! Component: 4050 - total flux cleaned = 0.00612498 Jy ! Total flux subtracted in 4096 components = 0.00614946 Jy ! Clean residual min=-0.001355 max=0.001484 Jy/beam ! Clean residual mean=0.000015 rms=0.000399 Jy/beam ! Combined flux in latest and established models = 2.08492 Jy keep ! Adding 399 model components to the UV plane model. ! The established model now contains 596 components and 2.08491 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv102_uvs/J1642+6856_X_map.fits ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5788 x 0.7186 at -18.72 degrees (North through East) ! Clean map min=-0.0019125 max=1.089 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv102_uvs/J1642+6856_X_map.fits wmodel /scr/pima/rv102_uvs/J1642+6856_X_map.mod ! Writing 596 model components to file: /scr/pima/rv102_uvs/J1642+6856_X_map.mod wobs /scr/pima/rv102_uvs/J1642+6856_X_uvs.fits ! Writing UV FITS file: /scr/pima/rv102_uvs/J1642+6856_X_uvs.fits wwins /scr/pima/rv102_uvs/J1642+6856_X_map.win ! wwins: Wrote 19 windows to /scr/pima/rv102_uvs/J1642+6856_X_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 1.08898 invert ! Inverting map ! Your choice of large map pixels excluded 1.32% of the data. ! The x-axis pixel size should ideally be below 0.1818 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000399462 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.110047 0.110047 0.220093 0.440187 0.880374 1.76075 3.52149 7.04299 14.086 28.172 56.3439 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0011 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.088 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2726.11 quit ! Quitting program ! Log file /scr/pima/rv102_uvs/J1642+6856_X_dfm.log closed on Sat Nov 9 10:20:44 2019