! Started logfile: /scr/pima/rv123_uvs/J0006-0623_X_dfm.log on Sun Dec 23 00:38:46 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/J0006-0623_X_uva.fits ! Reading UV FITS file: /scr/pima/rv123_uvs/J0006-0623_X_uva.fits ! AN table 1: 51 integrations on 78 of 78 possible baselines. ! AN table 2: 15 integrations on 45 of 45 possible baselines. ! AN table 3: 14 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.367691 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0006-0623 ! ! 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 7420 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 J0006-0623_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 318 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! A total of 58 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 318 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! A total of 58 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 318 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! A total of 58 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 318 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! A total of 58 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=2.662374Jy sigma=21.605784 ! Fit after self-cal, rms=1.630191Jy sigma=13.199858 ! 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 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5177 mas, bmaj=1.252 mas, bpa=-9.043 degrees ! Estimated noise=1.58386 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 = 1.61629 Jy ! Component: 100 - total flux cleaned = 2.22603 Jy ! Total flux subtracted in 100 components = 2.22603 Jy ! Clean residual min=-0.110142 max=0.309776 Jy/beam ! Clean residual mean=-0.000489 rms=0.024866 Jy/beam ! Combined flux in latest and established models = 2.22603 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 7 components and 2.22603 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.893067Jy sigma=7.360962 ! Fit after self-cal, rms=0.892911Jy sigma=7.360303 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.369 Jy ! Component: 100 - total flux cleaned = 0.606906 Jy ! Total flux subtracted in 100 components = 0.606906 Jy ! Clean residual min=-0.091464 max=0.152626 Jy/beam ! Clean residual mean=-0.000682 rms=0.019621 Jy/beam ! Combined flux in latest and established models = 2.83294 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 15 components and 2.83293 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.680587Jy sigma=5.656047 ! Fit after self-cal, rms=0.680480Jy sigma=5.655729 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (0.4, -1.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.170671 Jy ! Component: 100 - total flux cleaned = 0.284373 Jy ! Total flux subtracted in 100 components = 0.284373 Jy ! Clean residual min=-0.095544 max=0.113195 Jy/beam ! Clean residual mean=-0.000771 rms=0.018157 Jy/beam ! Combined flux in latest and established models = 3.11731 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 20 components and 3.11731 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.639052Jy sigma=5.258569 ! Fit after self-cal, rms=0.638237Jy sigma=5.252544 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-4.4, -2.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0986987 Jy ! Component: 100 - total flux cleaned = 0.151422 Jy ! Total flux subtracted in 100 components = 0.151422 Jy ! Clean residual min=-0.100007 max=0.108401 Jy/beam ! Clean residual mean=-0.000814 rms=0.017354 Jy/beam ! Combined flux in latest and established models = 3.26873 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 26 components and 3.26873 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.624443Jy sigma=5.097281 ! Fit after self-cal, rms=0.623631Jy sigma=5.087144 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-3.6, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.112423 Jy ! Component: 100 - total flux cleaned = 0.170994 Jy ! Total flux subtracted in 100 components = 0.170994 Jy ! Clean residual min=-0.094274 max=0.103501 Jy/beam ! Clean residual mean=-0.000863 rms=0.016533 Jy/beam ! Combined flux in latest and established models = 3.43973 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 35 components and 3.43972 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.602645Jy sigma=4.881003 ! Fit after self-cal, rms=0.601151Jy sigma=4.867849 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-1.6, -0.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0946016 Jy ! Component: 100 - total flux cleaned = 0.15728 Jy ! Total flux subtracted in 100 components = 0.15728 Jy ! Clean residual min=-0.092866 max=0.089214 Jy/beam ! Clean residual mean=-0.000909 rms=0.015841 Jy/beam ! Combined flux in latest and established models = 3.597 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 44 components and 3.597 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.602541Jy sigma=4.793375 ! Fit after self-cal, rms=0.601418Jy sigma=4.783874 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 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 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5384 mas, bmaj=1.5 mas, bpa=-8.481 degrees ! Estimated noise=1.11966 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.0325861 Jy ! Component: 100 - total flux cleaned = -0.0400482 Jy ! Component: 150 - total flux cleaned = -0.0486147 Jy ! Component: 200 - total flux cleaned = -0.059779 Jy ! Total flux subtracted in 200 components = -0.059779 Jy ! Clean residual min=-0.119601 max=0.111365 Jy/beam ! Clean residual mean=-0.028255 rms=0.021038 Jy/beam ! Combined flux in latest and established models = 3.53723 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 54 components and 3.53723 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.581262Jy sigma=4.630610 ! Fit after self-cal, rms=0.580403Jy sigma=4.620538 ! Inverting map ! Your choice of large map pixels excluded 5.02% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.5384 x 1.5 at -8.481 degrees (North through East) ! Clean map min=-0.11631 max=1.8987 Jy/beam gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.88 HH 0.84 HN 1.05 ! KP 1.09 LA 1.04 NL 1.01 ON 0.99 ! OV 1.02 PT 0.86 SC 1.10 WN 1.57 ! WZ 1.01 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 0.87 HN 1.00 KP 1.08 ! LA 0.99 MK 0.96 NL 1.04 OV 0.95 ! PT 0.88 SC 0.86 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 0.90 KP 1.06 LA 1.03 ! MK 0.92 OV 0.98 PT 0.84 YG 1.10 ! ! ! Correcting IF 2. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 0.89 HH 0.83 HN 1.03 ! KP 1.10 LA 1.04 NL 1.01 ON 0.96 ! OV 1.01 PT 0.85 SC 1.17 WN 1.51 ! WZ 0.99 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.84 HN 0.93 KP 1.08 ! LA 1.01 MK 0.97 NL 1.05 OV 0.95 ! PT 0.87 SC 0.90 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 0.90 KP 1.05 LA 1.00 ! MK 0.90 OV 0.96 PT 0.88 YG 0.99 ! ! ! Correcting IF 3. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.89 HH 0.88 HN 1.03 ! KP 1.09 LA 1.04 NL 1.00 ON 1.31 ! OV 1.03 PT 0.86 SC 1.09 WN 1.38 ! WZ 1.07 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 0.83 HN 0.99 KP 1.05 ! LA 1.00 MK 0.98 NL 1.04 OV 0.93 ! PT 0.87 SC 0.95 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 0.91 KP 1.06 LA 1.01 ! MK 0.95 OV 0.96 PT 0.86 YG 1.21 ! ! ! Correcting IF 4. ! A total of 27 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.89 HH 0.99 HN 1.03 ! KP 1.11 LA 1.02 NL 1.01 ON 1.27 ! OV 0.99 PT 0.85 SC 1.06 WN 1.30 ! WZ 1.08 ! ! A total of 24 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 0.85 HN 0.97 KP 1.07 ! LA 0.99 MK 0.97 NL 1.04 OV 0.94 ! PT 0.87 SC 0.89 ! ! A total of 30 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.90 FD 0.92 KP 1.03 LA 1.00 ! MK 0.94 OV 0.98 PT 0.86 YG 1.33 ! ! ! Fit before self-cal, rms=0.580403Jy sigma=4.620538 ! Fit after self-cal, rms=0.433921Jy sigma=3.375722 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.547 mas, bmaj=1.474 mas, bpa=-8.116 degrees ! Estimated noise=1.19274 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.0639941 Jy ! Component: 100 - total flux cleaned = -0.0684256 Jy ! Component: 150 - total flux cleaned = -0.0702893 Jy ! Component: 200 - total flux cleaned = -0.0735421 Jy ! Total flux subtracted in 200 components = -0.0735421 Jy ! Clean residual min=-0.086115 max=0.080994 Jy/beam ! Clean residual mean=-0.006524 rms=0.013586 Jy/beam ! Combined flux in latest and established models = 3.46368 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 64 components and 3.46368 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.408400Jy sigma=3.232343 ! Fit after self-cal, rms=0.407948Jy sigma=3.218312 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-5.4, 1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0485456 Jy ! Component: 100 - total flux cleaned = 0.046248 Jy ! Component: 150 - total flux cleaned = 0.042201 Jy ! Component: 200 - total flux cleaned = 0.0367412 Jy ! Total flux subtracted in 200 components = 0.0367412 Jy ! Clean residual min=-0.069618 max=0.070618 Jy/beam ! Clean residual mean=-0.006847 rms=0.012179 Jy/beam ! Combined flux in latest and established models = 3.50042 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 78 components and 3.50042 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.383856Jy sigma=2.997403 ! Fit after self-cal, rms=0.382001Jy sigma=2.973861 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.382001Jy sigma=2.973861 ! Fit after self-cal, rms=0.376747Jy sigma=2.780129 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5508 mas, bmaj=1.486 mas, bpa=-8.161 degrees ! Estimated noise=1.20555 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0132598 Jy ! Component: 100 - total flux cleaned = -0.0149332 Jy ! Component: 150 - total flux cleaned = -0.0150054 Jy ! Component: 200 - total flux cleaned = -0.0150906 Jy ! Total flux subtracted in 200 components = -0.0150906 Jy ! Clean residual min=-0.070385 max=0.055282 Jy/beam ! Clean residual mean=-0.004783 rms=0.010710 Jy/beam ! Combined flux in latest and established models = 3.48533 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 90 components and 3.48533 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.362027Jy sigma=2.695349 ! Fit after self-cal, rms=0.362222Jy sigma=2.690967 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.00299952 Jy ! Component: 100 - total flux cleaned = 0.00298456 Jy ! Component: 150 - total flux cleaned = 0.00299011 Jy ! Component: 200 - total flux cleaned = 0.0030125 Jy ! Total flux subtracted in 200 components = 0.0030125 Jy ! Clean residual min=-0.061556 max=0.051914 Jy/beam ! Clean residual mean=-0.004808 rms=0.010325 Jy/beam ! Combined flux in latest and established models = 3.48835 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 96 components and 3.48835 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.354009Jy sigma=2.640595 ! Fit after self-cal, rms=0.354169Jy sigma=2.638551 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.354169Jy sigma=2.638551 ! Fit after self-cal, rms=0.253424Jy sigma=1.594572 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.551 mas, bmaj=1.497 mas, bpa=-8.408 degrees ! Estimated noise=1.21368 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0336161 Jy ! Component: 100 - total flux cleaned = -0.0411655 Jy ! Component: 150 - total flux cleaned = -0.0390033 Jy ! Component: 200 - total flux cleaned = -0.0380445 Jy ! Total flux subtracted in 200 components = -0.0380445 Jy ! Clean residual min=-0.036319 max=0.035109 Jy/beam ! Clean residual mean=-0.004249 rms=0.005222 Jy/beam ! Combined flux in latest and established models = 3.4503 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 115 components and 3.4503 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.241168Jy sigma=1.507238 ! Fit after self-cal, rms=0.240793Jy sigma=1.500854 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.00583666 Jy ! Component: 100 - total flux cleaned = -0.00496463 Jy ! Component: 150 - total flux cleaned = -0.00582547 Jy ! Component: 200 - total flux cleaned = -0.00578589 Jy ! Total flux subtracted in 200 components = -0.00578589 Jy ! Clean residual min=-0.031898 max=0.032576 Jy/beam ! Clean residual mean=-0.004203 rms=0.004932 Jy/beam ! Combined flux in latest and established models = 3.44452 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 121 components and 3.44452 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.234711Jy sigma=1.453497 ! Fit after self-cal, rms=0.234362Jy sigma=1.450644 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-0.8, -3). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0161718 Jy ! Component: 100 - total flux cleaned = 0.0119288 Jy ! Component: 150 - total flux cleaned = 0.00960361 Jy ! Component: 200 - total flux cleaned = 0.00747242 Jy ! Total flux subtracted in 200 components = 0.00747242 Jy ! Clean residual min=-0.029432 max=0.019032 Jy/beam ! Clean residual mean=-0.004264 rms=0.004569 Jy/beam ! Combined flux in latest and established models = 3.45199 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 130 components and 3.45199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.227646Jy sigma=1.396010 ! Fit after self-cal, rms=0.227280Jy sigma=1.392656 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.227280Jy sigma=1.392656 ! Fit after self-cal, rms=0.224758Jy sigma=1.363066 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=2.241 mas, bmaj=3.156 mas, bpa=-12.81 degrees ! Estimated noise=1.76548 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0290472 Jy ! Component: 100 - total flux cleaned = -0.0317558 Jy ! Component: 150 - total flux cleaned = -0.0344368 Jy ! Component: 200 - total flux cleaned = -0.0358189 Jy ! Total flux subtracted in 200 components = -0.0358189 Jy ! Clean residual min=-0.029349 max=0.012400 Jy/beam ! Clean residual mean=-0.009542 rms=0.005848 Jy/beam ! Combined flux in latest and established models = 3.41617 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 134 components and 3.41617 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220609Jy sigma=1.338950 ! Fit after self-cal, rms=0.220192Jy sigma=1.331565 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.00295266 Jy ! Component: 100 - total flux cleaned = -0.00431551 Jy ! Component: 150 - total flux cleaned = -0.00513123 Jy ! Component: 200 - total flux cleaned = -0.00584768 Jy ! Total flux subtracted in 200 components = -0.00584768 Jy ! Clean residual min=-0.028915 max=0.012301 Jy/beam ! Clean residual mean=-0.009412 rms=0.005697 Jy/beam ! Combined flux in latest and established models = 3.41032 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 136 components and 3.41032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.218953Jy sigma=1.321707 ! Fit after self-cal, rms=0.218738Jy sigma=1.319649 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5525 mas, bmaj=1.503 mas, bpa=-8.46 degrees ! Estimated noise=1.21396 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5525 x 1.503 at -8.46 degrees (North through East) ! Clean map min=-0.052889 max=1.9765 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.218738Jy sigma=1.319649 ! Fit after self-cal, rms=0.202735Jy sigma=1.190105 ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5523 mas, bmaj=1.509 mas, bpa=-8.531 degrees ! Estimated noise=1.20634 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0242495 Jy ! Component: 100 - total flux cleaned = -0.0296199 Jy ! Component: 150 - total flux cleaned = -0.0315946 Jy ! Component: 200 - total flux cleaned = -0.0333784 Jy ! Total flux subtracted in 200 components = -0.0333784 Jy ! Clean residual min=-0.022521 max=0.014966 Jy/beam ! Clean residual mean=-0.003225 rms=0.003544 Jy/beam ! Combined flux in latest and established models = 3.37694 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 142 components and 3.37694 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.195818Jy sigma=1.145565 ! Fit after self-cal, rms=0.195301Jy sigma=1.142607 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00325067 Jy ! Component: 100 - total flux cleaned = -0.00436461 Jy ! Component: 150 - total flux cleaned = -0.00488132 Jy ! Component: 200 - total flux cleaned = -0.00535106 Jy ! Total flux subtracted in 200 components = -0.00535106 Jy ! Clean residual min=-0.020750 max=0.015811 Jy/beam ! Clean residual mean=-0.003182 rms=0.003433 Jy/beam ! Combined flux in latest and established models = 3.37159 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 146 components and 3.37159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192199Jy sigma=1.119663 ! Fit after self-cal, rms=0.191949Jy sigma=1.118444 ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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.191949Jy sigma=1.118444 ! Fit after self-cal, rms=0.190216Jy sigma=1.103440 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 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5355 mas, bmaj=1.268 mas, bpa=-9.4 degrees ! Estimated noise=1.67443 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.66969 Jy ! Component: 100 - total flux cleaned = 2.25783 Jy ! Component: 150 - total flux cleaned = 2.59363 Jy ! Component: 200 - total flux cleaned = 2.79024 Jy ! Total flux subtracted in 200 components = 2.79024 Jy ! Clean residual min=-0.054213 max=0.099545 Jy/beam ! Clean residual mean=-0.000072 rms=0.008915 Jy/beam ! Combined flux in latest and established models = 2.79024 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 11 model components to the UV plane model. ! The established model now contains 11 components and 2.79024 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5539 mas, bmaj=1.514 mas, bpa=-8.543 degrees ! Estimated noise=1.20639 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.166154 Jy ! Component: 100 - total flux cleaned = 0.279006 Jy ! Component: 150 - total flux cleaned = 0.360601 Jy ! Component: 200 - total flux cleaned = 0.421301 Jy ! Total flux subtracted in 200 components = 0.421301 Jy ! Clean residual min=-0.031084 max=0.049839 Jy/beam ! Clean residual mean=-0.001544 rms=0.006174 Jy/beam ! Combined flux in latest and established models = 3.21154 Jy ! Adding 30 model components to the UV plane model. ! The established model now contains 36 components and 3.21154 Jy ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Added new window around map position (-2.6, 1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0510077 Jy ! Component: 100 - total flux cleaned = 0.0758544 Jy ! Component: 150 - total flux cleaned = 0.0930052 Jy ! Component: 200 - total flux cleaned = 0.105713 Jy ! Total flux subtracted in 200 components = 0.105713 Jy ! Clean residual min=-0.025832 max=0.020651 Jy/beam ! Clean residual mean=-0.002380 rms=0.004279 Jy/beam ! Combined flux in latest and established models = 3.31726 Jy ! Adding 27 model components to the UV plane model. ! The established model now contains 50 components and 3.31725 Jy ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 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 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=2.244 mas, bmaj=3.163 mas, bpa=-12.78 degrees ! Estimated noise=1.74504 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0333185 Jy ! Component: 100 - total flux cleaned = 0.042815 Jy ! Component: 150 - total flux cleaned = 0.0450857 Jy ! Component: 200 - total flux cleaned = 0.0461267 Jy ! Total flux subtracted in 200 components = 0.0461267 Jy ! Clean residual min=-0.026369 max=0.012471 Jy/beam ! Clean residual mean=-0.007645 rms=0.004771 Jy/beam ! Combined flux in latest and established models = 3.36338 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.00094001 Jy ! Component: 100 - total flux cleaned = 0.00138642 Jy ! Component: 150 - total flux cleaned = 0.00140151 Jy ! Component: 200 - total flux cleaned = 0.00139335 Jy ! Total flux subtracted in 200 components = 0.00139335 Jy ! Clean residual min=-0.025631 max=0.010857 Jy/beam ! Clean residual mean=-0.007677 rms=0.004592 Jy/beam ! Combined flux in latest and established models = 3.36477 Jy ! Adding 31 model components to the UV plane model. ! The established model now contains 75 components and 3.36477 Jy ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! Estimated beam: bmin=0.5539 mas, bmaj=1.514 mas, bpa=-8.543 degrees ! Estimated noise=1.20639 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00227083 Jy ! Component: 100 - total flux cleaned = 0.00312649 Jy ! Component: 150 - total flux cleaned = 0.0038931 Jy ! Component: 200 - total flux cleaned = 0.00458035 Jy ! Total flux subtracted in 200 components = 0.00458035 Jy ! Clean residual min=-0.026941 max=0.015197 Jy/beam ! Clean residual mean=-0.002794 rms=0.003424 Jy/beam ! Combined flux in latest and established models = 3.36936 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 89 components and 3.36935 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.203316Jy sigma=1.193429 ! Fit after self-cal, rms=0.191078Jy sigma=1.101322 wmodel J0006-0623_X_map.mod ! Writing 89 model components to file: J0006-0623_X_map.mod wobs J0006-0623_X_uvs.fits ! Writing UV FITS file: J0006-0623_X_uvs.fits wwins J0006-0623_X_map.win ! wwins: Wrote 8 windows to J0006-0623_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 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.0100699 Jy ! Component: 100 - total flux cleaned = -0.0188623 Jy ! Component: 150 - total flux cleaned = -0.0267793 Jy ! Component: 200 - total flux cleaned = -0.0315312 Jy ! Component: 250 - total flux cleaned = -0.0360686 Jy ! Component: 300 - total flux cleaned = -0.039909 Jy ! Component: 350 - total flux cleaned = -0.0444413 Jy ! Component: 400 - total flux cleaned = -0.0491111 Jy ! Component: 450 - total flux cleaned = -0.0536696 Jy ! Component: 500 - total flux cleaned = -0.0568943 Jy ! Component: 550 - total flux cleaned = -0.0595829 Jy ! Component: 600 - total flux cleaned = -0.0622285 Jy ! Component: 650 - total flux cleaned = -0.065074 Jy ! Component: 700 - total flux cleaned = -0.0678763 Jy ! Component: 750 - total flux cleaned = -0.070865 Jy ! Component: 800 - total flux cleaned = -0.0738092 Jy ! Component: 850 - total flux cleaned = -0.0762673 Jy ! Component: 900 - total flux cleaned = -0.0793517 Jy ! Component: 950 - total flux cleaned = -0.0819567 Jy ! Component: 1000 - total flux cleaned = -0.0845344 Jy ! Component: 1050 - total flux cleaned = -0.0877124 Jy ! Component: 1100 - total flux cleaned = -0.0904321 Jy ! Component: 1150 - total flux cleaned = -0.0929128 Jy ! Component: 1200 - total flux cleaned = -0.0957719 Jy ! Component: 1250 - total flux cleaned = -0.0981967 Jy ! Component: 1300 - total flux cleaned = -0.100595 Jy ! Component: 1350 - total flux cleaned = -0.103164 Jy ! Component: 1400 - total flux cleaned = -0.105509 Jy ! Component: 1450 - total flux cleaned = -0.108024 Jy ! Component: 1500 - total flux cleaned = -0.110321 Jy ! Component: 1550 - total flux cleaned = -0.112974 Jy ! Component: 1600 - total flux cleaned = -0.115598 Jy ! Component: 1650 - total flux cleaned = -0.117454 Jy ! Component: 1700 - total flux cleaned = -0.120028 Jy ! Component: 1750 - total flux cleaned = -0.122213 Jy ! Component: 1800 - total flux cleaned = -0.124736 Jy ! Component: 1850 - total flux cleaned = -0.126699 Jy ! Component: 1900 - total flux cleaned = -0.128998 Jy ! Component: 1950 - total flux cleaned = -0.131275 Jy ! Component: 2000 - total flux cleaned = -0.133358 Jy ! Component: 2050 - total flux cleaned = -0.135423 Jy ! Component: 2100 - total flux cleaned = -0.137469 Jy ! Component: 2150 - total flux cleaned = -0.139161 Jy ! Component: 2200 - total flux cleaned = -0.141675 Jy ! Component: 2250 - total flux cleaned = -0.143505 Jy ! Component: 2300 - total flux cleaned = -0.146141 Jy ! Component: 2350 - total flux cleaned = -0.147776 Jy ! Component: 2400 - total flux cleaned = -0.150044 Jy ! Component: 2450 - total flux cleaned = -0.152292 Jy ! Component: 2500 - total flux cleaned = -0.154045 Jy ! Component: 2550 - total flux cleaned = -0.15673 Jy ! Component: 2600 - total flux cleaned = -0.158454 Jy ! Component: 2650 - total flux cleaned = -0.160319 Jy ! Component: 2700 - total flux cleaned = -0.162323 Jy ! Component: 2750 - total flux cleaned = -0.164464 Jy ! Component: 2800 - total flux cleaned = -0.166587 Jy ! Component: 2850 - total flux cleaned = -0.168242 Jy ! Component: 2900 - total flux cleaned = -0.170333 Jy ! Component: 2950 - total flux cleaned = -0.172109 Jy ! Component: 3000 - total flux cleaned = -0.173874 Jy ! Component: 3050 - total flux cleaned = -0.175917 Jy ! Component: 3100 - total flux cleaned = -0.177364 Jy ! Component: 3150 - total flux cleaned = -0.178946 Jy ! Component: 3200 - total flux cleaned = -0.180801 Jy ! Component: 3250 - total flux cleaned = -0.182642 Jy ! Component: 3300 - total flux cleaned = -0.184328 Jy ! Component: 3350 - total flux cleaned = -0.186143 Jy ! Component: 3400 - total flux cleaned = -0.187668 Jy ! Component: 3450 - total flux cleaned = -0.189593 Jy ! Component: 3500 - total flux cleaned = -0.191369 Jy ! Component: 3550 - total flux cleaned = -0.192861 Jy ! Component: 3600 - total flux cleaned = -0.194611 Jy ! Component: 3650 - total flux cleaned = -0.196215 Jy ! Component: 3700 - total flux cleaned = -0.19781 Jy ! Component: 3750 - total flux cleaned = -0.199261 Jy ! Component: 3800 - total flux cleaned = -0.200703 Jy ! Component: 3850 - total flux cleaned = -0.202266 Jy ! Component: 3900 - total flux cleaned = -0.203818 Jy ! Component: 3950 - total flux cleaned = -0.205231 Jy ! Component: 4000 - total flux cleaned = -0.207017 Jy ! Component: 4050 - total flux cleaned = -0.208157 Jy ! Total flux subtracted in 4096 components = -0.209669 Jy ! Clean residual min=-0.006757 max=0.006284 Jy/beam ! Clean residual mean=-0.001090 rms=0.001840 Jy/beam ! Combined flux in latest and established models = 3.15969 Jy keep ! Adding 290 model components to the UV plane model. ! The established model now contains 377 components and 3.15969 Jy ![Exited script file: /opt64/bin/pima_mupet_01.dfm] wmap /scr/pima/rv123_uvs/J0006-0623_X_map.fits ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5539 x 1.514 at -8.543 degrees (North through East) ! Clean map min=-0.026882 max=2.0152 Jy/beam ! Writing clean map to FITS file: /scr/pima/rv123_uvs/J0006-0623_X_map.fits wmodel /scr/pima/rv123_uvs/J0006-0623_X_map.mod ! Writing 377 model components to file: /scr/pima/rv123_uvs/J0006-0623_X_map.mod wobs /scr/pima/rv123_uvs/J0006-0623_X_uvs.fits ! Writing UV FITS file: /scr/pima/rv123_uvs/J0006-0623_X_uvs.fits wwins /scr/pima/rv123_uvs/J0006-0623_X_map.win ! wwins: Wrote 9 windows to /scr/pima/rv123_uvs/J0006-0623_X_map.win float peak_flux peak_flux = peak(flux,max) print "The peak flux is", peak_flux ! The peak flux is 2.01515 invert ! Inverting map ! Your choice of large map pixels excluded 4.21% of the data. ! The x-axis pixel size should ideally be below 0.1642 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00184052 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.274002 0.274002 0.548005 1.09601 2.19202 4.38404 8.76808 17.5362 35.0723 70.1446 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0055 print "peak= ", int(1000*peak_flux)/1000 ! peak= 2.015 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1094.88 quit ! Quitting program ! Log file /scr/pima/rv123_uvs/J0006-0623_X_dfm.log closed on Sun Dec 23 00:38:50 2018