! Started logfile: /scr/pima/rv123_uvs/J2225-0457_X_dfm.log on Sun Dec 23 00:45:29 2018 float field_size field_size = 1024 float field_cell field_cell = 0.2000 float taper_size taper_size = 70.0 float freq freq = 8.36575 obs /scr/pima/rv123_uvs/J2225-0457_X_uva.fits ! Reading UV FITS file: /scr/pima/rv123_uvs/J2225-0457_X_uva.fits ! AN table 1: 25 integrations on 15 of 15 possible baselines. ! AN table 2: 27 integrations on 21 of 21 possible baselines. ! AN table 3: 2 integrations on 21 of 21 possible baselines. ! AN table 4: 35 integrations on 45 of 45 possible baselines. ! AN table 5: 2 integrations on 6 of 6 possible baselines. ! Apparent sampling: 0.438351 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2225-0457 ! ! 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.36575e+09 8e+06 1 8e+06 ! 02 2 8.44575e+09 8e+06 1 8e+06 ! 03 3 8.80575e+09 8e+06 1 8e+06 ! 04 4 8.92575e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 89 lines of history. ! ! Reading 4508 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 J2225-0457_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 47 telescope corrections were flagged in sub-array 1. ! A total of 113 telescope corrections were flagged in sub-array 2. ! A total of 192 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 47 telescope corrections were flagged in sub-array 1. ! A total of 113 telescope corrections were flagged in sub-array 2. ! A total of 192 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 47 telescope corrections were flagged in sub-array 1. ! A total of 113 telescope corrections were flagged in sub-array 2. ! A total of 192 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 47 telescope corrections were flagged in sub-array 1. ! A total of 113 telescope corrections were flagged in sub-array 2. ! A total of 192 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.825022Jy sigma=20.266031 ! Fit after self-cal, rms=0.724408Jy sigma=7.539768 ! 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 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.6784 mas, bmaj=0.9044 mas, bpa=6.038 degrees ! Estimated noise=1.39027 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.982941 Jy ! Component: 100 - total flux cleaned = 1.24056 Jy ! Total flux subtracted in 100 components = 1.24056 Jy ! Clean residual min=-0.069822 max=0.238342 Jy/beam ! Clean residual mean=0.000193 rms=0.014044 Jy/beam ! Combined flux in latest and established models = 1.24056 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 1.24056 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.610959Jy sigma=6.452553 ! Fit after self-cal, rms=0.610598Jy sigma=6.451331 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (1.8, -1.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.232443 Jy ! Component: 100 - total flux cleaned = 0.378133 Jy ! Total flux subtracted in 100 components = 0.378133 Jy ! Clean residual min=-0.051013 max=0.102503 Jy/beam ! Clean residual mean=0.000065 rms=0.010256 Jy/beam ! Combined flux in latest and established models = 1.6187 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 14 components and 1.6187 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.457708Jy sigma=4.797130 ! Fit after self-cal, rms=0.447766Jy sigma=4.668365 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (2.6, -0.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.129774 Jy ! Component: 100 - total flux cleaned = 0.225392 Jy ! Total flux subtracted in 100 components = 0.225392 Jy ! Clean residual min=-0.042152 max=0.065657 Jy/beam ! Clean residual mean=-0.000013 rms=0.008174 Jy/beam ! Combined flux in latest and established models = 1.84409 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 26 components and 1.84409 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.388363Jy sigma=4.019789 ! Fit after self-cal, rms=0.385940Jy sigma=3.986348 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0813663 Jy ! Component: 100 - total flux cleaned = 0.144217 Jy ! Total flux subtracted in 100 components = 0.144217 Jy ! Clean residual min=-0.036406 max=0.055406 Jy/beam ! Clean residual mean=-0.000062 rms=0.007271 Jy/beam ! Combined flux in latest and established models = 1.98831 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 33 components and 1.98831 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.365743Jy sigma=3.750238 ! Fit after self-cal, rms=0.364922Jy sigma=3.738501 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-0.8, -2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0580393 Jy ! Component: 100 - total flux cleaned = 0.102312 Jy ! Total flux subtracted in 100 components = 0.102312 Jy ! Clean residual min=-0.031579 max=0.038551 Jy/beam ! Clean residual mean=-0.000097 rms=0.006781 Jy/beam ! Combined flux in latest and established models = 2.09062 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 38 components and 2.09062 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.359056Jy sigma=3.642843 ! Fit after self-cal, rms=0.358634Jy sigma=3.637640 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 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 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7552 mas, bmaj=1.024 mas, bpa=19.36 degrees ! Estimated noise=1.15459 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.0168368 Jy ! Component: 100 - total flux cleaned = 0.0235004 Jy ! Component: 150 - total flux cleaned = 0.0286945 Jy ! Component: 200 - total flux cleaned = 0.031091 Jy ! Total flux subtracted in 200 components = 0.031091 Jy ! Clean residual min=-0.057150 max=0.056727 Jy/beam ! Clean residual mean=-0.002212 rms=0.013626 Jy/beam ! Combined flux in latest and established models = 2.12171 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 48 components and 2.12171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.354973Jy sigma=3.581107 ! Fit after self-cal, rms=0.354769Jy sigma=3.577586 ! Inverting map ! Your choice of large map pixels excluded 0.352% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.7552 x 1.024 at 19.36 degrees (North through East) ! Clean map min=-0.055141 max=1.2795 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.81 MC 0.91 ON 1.01 WN 1.40 ! WZ 0.92 YG 0.99 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.89 HN 1.17 MC 0.88 ON 0.94 ! SC 0.90 WN 1.48 WZ 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.81 HN 2.90 KP 1.10 LA 0.96 ! NL 1.03 ON 1.05 PT 0.81 ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 FD 0.86 HN 1.03 KP 1.05 ! LA 0.97 MK 0.99 NL 1.01 OV 0.98 ! PT 0.84 SC 0.98 ! ! Telescope amplitude corrections in sub-array 5: ! KP 1.05 MK 0.91 OV 0.99 YG 1.13 ! ! ! Correcting IF 2. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.78 MC 0.97 ON 0.94 WN 1.46 ! WZ 0.94 YG 0.84 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.81 HN 1.16 MC 0.92 ON 0.92 ! SC 0.90 WN 1.48 WZ 0.97 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.81 HN 2.98 KP 1.09 LA 0.97 ! NL 1.04 ON 1.02 PT 0.82 ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 FD 0.86 HN 1.03 KP 1.06 ! LA 0.96 MK 0.97 NL 1.00 OV 0.97 ! PT 0.84 SC 0.99 ! ! Telescope amplitude corrections in sub-array 5: ! KP 1.05 MK 0.89 OV 1.01 YG 1.03 ! ! ! Correcting IF 3. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.89 MC 1.42 ON 1.32 WN 1.29 ! WZ 0.97 YG 0.86 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.87 HN 1.12 MC 1.32 ON 1.25 ! SC 0.90 WN 1.37 WZ 1.06 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.81 HN 2.96 KP 1.11 LA 0.97 ! NL 1.04 ON 1.35 PT 0.80 ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.99 FD 0.86 HN 1.02 KP 1.07 ! LA 0.97 MK 0.98 NL 0.99 OV 0.96 ! PT 0.84 SC 0.99 ! ! Telescope amplitude corrections in sub-array 5: ! KP 1.05 MK 0.92 OV 0.98 YG 1.04 ! ! ! Correcting IF 4. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 1.02 MC 1.41 ON 1.28 WN 1.24 ! WZ 1.00 YG 1.03 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! HH 0.95 HN 1.21 MC 1.30 ON 1.22 ! SC 0.92 WN 1.27 WZ 1.06 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.81 HN 2.90 KP 1.11 LA 0.96 ! NL 1.04 ON 1.42 PT 0.81 ! ! A total of 36 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.98 FD 0.86 HN 1.02 KP 1.05 ! LA 0.96 MK 0.98 NL 1.01 OV 0.96 ! PT 0.84 SC 1.00 ! ! Telescope amplitude corrections in sub-array 5: ! KP 1.08 MK 0.92 OV 0.95 YG 1.13 ! ! ! Fit before self-cal, rms=0.354769Jy sigma=3.577586 ! Fit after self-cal, rms=0.226342Jy sigma=1.640056 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7724 mas, bmaj=1.049 mas, bpa=14.27 degrees ! Estimated noise=1.18364 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.0307727 Jy ! Component: 100 - total flux cleaned = 0.03507 Jy ! Component: 150 - total flux cleaned = 0.0356706 Jy ! Component: 200 - total flux cleaned = 0.035274 Jy ! Total flux subtracted in 200 components = 0.035274 Jy ! Clean residual min=-0.011279 max=0.026617 Jy/beam ! Clean residual mean=0.005235 rms=0.003821 Jy/beam ! Combined flux in latest and established models = 2.15698 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 62 components and 2.15698 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220833Jy sigma=1.521395 ! Fit after self-cal, rms=0.220088Jy sigma=1.503441 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-1.6, 9.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0261011 Jy ! Component: 100 - total flux cleaned = 0.0407208 Jy ! Component: 150 - total flux cleaned = 0.0503752 Jy ! Component: 200 - total flux cleaned = 0.0542072 Jy ! Total flux subtracted in 200 components = 0.0542072 Jy ! Clean residual min=-0.009524 max=0.022466 Jy/beam ! Clean residual mean=0.004666 rms=0.003453 Jy/beam ! Combined flux in latest and established models = 2.21119 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 82 components and 2.21119 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.214224Jy sigma=1.460136 ! Fit after self-cal, rms=0.213840Jy sigma=1.454202 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-0.2, -7.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0211158 Jy ! Component: 100 - total flux cleaned = 0.032529 Jy ! Component: 150 - total flux cleaned = 0.0383278 Jy ! Component: 200 - total flux cleaned = 0.0429326 Jy ! Total flux subtracted in 200 components = 0.0429326 Jy ! Clean residual min=-0.009251 max=0.019137 Jy/beam ! Clean residual mean=0.004218 rms=0.003187 Jy/beam ! Combined flux in latest and established models = 2.25412 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 100 components and 2.25412 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.209272Jy sigma=1.431193 ! Fit after self-cal, rms=0.209091Jy sigma=1.426489 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-0.4, 34.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0188268 Jy ! Component: 100 - total flux cleaned = 0.0290184 Jy ! Component: 150 - total flux cleaned = 0.0372302 Jy ! Component: 200 - total flux cleaned = 0.0443097 Jy ! Total flux subtracted in 200 components = 0.0443097 Jy ! Clean residual min=-0.009578 max=0.017266 Jy/beam ! Clean residual mean=0.003711 rms=0.003019 Jy/beam ! Combined flux in latest and established models = 2.29843 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 123 components and 2.29843 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.205546Jy sigma=1.403832 ! Fit after self-cal, rms=0.205405Jy sigma=1.401122 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (50, -18.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0172003 Jy ! Component: 100 - total flux cleaned = 0.0253558 Jy ! Component: 150 - total flux cleaned = 0.0299431 Jy ! Component: 200 - total flux cleaned = 0.033737 Jy ! Total flux subtracted in 200 components = 0.033737 Jy ! Clean residual min=-0.009940 max=0.015514 Jy/beam ! Clean residual mean=0.003296 rms=0.002846 Jy/beam ! Combined flux in latest and established models = 2.33217 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 152 components and 2.33217 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.202861Jy sigma=1.379667 ! Fit after self-cal, rms=0.202625Jy sigma=1.376888 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (9.8, -0.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0152753 Jy ! Component: 100 - total flux cleaned = 0.0234535 Jy ! Component: 150 - total flux cleaned = 0.0288943 Jy ! Component: 200 - total flux cleaned = 0.0316923 Jy ! Total flux subtracted in 200 components = 0.0316923 Jy ! Clean residual min=-0.009093 max=0.013515 Jy/beam ! Clean residual mean=0.002971 rms=0.002740 Jy/beam ! Combined flux in latest and established models = 2.36386 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 175 components and 2.36386 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.200112Jy sigma=1.364735 ! Fit after self-cal, rms=0.200069Jy sigma=1.362881 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (7.4, -45.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0130556 Jy ! Component: 100 - total flux cleaned = 0.0196523 Jy ! Component: 150 - total flux cleaned = 0.0223754 Jy ! Component: 200 - total flux cleaned = 0.024466 Jy ! Total flux subtracted in 200 components = 0.024466 Jy ! Clean residual min=-0.008886 max=0.013899 Jy/beam ! Clean residual mean=0.002671 rms=0.002744 Jy/beam ! Combined flux in latest and established models = 2.38833 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 196 components and 2.38833 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.198535Jy sigma=1.350992 ! Fit after self-cal, rms=0.198434Jy sigma=1.349206 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.198434Jy sigma=1.349206 ! Fit after self-cal, rms=0.190896Jy sigma=1.173246 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7747 mas, bmaj=1.051 mas, bpa=14.08 degrees ! Estimated noise=1.1856 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00677928 Jy ! Component: 100 - total flux cleaned = 0.0110621 Jy ! Component: 150 - total flux cleaned = 0.0134237 Jy ! Component: 200 - total flux cleaned = 0.0152163 Jy ! Total flux subtracted in 200 components = 0.0152163 Jy ! Clean residual min=-0.008440 max=0.011788 Jy/beam ! Clean residual mean=0.002484 rms=0.001897 Jy/beam ! Combined flux in latest and established models = 2.40354 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 212 components and 2.40354 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190144Jy sigma=1.166857 ! Fit after self-cal, rms=0.190162Jy sigma=1.166340 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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 (0.4, -18.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0117151 Jy ! Component: 100 - total flux cleaned = 0.0178219 Jy ! Component: 150 - total flux cleaned = 0.0212631 Jy ! Component: 200 - total flux cleaned = 0.0241694 Jy ! Total flux subtracted in 200 components = 0.0241694 Jy ! Clean residual min=-0.007974 max=0.011067 Jy/beam ! Clean residual mean=0.002223 rms=0.001836 Jy/beam ! Combined flux in latest and established models = 2.42771 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 232 components and 2.42771 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.188714Jy sigma=1.155357 ! Fit after self-cal, rms=0.188577Jy sigma=1.153620 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (2.2, -36.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0100329 Jy ! Component: 100 - total flux cleaned = 0.0158156 Jy ! Component: 150 - total flux cleaned = 0.0204443 Jy ! Component: 200 - total flux cleaned = 0.0227303 Jy ! Total flux subtracted in 200 components = 0.0227303 Jy ! Clean residual min=-0.007860 max=0.011175 Jy/beam ! Clean residual mean=0.001963 rms=0.001767 Jy/beam ! Combined flux in latest and established models = 2.45044 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 247 components and 2.45044 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187323Jy sigma=1.146236 ! Fit after self-cal, rms=0.187224Jy sigma=1.145250 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (2.8, -7.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0100414 Jy ! Component: 100 - total flux cleaned = 0.0138175 Jy ! Component: 150 - total flux cleaned = 0.0171151 Jy ! Component: 200 - total flux cleaned = 0.0200871 Jy ! Total flux subtracted in 200 components = 0.0200871 Jy ! Clean residual min=-0.007402 max=0.009550 Jy/beam ! Clean residual mean=0.001755 rms=0.001725 Jy/beam ! Combined flux in latest and established models = 2.47053 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 256 components and 2.47053 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.186194Jy sigma=1.134163 ! Fit after self-cal, rms=0.186151Jy sigma=1.133507 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-3.4, -5.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00773478 Jy ! Component: 100 - total flux cleaned = 0.00923816 Jy ! Component: 150 - total flux cleaned = 0.0109854 Jy ! Component: 200 - total flux cleaned = 0.0124384 Jy ! Total flux subtracted in 200 components = 0.0124384 Jy ! Clean residual min=-0.006924 max=0.008586 Jy/beam ! Clean residual mean=0.001621 rms=0.001694 Jy/beam ! Combined flux in latest and established models = 2.48297 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 268 components and 2.48297 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185354Jy sigma=1.126473 ! Fit after self-cal, rms=0.185331Jy sigma=1.125928 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-8.2, -10.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00833966 Jy ! Component: 100 - total flux cleaned = 0.0117825 Jy ! Component: 150 - total flux cleaned = 0.0128726 Jy ! Component: 200 - total flux cleaned = 0.0138966 Jy ! Total flux subtracted in 200 components = 0.0138966 Jy ! Clean residual min=-0.006392 max=0.008361 Jy/beam ! Clean residual mean=0.001476 rms=0.001652 Jy/beam ! Combined flux in latest and established models = 2.49687 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 284 components and 2.49687 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.184643Jy sigma=1.118966 ! Fit after self-cal, rms=0.184589Jy sigma=1.118148 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-12.2, 26.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00755821 Jy ! Component: 100 - total flux cleaned = 0.0117488 Jy ! Component: 150 - total flux cleaned = 0.0137489 Jy ! Component: 200 - total flux cleaned = 0.015022 Jy ! Total flux subtracted in 200 components = 0.015022 Jy ! Clean residual min=-0.006330 max=0.007568 Jy/beam ! Clean residual mean=0.001308 rms=0.001588 Jy/beam ! Combined flux in latest and established models = 2.51189 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 304 components and 2.51189 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183990Jy sigma=1.112054 ! Fit after self-cal, rms=0.183987Jy sigma=1.111427 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-9.6, 22.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00688541 Jy ! Component: 100 - total flux cleaned = 0.00857412 Jy ! Component: 150 - total flux cleaned = 0.00953944 Jy ! Component: 200 - total flux cleaned = 0.0101903 Jy ! Total flux subtracted in 200 components = 0.0101903 Jy ! Clean residual min=-0.006367 max=0.007427 Jy/beam ! Clean residual mean=0.001194 rms=0.001553 Jy/beam ! Combined flux in latest and established models = 2.52208 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 320 components and 2.52208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183593Jy sigma=1.105433 ! Fit after self-cal, rms=0.183529Jy sigma=1.104897 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-9.4, -0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00668605 Jy ! Component: 100 - total flux cleaned = 0.0102037 Jy ! Component: 150 - total flux cleaned = 0.0113355 Jy ! Component: 200 - total flux cleaned = 0.0119975 Jy ! Total flux subtracted in 200 components = 0.0119975 Jy ! Clean residual min=-0.006505 max=0.007473 Jy/beam ! Clean residual mean=0.001060 rms=0.001514 Jy/beam ! Combined flux in latest and established models = 2.53407 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 336 components and 2.53408 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183034Jy sigma=1.099911 ! Fit after self-cal, rms=0.182987Jy sigma=1.099495 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-17.8, -14.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00658245 Jy ! Component: 100 - total flux cleaned = 0.00940043 Jy ! Component: 150 - total flux cleaned = 0.0113018 Jy ! Component: 200 - total flux cleaned = 0.0121918 Jy ! Total flux subtracted in 200 components = 0.0121918 Jy ! Clean residual min=-0.006464 max=0.006848 Jy/beam ! Clean residual mean=0.000931 rms=0.001476 Jy/beam ! Combined flux in latest and established models = 2.54627 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 350 components and 2.54627 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182549Jy sigma=1.094571 ! Fit after self-cal, rms=0.182528Jy sigma=1.094207 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (5.4, -0.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00648763 Jy ! Component: 100 - total flux cleaned = 0.0068067 Jy ! Component: 150 - total flux cleaned = 0.00653666 Jy ! Component: 200 - total flux cleaned = 0.00653186 Jy ! Total flux subtracted in 200 components = 0.00653186 Jy ! Clean residual min=-0.005696 max=0.006350 Jy/beam ! Clean residual mean=0.000859 rms=0.001452 Jy/beam ! Combined flux in latest and established models = 2.5528 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 365 components and 2.5528 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.182191Jy sigma=1.088323 ! Fit after self-cal, rms=0.182179Jy sigma=1.087856 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (1.4, 48.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00605528 Jy ! Component: 100 - total flux cleaned = 0.00801363 Jy ! Component: 150 - total flux cleaned = 0.0083924 Jy ! Component: 200 - total flux cleaned = 0.00850758 Jy ! Total flux subtracted in 200 components = 0.00850758 Jy ! Clean residual min=-0.005144 max=0.005930 Jy/beam ! Clean residual mean=0.000751 rms=0.001430 Jy/beam ! Combined flux in latest and established models = 2.56131 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 381 components and 2.56131 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181764Jy sigma=1.083722 ! Fit after self-cal, rms=0.181736Jy sigma=1.083329 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.181736Jy sigma=1.083329 ! Fit after self-cal, rms=0.178928Jy sigma=1.047670 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7751 mas, bmaj=1.051 mas, bpa=14.2 degrees ! Estimated noise=1.18569 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00465345 Jy ! Component: 100 - total flux cleaned = 0.00772842 Jy ! Component: 150 - total flux cleaned = 0.00898749 Jy ! Component: 200 - total flux cleaned = 0.0101719 Jy ! Total flux subtracted in 200 components = 0.0101719 Jy ! Clean residual min=-0.004679 max=0.006100 Jy/beam ! Clean residual mean=0.001081 rms=0.001294 Jy/beam ! Combined flux in latest and established models = 2.57148 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 401 components and 2.57148 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.178573Jy sigma=1.043988 ! Fit after self-cal, rms=0.178576Jy sigma=1.043786 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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 (40.4, -25). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00568418 Jy ! Component: 100 - total flux cleaned = 0.00765802 Jy ! Component: 150 - total flux cleaned = 0.00893806 Jy ! Component: 200 - total flux cleaned = 0.0102715 Jy ! Total flux subtracted in 200 components = 0.0102715 Jy ! Clean residual min=-0.005007 max=0.005908 Jy/beam ! Clean residual mean=0.000965 rms=0.001264 Jy/beam ! Combined flux in latest and established models = 2.58175 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 418 components and 2.58175 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.178145Jy sigma=1.039695 ! Fit after self-cal, rms=0.178127Jy sigma=1.039371 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (32.4, 20.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00530785 Jy ! Component: 100 - total flux cleaned = 0.00720137 Jy ! Component: 150 - total flux cleaned = 0.00826503 Jy ! Component: 200 - total flux cleaned = 0.00916522 Jy ! Total flux subtracted in 200 components = 0.00916522 Jy ! Clean residual min=-0.005011 max=0.005773 Jy/beam ! Clean residual mean=0.000868 rms=0.001236 Jy/beam ! Combined flux in latest and established models = 2.59092 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 432 components and 2.59092 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177781Jy sigma=1.036020 ! Fit after self-cal, rms=0.177774Jy sigma=1.035646 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (1.2, 22.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00525618 Jy ! Component: 100 - total flux cleaned = 0.00809033 Jy ! Component: 150 - total flux cleaned = 0.00991077 Jy ! Component: 200 - total flux cleaned = 0.0114052 Jy ! Total flux subtracted in 200 components = 0.0114052 Jy ! Clean residual min=-0.004537 max=0.005376 Jy/beam ! Clean residual mean=0.000743 rms=0.001213 Jy/beam ! Combined flux in latest and established models = 2.60232 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 446 components and 2.60232 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177422Jy sigma=1.032475 ! Fit after self-cal, rms=0.177415Jy sigma=1.032241 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (9.2, 6.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.005341 Jy ! Component: 100 - total flux cleaned = 0.00798459 Jy ! Component: 150 - total flux cleaned = 0.00985917 Jy ! Component: 200 - total flux cleaned = 0.0110494 Jy ! Total flux subtracted in 200 components = 0.0110494 Jy ! Clean residual min=-0.004986 max=0.005177 Jy/beam ! Clean residual mean=0.000625 rms=0.001173 Jy/beam ! Combined flux in latest and established models = 2.61337 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 462 components and 2.61337 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177097Jy sigma=1.029129 ! Fit after self-cal, rms=0.177093Jy sigma=1.028751 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (29.6, 22.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00458531 Jy ! Component: 100 - total flux cleaned = 0.00621722 Jy ! Component: 150 - total flux cleaned = 0.00672909 Jy ! Component: 200 - total flux cleaned = 0.00730379 Jy ! Total flux subtracted in 200 components = 0.00730379 Jy ! Clean residual min=-0.004889 max=0.005130 Jy/beam ! Clean residual mean=0.000543 rms=0.001151 Jy/beam ! Combined flux in latest and established models = 2.62067 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 477 components and 2.62068 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176851Jy sigma=1.026073 ! Fit after self-cal, rms=0.176833Jy sigma=1.025777 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-4.2, -15.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00469859 Jy ! Component: 100 - total flux cleaned = 0.00702822 Jy ! Component: 150 - total flux cleaned = 0.00794075 Jy ! Component: 200 - total flux cleaned = 0.00851621 Jy ! Total flux subtracted in 200 components = 0.00851621 Jy ! Clean residual min=-0.004996 max=0.004728 Jy/beam ! Clean residual mean=0.000455 rms=0.001126 Jy/beam ! Combined flux in latest and established models = 2.62919 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 497 components and 2.62919 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176599Jy sigma=1.022589 ! Fit after self-cal, rms=0.176583Jy sigma=1.022422 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.176583Jy sigma=1.022422 ! Fit after self-cal, rms=0.176453Jy sigma=1.018995 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=2.417 mas, bmaj=3.634 mas, bpa=-2.751 degrees ! Estimated noise=1.5733 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00550019 Jy ! Component: 100 - total flux cleaned = 0.010161 Jy ! Component: 150 - total flux cleaned = 0.0142708 Jy ! Component: 200 - total flux cleaned = 0.0179948 Jy ! Total flux subtracted in 200 components = 0.0179948 Jy ! Clean residual min=-0.004478 max=0.006211 Jy/beam ! Clean residual mean=0.001326 rms=0.001559 Jy/beam ! Combined flux in latest and established models = 2.64719 Jy selfcal ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 524 components and 2.64719 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176158Jy sigma=1.017600 ! Fit after self-cal, rms=0.176139Jy sigma=1.017394 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.00372982 Jy ! Component: 100 - total flux cleaned = 0.00693153 Jy ! Component: 150 - total flux cleaned = 0.00954373 Jy ! Component: 200 - total flux cleaned = 0.0114236 Jy ! Total flux subtracted in 200 components = 0.0114236 Jy ! Clean residual min=-0.004436 max=0.005960 Jy/beam ! Clean residual mean=0.001063 rms=0.001524 Jy/beam ! Combined flux in latest and established models = 2.65861 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 540 components and 2.65861 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176002Jy sigma=1.016958 ! Fit after self-cal, rms=0.176001Jy sigma=1.016889 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.775 mas, bmaj=1.05 mas, bpa=14.29 degrees ! Estimated noise=1.18552 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.775 x 1.05 at 14.29 degrees (North through East) ! Clean map min=-0.0046877 max=1.2758 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.176001Jy sigma=1.016889 ! Fit after self-cal, rms=0.831807Jy sigma=0.759341 ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7754 mas, bmaj=1.051 mas, bpa=14.18 degrees ! Estimated noise=1.18567 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00219078 Jy ! Component: 100 - total flux cleaned = 0.00384949 Jy ! Component: 150 - total flux cleaned = 0.00536741 Jy ! Component: 200 - total flux cleaned = 0.00632398 Jy ! Total flux subtracted in 200 components = 0.00632398 Jy ! Clean residual min=-0.004102 max=0.004703 Jy/beam ! Clean residual mean=0.000557 rms=0.001051 Jy/beam ! Combined flux in latest and established models = 2.66493 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 551 components and 2.66493 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.831768Jy sigma=0.756495 ! Fit after self-cal, rms=0.831742Jy sigma=0.756410 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-36.8, 32.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00349889 Jy ! Component: 100 - total flux cleaned = 0.00419866 Jy ! Component: 150 - total flux cleaned = 0.00484333 Jy ! Component: 200 - total flux cleaned = 0.00519609 Jy ! Total flux subtracted in 200 components = 0.00519609 Jy ! Clean residual min=-0.004136 max=0.004581 Jy/beam ! Clean residual mean=0.000497 rms=0.001040 Jy/beam ! Combined flux in latest and established models = 2.67013 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 566 components and 2.67013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.831770Jy sigma=0.753057 ! Fit after self-cal, rms=0.831829Jy sigma=0.752915 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-27.2, -9.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00451742 Jy ! Component: 100 - total flux cleaned = 0.00579105 Jy ! Component: 150 - total flux cleaned = 0.00569945 Jy ! Component: 200 - total flux cleaned = 0.00595594 Jy ! Total flux subtracted in 200 components = 0.00595594 Jy ! Clean residual min=-0.004448 max=0.004288 Jy/beam ! Clean residual mean=0.000440 rms=0.001023 Jy/beam ! Combined flux in latest and established models = 2.67609 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 581 components and 2.67609 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.831869Jy sigma=0.749151 ! Fit after self-cal, rms=0.831828Jy sigma=0.748847 ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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.831828Jy sigma=0.748847 ! Fit after self-cal, rms=0.862120Jy sigma=0.747815 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 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.688 mas, bmaj=0.9433 mas, bpa=2.81 degrees ! Estimated noise=1.38715 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.947201 Jy ! Component: 100 - total flux cleaned = 1.25033 Jy ! Component: 150 - total flux cleaned = 1.46185 Jy ! Component: 200 - total flux cleaned = 1.61813 Jy ! Total flux subtracted in 200 components = 1.61813 Jy ! Clean residual min=-0.030248 max=0.089328 Jy/beam ! Clean residual mean=0.000253 rms=0.006800 Jy/beam ! Combined flux in latest and established models = 1.61813 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 21 model components to the UV plane model. ! The established model now contains 21 components and 1.61813 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7753 mas, bmaj=1.05 mas, bpa=14.2 degrees ! Estimated noise=1.18562 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.198469 Jy ! Component: 100 - total flux cleaned = 0.320183 Jy ! Component: 150 - total flux cleaned = 0.397974 Jy ! Component: 200 - total flux cleaned = 0.448921 Jy ! Total flux subtracted in 200 components = 0.448921 Jy ! Clean residual min=-0.006629 max=0.031961 Jy/beam ! Clean residual mean=0.007070 rms=0.003706 Jy/beam ! Combined flux in latest and established models = 2.06706 Jy ! Adding 38 model components to the UV plane model. ! The established model now contains 45 components and 2.06706 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (1.4, 0). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0374884 Jy ! Component: 100 - total flux cleaned = 0.0679758 Jy ! Component: 150 - total flux cleaned = 0.0942379 Jy ! Component: 200 - total flux cleaned = 0.117828 Jy ! Total flux subtracted in 200 components = 0.117828 Jy ! Clean residual min=-0.004475 max=0.017958 Jy/beam ! Clean residual mean=0.005776 rms=0.002571 Jy/beam ! Combined flux in latest and established models = 2.18488 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 91 components and 2.18488 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (3.6, -1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0219883 Jy ! Component: 100 - total flux cleaned = 0.0418764 Jy ! Component: 150 - total flux cleaned = 0.0603847 Jy ! Component: 200 - total flux cleaned = 0.0777273 Jy ! Total flux subtracted in 200 components = 0.0777273 Jy ! Clean residual min=-0.003981 max=0.014066 Jy/beam ! Clean residual mean=0.004956 rms=0.002094 Jy/beam ! Combined flux in latest and established models = 2.26261 Jy ! Adding 74 model components to the UV plane model. ! The established model now contains 133 components and 2.26261 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (43, -7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0170411 Jy ! Component: 100 - total flux cleaned = 0.0325236 Jy ! Component: 150 - total flux cleaned = 0.0471862 Jy ! Component: 200 - total flux cleaned = 0.0610974 Jy ! Total flux subtracted in 200 components = 0.0610974 Jy ! Clean residual min=-0.003950 max=0.011590 Jy/beam ! Clean residual mean=0.004306 rms=0.001781 Jy/beam ! Combined flux in latest and established models = 2.32371 Jy ! Adding 86 model components to the UV plane model. ! The established model now contains 176 components and 2.32371 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (1.2, 51). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0136544 Jy ! Component: 100 - total flux cleaned = 0.0263156 Jy ! Component: 150 - total flux cleaned = 0.038415 Jy ! Component: 200 - total flux cleaned = 0.0499886 Jy ! Total flux subtracted in 200 components = 0.0499886 Jy ! Clean residual min=-0.003805 max=0.009865 Jy/beam ! Clean residual mean=0.003776 rms=0.001577 Jy/beam ! Combined flux in latest and established models = 2.3737 Jy ! Adding 100 model components to the UV plane model. ! The established model now contains 219 components and 2.3737 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-6.8, -10.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0116045 Jy ! Component: 100 - total flux cleaned = 0.0223314 Jy ! Component: 150 - total flux cleaned = 0.0326226 Jy ! Component: 200 - total flux cleaned = 0.0425222 Jy ! Total flux subtracted in 200 components = 0.0425222 Jy ! Clean residual min=-0.003753 max=0.008762 Jy/beam ! Clean residual mean=0.003332 rms=0.001440 Jy/beam ! Combined flux in latest and established models = 2.41622 Jy ! Adding 109 model components to the UV plane model. ! The established model now contains 255 components and 2.41622 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (11.6, 25.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00996658 Jy ! Component: 100 - total flux cleaned = 0.0191684 Jy ! Component: 150 - total flux cleaned = 0.0280294 Jy ! Component: 200 - total flux cleaned = 0.0365788 Jy ! Total flux subtracted in 200 components = 0.0365788 Jy ! Clean residual min=-0.003733 max=0.007779 Jy/beam ! Clean residual mean=0.002943 rms=0.001329 Jy/beam ! Combined flux in latest and established models = 2.4528 Jy ! Adding 112 model components to the UV plane model. ! The established model now contains 281 components and 2.4528 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-24.2, -50.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00876872 Jy ! Component: 100 - total flux cleaned = 0.0167556 Jy ! Component: 150 - total flux cleaned = 0.0244574 Jy ! Component: 200 - total flux cleaned = 0.0319067 Jy ! Total flux subtracted in 200 components = 0.0319067 Jy ! Clean residual min=-0.003681 max=0.007136 Jy/beam ! Clean residual mean=0.002600 rms=0.001245 Jy/beam ! Combined flux in latest and established models = 2.4847 Jy ! Adding 117 model components to the UV plane model. ! The established model now contains 311 components and 2.4847 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (50.6, 50.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00776348 Jy ! Component: 100 - total flux cleaned = 0.0147481 Jy ! Component: 150 - total flux cleaned = 0.02149 Jy ! Component: 200 - total flux cleaned = 0.0280146 Jy ! Total flux subtracted in 200 components = 0.0280146 Jy ! Clean residual min=-0.003798 max=0.006697 Jy/beam ! Clean residual mean=0.002299 rms=0.001179 Jy/beam ! Combined flux in latest and established models = 2.51272 Jy ! Adding 116 model components to the UV plane model. ! The established model now contains 340 components and 2.51272 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-11.4, -21.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00714072 Jy ! Component: 100 - total flux cleaned = 0.0133351 Jy ! Component: 150 - total flux cleaned = 0.0192991 Jy ! Component: 200 - total flux cleaned = 0.0250744 Jy ! Total flux subtracted in 200 components = 0.0250744 Jy ! Clean residual min=-0.003786 max=0.006217 Jy/beam ! Clean residual mean=0.002042 rms=0.001118 Jy/beam ! Combined flux in latest and established models = 2.53779 Jy ! Adding 123 model components to the UV plane model. ! The established model now contains 364 components and 2.53779 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-42.8, -21.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00641942 Jy ! Component: 100 - total flux cleaned = 0.0118791 Jy ! Component: 150 - total flux cleaned = 0.0171524 Jy ! Component: 200 - total flux cleaned = 0.022264 Jy ! Total flux subtracted in 200 components = 0.022264 Jy ! Clean residual min=-0.003764 max=0.005777 Jy/beam ! Clean residual mean=0.001804 rms=0.001076 Jy/beam ! Combined flux in latest and established models = 2.56006 Jy ! Adding 117 model components to the UV plane model. ! The established model now contains 384 components and 2.56006 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (36.8, -23.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00595279 Jy ! Component: 100 - total flux cleaned = 0.0108718 Jy ! Component: 150 - total flux cleaned = 0.0156116 Jy ! Component: 200 - total flux cleaned = 0.0201979 Jy ! Total flux subtracted in 200 components = 0.0201979 Jy ! Clean residual min=-0.003817 max=0.005498 Jy/beam ! Clean residual mean=0.001589 rms=0.001039 Jy/beam ! Combined flux in latest and established models = 2.58026 Jy ! Adding 110 model components to the UV plane model. ! The established model now contains 409 components and 2.58026 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (23.4, 48.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00536631 Jy ! Component: 100 - total flux cleaned = 0.00975523 Jy ! Component: 150 - total flux cleaned = 0.0139861 Jy ! Component: 200 - total flux cleaned = 0.018084 Jy ! Total flux subtracted in 200 components = 0.018084 Jy ! Clean residual min=-0.003768 max=0.005435 Jy/beam ! Clean residual mean=0.001393 rms=0.001013 Jy/beam ! Combined flux in latest and established models = 2.59834 Jy ! Adding 111 model components to the UV plane model. ! The established model now contains 426 components and 2.59834 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (27.2, -50.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00498047 Jy ! Component: 100 - total flux cleaned = 0.00893709 Jy ! Component: 150 - total flux cleaned = 0.0127498 Jy ! Component: 200 - total flux cleaned = 0.0161518 Jy ! Total flux subtracted in 200 components = 0.0161518 Jy ! Clean residual min=-0.003763 max=0.005091 Jy/beam ! Clean residual mean=0.001222 rms=0.000996 Jy/beam ! Combined flux in latest and established models = 2.61449 Jy ! Adding 107 model components to the UV plane model. ! The established model now contains 446 components and 2.61449 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-10.2, 32). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00469375 Jy ! Component: 100 - total flux cleaned = 0.00831647 Jy ! Component: 150 - total flux cleaned = 0.0112444 Jy ! Component: 200 - total flux cleaned = 0.0136822 Jy ! Total flux subtracted in 200 components = 0.0136822 Jy ! Clean residual min=-0.003757 max=0.004765 Jy/beam ! Clean residual mean=0.001081 rms=0.000977 Jy/beam ! Combined flux in latest and established models = 2.62817 Jy ! Adding 88 model components to the UV plane model. ! The established model now contains 456 components and 2.62817 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-23.4, 32.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00395655 Jy ! Component: 100 - total flux cleaned = 0.00648827 Jy ! Component: 150 - total flux cleaned = 0.00879145 Jy ! Component: 200 - total flux cleaned = 0.0110319 Jy ! Total flux subtracted in 200 components = 0.0110319 Jy ! Clean residual min=-0.003822 max=0.004667 Jy/beam ! Clean residual mean=0.000963 rms=0.000955 Jy/beam ! Combined flux in latest and established models = 2.63921 Jy ! Adding 82 model components to the UV plane model. ! The established model now contains 472 components and 2.6392 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (0.2, 26). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00430978 Jy ! Component: 100 - total flux cleaned = 0.00666938 Jy ! Component: 150 - total flux cleaned = 0.00858054 Jy ! Component: 200 - total flux cleaned = 0.0106681 Jy ! Total flux subtracted in 200 components = 0.0106681 Jy ! Clean residual min=-0.003723 max=0.004613 Jy/beam ! Clean residual mean=0.000848 rms=0.000945 Jy/beam ! Combined flux in latest and established models = 2.64987 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 491 components and 2.64987 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (40.6, 0.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00334724 Jy ! Component: 100 - total flux cleaned = 0.00527278 Jy ! Component: 150 - total flux cleaned = 0.00702156 Jy ! Component: 200 - total flux cleaned = 0.00872792 Jy ! Total flux subtracted in 200 components = 0.00872792 Jy ! Clean residual min=-0.003668 max=0.004377 Jy/beam ! Clean residual mean=0.000753 rms=0.000934 Jy/beam ! Combined flux in latest and established models = 2.6586 Jy ! Adding 76 model components to the UV plane model. ! The established model now contains 502 components and 2.6586 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-37, 7.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00429348 Jy ! Component: 100 - total flux cleaned = 0.00646433 Jy ! Component: 150 - total flux cleaned = 0.00811193 Jy ! Component: 200 - total flux cleaned = 0.00921372 Jy ! Total flux subtracted in 200 components = 0.00921372 Jy ! Clean residual min=-0.003553 max=0.004182 Jy/beam ! Clean residual mean=0.000657 rms=0.000915 Jy/beam ! Combined flux in latest and established models = 2.66781 Jy ! Adding 67 model components to the UV plane model. ! The established model now contains 521 components and 2.66781 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-38.4, -16). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00353203 Jy ! Component: 100 - total flux cleaned = 0.00510766 Jy ! Component: 150 - total flux cleaned = 0.0064477 Jy ! Component: 200 - total flux cleaned = 0.00748125 Jy ! Total flux subtracted in 200 components = 0.00748125 Jy ! Clean residual min=-0.003449 max=0.004094 Jy/beam ! Clean residual mean=0.000579 rms=0.000907 Jy/beam ! Combined flux in latest and established models = 2.67529 Jy ! Adding 83 model components to the UV plane model. ! The established model now contains 539 components and 2.6753 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (-22.2, 3.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00238607 Jy ! Component: 100 - total flux cleaned = 0.00305763 Jy ! Component: 150 - total flux cleaned = 0.00407679 Jy ! Component: 200 - total flux cleaned = 0.00497838 Jy ! Total flux subtracted in 200 components = 0.00497838 Jy ! Clean residual min=-0.003453 max=0.003919 Jy/beam ! Clean residual mean=0.000528 rms=0.000898 Jy/beam ! Combined flux in latest and established models = 2.68027 Jy ! Adding 73 model components to the UV plane model. ! The established model now contains 551 components and 2.68027 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (22.2, -51.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0032401 Jy ! Component: 100 - total flux cleaned = 0.00434376 Jy ! Component: 150 - total flux cleaned = 0.00531516 Jy ! Component: 200 - total flux cleaned = 0.00617724 Jy ! Total flux subtracted in 200 components = 0.00617724 Jy ! Clean residual min=-0.003592 max=0.003857 Jy/beam ! Clean residual mean=0.000459 rms=0.000885 Jy/beam ! Combined flux in latest and established models = 2.68645 Jy ! Adding 71 model components to the UV plane model. ! The established model now contains 561 components and 2.68645 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (9, 49.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00312726 Jy ! Component: 100 - total flux cleaned = 0.00348378 Jy ! Component: 150 - total flux cleaned = 0.00425016 Jy ! Component: 200 - total flux cleaned = 0.0047501 Jy ! Total flux subtracted in 200 components = 0.0047501 Jy ! Clean residual min=-0.003627 max=0.003821 Jy/beam ! Clean residual mean=0.000408 rms=0.000869 Jy/beam ! Combined flux in latest and established models = 2.6912 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 575 components and 2.6912 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Added new window around map position (22.8, -25). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.002881 Jy ! Component: 100 - total flux cleaned = 0.00338738 Jy ! Component: 150 - total flux cleaned = 0.00395341 Jy ! Component: 200 - total flux cleaned = 0.00435106 Jy ! Total flux subtracted in 200 components = 0.00435106 Jy ! Clean residual min=-0.003759 max=0.003738 Jy/beam ! Clean residual mean=0.000366 rms=0.000854 Jy/beam ! Combined flux in latest and established models = 2.69555 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 593 components and 2.69555 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 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 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=2.416 mas, bmaj=3.631 mas, bpa=-2.745 degrees ! Estimated noise=1.57413 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0030657 Jy ! Component: 100 - total flux cleaned = 0.00530422 Jy ! Component: 150 - total flux cleaned = 0.00751102 Jy ! Component: 200 - total flux cleaned = 0.00944259 Jy ! Total flux subtracted in 200 components = 0.00944259 Jy ! Clean residual min=-0.003455 max=0.004353 Jy/beam ! Clean residual mean=0.000573 rms=0.001097 Jy/beam ! Combined flux in latest and established models = 2.70499 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.00188621 Jy ! Component: 100 - total flux cleaned = 0.00357553 Jy ! Component: 150 - total flux cleaned = 0.00516715 Jy ! Component: 200 - total flux cleaned = 0.00603317 Jy ! Total flux subtracted in 200 components = 0.00603317 Jy ! Clean residual min=-0.003580 max=0.004091 Jy/beam ! Clean residual mean=0.000425 rms=0.001076 Jy/beam ! Combined flux in latest and established models = 2.71103 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 644 components and 2.71103 Jy ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! Estimated beam: bmin=0.7753 mas, bmaj=1.05 mas, bpa=14.2 degrees ! Estimated noise=1.18562 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00145541 Jy ! Component: 100 - total flux cleaned = -0.00218288 Jy ! Component: 150 - total flux cleaned = -0.00249377 Jy ! Component: 200 - total flux cleaned = -0.00279561 Jy ! Total flux subtracted in 200 components = -0.00279561 Jy ! Clean residual min=-0.003800 max=0.003482 Jy/beam ! Clean residual mean=0.000234 rms=0.000826 Jy/beam ! Combined flux in latest and established models = 2.70823 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 657 components and 2.70823 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.861907Jy sigma=0.710887 ! Fit after self-cal, rms=0.861745Jy sigma=0.703935 wmodel J2225-0457_X_map.mod ! Writing 657 model components to file: J2225-0457_X_map.mod wobs J2225-0457_X_uvs.fits ! Writing UV FITS file: J2225-0457_X_uvs.fits wwins J2225-0457_X_map.win ! wwins: Wrote 52 windows to J2225-0457_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 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000540543 Jy ! Component: 100 - total flux cleaned = 0.000912807 Jy ! Component: 150 - total flux cleaned = 0.00132857 Jy ! Component: 200 - total flux cleaned = 0.00190688 Jy ! Component: 250 - total flux cleaned = 0.00236147 Jy ! Component: 300 - total flux cleaned = 0.00286388 Jy ! Component: 350 - total flux cleaned = 0.00319389 Jy ! Component: 400 - total flux cleaned = 0.003628 Jy ! Component: 450 - total flux cleaned = 0.00416282 Jy ! Component: 500 - total flux cleaned = 0.00458537 Jy ! Component: 550 - total flux cleaned = 0.00505527 Jy ! Component: 600 - total flux cleaned = 0.00552015 Jy ! Component: 650 - total flux cleaned = 0.00582678 Jy ! Component: 700 - total flux cleaned = 0.006182 Jy ! Component: 750 - total flux cleaned = 0.0065838 Jy ! Component: 800 - total flux cleaned = 0.00698175 Jy ! Component: 850 - total flux cleaned = 0.00727756 Jy ! Component: 900 - total flux cleaned = 0.00752236 Jy ! Component: 950 - total flux cleaned = 0.00786202 Jy ! Component: 1000 - total flux cleaned = 0.00815066 Jy ! Component: 1050 - total flux cleaned = 0.00848477 Jy ! Component: 1100 - total flux cleaned = 0.00867413 Jy ! Component: 1150 - total flux cleaned = 0.00895624 Jy ! Component: 1200 - total flux cleaned = 0.00900298 Jy ! Component: 1250 - total flux cleaned = 0.00928164 Jy ! Component: 1300 - total flux cleaned = 0.00946613 Jy ! Component: 1350 - total flux cleaned = 0.00974119 Jy ! Component: 1400 - total flux cleaned = 0.00992316 Jy ! Component: 1450 - total flux cleaned = 0.010149 Jy ! Component: 1500 - total flux cleaned = 0.0102839 Jy ! Component: 1550 - total flux cleaned = 0.0103733 Jy ! Component: 1600 - total flux cleaned = 0.0105955 Jy ! Component: 1650 - total flux cleaned = 0.010772 Jy ! Component: 1700 - total flux cleaned = 0.010816 Jy ! Component: 1750 - total flux cleaned = 0.0109471 Jy ! Component: 1800 - total flux cleaned = 0.0111208 Jy ! Component: 1850 - total flux cleaned = 0.0112072 Jy ! Component: 1900 - total flux cleaned = 0.0112932 Jy ! Component: 1950 - total flux cleaned = 0.0113788 Jy ! Component: 2000 - total flux cleaned = 0.0115065 Jy ! Component: 2050 - total flux cleaned = 0.0116336 Jy ! Component: 2100 - total flux cleaned = 0.0116337 Jy ! Component: 2150 - total flux cleaned = 0.0118017 Jy ! Component: 2200 - total flux cleaned = 0.0118018 Jy ! Component: 2250 - total flux cleaned = 0.011885 Jy ! Component: 2300 - total flux cleaned = 0.0119265 Jy ! Component: 2350 - total flux cleaned = 0.0120501 Jy ! Component: 2400 - total flux cleaned = 0.0121322 Jy ! Component: 2450 - total flux cleaned = 0.012214 Jy ! Component: 2500 - total flux cleaned = 0.012336 Jy ! Component: 2550 - total flux cleaned = 0.0123764 Jy ! Component: 2600 - total flux cleaned = 0.0124167 Jy ! Component: 2650 - total flux cleaned = 0.0124972 Jy ! Component: 2700 - total flux cleaned = 0.0124571 Jy ! Component: 2750 - total flux cleaned = 0.0126962 Jy ! Component: 2800 - total flux cleaned = 0.0126961 Jy ! Component: 2850 - total flux cleaned = 0.012775 Jy ! Component: 2900 - total flux cleaned = 0.0127751 Jy ! Component: 2950 - total flux cleaned = 0.0129318 Jy ! Component: 3000 - total flux cleaned = 0.0128927 Jy ! Component: 3050 - total flux cleaned = 0.0129316 Jy ! Component: 3100 - total flux cleaned = 0.0130093 Jy ! Component: 3150 - total flux cleaned = 0.0131637 Jy ! Component: 3200 - total flux cleaned = 0.0132406 Jy ! Component: 3250 - total flux cleaned = 0.0132406 Jy ! Component: 3300 - total flux cleaned = 0.0133169 Jy ! Component: 3350 - total flux cleaned = 0.013431 Jy ! Component: 3400 - total flux cleaned = 0.0135067 Jy ! Component: 3450 - total flux cleaned = 0.0135068 Jy ! Component: 3500 - total flux cleaned = 0.0136572 Jy ! Component: 3550 - total flux cleaned = 0.0136945 Jy ! Component: 3600 - total flux cleaned = 0.0136946 Jy ! Component: 3650 - total flux cleaned = 0.0138433 Jy ! Component: 3700 - total flux cleaned = 0.0139172 Jy ! Component: 3750 - total flux cleaned = 0.0139172 Jy ! Component: 3800 - total flux cleaned = 0.0139907 Jy ! Component: 3850 - total flux cleaned = 0.0140641 Jy ! Component: 3900 - total flux cleaned = 0.0141371 Jy ! Component: 3950 - total flux cleaned = 0.014137 Jy ! Component: 4000 - total flux cleaned = 0.0142095 Jy ! Component: 4050 - total flux cleaned = 0.0142456 Jy ! Total flux subtracted in 4096 components = 0.0143177 Jy ! Clean residual min=-0.001936 max=0.002015 Jy/beam ! Clean residual mean=0.000068 rms=0.000607 Jy/beam ! Combined flux in latest and established models = 2.72255 Jy keep ! Adding 338 model components to the UV plane model. ! The established model now contains 995 components and 2.72255 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv123_uvs/J2225-0457_X_map.fits ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7753 x 1.05 at 14.2 degrees (North through East) ! Clean map min=-0.0037182 max=1.2576 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv123_uvs/J2225-0457_X_map.fits wmodel /scr/pima/rv123_uvs/J2225-0457_X_map.mod ! Writing 995 model components to file: /scr/pima/rv123_uvs/J2225-0457_X_map.mod wobs /scr/pima/rv123_uvs/J2225-0457_X_uvs.fits ! Writing UV FITS file: /scr/pima/rv123_uvs/J2225-0457_X_uvs.fits wwins /scr/pima/rv123_uvs/J2225-0457_X_map.win ! wwins: Wrote 53 windows to /scr/pima/rv123_uvs/J2225-0457_X_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 1.25765 invert ! Inverting map ! Your choice of large map pixels excluded 0.392% of the data. ! The x-axis pixel size should ideally be below 0.1771 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000606719 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.144728 0.144728 0.289455 0.57891 1.15782 2.31564 4.63128 9.26256 18.5251 37.0502 74.1005 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0018 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.257 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2072.86 quit ! Quitting program ! Log file /scr/pima/rv123_uvs/J2225-0457_X_dfm.log closed on Sun Dec 23 00:45:37 2018