! Started logfile: J0118-2141_S_2013_01_09_pus_uvs.log on Tue Nov 4 18:58:59 2014 obs J0118-2141_S_2013_01_09_pus_uva.fits ! Reading UV FITS file: J0118-2141_S_2013_01_09_pus_uva.fits ! AN table 1: 55 integrations on 120 of 120 possible baselines. ! AN table 2: 78 integrations on 120 of 120 possible baselines. ! Apparent sampling: 0.099812 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0118-2141 ! ! 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 2.23697e+09 8e+06 1 8e+06 ! 02 2 2.26697e+09 8e+06 1 8e+06 ! 03 3 2.35697e+09 8e+06 1 8e+06 ! 04 4 2.37697e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 3254 lines of history. ! ! Reading 6372 visibilities. select I ! Polarization I is unavailable. ! 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 = 512 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J0118-2141_S_2013_01_09_pus] 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 mapsize field_size, field_cell ! Map grid = 512x512 pixels with 0.500x0.500 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 711 telescope corrections were flagged in sub-array 1. ! A total of 991 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 711 telescope corrections were flagged in sub-array 1. ! A total of 991 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 711 telescope corrections were flagged in sub-array 1. ! A total of 991 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 711 telescope corrections were flagged in sub-array 1. ! A total of 991 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.465183Jy sigma=0.003968 ! Fit after self-cal, rms=0.438271Jy sigma=0.003819 ! 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 ! Estimated beam: bmin=2.897 mas, bmaj=7.889 mas, bpa=-1.931 degrees ! Estimated noise=1310.42 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.489572 Jy ! Component: 100 - total flux cleaned = 0.596609 Jy ! Total flux subtracted in 100 components = 0.596609 Jy ! Clean residual min=-0.025469 max=0.036958 Jy/beam ! Clean residual mean=0.000412 rms=0.008060 Jy/beam ! Combined flux in latest and established models = 0.596609 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.596609 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.198055Jy sigma=0.001623 ! Fit after self-cal, rms=0.198049Jy sigma=0.001623 ! Inverting map 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 ! Estimated beam: bmin=3.617 mas, bmaj=9.601 mas, bpa=-3.757 degrees ! Estimated noise=1051.62 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.0295479 Jy ! Component: 100 - total flux cleaned = 0.0463017 Jy ! Component: 150 - total flux cleaned = 0.0540845 Jy ! Component: 200 - total flux cleaned = 0.0553027 Jy ! Total flux subtracted in 200 components = 0.0553027 Jy ! Clean residual min=-0.022438 max=0.029798 Jy/beam ! Clean residual mean=-0.003409 rms=0.007279 Jy/beam ! Combined flux in latest and established models = 0.651912 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 14 components and 0.651912 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.193766Jy sigma=0.001580 ! Fit after self-cal, rms=0.193743Jy sigma=0.001580 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 1.20 HH 1.00* HN 1.08 ! KP 0.97 LA 1.03 MC 1.00* MK 1.00* ! NL 0.86 NY 1.00* OV 0.93 PT 1.02 ! SC 0.90 TS 1.00* WF 1.02 WZ 1.19 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.91 FD 1.18 HH 1.00* HN 1.20 ! KP 0.98 LA 1.02 MC 1.00* MK 0.77 ! NL 0.91 NY 1.00* OV 0.89 PT 1.04 ! SC 0.96 TS 1.00* WF 1.03 WZ 1.00* ! ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 1.32 HH 1.00* HN 1.17 ! KP 1.02 LA 0.95 MC 1.00* MK 1.00* ! NL 0.95 NY 1.00* OV 0.96 PT 0.97 ! SC 1.09 TS 1.00* WF 1.10 WZ 1.06 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.31 HH 1.00* HN 1.29 ! KP 0.98 LA 0.95 MC 1.00* MK 0.90 ! NL 0.98 NY 1.00* OV 0.97 PT 0.99 ! SC 1.09 TS 1.00* WF 1.12 WZ 1.00* ! ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 1.42 HH 1.00* HN 1.13 ! KP 0.81 LA 0.89 MC 1.00* MK 1.00* ! NL 0.94 NY 1.00* OV 0.94 PT 0.97 ! SC 0.86 TS 1.00* WF 1.07 WZ 1.16 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 1.32 HH 1.00* HN 1.25 ! KP 0.81 LA 0.92 MC 1.00* MK 0.87 ! NL 0.98 NY 1.00* OV 0.94 PT 0.98 ! SC 0.95 TS 1.00* WF 1.02 WZ 1.00* ! ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.33 HH 1.00* HN 1.05 ! KP 0.80 LA 0.75 MC 1.00* MK 1.00* ! NL 0.88 NY 1.00* OV 0.82 PT 0.85 ! SC 0.89 TS 1.00* WF 1.15 WZ 1.01 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 1.35 HH 1.00* HN 1.15 ! KP 0.79 LA 0.74 MC 1.00* MK 1.00 ! NL 0.91 NY 1.00* OV 0.84 PT 0.87 ! SC 0.97 TS 1.00* WF 1.13 WZ 1.00* ! ! ! Fit before self-cal, rms=0.193743Jy sigma=0.001580 ! Fit after self-cal, rms=0.158269Jy sigma=0.001161 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.591 mas, bmaj=9.504 mas, bpa=-3.671 degrees ! Estimated noise=1059.93 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.00542183 Jy ! Component: 100 - total flux cleaned = -0.00662785 Jy ! Component: 150 - total flux cleaned = -0.00706488 Jy ! Component: 200 - total flux cleaned = -0.00758638 Jy ! Total flux subtracted in 200 components = -0.00758638 Jy ! Clean residual min=-0.006631 max=0.006394 Jy/beam ! Clean residual mean=-0.000216 rms=0.001689 Jy/beam ! Combined flux in latest and established models = 0.644325 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 21 components and 0.644325 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.157781Jy sigma=0.001158 ! Fit after self-cal, rms=0.158105Jy sigma=0.001157 ! Inverting map 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.158105Jy sigma=0.001157 ! Fit after self-cal, rms=0.156522Jy sigma=0.001147 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.595 mas, bmaj=9.527 mas, bpa=-3.729 degrees ! Estimated noise=1057.56 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00412683 Jy ! Component: 100 - total flux cleaned = -0.00494375 Jy ! Component: 150 - total flux cleaned = -0.00512242 Jy ! Component: 200 - total flux cleaned = -0.00535999 Jy ! Total flux subtracted in 200 components = -0.00535999 Jy ! Clean residual min=-0.004596 max=0.004634 Jy/beam ! Clean residual mean=-0.000110 rms=0.001312 Jy/beam ! Combined flux in latest and established models = 0.638965 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 24 components and 0.638965 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156307Jy sigma=0.001145 ! Fit after self-cal, rms=0.156305Jy sigma=0.001145 flux_cutoff = imstat(rms) * dynam ! Inverting map 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.000179319 Jy ! Component: 100 - total flux cleaned = -0.000437509 Jy ! Component: 150 - total flux cleaned = -0.000611173 Jy ! Component: 200 - total flux cleaned = -0.000770211 Jy ! Total flux subtracted in 200 components = -0.000770211 Jy ! Clean residual min=-0.004396 max=0.004687 Jy/beam ! Clean residual mean=-0.000093 rms=0.001293 Jy/beam ! Combined flux in latest and established models = 0.638195 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 28 components and 0.638195 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156247Jy sigma=0.001145 ! Fit after self-cal, rms=0.156238Jy sigma=0.001145 ! Inverting map 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.156238Jy sigma=0.001145 ! Fit after self-cal, rms=0.156345Jy sigma=0.001144 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.592 mas, bmaj=9.53 mas, bpa=-3.716 degrees ! Estimated noise=1056.4 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00213783 Jy ! Component: 100 - total flux cleaned = -0.00268873 Jy ! Component: 150 - total flux cleaned = -0.00284616 Jy ! Component: 200 - total flux cleaned = -0.00298629 Jy ! Total flux subtracted in 200 components = -0.00298629 Jy ! Clean residual min=-0.004108 max=0.004765 Jy/beam ! Clean residual mean=0.000042 rms=0.001216 Jy/beam ! Combined flux in latest and established models = 0.635209 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 29 components and 0.635209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156281Jy sigma=0.001143 ! Fit after self-cal, rms=0.156278Jy sigma=0.001143 flux_cutoff = imstat(rms) * dynam ! Inverting map 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 = -7.184e-05 Jy ! Component: 100 - total flux cleaned = -0.000230541 Jy ! Component: 150 - total flux cleaned = -0.000265139 Jy ! Component: 200 - total flux cleaned = -0.000298997 Jy ! Total flux subtracted in 200 components = -0.000298997 Jy ! Clean residual min=-0.004016 max=0.004772 Jy/beam ! Clean residual mean=0.000047 rms=0.001207 Jy/beam ! Combined flux in latest and established models = 0.63491 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 32 components and 0.63491 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156258Jy sigma=0.001143 ! Fit after self-cal, rms=0.156253Jy sigma=0.001143 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 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.156253Jy sigma=0.001143 ! Fit after self-cal, rms=0.156455Jy sigma=0.001143 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.575 mas, bmaj=15.83 mas, bpa=-7.035 degrees ! Estimated noise=1226.05 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.48729e-05 Jy ! Component: 100 - total flux cleaned = -0.000126477 Jy ! Component: 150 - total flux cleaned = -0.000282323 Jy ! Component: 200 - total flux cleaned = -0.000404607 Jy ! Total flux subtracted in 200 components = -0.000404607 Jy ! Clean residual min=-0.004336 max=0.005089 Jy/beam ! Clean residual mean=0.000228 rms=0.001343 Jy/beam ! Combined flux in latest and established models = 0.634505 Jy selfcal ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 44 components and 0.634505 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156463Jy sigma=0.001143 ! Fit after self-cal, rms=0.156444Jy sigma=0.001142 flux_cutoff = imstat(rms) * dynam ! Inverting map 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 = 7.28462e-05 Jy ! Component: 100 - total flux cleaned = -1.32453e-05 Jy ! Component: 150 - total flux cleaned = -7.96967e-05 Jy ! Component: 200 - total flux cleaned = -0.0001068 Jy ! Total flux subtracted in 200 components = -0.0001068 Jy ! Clean residual min=-0.004312 max=0.005030 Jy/beam ! Clean residual mean=0.000231 rms=0.001330 Jy/beam ! Combined flux in latest and established models = 0.634399 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 53 components and 0.634399 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.156458Jy sigma=0.001142 ! Fit after self-cal, rms=0.156439Jy sigma=0.001142 ! Inverting map 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 ! Estimated beam: bmin=3.589 mas, bmaj=9.531 mas, bpa=-3.719 degrees ! Estimated noise=1055.58 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.589 x 9.531 at -3.719 degrees (North through East) ! Clean map min=-0.0040217 max=0.59991 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.156439Jy sigma=0.001142 ! Fit after self-cal, rms=0.148407Jy sigma=0.001056 ! Inverting map and beam ! Estimated beam: bmin=3.586 mas, bmaj=9.533 mas, bpa=-3.702 degrees ! Estimated noise=1055.87 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00103041 Jy ! Component: 100 - total flux cleaned = -0.00109997 Jy ! Component: 150 - total flux cleaned = -0.00116252 Jy ! Component: 200 - total flux cleaned = -0.00121955 Jy ! Total flux subtracted in 200 components = -0.00121955 Jy ! Clean residual min=-0.003875 max=0.004969 Jy/beam ! Clean residual mean=0.000213 rms=0.001156 Jy/beam ! Combined flux in latest and established models = 0.633179 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 57 components and 0.633179 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148388Jy sigma=0.001056 ! Fit after self-cal, rms=0.148386Jy sigma=0.001056 ! Inverting map ! Added new window around map position (6.5, -10.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00488364 Jy ! Component: 100 - total flux cleaned = 0.00589947 Jy ! Component: 150 - total flux cleaned = 0.00575037 Jy ! Component: 200 - total flux cleaned = 0.0054739 Jy ! Total flux subtracted in 200 components = 0.0054739 Jy ! Clean residual min=-0.003369 max=0.004025 Jy/beam ! Clean residual mean=0.000079 rms=0.001022 Jy/beam ! Combined flux in latest and established models = 0.638653 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 72 components and 0.638653 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148191Jy sigma=0.001054 ! Fit after self-cal, rms=0.148171Jy sigma=0.001054 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148171Jy sigma=0.001054 ! Fit after self-cal, rms=0.148148Jy sigma=0.001054 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 ! Estimated beam: bmin=2.864 mas, bmaj=7.782 mas, bpa=-1.932 degrees ! Estimated noise=1280.79 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.476955 Jy ! Component: 100 - total flux cleaned = 0.580962 Jy ! Component: 150 - total flux cleaned = 0.607398 Jy ! Component: 200 - total flux cleaned = 0.61858 Jy ! Total flux subtracted in 200 components = 0.61858 Jy ! Clean residual min=-0.005604 max=0.005860 Jy/beam ! Clean residual mean=0.000171 rms=0.001497 Jy/beam ! Combined flux in latest and established models = 0.61858 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 15 model components to the UV plane model. ! The established model now contains 15 components and 0.61858 Jy ! Inverting map and beam ! Estimated beam: bmin=3.584 mas, bmaj=9.533 mas, bpa=-3.724 degrees ! Estimated noise=1055.7 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.00789646 Jy ! Component: 100 - total flux cleaned = 0.0128406 Jy ! Component: 150 - total flux cleaned = 0.0160579 Jy ! Component: 200 - total flux cleaned = 0.0174441 Jy ! Total flux subtracted in 200 components = 0.0174441 Jy ! Clean residual min=-0.003411 max=0.004282 Jy/beam ! Clean residual mean=0.000166 rms=0.001031 Jy/beam ! Combined flux in latest and established models = 0.636024 Jy ! Adding 30 model components to the UV plane model. ! The established model now contains 43 components and 0.636024 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.567 mas, bmaj=15.83 mas, bpa=-6.935 degrees ! Estimated noise=1227.74 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00215329 Jy ! Component: 100 - total flux cleaned = 0.0029668 Jy ! Component: 150 - total flux cleaned = 0.00315778 Jy ! Component: 200 - total flux cleaned = 0.00329522 Jy ! Total flux subtracted in 200 components = 0.00329522 Jy ! Clean residual min=-0.003879 max=0.004240 Jy/beam ! Clean residual mean=0.000160 rms=0.001185 Jy/beam ! Combined flux in latest and established models = 0.639319 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.000126473 Jy ! Component: 100 - total flux cleaned = 0.000241179 Jy ! Component: 150 - total flux cleaned = 0.000320993 Jy ! Component: 200 - total flux cleaned = 0.000370909 Jy ! Total flux subtracted in 200 components = 0.000370909 Jy ! Clean residual min=-0.003895 max=0.004212 Jy/beam ! Clean residual mean=0.000142 rms=0.001174 Jy/beam ! Combined flux in latest and established models = 0.63969 Jy ! Adding 25 model components to the UV plane model. ! The established model now contains 64 components and 0.63969 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.584 mas, bmaj=9.533 mas, bpa=-3.724 degrees ! Estimated noise=1055.7 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000991907 Jy ! Component: 100 - total flux cleaned = -0.00122659 Jy ! Component: 150 - total flux cleaned = -0.00122606 Jy ! Component: 200 - total flux cleaned = -0.00113028 Jy ! Total flux subtracted in 200 components = -0.00113028 Jy ! Clean residual min=-0.003395 max=0.004041 Jy/beam ! Clean residual mean=0.000102 rms=0.000972 Jy/beam ! Combined flux in latest and established models = 0.63856 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 84 components and 0.63856 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148078Jy sigma=0.001053 ! Fit after self-cal, rms=0.148050Jy sigma=0.001053 wmodel J0118-2141_S_2013_01_09_pus_map.mod ! Writing 84 model components to file: J0118-2141_S_2013_01_09_pus_map.mod wobs J0118-2141_S_2013_01_09_pus_uvs.fits ! Writing UV FITS file: J0118-2141_S_2013_01_09_pus_uvs.fits wwins J0118-2141_S_2013_01_09_pus_map.win ! wwins: Wrote 2 windows to J0118-2141_S_2013_01_09_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000872823 Jy ! Component: 100 - total flux cleaned = 0.00135531 Jy ! Component: 150 - total flux cleaned = 0.00197135 Jy ! Component: 200 - total flux cleaned = 0.00218274 Jy ! Component: 250 - total flux cleaned = 0.00233216 Jy ! Component: 300 - total flux cleaned = 0.00238097 Jy ! Component: 350 - total flux cleaned = 0.00247501 Jy ! Component: 400 - total flux cleaned = 0.00256506 Jy ! Component: 450 - total flux cleaned = 0.00256534 Jy ! Component: 500 - total flux cleaned = 0.0026513 Jy ! Component: 550 - total flux cleaned = 0.00269297 Jy ! Component: 600 - total flux cleaned = 0.00273404 Jy ! Component: 650 - total flux cleaned = 0.00281357 Jy ! Component: 700 - total flux cleaned = 0.0028522 Jy ! Component: 750 - total flux cleaned = 0.00285211 Jy ! Component: 800 - total flux cleaned = 0.00285222 Jy ! Component: 850 - total flux cleaned = 0.00285238 Jy ! Component: 900 - total flux cleaned = 0.00285241 Jy ! Component: 950 - total flux cleaned = 0.00278239 Jy ! Component: 1000 - total flux cleaned = 0.0028172 Jy ! Component: 1050 - total flux cleaned = 0.00285114 Jy ! Component: 1100 - total flux cleaned = 0.00285109 Jy ! Component: 1150 - total flux cleaned = 0.0029171 Jy ! Component: 1200 - total flux cleaned = 0.00298171 Jy ! Component: 1250 - total flux cleaned = 0.00301348 Jy ! Component: 1300 - total flux cleaned = 0.00304532 Jy ! Component: 1350 - total flux cleaned = 0.00310753 Jy ! Component: 1400 - total flux cleaned = 0.00319959 Jy ! Component: 1450 - total flux cleaned = 0.0032299 Jy ! Component: 1500 - total flux cleaned = 0.00331994 Jy ! Component: 1550 - total flux cleaned = 0.00340888 Jy ! Component: 1600 - total flux cleaned = 0.0034965 Jy ! Component: 1650 - total flux cleaned = 0.00352521 Jy ! Component: 1700 - total flux cleaned = 0.00358237 Jy ! Component: 1750 - total flux cleaned = 0.00363885 Jy ! Component: 1800 - total flux cleaned = 0.00366693 Jy ! Component: 1850 - total flux cleaned = 0.00377755 Jy ! Component: 1900 - total flux cleaned = 0.003805 Jy ! Component: 1950 - total flux cleaned = 0.00383236 Jy ! Component: 2000 - total flux cleaned = 0.00391306 Jy ! Total flux subtracted in 2048 components = 0.00399272 Jy ! Clean residual min=-0.001460 max=0.001661 Jy/beam ! Clean residual mean=0.000002 rms=0.000562 Jy/beam ! Combined flux in latest and established models = 0.642553 Jy keep ! Adding 148 model components to the UV plane model. ! The established model now contains 232 components and 0.642553 Jy ![Exited script file: muppet] wmap "J0118-2141_S_2013_01_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.584 x 9.533 at -3.724 degrees (North through East) ! Clean map min=-0.0027321 max=0.60322 Jy/beam ! Writing clean map to FITS file: J0118-2141_S_2013_01_09_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.603217 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000560898 device "J0118-2141_S_2013_01_09_pus\_map.eps/VPS" ! Attempting to open device: 'J0118-2141_S_2013_01_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.278953 0.278953 0.557907 1.11581 2.23163 4.46325 8.92651 17.853 35.706 71.4121 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0016 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.603 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.584 x 9.533 at -3.724 degrees (North through East) ! Clean map min=-0.0027321 max=0.60322 Jy/beam quit ! Quitting program ! Log file J0118-2141_S_2013_01_09_pus_uvs.log closed on Tue Nov 4 18:59:02 2014