! Started logfile: /scr/pima/rv124_uvs/J2358-1020_X_dfm.log on Mon Dec 24 00:22:19 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/rv124_uvs/J2358-1020_X_uva.fits ! Reading UV FITS file: /scr/pima/rv124_uvs/J2358-1020_X_uva.fits ! AN table 1: 113 integrations on 6 of 6 possible baselines. ! AN table 2: 8 integrations on 15 of 15 possible baselines. ! AN table 3: 11 integrations on 21 of 21 possible baselines. ! AN table 4: 31 integrations on 55 of 55 possible baselines. ! AN table 5: 4 integrations on 21 of 21 possible baselines. ! AN table 6: 41 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.518619 visibilities/baseline/integration-bin. ! Found source: J2358-1020 ! ! 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 79 lines of history. ! ! Reading 7632 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 J2358-1020_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 167 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 214 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 113 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 167 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 214 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 113 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 167 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 214 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 113 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 167 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 13 telescope corrections were flagged in sub-array 3. ! A total of 214 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 113 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.275773Jy sigma=73.292976 ! Fit after self-cal, rms=0.546091Jy sigma=70.726673 ! 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 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.3994 mas, bmaj=0.4226 mas, bpa=1.461 degrees ! Estimated noise=0.509889 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.355236 Jy ! Component: 100 - total flux cleaned = 0.560438 Jy ! Total flux subtracted in 100 components = 0.560438 Jy ! Clean residual min=-0.066128 max=0.181361 Jy/beam ! Clean residual mean=0.000070 rms=0.011794 Jy/beam ! Combined flux in latest and established models = 0.560438 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 0.560438 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.347756Jy sigma=9.877227 ! Fit after self-cal, rms=0.335207Jy sigma=9.780280 ! Inverting map ! Your choice of large map pixels excluded 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (0, -0.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.166242 Jy ! Component: 100 - total flux cleaned = 0.25608 Jy ! Total flux subtracted in 100 components = 0.25608 Jy ! Clean residual min=-0.038003 max=0.077584 Jy/beam ! Clean residual mean=-0.000028 rms=0.005788 Jy/beam ! Combined flux in latest and established models = 0.816518 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 11 components and 0.816518 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.195420Jy sigma=4.715491 ! Fit after self-cal, rms=0.185338Jy sigma=4.375202 ! Inverting map ! Your choice of large map pixels excluded 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (0, 0.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0681351 Jy ! Component: 100 - total flux cleaned = 0.109554 Jy ! Total flux subtracted in 100 components = 0.109554 Jy ! Clean residual min=-0.025306 max=0.027516 Jy/beam ! Clean residual mean=-0.000073 rms=0.003744 Jy/beam ! Combined flux in latest and established models = 0.926072 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 15 components and 0.926072 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166538Jy sigma=3.364613 ! Fit after self-cal, rms=0.164904Jy sigma=3.170094 ! Inverting map ! Your choice of large map pixels excluded 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0297825 Jy ! Component: 100 - total flux cleaned = 0.0484337 Jy ! Total flux subtracted in 100 components = 0.0484337 Jy ! Clean residual min=-0.019270 max=0.018934 Jy/beam ! Clean residual mean=-0.000094 rms=0.003056 Jy/beam ! Combined flux in latest and established models = 0.974505 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 16 components and 0.974505 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168754Jy sigma=2.996590 ! Fit after self-cal, rms=0.168368Jy sigma=2.953288 ! Inverting map ! Your choice of large map pixels excluded 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 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 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.3188 mas, bmaj=0.5573 mas, bpa=45.05 degrees ! Estimated noise=0.281239 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.00653931 Jy ! Component: 100 - total flux cleaned = -0.00824879 Jy ! Component: 150 - total flux cleaned = -0.0087656 Jy ! Component: 200 - total flux cleaned = -0.00925765 Jy ! Total flux subtracted in 200 components = -0.00925765 Jy ! Clean residual min=-0.021862 max=0.015756 Jy/beam ! Clean residual mean=-0.006408 rms=0.002971 Jy/beam ! Combined flux in latest and established models = 0.965248 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 18 components and 0.965248 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166877Jy sigma=3.019920 ! Fit after self-cal, rms=0.166907Jy sigma=3.013386 ! Inverting map ! Your choice of large map pixels excluded 6.88% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.3188 x 0.5573 at 45.05 degrees (North through East) ! Clean map min=-0.016425 max=0.33613 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 285 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* WN 0.00* WZ 0.00* YG 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.94 HH 1.00* HN 1.05 LA 1.16 ! NL 1.05 SC 0.89 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.90 HN 1.01 KP 1.11 LA 1.20 ! NL 1.05 PT 0.97 SC 0.91 ! ! A total of 27 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.05 FD 0.88 HH 1.00* HN 1.01 ! KP 1.06 NL 1.04 OV 0.97 PT 0.89 ! SC 0.91 WN 1.71 WZ 1.14 ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 0.90 KP 1.04 MK 0.89 ! NL 1.06 OV 0.96 PT 0.88 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 0.90 KP 1.07 MK 0.92 ! OV 0.95 PT 0.91 YG 4.43 ! ! ! Correcting IF 2. ! A total of 285 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* WN 0.00* WZ 0.00* YG 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.94 HH 1.00* HN 1.05 LA 1.21 ! NL 1.06 SC 0.90 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.92 HN 1.02 KP 1.09 LA 1.20 ! NL 1.03 PT 0.97 SC 0.89 ! ! A total of 27 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.05 FD 0.90 HH 1.00* HN 1.02 ! KP 1.06 NL 1.03 OV 0.97 PT 0.89 ! SC 0.89 WN 1.75 WZ 1.09 ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 0.89 KP 1.03 MK 0.93 ! NL 1.05 OV 0.95 PT 0.90 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 0.91 KP 1.05 MK 0.92 ! OV 0.95 PT 0.91 YG 4.21 ! ! ! Correcting IF 3. ! A total of 285 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* WN 0.00* WZ 0.00* YG 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.93 HH 1.00* HN 1.04 LA 1.19 ! NL 1.06 SC 0.89 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.92 HN 1.02 KP 1.07 LA 1.21 ! NL 1.03 PT 1.04 SC 0.90 ! ! A total of 27 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.08 FD 0.89 HH 1.00* HN 1.02 ! KP 1.04 NL 1.04 OV 0.97 PT 0.92 ! SC 0.89 WN 1.63 WZ 1.18 ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.99 FD 0.89 KP 1.03 MK 0.90 ! NL 1.06 OV 0.97 PT 0.93 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.00 FD 0.92 KP 1.05 MK 0.93 ! OV 0.95 PT 0.95 YG 2.96 ! ! ! Correcting IF 4. ! A total of 285 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! HH 0.00* WN 0.00* WZ 0.00* YG 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.96 HH 1.00* HN 1.01 LA 1.20 ! NL 1.05 SC 0.91 ! ! Telescope amplitude corrections in sub-array 3: ! FD 0.92 HN 1.01 KP 1.08 LA 1.21 ! NL 1.04 PT 1.02 SC 0.92 ! ! A total of 27 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.08 FD 0.89 HH 1.00* HN 1.02 ! KP 1.06 NL 1.04 OV 0.97 PT 0.92 ! SC 0.92 WN 1.51 WZ 1.22 ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.99 FD 0.90 KP 1.05 MK 0.88 ! NL 1.08 OV 0.97 PT 0.93 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.02 FD 0.92 KP 1.07 MK 0.91 ! OV 0.96 PT 0.93 YG 3.34 ! ! ! Fit before self-cal, rms=0.166907Jy sigma=3.013386 ! Fit after self-cal, rms=0.085098Jy sigma=1.255421 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.5612 mas, bmaj=1.655 mas, bpa=-35.21 degrees ! Estimated noise=0.760406 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.000765063 Jy ! Component: 100 - total flux cleaned = -0.000809061 Jy ! Component: 150 - total flux cleaned = -0.000844533 Jy ! Component: 200 - total flux cleaned = -0.000555909 Jy ! Total flux subtracted in 200 components = -0.000555909 Jy ! Clean residual min=-0.012550 max=0.013631 Jy/beam ! Clean residual mean=0.000265 rms=0.002179 Jy/beam ! Combined flux in latest and established models = 0.964692 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 20 components and 0.964692 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.082789Jy sigma=1.145360 ! Fit after self-cal, rms=0.082151Jy sigma=1.131806 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.082151Jy sigma=1.131806 ! Fit after self-cal, rms=0.081376Jy sigma=1.093264 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.5991 mas, bmaj=1.651 mas, bpa=-33.76 degrees ! Estimated noise=0.782965 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000855963 Jy ! Component: 100 - total flux cleaned = 0.000894835 Jy ! Component: 150 - total flux cleaned = 0.000636201 Jy ! Component: 200 - total flux cleaned = 0.000626615 Jy ! Total flux subtracted in 200 components = 0.000626615 Jy ! Clean residual min=-0.011209 max=0.011941 Jy/beam ! Clean residual mean=0.000250 rms=0.002017 Jy/beam ! Combined flux in latest and established models = 0.965318 Jy selfcal ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 20 components and 0.965318 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.080385Jy sigma=1.070570 ! Fit after self-cal, rms=0.080315Jy sigma=1.068957 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.00126367 Jy ! Component: 100 - total flux cleaned = -0.00105828 Jy ! Component: 150 - total flux cleaned = -0.000876972 Jy ! Component: 200 - total flux cleaned = -0.000543299 Jy ! Total flux subtracted in 200 components = -0.000543299 Jy ! Clean residual min=-0.010496 max=0.011488 Jy/beam ! Clean residual mean=0.000252 rms=0.001936 Jy/beam ! Combined flux in latest and established models = 0.964775 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 21 components and 0.964775 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079787Jy sigma=1.078223 ! Fit after self-cal, rms=0.079720Jy sigma=1.076753 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.079720Jy sigma=1.076753 ! Fit after self-cal, rms=0.077552Jy sigma=0.990026 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.647 mas, bmaj=1.658 mas, bpa=-31.75 degrees ! Estimated noise=0.80609 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000630678 Jy ! Component: 100 - total flux cleaned = 0.000880061 Jy ! Component: 150 - total flux cleaned = 0.000696906 Jy ! Component: 200 - total flux cleaned = 0.00033117 Jy ! Total flux subtracted in 200 components = 0.00033117 Jy ! Clean residual min=-0.009873 max=0.010930 Jy/beam ! Clean residual mean=0.000277 rms=0.001625 Jy/beam ! Combined flux in latest and established models = 0.965106 Jy selfcal ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 21 components and 0.965106 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.077117Jy sigma=0.986412 ! Fit after self-cal, rms=0.077075Jy sigma=0.985367 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.000713098 Jy ! Component: 100 - total flux cleaned = -0.000705526 Jy ! Component: 150 - total flux cleaned = -0.000703953 Jy ! Component: 200 - total flux cleaned = -0.000576333 Jy ! Total flux subtracted in 200 components = -0.000576333 Jy ! Clean residual min=-0.008927 max=0.010276 Jy/beam ! Clean residual mean=0.000279 rms=0.001568 Jy/beam ! Combined flux in latest and established models = 0.96453 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 21 components and 0.96453 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.076851Jy sigma=0.993483 ! Fit after self-cal, rms=0.076804Jy sigma=0.992255 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.076804Jy sigma=0.992255 ! Fit after self-cal, rms=0.076880Jy sigma=0.957303 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=2.176 mas, bmaj=3.068 mas, bpa=8.547 degrees ! Estimated noise=1.04134 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00372234 Jy ! Component: 100 - total flux cleaned = 0.00372538 Jy ! Component: 150 - total flux cleaned = 0.00372442 Jy ! Component: 200 - total flux cleaned = 0.00372291 Jy ! Total flux subtracted in 200 components = 0.00372291 Jy ! Clean residual min=-0.004981 max=0.008230 Jy/beam ! Clean residual mean=0.000543 rms=0.001387 Jy/beam ! Combined flux in latest and established models = 0.968253 Jy selfcal ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 21 components and 0.968253 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.076805Jy sigma=0.958458 ! Fit after self-cal, rms=0.076691Jy sigma=0.956925 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.2, -3.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00618112 Jy ! Component: 100 - total flux cleaned = 0.00761364 Jy ! Component: 150 - total flux cleaned = 0.00845298 Jy ! Component: 200 - total flux cleaned = 0.00869083 Jy ! Total flux subtracted in 200 components = 0.00869083 Jy ! Clean residual min=-0.003928 max=0.005725 Jy/beam ! Clean residual mean=0.000483 rms=0.001180 Jy/beam ! Combined flux in latest and established models = 0.976944 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 44 components and 0.976944 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.075358Jy sigma=0.932924 ! Fit after self-cal, rms=0.075303Jy sigma=0.931967 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (-18.6, -38). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00291887 Jy ! Component: 100 - total flux cleaned = 0.00422014 Jy ! Component: 150 - total flux cleaned = 0.00502065 Jy ! Component: 200 - total flux cleaned = 0.00529322 Jy ! Total flux subtracted in 200 components = 0.00529322 Jy ! Clean residual min=-0.003639 max=0.004812 Jy/beam ! Clean residual mean=0.000458 rms=0.001067 Jy/beam ! Combined flux in latest and established models = 0.982237 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 71 components and 0.982237 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074809Jy sigma=0.923856 ! Fit after self-cal, rms=0.074766Jy sigma=0.923269 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (-9, -21.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00467919 Jy ! Component: 100 - total flux cleaned = 0.00585157 Jy ! Component: 150 - total flux cleaned = 0.00651246 Jy ! Component: 200 - total flux cleaned = 0.00731915 Jy ! Total flux subtracted in 200 components = 0.00731915 Jy ! Clean residual min=-0.003383 max=0.004493 Jy/beam ! Clean residual mean=0.000426 rms=0.000957 Jy/beam ! Combined flux in latest and established models = 0.989556 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 98 components and 0.989556 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074349Jy sigma=0.918806 ! Fit after self-cal, rms=0.074310Jy sigma=0.918484 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (7, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00124206 Jy ! Component: 100 - total flux cleaned = 0.0023121 Jy ! Component: 150 - total flux cleaned = 0.00303131 Jy ! Component: 200 - total flux cleaned = 0.00368223 Jy ! Total flux subtracted in 200 components = 0.00368223 Jy ! Clean residual min=-0.003215 max=0.003904 Jy/beam ! Clean residual mean=0.000411 rms=0.000828 Jy/beam ! Combined flux in latest and established models = 0.993238 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 125 components and 0.993238 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074068Jy sigma=0.915916 ! Fit after self-cal, rms=0.074040Jy sigma=0.915657 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (15.6, 1.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0032181 Jy ! Component: 100 - total flux cleaned = 0.00427213 Jy ! Component: 150 - total flux cleaned = 0.00505018 Jy ! Component: 200 - total flux cleaned = 0.00569337 Jy ! Total flux subtracted in 200 components = 0.00569337 Jy ! Clean residual min=-0.003100 max=0.003107 Jy/beam ! Clean residual mean=0.000382 rms=0.000753 Jy/beam ! Combined flux in latest and established models = 0.998932 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 156 components and 0.998932 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.073774Jy sigma=0.912920 ! Fit after self-cal, rms=0.073735Jy sigma=0.912718 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.6857 mas, bmaj=1.671 mas, bpa=-29.84 degrees ! Estimated noise=0.821575 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6857 x 1.671 at -29.84 degrees (North through East) ! Clean map min=-0.010712 max=0.73783 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.073735Jy sigma=0.912718 ! Fit after self-cal, rms=0.067961Jy sigma=0.815810 ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.6872 mas, bmaj=1.673 mas, bpa=-29.83 degrees ! Estimated noise=0.823233 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00234666 Jy ! Component: 100 - total flux cleaned = 0.00431064 Jy ! Component: 150 - total flux cleaned = 0.00505117 Jy ! Component: 200 - total flux cleaned = 0.00594291 Jy ! Total flux subtracted in 200 components = 0.00594291 Jy ! Clean residual min=-0.005027 max=0.005864 Jy/beam ! Clean residual mean=0.000208 rms=0.000993 Jy/beam ! Combined flux in latest and established models = 1.00487 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 176 components and 1.00487 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065553Jy sigma=0.784188 ! Fit after self-cal, rms=0.065429Jy sigma=0.780477 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (-3.4, 22.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00242621 Jy ! Component: 100 - total flux cleaned = 0.00342338 Jy ! Component: 150 - total flux cleaned = 0.00414789 Jy ! Component: 200 - total flux cleaned = 0.00447984 Jy ! Total flux subtracted in 200 components = 0.00447984 Jy ! Clean residual min=-0.004113 max=0.004311 Jy/beam ! Clean residual mean=0.000197 rms=0.000864 Jy/beam ! Combined flux in latest and established models = 1.00935 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 196 components and 1.00935 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.064176Jy sigma=0.765370 ! Fit after self-cal, rms=0.064124Jy sigma=0.763829 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (-3.2, 18.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00150321 Jy ! Component: 100 - total flux cleaned = 0.00275948 Jy ! Component: 150 - total flux cleaned = 0.00378657 Jy ! Component: 200 - total flux cleaned = 0.004748 Jy ! Total flux subtracted in 200 components = 0.004748 Jy ! Clean residual min=-0.003770 max=0.003823 Jy/beam ! Clean residual mean=0.000186 rms=0.000779 Jy/beam ! Combined flux in latest and established models = 1.0141 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 211 components and 1.0141 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063325Jy sigma=0.755107 ! Fit after self-cal, rms=0.063297Jy sigma=0.754499 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Added new window around map position (18.6, 12). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00228127 Jy ! Component: 100 - total flux cleaned = 0.0034559 Jy ! Component: 150 - total flux cleaned = 0.0040685 Jy ! Component: 200 - total flux cleaned = 0.00476611 Jy ! Total flux subtracted in 200 components = 0.00476611 Jy ! Clean residual min=-0.003341 max=0.003223 Jy/beam ! Clean residual mean=0.000177 rms=0.000712 Jy/beam ! Combined flux in latest and established models = 1.01887 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 228 components and 1.01887 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.062729Jy sigma=0.748782 ! Fit after self-cal, rms=0.062710Jy sigma=0.748427 ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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.062710Jy sigma=0.748427 ! Fit after self-cal, rms=0.061825Jy sigma=0.743955 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 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.6015 mas, bmaj=1.515 mas, bpa=-29.23 degrees ! Estimated noise=1.23165 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.646361 Jy ! Component: 100 - total flux cleaned = 0.833391 Jy ! Component: 150 - total flux cleaned = 0.903666 Jy ! Component: 200 - total flux cleaned = 0.939654 Jy ! Total flux subtracted in 200 components = 0.939654 Jy ! Clean residual min=-0.011665 max=0.017605 Jy/beam ! Clean residual mean=0.000042 rms=0.001971 Jy/beam ! Combined flux in latest and established models = 0.939654 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 19 model components to the UV plane model. ! The established model now contains 19 components and 0.939654 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.6929 mas, bmaj=1.671 mas, bpa=-29.36 degrees ! Estimated noise=0.82548 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.0227568 Jy ! Component: 100 - total flux cleaned = 0.0357785 Jy ! Component: 150 - total flux cleaned = 0.041783 Jy ! Component: 200 - total flux cleaned = 0.0465269 Jy ! Total flux subtracted in 200 components = 0.0465269 Jy ! Clean residual min=-0.008781 max=0.005950 Jy/beam ! Clean residual mean=0.000186 rms=0.001283 Jy/beam ! Combined flux in latest and established models = 0.986181 Jy ! Adding 27 model components to the UV plane model. ! The established model now contains 40 components and 0.986181 Jy ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 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 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=2.174 mas, bmaj=3.067 mas, bpa=8.665 degrees ! Estimated noise=1.04432 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00480518 Jy ! Component: 100 - total flux cleaned = 0.00851212 Jy ! Component: 150 - total flux cleaned = 0.0117885 Jy ! Component: 200 - total flux cleaned = 0.0145064 Jy ! Total flux subtracted in 200 components = 0.0145064 Jy ! Clean residual min=-0.003537 max=0.003217 Jy/beam ! Clean residual mean=0.000323 rms=0.000877 Jy/beam ! Combined flux in latest and established models = 1.00069 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.00203573 Jy ! Component: 100 - total flux cleaned = 0.00391003 Jy ! Component: 150 - total flux cleaned = 0.0053183 Jy ! Component: 200 - total flux cleaned = 0.00664927 Jy ! Total flux subtracted in 200 components = 0.00664927 Jy ! Clean residual min=-0.003089 max=0.002854 Jy/beam ! Clean residual mean=0.000289 rms=0.000771 Jy/beam ! Combined flux in latest and established models = 1.00734 Jy ! Adding 58 model components to the UV plane model. ! The established model now contains 95 components and 1.00734 Jy ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! Estimated beam: bmin=0.6929 mas, bmaj=1.671 mas, bpa=-29.36 degrees ! Estimated noise=0.82548 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000101527 Jy ! Component: 100 - total flux cleaned = 0.000638562 Jy ! Component: 150 - total flux cleaned = 0.00132414 Jy ! Component: 200 - total flux cleaned = 0.00216601 Jy ! Total flux subtracted in 200 components = 0.00216601 Jy ! Clean residual min=-0.006496 max=0.004670 Jy/beam ! Clean residual mean=0.000129 rms=0.000940 Jy/beam ! Combined flux in latest and established models = 1.0095 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 112 components and 1.0095 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065782Jy sigma=0.858164 ! Fit after self-cal, rms=0.065167Jy sigma=0.841269 wmodel J2358-1020_X_map.mod ! Writing 112 model components to file: J2358-1020_X_map.mod wobs J2358-1020_X_uvs.fits ! Writing UV FITS file: J2358-1020_X_uvs.fits wwins J2358-1020_X_map.win ! wwins: Wrote 11 windows to J2358-1020_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 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00109716 Jy ! Component: 100 - total flux cleaned = -0.00128542 Jy ! Component: 150 - total flux cleaned = -0.00112939 Jy ! Component: 200 - total flux cleaned = -0.000833019 Jy ! Component: 250 - total flux cleaned = -0.000550341 Jy ! Component: 300 - total flux cleaned = 5.9112e-05 Jy ! Component: 350 - total flux cleaned = 0.000192702 Jy ! Component: 400 - total flux cleaned = 0.000896495 Jy ! Component: 450 - total flux cleaned = 0.00133221 Jy ! Component: 500 - total flux cleaned = 0.00194149 Jy ! Component: 550 - total flux cleaned = 0.00247653 Jy ! Component: 600 - total flux cleaned = 0.00305866 Jy ! Component: 650 - total flux cleaned = 0.00340211 Jy ! Component: 700 - total flux cleaned = 0.00402057 Jy ! Component: 750 - total flux cleaned = 0.00462633 Jy ! Component: 800 - total flux cleaned = 0.00500415 Jy ! Component: 850 - total flux cleaned = 0.00564191 Jy ! Component: 900 - total flux cleaned = 0.00605971 Jy ! Component: 950 - total flux cleaned = 0.00652252 Jy ! Component: 1000 - total flux cleaned = 0.00707914 Jy ! Component: 1050 - total flux cleaned = 0.00757708 Jy ! Component: 1100 - total flux cleaned = 0.00801654 Jy ! Component: 1150 - total flux cleaned = 0.00859454 Jy ! Component: 1200 - total flux cleaned = 0.00897471 Jy ! Component: 1250 - total flux cleaned = 0.00958225 Jy ! Component: 1300 - total flux cleaned = 0.0100875 Jy ! Component: 1350 - total flux cleaned = 0.0105863 Jy ! Component: 1400 - total flux cleaned = 0.0110331 Jy ! Component: 1450 - total flux cleaned = 0.0115623 Jy ! Component: 1500 - total flux cleaned = 0.0119974 Jy ! Component: 1550 - total flux cleaned = 0.0125553 Jy ! Component: 1600 - total flux cleaned = 0.0129363 Jy ! Component: 1650 - total flux cleaned = 0.0133957 Jy ! Component: 1700 - total flux cleaned = 0.0140133 Jy ! Component: 1750 - total flux cleaned = 0.0143793 Jy ! Component: 1800 - total flux cleaned = 0.0148203 Jy ! Component: 1850 - total flux cleaned = 0.0152959 Jy ! Component: 1900 - total flux cleaned = 0.0155309 Jy ! Component: 1950 - total flux cleaned = 0.0159188 Jy ! Component: 2000 - total flux cleaned = 0.0162254 Jy ! Component: 2050 - total flux cleaned = 0.0167551 Jy ! Component: 2100 - total flux cleaned = 0.017129 Jy ! Component: 2150 - total flux cleaned = 0.0174254 Jy ! Component: 2200 - total flux cleaned = 0.0178649 Jy ! Component: 2250 - total flux cleaned = 0.0182271 Jy ! Component: 2300 - total flux cleaned = 0.0186929 Jy ! Component: 2350 - total flux cleaned = 0.0190123 Jy ! Component: 2400 - total flux cleaned = 0.0194687 Jy ! Component: 2450 - total flux cleaned = 0.019712 Jy ! Component: 2500 - total flux cleaned = 0.0201592 Jy ! Component: 2550 - total flux cleaned = 0.0204663 Jy ! Component: 2600 - total flux cleaned = 0.0208374 Jy ! Component: 2650 - total flux cleaned = 0.0211047 Jy ! Component: 2700 - total flux cleaned = 0.0215352 Jy ! Component: 2750 - total flux cleaned = 0.0217646 Jy ! Component: 2800 - total flux cleaned = 0.0220895 Jy ! Component: 2850 - total flux cleaned = 0.022507 Jy ! Component: 2900 - total flux cleaned = 0.0227937 Jy ! Component: 2950 - total flux cleaned = 0.0230777 Jy ! Component: 3000 - total flux cleaned = 0.0234531 Jy ! Component: 3050 - total flux cleaned = 0.0236704 Jy ! Component: 3100 - total flux cleaned = 0.0239164 Jy ! Component: 3150 - total flux cleaned = 0.0242518 Jy ! Component: 3200 - total flux cleaned = 0.0245236 Jy ! Component: 3250 - total flux cleaned = 0.0247334 Jy ! Component: 3300 - total flux cleaned = 0.0249116 Jy ! Component: 3350 - total flux cleaned = 0.025177 Jy ! Component: 3400 - total flux cleaned = 0.0254988 Jy ! Component: 3450 - total flux cleaned = 0.0257307 Jy ! Component: 3500 - total flux cleaned = 0.0259899 Jy ! Component: 3550 - total flux cleaned = 0.0262754 Jy ! Component: 3600 - total flux cleaned = 0.026587 Jy ! Component: 3650 - total flux cleaned = 0.0268397 Jy ! Component: 3700 - total flux cleaned = 0.0270068 Jy ! Component: 3750 - total flux cleaned = 0.0272559 Jy ! Component: 3800 - total flux cleaned = 0.0274206 Jy ! Component: 3850 - total flux cleaned = 0.0276387 Jy ! Component: 3900 - total flux cleaned = 0.0278822 Jy ! Component: 3950 - total flux cleaned = 0.0280436 Jy ! Component: 4000 - total flux cleaned = 0.0283369 Jy ! Component: 4050 - total flux cleaned = 0.0284693 Jy ! Total flux subtracted in 4096 components = 0.0286534 Jy ! Clean residual min=-0.001702 max=0.001513 Jy/beam ! Clean residual mean=0.000066 rms=0.000433 Jy/beam ! Combined flux in latest and established models = 1.03816 Jy keep ! Adding 297 model components to the UV plane model. ! The established model now contains 401 components and 1.03816 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv124_uvs/J2358-1020_X_map.fits ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.6929 x 1.671 at -29.36 degrees (North through East) ! Clean map min=-0.0065245 max=0.76688 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv124_uvs/J2358-1020_X_map.fits wmodel /scr/pima/rv124_uvs/J2358-1020_X_map.mod ! Writing 401 model components to file: /scr/pima/rv124_uvs/J2358-1020_X_map.mod wobs /scr/pima/rv124_uvs/J2358-1020_X_uvs.fits ! Writing UV FITS file: /scr/pima/rv124_uvs/J2358-1020_X_uvs.fits wwins /scr/pima/rv124_uvs/J2358-1020_X_map.win ! wwins: Wrote 12 windows to /scr/pima/rv124_uvs/J2358-1020_X_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 0.766876 invert ! Inverting map ! Your choice of large map pixels excluded 8.82% of the data. ! The x-axis pixel size should ideally be below 0.1648 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00043155 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.168821 0.168821 0.337642 0.675285 1.35057 2.70114 5.40228 10.8046 21.6091 43.2182 86.4364 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0012 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.766 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1777.03 quit ! Quitting program ! Log file /scr/pima/rv124_uvs/J2358-1020_X_dfm.log closed on Mon Dec 24 00:22:24 2018