! Started logfile: /scr/pima/rv123_uvs/J0405-1308_X_dfm.log on Sun Dec 23 00:39:53 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/J0405-1308_X_uva.fits ! Reading UV FITS file: /scr/pima/rv123_uvs/J0405-1308_X_uva.fits ! AN table 1: 22 integrations on 45 of 45 possible baselines. ! AN table 2: 3 integrations on 6 of 6 possible baselines. ! AN table 3: 3 integrations on 6 of 6 possible baselines. ! AN table 4: 125 integrations on 15 of 15 possible baselines. ! AN table 5: 44 integrations on 36 of 36 possible baselines. ! AN table 6: 12 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.427505 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0405-1308 ! ! 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 8244 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 J0405-1308_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 94 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 3. ! A total of 246 telescope corrections were flagged in sub-array 4. ! A total of 228 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 94 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 3. ! A total of 246 telescope corrections were flagged in sub-array 4. ! A total of 228 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 94 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 3. ! A total of 246 telescope corrections were flagged in sub-array 4. ! A total of 228 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 94 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 3. ! A total of 246 telescope corrections were flagged in sub-array 4. ! A total of 228 telescope corrections were flagged in sub-array 5. ! A total of 28 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.204992Jy sigma=13.096761 ! Fit after self-cal, rms=0.413446Jy sigma=3.745712 ! 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 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! Estimated beam: bmin=0.6099 mas, bmaj=0.8458 mas, bpa=21.18 degrees ! Estimated noise=1.18033 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.566158 Jy ! Component: 100 - total flux cleaned = 0.707795 Jy ! Total flux subtracted in 100 components = 0.707795 Jy ! Clean residual min=-0.024116 max=0.050461 Jy/beam ! Clean residual mean=0.000137 rms=0.004659 Jy/beam ! Combined flux in latest and established models = 0.707795 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.707795 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227602Jy sigma=1.988376 ! Fit after self-cal, rms=0.227555Jy sigma=1.988104 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0484442 Jy ! Component: 100 - total flux cleaned = 0.0689941 Jy ! Total flux subtracted in 100 components = 0.0689941 Jy ! Clean residual min=-0.024876 max=0.028287 Jy/beam ! Clean residual mean=0.000101 rms=0.004105 Jy/beam ! Combined flux in latest and established models = 0.77679 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 14 components and 0.77679 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.224268Jy sigma=1.954460 ! Fit after self-cal, rms=0.224257Jy sigma=1.954353 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! Added new window around map position (3.2, 4.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0180361 Jy ! Component: 100 - total flux cleaned = 0.0196896 Jy ! Total flux subtracted in 100 components = 0.0196896 Jy ! Clean residual min=-0.023923 max=0.027753 Jy/beam ! Clean residual mean=0.000091 rms=0.003991 Jy/beam ! Combined flux in latest and established models = 0.796479 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 17 components and 0.796479 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220948Jy sigma=1.925198 ! Fit after self-cal, rms=0.220735Jy sigma=1.921747 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! Added new window around map position (0.8, -2.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0172698 Jy ! Component: 100 - total flux cleaned = 0.0184868 Jy ! Total flux subtracted in 100 components = 0.0184868 Jy ! Clean residual min=-0.022321 max=0.024775 Jy/beam ! Clean residual mean=0.000081 rms=0.003858 Jy/beam ! Combined flux in latest and established models = 0.814966 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 21 components and 0.814966 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.217459Jy sigma=1.907109 ! Fit after self-cal, rms=0.217108Jy sigma=1.904304 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! Added new window around map position (1, 2.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.015625 Jy ! Component: 100 - total flux cleaned = 0.0157042 Jy ! Total flux subtracted in 100 components = 0.0157042 Jy ! Clean residual min=-0.021596 max=0.019322 Jy/beam ! Clean residual mean=0.000073 rms=0.003719 Jy/beam ! Combined flux in latest and established models = 0.83067 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 25 components and 0.83067 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.215021Jy sigma=1.896342 ! Fit after self-cal, rms=0.214871Jy sigma=1.894911 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 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 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 milli-arcsec ! Estimated beam: bmin=0.6131 mas, bmaj=0.9353 mas, bpa=44.58 degrees ! Estimated noise=0.833469 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.0349679 Jy ! Component: 100 - total flux cleaned = -0.0366214 Jy ! Component: 150 - total flux cleaned = -0.0367402 Jy ! Component: 200 - total flux cleaned = -0.0346338 Jy ! Total flux subtracted in 200 components = -0.0346338 Jy ! Clean residual min=-0.025812 max=0.038127 Jy/beam ! Clean residual mean=0.007394 rms=0.007945 Jy/beam ! Combined flux in latest and established models = 0.796036 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 37 components and 0.796036 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.205796Jy sigma=1.836962 ! Fit after self-cal, rms=0.205511Jy sigma=1.833768 ! Inverting map ! Your choice of large map pixels excluded 3.34% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1971 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.6131 x 0.9353 at 44.58 degrees (North through East) ! Clean map min=-0.022501 max=0.69641 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.85 HN 1.07 KP 1.01 ! LA 0.97 MK 0.93 NL 1.01 OV 0.94 ! PT 0.86 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! KP 0.97 MK 0.93 OV 0.97 YG 1.07 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 MK 0.83 OV 0.98 YG 1.20 ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! HH 0.87 MC 0.95 ON 1.06 WN 1.20 ! WZ 0.92 YG 0.97 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.91 HH 0.73 HN 1.11 LA 1.02 ! NL 1.04 ON 0.99 PT 0.87 SC 1.25 ! WZ 1.08 ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.85 HH 0.72 HN 1.02 KP 1.02 ! LA 0.96 NL 1.05 OV 0.95 SC 0.89 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.86 HN 1.08 KP 1.01 ! LA 0.98 MK 0.94 NL 1.01 OV 0.93 ! PT 0.85 SC 1.03 ! ! Telescope amplitude corrections in sub-array 2: ! KP 0.98 MK 0.93 OV 0.98 YG 0.97 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.02 MK 0.88 OV 0.96 YG 0.90 ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! HH 0.83 MC 1.03 ON 0.99 WN 1.23 ! WZ 0.92 YG 0.88 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.89 HH 0.72 HN 1.10 LA 1.01 ! NL 1.03 ON 0.97 PT 0.86 SC 1.35 ! WZ 1.05 ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.84 HH 0.69 HN 1.04 KP 1.01 ! LA 0.96 NL 1.05 OV 0.96 SC 0.88 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.84 HN 1.03 KP 1.01 ! LA 0.98 MK 0.95 NL 1.01 OV 0.93 ! PT 0.86 SC 0.99 ! ! Telescope amplitude corrections in sub-array 2: ! KP 1.03 MK 0.87 OV 0.99 YG 0.99 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.05 MK 0.91 OV 0.93 YG 1.20 ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! HH 0.95 MC 1.55 ON 1.44 WN 1.02 ! WZ 0.92 YG 0.97 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.89 HH 0.73 HN 1.10 LA 1.01 ! NL 1.05 ON 1.30 PT 0.86 SC 1.33 ! WZ 1.17 ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.84 HH 0.72 HN 1.02 KP 1.05 ! LA 0.97 NL 1.05 OV 0.95 SC 0.88 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.85 HN 1.07 KP 1.03 ! LA 0.97 MK 0.92 NL 1.03 OV 0.93 ! PT 0.86 SC 1.02 ! ! Telescope amplitude corrections in sub-array 2: ! KP 1.04 MK 0.95 OV 0.96 YG 1.01 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 MK 0.83 OV 1.01 YG 1.20 ! ! A total of 72 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! HH 1.07 MC 1.49 ON 1.40 WN 0.93 ! WZ 0.94 YG 1.06 ! ! A total of 54 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.90 HH 0.85 HN 1.11 LA 1.02 ! NL 1.03 ON 1.29 PT 0.86 SC 1.30 ! WZ 1.19 ! ! Telescope amplitude corrections in sub-array 6: ! FD 0.83 HH 0.88 HN 1.03 KP 1.01 ! LA 0.97 NL 1.06 OV 0.96 SC 0.90 ! ! ! Fit before self-cal, rms=0.205511Jy sigma=1.833768 ! Fit after self-cal, rms=0.162260Jy sigma=1.170536 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6237 mas, bmaj=0.9383 mas, bpa=42.73 degrees ! Estimated noise=0.845107 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.2, -6.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0201211 Jy ! Component: 100 - total flux cleaned = 0.0288115 Jy ! Component: 150 - total flux cleaned = 0.030643 Jy ! Component: 200 - total flux cleaned = 0.0325718 Jy ! Total flux subtracted in 200 components = 0.0325718 Jy ! Clean residual min=-0.008237 max=0.017526 Jy/beam ! Clean residual mean=0.005697 rms=0.002781 Jy/beam ! Combined flux in latest and established models = 0.828608 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 60 components and 0.828608 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.157384Jy sigma=1.148357 ! Fit after self-cal, rms=0.157044Jy sigma=1.140563 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-9.8, 0). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0174998 Jy ! Component: 100 - total flux cleaned = 0.025908 Jy ! Component: 150 - total flux cleaned = 0.0306666 Jy ! Component: 200 - total flux cleaned = 0.0329499 Jy ! Total flux subtracted in 200 components = 0.0329499 Jy ! Clean residual min=-0.006128 max=0.016659 Jy/beam ! Clean residual mean=0.004828 rms=0.002461 Jy/beam ! Combined flux in latest and established models = 0.861558 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 80 components and 0.861558 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.153561Jy sigma=1.122222 ! Fit after self-cal, rms=0.153328Jy sigma=1.118608 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (28.2, -15.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0140305 Jy ! Component: 100 - total flux cleaned = 0.0189673 Jy ! Component: 150 - total flux cleaned = 0.0230082 Jy ! Component: 200 - total flux cleaned = 0.024466 Jy ! Total flux subtracted in 200 components = 0.024466 Jy ! Clean residual min=-0.005025 max=0.013370 Jy/beam ! Clean residual mean=0.004164 rms=0.002180 Jy/beam ! Combined flux in latest and established models = 0.886024 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 95 components and 0.886024 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.150964Jy sigma=1.107776 ! Fit after self-cal, rms=0.150902Jy sigma=1.105317 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-6.2, -17.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0130382 Jy ! Component: 100 - total flux cleaned = 0.0195567 Jy ! Component: 150 - total flux cleaned = 0.0229226 Jy ! Component: 200 - total flux cleaned = 0.0252151 Jy ! Total flux subtracted in 200 components = 0.0252151 Jy ! Clean residual min=-0.004972 max=0.011977 Jy/beam ! Clean residual mean=0.003497 rms=0.002135 Jy/beam ! Combined flux in latest and established models = 0.911239 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 110 components and 0.911239 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148771Jy sigma=1.096325 ! Fit after self-cal, rms=0.148648Jy sigma=1.094465 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (9.8, 1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0111938 Jy ! Component: 100 - total flux cleaned = 0.0157734 Jy ! Component: 150 - total flux cleaned = 0.0183587 Jy ! Component: 200 - total flux cleaned = 0.0197266 Jy ! Total flux subtracted in 200 components = 0.0197266 Jy ! Clean residual min=-0.004872 max=0.011157 Jy/beam ! Clean residual mean=0.002973 rms=0.001978 Jy/beam ! Combined flux in latest and established models = 0.930966 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 125 components and 0.930966 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146961Jy sigma=1.088497 ! Fit after self-cal, rms=0.146854Jy sigma=1.087076 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (20.8, -31.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0108478 Jy ! Component: 100 - total flux cleaned = 0.0162281 Jy ! Component: 150 - total flux cleaned = 0.0173635 Jy ! Component: 200 - total flux cleaned = 0.0183776 Jy ! Total flux subtracted in 200 components = 0.0183776 Jy ! Clean residual min=-0.004739 max=0.009578 Jy/beam ! Clean residual mean=0.002473 rms=0.001835 Jy/beam ! Combined flux in latest and established models = 0.949343 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 143 components and 0.949343 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.145516Jy sigma=1.078542 ! Fit after self-cal, rms=0.145443Jy sigma=1.077111 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (21.4, 13.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00848128 Jy ! Component: 100 - total flux cleaned = 0.011471 Jy ! Component: 150 - total flux cleaned = 0.0123647 Jy ! Component: 200 - total flux cleaned = 0.0128543 Jy ! Total flux subtracted in 200 components = 0.0128543 Jy ! Clean residual min=-0.004591 max=0.009479 Jy/beam ! Clean residual mean=0.002145 rms=0.001716 Jy/beam ! Combined flux in latest and established models = 0.962198 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 161 components and 0.962198 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144465Jy sigma=1.071712 ! Fit after self-cal, rms=0.144414Jy sigma=1.070888 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (7.6, 14). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00928651 Jy ! Component: 100 - total flux cleaned = 0.0133835 Jy ! Component: 150 - total flux cleaned = 0.0138864 Jy ! Component: 200 - total flux cleaned = 0.014822 Jy ! Total flux subtracted in 200 components = 0.014822 Jy ! Clean residual min=-0.004768 max=0.008514 Jy/beam ! Clean residual mean=0.001757 rms=0.001647 Jy/beam ! Combined flux in latest and established models = 0.97702 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 177 components and 0.97702 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143316Jy sigma=1.064319 ! Fit after self-cal, rms=0.143294Jy sigma=1.063512 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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.143294Jy sigma=1.063512 ! Fit after self-cal, rms=0.136874Jy sigma=1.014748 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6197 mas, bmaj=0.9369 mas, bpa=43.24 degrees ! Estimated noise=0.841245 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00191953 Jy ! Component: 100 - total flux cleaned = 0.00297007 Jy ! Component: 150 - total flux cleaned = 0.00472658 Jy ! Component: 200 - total flux cleaned = 0.00604652 Jy ! Total flux subtracted in 200 components = 0.00604652 Jy ! Clean residual min=-0.003826 max=0.007853 Jy/beam ! Clean residual mean=0.001985 rms=0.001436 Jy/beam ! Combined flux in latest and established models = 0.983066 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 185 components and 0.983066 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136245Jy sigma=1.011993 ! Fit after self-cal, rms=0.136315Jy sigma=1.011780 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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 (42.6, 15). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0070643 Jy ! Component: 100 - total flux cleaned = 0.00722866 Jy ! Component: 150 - total flux cleaned = 0.00777981 Jy ! Component: 200 - total flux cleaned = 0.00816034 Jy ! Total flux subtracted in 200 components = 0.00816034 Jy ! Clean residual min=-0.003869 max=0.007355 Jy/beam ! Clean residual mean=0.001758 rms=0.001371 Jy/beam ! Combined flux in latest and established models = 0.991227 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 197 components and 0.991227 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135573Jy sigma=1.007760 ! Fit after self-cal, rms=0.135536Jy sigma=1.007118 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-20.2, -10.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00625852 Jy ! Component: 100 - total flux cleaned = 0.00773923 Jy ! Component: 150 - total flux cleaned = 0.00798269 Jy ! Component: 200 - total flux cleaned = 0.00809383 Jy ! Total flux subtracted in 200 components = 0.00809383 Jy ! Clean residual min=-0.004072 max=0.007530 Jy/beam ! Clean residual mean=0.001561 rms=0.001348 Jy/beam ! Combined flux in latest and established models = 0.999321 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 208 components and 0.999321 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134986Jy sigma=1.003461 ! Fit after self-cal, rms=0.134959Jy sigma=1.002936 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-16.2, -1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00625441 Jy ! Component: 100 - total flux cleaned = 0.00916612 Jy ! Component: 150 - total flux cleaned = 0.0100088 Jy ! Component: 200 - total flux cleaned = 0.0104468 Jy ! Total flux subtracted in 200 components = 0.0104468 Jy ! Clean residual min=-0.003998 max=0.006497 Jy/beam ! Clean residual mean=0.001294 rms=0.001308 Jy/beam ! Combined flux in latest and established models = 1.00977 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 221 components and 1.00977 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134407Jy sigma=1.000191 ! Fit after self-cal, rms=0.134374Jy sigma=0.999658 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-32.4, 46). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0055732 Jy ! Component: 100 - total flux cleaned = 0.00708492 Jy ! Component: 150 - total flux cleaned = 0.00708231 Jy ! Component: 200 - total flux cleaned = 0.00718907 Jy ! Total flux subtracted in 200 components = 0.00718907 Jy ! Clean residual min=-0.003956 max=0.005974 Jy/beam ! Clean residual mean=0.001086 rms=0.001242 Jy/beam ! Combined flux in latest and established models = 1.01696 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 237 components and 1.01696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133968Jy sigma=0.996606 ! Fit after self-cal, rms=0.133942Jy sigma=0.996236 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-5.2, -5.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0058228 Jy ! Component: 100 - total flux cleaned = 0.00640704 Jy ! Component: 150 - total flux cleaned = 0.00677063 Jy ! Component: 200 - total flux cleaned = 0.00720931 Jy ! Total flux subtracted in 200 components = 0.00720931 Jy ! Clean residual min=-0.004012 max=0.005688 Jy/beam ! Clean residual mean=0.000917 rms=0.001201 Jy/beam ! Combined flux in latest and established models = 1.02417 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 256 components and 1.02417 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133548Jy sigma=0.991834 ! Fit after self-cal, rms=0.133532Jy sigma=0.991556 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-37.8, 25). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00488428 Jy ! Component: 100 - total flux cleaned = 0.00466434 Jy ! Component: 150 - total flux cleaned = 0.00445032 Jy ! Component: 200 - total flux cleaned = 0.00435022 Jy ! Total flux subtracted in 200 components = 0.00435022 Jy ! Clean residual min=-0.004246 max=0.005418 Jy/beam ! Clean residual mean=0.000796 rms=0.001159 Jy/beam ! Combined flux in latest and established models = 1.02852 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 271 components and 1.02852 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133174Jy sigma=0.988538 ! Fit after self-cal, rms=0.133150Jy sigma=0.988211 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (10.4, 39.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00508132 Jy ! Component: 100 - total flux cleaned = 0.00594346 Jy ! Component: 150 - total flux cleaned = 0.00553679 Jy ! Component: 200 - total flux cleaned = 0.00544163 Jy ! Total flux subtracted in 200 components = 0.00544163 Jy ! Clean residual min=-0.004066 max=0.005007 Jy/beam ! Clean residual mean=0.000654 rms=0.001110 Jy/beam ! Combined flux in latest and established models = 1.03396 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 285 components and 1.03396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132881Jy sigma=0.984777 ! Fit after self-cal, rms=0.132856Jy sigma=0.984551 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-49.4, -21.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00372246 Jy ! Component: 100 - total flux cleaned = 0.0033235 Jy ! Component: 150 - total flux cleaned = 0.00285428 Jy ! Component: 200 - total flux cleaned = 0.00267349 Jy ! Total flux subtracted in 200 components = 0.00267349 Jy ! Clean residual min=-0.004014 max=0.004952 Jy/beam ! Clean residual mean=0.000577 rms=0.001092 Jy/beam ! Combined flux in latest and established models = 1.03663 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 296 components and 1.03663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132611Jy sigma=0.982151 ! Fit after self-cal, rms=0.132598Jy sigma=0.981976 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-24.6, -6.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00432048 Jy ! Component: 100 - total flux cleaned = 0.00455299 Jy ! Component: 150 - total flux cleaned = 0.00395997 Jy ! Component: 200 - total flux cleaned = 0.0035006 Jy ! Total flux subtracted in 200 components = 0.0035006 Jy ! Clean residual min=-0.003984 max=0.004854 Jy/beam ! Clean residual mean=0.000498 rms=0.001081 Jy/beam ! Combined flux in latest and established models = 1.04013 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 312 components and 1.04013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132256Jy sigma=0.979637 ! Fit after self-cal, rms=0.132245Jy sigma=0.979437 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (17.2, 14.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0033295 Jy ! Component: 100 - total flux cleaned = 0.00285371 Jy ! Component: 150 - total flux cleaned = 0.00221953 Jy ! Component: 200 - total flux cleaned = 0.00178528 Jy ! Total flux subtracted in 200 components = 0.00178528 Jy ! Clean residual min=-0.004091 max=0.004621 Jy/beam ! Clean residual mean=0.000455 rms=0.001069 Jy/beam ! Combined flux in latest and established models = 1.04192 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 319 components and 1.04192 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132027Jy sigma=0.977114 ! Fit after self-cal, rms=0.132030Jy sigma=0.976949 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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.132030Jy sigma=0.976949 ! Fit after self-cal, rms=0.129182Jy sigma=0.953306 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.621 mas, bmaj=0.9356 mas, bpa=43.38 degrees ! Estimated noise=0.841039 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00446344 Jy ! Component: 100 - total flux cleaned = 0.00810475 Jy ! Component: 150 - total flux cleaned = 0.0112933 Jy ! Component: 200 - total flux cleaned = 0.0132421 Jy ! Total flux subtracted in 200 components = 0.0132421 Jy ! Clean residual min=-0.002904 max=0.005802 Jy/beam ! Clean residual mean=0.000806 rms=0.000936 Jy/beam ! Combined flux in latest and established models = 1.05516 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 345 components and 1.05516 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128752Jy sigma=0.951052 ! Fit after self-cal, rms=0.128755Jy sigma=0.950838 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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 (12.6, -2.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00508273 Jy ! Component: 100 - total flux cleaned = 0.00658215 Jy ! Component: 150 - total flux cleaned = 0.00770271 Jy ! Component: 200 - total flux cleaned = 0.00855041 Jy ! Total flux subtracted in 200 components = 0.00855041 Jy ! Clean residual min=-0.002991 max=0.004273 Jy/beam ! Clean residual mean=0.000600 rms=0.000886 Jy/beam ! Combined flux in latest and established models = 1.06371 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 358 components and 1.06371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128414Jy sigma=0.948203 ! Fit after self-cal, rms=0.128394Jy sigma=0.947952 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (17.4, 8.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00375129 Jy ! Component: 100 - total flux cleaned = 0.0043412 Jy ! Component: 150 - total flux cleaned = 0.00478114 Jy ! Component: 200 - total flux cleaned = 0.00520228 Jy ! Total flux subtracted in 200 components = 0.00520228 Jy ! Clean residual min=-0.003208 max=0.003838 Jy/beam ! Clean residual mean=0.000467 rms=0.000862 Jy/beam ! Combined flux in latest and established models = 1.06891 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 373 components and 1.06891 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128189Jy sigma=0.946097 ! Fit after self-cal, rms=0.128174Jy sigma=0.945974 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-30.6, -8.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00289809 Jy ! Component: 100 - total flux cleaned = 0.00281385 Jy ! Component: 150 - total flux cleaned = 0.00272976 Jy ! Component: 200 - total flux cleaned = 0.00272664 Jy ! Total flux subtracted in 200 components = 0.00272664 Jy ! Clean residual min=-0.003100 max=0.003657 Jy/beam ! Clean residual mean=0.000399 rms=0.000837 Jy/beam ! Combined flux in latest and established models = 1.07164 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 390 components and 1.07164 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127986Jy sigma=0.944220 ! Fit after self-cal, rms=0.127975Jy sigma=0.944093 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-2.8, 7.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00279671 Jy ! Component: 100 - total flux cleaned = 0.00322938 Jy ! Component: 150 - total flux cleaned = 0.00338527 Jy ! Component: 200 - total flux cleaned = 0.00345721 Jy ! Total flux subtracted in 200 components = 0.00345721 Jy ! Clean residual min=-0.003135 max=0.003617 Jy/beam ! Clean residual mean=0.000314 rms=0.000827 Jy/beam ! Combined flux in latest and established models = 1.0751 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 403 components and 1.0751 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127804Jy sigma=0.942092 ! Fit after self-cal, rms=0.127796Jy sigma=0.942020 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-27.2, 16.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00311175 Jy ! Component: 100 - total flux cleaned = 0.00279416 Jy ! Component: 150 - total flux cleaned = 0.00271977 Jy ! Component: 200 - total flux cleaned = 0.00243367 Jy ! Total flux subtracted in 200 components = 0.00243367 Jy ! Clean residual min=-0.003116 max=0.003372 Jy/beam ! Clean residual mean=0.000251 rms=0.000787 Jy/beam ! Combined flux in latest and established models = 1.07753 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 420 components and 1.07753 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127639Jy sigma=0.940150 ! Fit after self-cal, rms=0.127627Jy sigma=0.940048 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (4, 11). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00184578 Jy ! Component: 100 - total flux cleaned = 0.00200176 Jy ! Component: 150 - total flux cleaned = 0.00171911 Jy ! Component: 200 - total flux cleaned = 0.00145165 Jy ! Total flux subtracted in 200 components = 0.00145165 Jy ! Clean residual min=-0.003064 max=0.003258 Jy/beam ! Clean residual mean=0.000218 rms=0.000769 Jy/beam ! Combined flux in latest and established models = 1.07898 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 428 components and 1.07898 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127520Jy sigma=0.938724 ! Fit after self-cal, rms=0.127518Jy sigma=0.938661 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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.127518Jy sigma=0.938661 ! Fit after self-cal, rms=0.127602Jy sigma=0.936078 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=2.503 mas, bmaj=4.212 mas, bpa=-7.739 degrees ! Estimated noise=1.05966 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0043617 Jy ! Component: 100 - total flux cleaned = 0.00769657 Jy ! Component: 150 - total flux cleaned = 0.01046 Jy ! Component: 200 - total flux cleaned = 0.0127985 Jy ! Total flux subtracted in 200 components = 0.0127985 Jy ! Clean residual min=-0.002660 max=0.004183 Jy/beam ! Clean residual mean=0.000804 rms=0.000958 Jy/beam ! Combined flux in latest and established models = 1.09178 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 452 components and 1.09178 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127470Jy sigma=0.935094 ! Fit after self-cal, rms=0.127484Jy sigma=0.934978 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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 (37.8, -22.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00376066 Jy ! Component: 100 - total flux cleaned = 0.00498369 Jy ! Component: 150 - total flux cleaned = 0.0058897 Jy ! Component: 200 - total flux cleaned = 0.00643918 Jy ! Total flux subtracted in 200 components = 0.00643918 Jy ! Clean residual min=-0.002837 max=0.003425 Jy/beam ! Clean residual mean=0.000442 rms=0.000902 Jy/beam ! Combined flux in latest and established models = 1.09822 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 473 components and 1.09822 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127444Jy sigma=0.933608 ! Fit after self-cal, rms=0.127442Jy sigma=0.933437 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6211 mas, bmaj=0.9352 mas, bpa=43.47 degrees ! Estimated noise=0.840491 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6211 x 0.9352 at 43.47 degrees (North through East) ! Clean map min=-0.011046 max=0.69502 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.127442Jy sigma=0.933437 ! Fit after self-cal, rms=0.116268Jy sigma=0.784897 ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6206 mas, bmaj=0.9352 mas, bpa=43.48 degrees ! Estimated noise=0.839972 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00283618 Jy ! Component: 100 - total flux cleaned = 0.00437253 Jy ! Component: 150 - total flux cleaned = 0.0057718 Jy ! Component: 200 - total flux cleaned = 0.00689945 Jy ! Total flux subtracted in 200 components = 0.00689945 Jy ! Clean residual min=-0.002763 max=0.003327 Jy/beam ! Clean residual mean=0.000348 rms=0.000701 Jy/beam ! Combined flux in latest and established models = 1.10512 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 500 components and 1.10512 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.116056Jy sigma=0.782812 ! Fit after self-cal, rms=0.116053Jy sigma=0.782734 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (9.6, 17.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.002371 Jy ! Component: 100 - total flux cleaned = 0.00271188 Jy ! Component: 150 - total flux cleaned = 0.00303098 Jy ! Component: 200 - total flux cleaned = 0.00326139 Jy ! Total flux subtracted in 200 components = 0.00326139 Jy ! Clean residual min=-0.002684 max=0.003085 Jy/beam ! Clean residual mean=0.000273 rms=0.000685 Jy/beam ! Combined flux in latest and established models = 1.10838 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 520 components and 1.10838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115888Jy sigma=0.780823 ! Fit after self-cal, rms=0.115875Jy sigma=0.780752 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-40, -17.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00259786 Jy ! Component: 100 - total flux cleaned = 0.00244301 Jy ! Component: 150 - total flux cleaned = 0.00222087 Jy ! Component: 200 - total flux cleaned = 0.00236571 Jy ! Total flux subtracted in 200 components = 0.00236571 Jy ! Clean residual min=-0.002470 max=0.002878 Jy/beam ! Clean residual mean=0.000212 rms=0.000664 Jy/beam ! Combined flux in latest and established models = 1.11074 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 538 components and 1.11075 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115721Jy sigma=0.779215 ! Fit after self-cal, rms=0.115722Jy sigma=0.779111 ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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.115722Jy sigma=0.779111 ! Fit after self-cal, rms=0.115782Jy sigma=0.778582 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 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6206 mas, bmaj=0.8618 mas, bpa=21.12 degrees ! Estimated noise=1.14333 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.54467 Jy ! Component: 100 - total flux cleaned = 0.679845 Jy ! Component: 150 - total flux cleaned = 0.727167 Jy ! Component: 200 - total flux cleaned = 0.749008 Jy ! Total flux subtracted in 200 components = 0.749008 Jy ! Clean residual min=-0.012854 max=0.012093 Jy/beam ! Clean residual mean=0.000147 rms=0.001899 Jy/beam ! Combined flux in latest and established models = 0.749008 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 16 model components to the UV plane model. ! The established model now contains 16 components and 0.749008 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6203 mas, bmaj=0.9351 mas, bpa=43.51 degrees ! Estimated noise=0.83976 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) ! Added new window around map position (-2.2, -8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0299408 Jy ! Component: 100 - total flux cleaned = 0.0546491 Jy ! Component: 150 - total flux cleaned = 0.0765677 Jy ! Component: 200 - total flux cleaned = 0.0963982 Jy ! Total flux subtracted in 200 components = 0.0963982 Jy ! Clean residual min=-0.012016 max=0.016284 Jy/beam ! Clean residual mean=0.006887 rms=0.002298 Jy/beam ! Combined flux in latest and established models = 0.845407 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 90 components and 0.845407 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-14.6, 2.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0189265 Jy ! Component: 100 - total flux cleaned = 0.033796 Jy ! Component: 150 - total flux cleaned = 0.0470531 Jy ! Component: 200 - total flux cleaned = 0.0565753 Jy ! Total flux subtracted in 200 components = 0.0565753 Jy ! Clean residual min=-0.010611 max=0.012496 Jy/beam ! Clean residual mean=0.005517 rms=0.001863 Jy/beam ! Combined flux in latest and established models = 0.901982 Jy ! Adding 82 model components to the UV plane model. ! The established model now contains 129 components and 0.901982 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-24.6, -38.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0111873 Jy ! Component: 100 - total flux cleaned = 0.0203078 Jy ! Component: 150 - total flux cleaned = 0.0289975 Jy ! Component: 200 - total flux cleaned = 0.0367972 Jy ! Total flux subtracted in 200 components = 0.0367972 Jy ! Clean residual min=-0.008964 max=0.010623 Jy/beam ! Clean residual mean=0.004618 rms=0.001580 Jy/beam ! Combined flux in latest and established models = 0.93878 Jy ! Adding 86 model components to the UV plane model. ! The established model now contains 158 components and 0.93878 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-12.4, -16.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00902039 Jy ! Component: 100 - total flux cleaned = 0.0166028 Jy ! Component: 150 - total flux cleaned = 0.0234111 Jy ! Component: 200 - total flux cleaned = 0.0299238 Jy ! Total flux subtracted in 200 components = 0.0299238 Jy ! Clean residual min=-0.007441 max=0.009316 Jy/beam ! Clean residual mean=0.003891 rms=0.001379 Jy/beam ! Combined flux in latest and established models = 0.968703 Jy ! Adding 90 model components to the UV plane model. ! The established model now contains 188 components and 0.968704 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (1.4, -11.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00812867 Jy ! Component: 100 - total flux cleaned = 0.0141242 Jy ! Component: 150 - total flux cleaned = 0.0198554 Jy ! Component: 200 - total flux cleaned = 0.0250115 Jy ! Total flux subtracted in 200 components = 0.0250115 Jy ! Clean residual min=-0.006397 max=0.007792 Jy/beam ! Clean residual mean=0.003283 rms=0.001219 Jy/beam ! Combined flux in latest and established models = 0.993715 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 208 components and 0.993715 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-0.6, 1.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00187448 Jy ! Component: 100 - total flux cleaned = 0.00743105 Jy ! Component: 150 - total flux cleaned = 0.0133312 Jy ! Component: 200 - total flux cleaned = 0.0189582 Jy ! Total flux subtracted in 200 components = 0.0189582 Jy ! Clean residual min=-0.004823 max=0.006964 Jy/beam ! Clean residual mean=0.002820 rms=0.001100 Jy/beam ! Combined flux in latest and established models = 1.01267 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 236 components and 1.01267 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (29.4, 4.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00633533 Jy ! Component: 100 - total flux cleaned = 0.0104362 Jy ! Component: 150 - total flux cleaned = 0.01463 Jy ! Component: 200 - total flux cleaned = 0.018408 Jy ! Total flux subtracted in 200 components = 0.018408 Jy ! Clean residual min=-0.004196 max=0.006222 Jy/beam ! Clean residual mean=0.002371 rms=0.000997 Jy/beam ! Combined flux in latest and established models = 1.03108 Jy ! Adding 86 model components to the UV plane model. ! The established model now contains 259 components and 1.03108 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (3.4, 27). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00540082 Jy ! Component: 100 - total flux cleaned = 0.00891789 Jy ! Component: 150 - total flux cleaned = 0.0122985 Jy ! Component: 200 - total flux cleaned = 0.0157593 Jy ! Total flux subtracted in 200 components = 0.0157593 Jy ! Clean residual min=-0.003591 max=0.005822 Jy/beam ! Clean residual mean=0.001985 rms=0.000936 Jy/beam ! Combined flux in latest and established models = 1.04684 Jy ! Adding 88 model components to the UV plane model. ! The established model now contains 284 components and 1.04684 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-7.2, -12.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00461933 Jy ! Component: 100 - total flux cleaned = 0.00767611 Jy ! Component: 150 - total flux cleaned = 0.0106031 Jy ! Component: 200 - total flux cleaned = 0.0134254 Jy ! Total flux subtracted in 200 components = 0.0134254 Jy ! Clean residual min=-0.003251 max=0.005086 Jy/beam ! Clean residual mean=0.001660 rms=0.000882 Jy/beam ! Combined flux in latest and established models = 1.06026 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 299 components and 1.06026 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-21.8, -8.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00401206 Jy ! Component: 100 - total flux cleaned = 0.00648068 Jy ! Component: 150 - total flux cleaned = 0.0088544 Jy ! Component: 200 - total flux cleaned = 0.0111484 Jy ! Total flux subtracted in 200 components = 0.0111484 Jy ! Clean residual min=-0.002957 max=0.004632 Jy/beam ! Clean residual mean=0.001388 rms=0.000839 Jy/beam ! Combined flux in latest and established models = 1.07141 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 316 components and 1.07141 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (5.8, 16.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00430773 Jy ! Component: 100 - total flux cleaned = 0.00632189 Jy ! Component: 150 - total flux cleaned = 0.00840268 Jy ! Component: 200 - total flux cleaned = 0.0102666 Jy ! Total flux subtracted in 200 components = 0.0102666 Jy ! Clean residual min=-0.002914 max=0.004361 Jy/beam ! Clean residual mean=0.001136 rms=0.000808 Jy/beam ! Combined flux in latest and established models = 1.08168 Jy ! Adding 72 model components to the UV plane model. ! The established model now contains 327 components and 1.08168 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (13.2, -5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00343548 Jy ! Component: 100 - total flux cleaned = 0.00520974 Jy ! Component: 150 - total flux cleaned = 0.00691385 Jy ! Component: 200 - total flux cleaned = 0.0084347 Jy ! Total flux subtracted in 200 components = 0.0084347 Jy ! Clean residual min=-0.002719 max=0.003973 Jy/beam ! Clean residual mean=0.000929 rms=0.000774 Jy/beam ! Combined flux in latest and established models = 1.09011 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 342 components and 1.09011 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (49.2, -8.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00331851 Jy ! Component: 100 - total flux cleaned = 0.00464018 Jy ! Component: 150 - total flux cleaned = 0.00602674 Jy ! Component: 200 - total flux cleaned = 0.00725983 Jy ! Total flux subtracted in 200 components = 0.00725983 Jy ! Clean residual min=-0.002552 max=0.003625 Jy/beam ! Clean residual mean=0.000744 rms=0.000741 Jy/beam ! Combined flux in latest and established models = 1.09737 Jy ! Adding 70 model components to the UV plane model. ! The established model now contains 357 components and 1.09737 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-43.2, 30.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00299046 Jy ! Component: 100 - total flux cleaned = 0.00385944 Jy ! Component: 150 - total flux cleaned = 0.00480056 Jy ! Component: 200 - total flux cleaned = 0.00561108 Jy ! Total flux subtracted in 200 components = 0.00561108 Jy ! Clean residual min=-0.002762 max=0.003436 Jy/beam ! Clean residual mean=0.000607 rms=0.000718 Jy/beam ! Combined flux in latest and established models = 1.10298 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 369 components and 1.10298 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-17.4, 36.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00230094 Jy ! Component: 100 - total flux cleaned = 0.00317435 Jy ! Component: 150 - total flux cleaned = 0.00392648 Jy ! Component: 200 - total flux cleaned = 0.00466131 Jy ! Total flux subtracted in 200 components = 0.00466131 Jy ! Clean residual min=-0.002982 max=0.003235 Jy/beam ! Clean residual mean=0.000498 rms=0.000699 Jy/beam ! Combined flux in latest and established models = 1.10765 Jy ! Adding 70 model components to the UV plane model. ! The established model now contains 382 components and 1.10765 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (-17.6, -24.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00195366 Jy ! Component: 100 - total flux cleaned = 0.00240009 Jy ! Component: 150 - total flux cleaned = 0.00283232 Jy ! Component: 200 - total flux cleaned = 0.00325364 Jy ! Total flux subtracted in 200 components = 0.00325364 Jy ! Clean residual min=-0.003021 max=0.003092 Jy/beam ! Clean residual mean=0.000422 rms=0.000679 Jy/beam ! Combined flux in latest and established models = 1.1109 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 403 components and 1.1109 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Added new window around map position (21.2, 19.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0012057 Jy ! Component: 100 - total flux cleaned = 0.00146338 Jy ! Component: 150 - total flux cleaned = 0.00179403 Jy ! Component: 200 - total flux cleaned = 0.00203515 Jy ! Total flux subtracted in 200 components = 0.00203515 Jy ! Clean residual min=-0.003015 max=0.002966 Jy/beam ! Clean residual mean=0.000373 rms=0.000666 Jy/beam ! Combined flux in latest and established models = 1.11293 Jy ! Adding 59 model components to the UV plane model. ! The established model now contains 414 components and 1.11293 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 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 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=2.502 mas, bmaj=4.212 mas, bpa=-7.766 degrees ! Estimated noise=1.05939 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00277368 Jy ! Component: 100 - total flux cleaned = 0.0049798 Jy ! Component: 150 - total flux cleaned = 0.00674746 Jy ! Component: 200 - total flux cleaned = 0.00764991 Jy ! Total flux subtracted in 200 components = 0.00764991 Jy ! Clean residual min=-0.002907 max=0.003031 Jy/beam ! Clean residual mean=0.000400 rms=0.000754 Jy/beam ! Combined flux in latest and established models = 1.12058 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.000415735 Jy ! Component: 100 - total flux cleaned = 0.000754043 Jy ! Component: 150 - total flux cleaned = 0.00107643 Jy ! Component: 200 - total flux cleaned = 0.00128233 Jy ! Total flux subtracted in 200 components = 0.00128233 Jy ! Clean residual min=-0.002902 max=0.002942 Jy/beam ! Clean residual mean=0.000326 rms=0.000730 Jy/beam ! Combined flux in latest and established models = 1.12187 Jy ! Adding 54 model components to the UV plane model. ! The established model now contains 458 components and 1.12187 Jy ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! Estimated beam: bmin=0.6203 mas, bmaj=0.9351 mas, bpa=43.51 degrees ! Estimated noise=0.83976 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00221246 Jy ! Component: 100 - total flux cleaned = -0.00309808 Jy ! Component: 150 - total flux cleaned = -0.00348229 Jy ! Component: 200 - total flux cleaned = -0.00363183 Jy ! Total flux subtracted in 200 components = -0.00363183 Jy ! Clean residual min=-0.002967 max=0.002761 Jy/beam ! Clean residual mean=0.000246 rms=0.000631 Jy/beam ! Combined flux in latest and established models = 1.11823 Jy selfcal ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 473 components and 1.11823 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.116284Jy sigma=0.780940 ! Fit after self-cal, rms=0.116145Jy sigma=0.776901 wmodel J0405-1308_X_map.mod ! Writing 473 model components to file: J0405-1308_X_map.mod wobs J0405-1308_X_uvs.fits ! Writing UV FITS file: J0405-1308_X_uvs.fits wwins J0405-1308_X_map.win ! wwins: Wrote 48 windows to J0405-1308_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 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000967075 Jy ! Component: 100 - total flux cleaned = 0.00159435 Jy ! Component: 150 - total flux cleaned = 0.00220318 Jy ! Component: 200 - total flux cleaned = 0.00274985 Jy ! Component: 250 - total flux cleaned = 0.0033703 Jy ! Component: 300 - total flux cleaned = 0.00380314 Jy ! Component: 350 - total flux cleaned = 0.00410047 Jy ! Component: 400 - total flux cleaned = 0.00435092 Jy ! Component: 450 - total flux cleaned = 0.00467918 Jy ! Component: 500 - total flux cleaned = 0.00488136 Jy ! Component: 550 - total flux cleaned = 0.00516028 Jy ! Component: 600 - total flux cleaned = 0.00539595 Jy ! Component: 650 - total flux cleaned = 0.00558993 Jy ! Component: 700 - total flux cleaned = 0.00570535 Jy ! Component: 750 - total flux cleaned = 0.00581967 Jy ! Component: 800 - total flux cleaned = 0.00600827 Jy ! Component: 850 - total flux cleaned = 0.00612027 Jy ! Component: 900 - total flux cleaned = 0.00623145 Jy ! Component: 950 - total flux cleaned = 0.00634172 Jy ! Component: 1000 - total flux cleaned = 0.0064509 Jy ! Component: 1050 - total flux cleaned = 0.00648709 Jy ! Component: 1100 - total flux cleaned = 0.00655911 Jy ! Component: 1150 - total flux cleaned = 0.00666599 Jy ! Component: 1200 - total flux cleaned = 0.00673661 Jy ! Component: 1250 - total flux cleaned = 0.00677181 Jy ! Component: 1300 - total flux cleaned = 0.00680676 Jy ! Component: 1350 - total flux cleaned = 0.0068762 Jy ! Component: 1400 - total flux cleaned = 0.00691061 Jy ! Component: 1450 - total flux cleaned = 0.00694481 Jy ! Component: 1500 - total flux cleaned = 0.00697886 Jy ! Component: 1550 - total flux cleaned = 0.0070127 Jy ! Component: 1600 - total flux cleaned = 0.00704629 Jy ! Component: 1650 - total flux cleaned = 0.00707974 Jy ! Component: 1700 - total flux cleaned = 0.00711293 Jy ! Component: 1750 - total flux cleaned = 0.00711287 Jy ! Component: 1800 - total flux cleaned = 0.00721149 Jy ! Component: 1850 - total flux cleaned = 0.00721127 Jy ! Component: 1900 - total flux cleaned = 0.0071787 Jy ! Component: 1950 - total flux cleaned = 0.0072436 Jy ! Component: 2000 - total flux cleaned = 0.00724368 Jy ! Component: 2050 - total flux cleaned = 0.00727576 Jy ! Component: 2100 - total flux cleaned = 0.00733956 Jy ! Component: 2150 - total flux cleaned = 0.00733955 Jy ! Component: 2200 - total flux cleaned = 0.00737116 Jy ! Component: 2250 - total flux cleaned = 0.00740264 Jy ! Component: 2300 - total flux cleaned = 0.00743397 Jy ! Component: 2350 - total flux cleaned = 0.00743394 Jy ! Component: 2400 - total flux cleaned = 0.00749599 Jy ! Component: 2450 - total flux cleaned = 0.00746503 Jy ! Component: 2500 - total flux cleaned = 0.0075266 Jy ! Component: 2550 - total flux cleaned = 0.00752667 Jy ! Component: 2600 - total flux cleaned = 0.00755713 Jy ! Component: 2650 - total flux cleaned = 0.00755708 Jy ! Component: 2700 - total flux cleaned = 0.00758737 Jy ! Component: 2750 - total flux cleaned = 0.00758734 Jy ! Component: 2800 - total flux cleaned = 0.00761742 Jy ! Component: 2850 - total flux cleaned = 0.00764727 Jy ! Component: 2900 - total flux cleaned = 0.00764729 Jy ! Component: 2950 - total flux cleaned = 0.00767696 Jy ! Component: 3000 - total flux cleaned = 0.00767693 Jy ! Component: 3050 - total flux cleaned = 0.00770639 Jy ! Component: 3100 - total flux cleaned = 0.00770638 Jy ! Component: 3150 - total flux cleaned = 0.00773565 Jy ! Component: 3200 - total flux cleaned = 0.0077356 Jy ! Component: 3250 - total flux cleaned = 0.00773565 Jy ! Component: 3300 - total flux cleaned = 0.00776458 Jy ! Component: 3350 - total flux cleaned = 0.00779335 Jy ! Component: 3400 - total flux cleaned = 0.00779337 Jy ! Component: 3450 - total flux cleaned = 0.00782195 Jy ! Component: 3500 - total flux cleaned = 0.00785045 Jy ! Component: 3550 - total flux cleaned = 0.00785047 Jy ! Component: 3600 - total flux cleaned = 0.00785051 Jy ! Component: 3650 - total flux cleaned = 0.00790686 Jy ! Component: 3700 - total flux cleaned = 0.00787884 Jy ! Component: 3750 - total flux cleaned = 0.00787886 Jy ! Component: 3800 - total flux cleaned = 0.00787898 Jy ! Component: 3850 - total flux cleaned = 0.00793465 Jy ! Component: 3900 - total flux cleaned = 0.00793459 Jy ! Component: 3950 - total flux cleaned = 0.00793465 Jy ! Component: 4000 - total flux cleaned = 0.00796219 Jy ! Component: 4050 - total flux cleaned = 0.00798965 Jy ! Total flux subtracted in 4096 components = 0.00798965 Jy ! Clean residual min=-0.002026 max=0.001633 Jy/beam ! Clean residual mean=0.000038 rms=0.000452 Jy/beam ! Combined flux in latest and established models = 1.12622 Jy keep ! Adding 352 model components to the UV plane model. ! The established model now contains 824 components and 1.12622 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv123_uvs/J0405-1308_X_map.fits ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6203 x 0.9351 at 43.51 degrees (North through East) ! Clean map min=-0.0086698 max=0.68118 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv123_uvs/J0405-1308_X_map.fits wmodel /scr/pima/rv123_uvs/J0405-1308_X_map.mod ! Writing 824 model components to file: /scr/pima/rv123_uvs/J0405-1308_X_map.mod wobs /scr/pima/rv123_uvs/J0405-1308_X_uvs.fits ! Writing UV FITS file: /scr/pima/rv123_uvs/J0405-1308_X_uvs.fits wwins /scr/pima/rv123_uvs/J0405-1308_X_map.win ! wwins: Wrote 49 windows to /scr/pima/rv123_uvs/J0405-1308_X_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.681184 invert ! Inverting map ! Your choice of large map pixels excluded 3.39% of the data. ! The x-axis pixel size should ideally be below 0.1677 milli-arcsec ! The y-axis pixel size should ideally be below 0.1975 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000449053 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.197767 0.197767 0.395534 0.791069 1.58214 3.16428 6.32855 12.6571 25.3142 50.6284 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0013 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.681 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1516.93 quit ! Quitting program ! Log file /scr/pima/rv123_uvs/J0405-1308_X_dfm.log closed on Sun Dec 23 00:40:01 2018