! Started logfile: J1316-3338_S_2013_01_09_pus_uvs.log on Tue Nov 4 19:01:41 2014 obs J1316-3338_S_2013_01_09_pus_uva.fits ! Reading UV FITS file: J1316-3338_S_2013_01_09_pus_uva.fits ! AN table 1: 10 integrations on 120 of 120 possible baselines. ! AN table 2: 71 integrations on 120 of 120 possible baselines. ! AN table 3: 10 integrations on 120 of 120 possible baselines. ! Apparent sampling: 0.191209 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1316-3338 ! ! 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 3255 lines of history. ! ! Reading 8352 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 J1316-3338_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 102 telescope corrections were flagged in sub-array 1. ! A total of 704 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 102 telescope corrections were flagged in sub-array 1. ! A total of 704 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 102 telescope corrections were flagged in sub-array 1. ! A total of 704 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 102 telescope corrections were flagged in sub-array 1. ! A total of 704 telescope corrections were flagged in sub-array 2. ! A total of 75 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.354763Jy sigma=0.002888 ! Fit after self-cal, rms=0.337936Jy sigma=0.002760 ! 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.304 mas, bmaj=5.286 mas, bpa=3.312 degrees ! Estimated noise=1434.74 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.591905 Jy ! Component: 100 - total flux cleaned = 0.723098 Jy ! Total flux subtracted in 100 components = 0.723098 Jy ! Clean residual min=-0.037684 max=0.077653 Jy/beam ! Clean residual mean=-0.000069 rms=0.009470 Jy/beam ! Combined flux in latest and established models = 0.723098 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.723098 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.281673Jy sigma=0.002296 ! Fit after self-cal, rms=0.281654Jy sigma=0.002296 ! Inverting map ! Added new window around map position (-5.5, -0.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0720619 Jy ! Component: 100 - total flux cleaned = 0.113732 Jy ! Total flux subtracted in 100 components = 0.113732 Jy ! Clean residual min=-0.036057 max=0.069533 Jy/beam ! Clean residual mean=-0.000157 rms=0.008047 Jy/beam ! Combined flux in latest and established models = 0.836829 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 22 components and 0.836829 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.249355Jy sigma=0.001966 ! Fit after self-cal, rms=0.245543Jy sigma=0.001933 ! Inverting map ! Added new window around map position (4, -6.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0517151 Jy ! Component: 100 - total flux cleaned = 0.0853423 Jy ! Total flux subtracted in 100 components = 0.0853423 Jy ! Clean residual min=-0.026744 max=0.035125 Jy/beam ! Clean residual mean=-0.000235 rms=0.006746 Jy/beam ! Combined flux in latest and established models = 0.922172 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 31 components and 0.922172 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.228336Jy sigma=0.001780 ! Fit after self-cal, rms=0.228153Jy sigma=0.001777 ! 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.217 mas, bmaj=8.361 mas, bpa=-0.9742 degrees ! Estimated noise=906.17 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.0364262 Jy ! Component: 100 - total flux cleaned = 0.0572544 Jy ! Component: 150 - total flux cleaned = 0.0687351 Jy ! Component: 200 - total flux cleaned = 0.072345 Jy ! Total flux subtracted in 200 components = 0.072345 Jy ! Clean residual min=-0.026285 max=0.026454 Jy/beam ! Clean residual mean=-0.002485 rms=0.006999 Jy/beam ! Combined flux in latest and established models = 0.994517 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 45 components and 0.994517 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.220532Jy sigma=0.001737 ! Fit after self-cal, rms=0.220398Jy sigma=0.001736 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.16 HH 1.00* HN 1.43 ! KP 1.03 LA 1.02 MC 1.00* MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.03 ! SC 0.74 TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 42 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.93 FD 1.11 HH 1.39 HN 1.22 ! KP 1.00 LA 1.03 MC 1.00* MK 0.81 ! NL 0.96 NY 1.00* OV 0.93 PT 1.03 ! SC 0.93 TS 1.00* WF 1.06 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.11 HH 1.00* HN 1.19 ! KP 0.98 LA 1.01 MC 1.00* MK 0.83 ! NL 0.91 NY 1.00* OV 0.93 PT 1.02 ! SC 0.96 TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.25 HH 1.00* HN 1.68 ! KP 1.08 LA 0.97 MC 1.00* MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.01 ! SC 0.83 TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 42 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.22 HH 1.04 HN 1.35 ! KP 1.01 LA 0.96 MC 1.00* MK 0.96 ! NL 1.04 NY 1.00* OV 0.99 PT 0.99 ! SC 1.04 TS 1.00* WF 1.10 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.24 HH 1.00* HN 1.28 ! KP 0.99 LA 0.94 MC 1.00* MK 0.96 ! NL 0.97 NY 1.00* OV 0.97 PT 1.00 ! SC 1.13 TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.23 HH 1.00* HN 1.50 ! KP 0.87 LA 0.94 MC 1.00* MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 1.01 ! SC 0.71 TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 42 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 1.22 HH 1.00 HN 1.22 ! KP 0.84 LA 0.93 MC 1.00* MK 0.93 ! NL 1.03 NY 1.00* OV 0.98 PT 0.98 ! SC 0.95 TS 1.00* WF 1.08 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.20 HH 1.00* HN 1.21 ! KP 0.83 LA 0.93 MC 1.00* MK 0.92 ! NL 0.90 NY 1.00* OV 0.96 PT 1.00 ! SC 0.99 TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.22 HH 1.00* HN 1.40 ! KP 0.89 LA 0.76 MC 1.00* MK 1.00* ! NL 1.00* NY 1.00* OV 1.00* PT 0.91 ! SC 0.79 TS 1.00* WF 1.00* WZ 1.00* ! ! A total of 42 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 1.23 HH 0.94 HN 1.16 ! KP 0.83 LA 0.77 MC 1.00* MK 1.02 ! NL 0.94 NY 1.00* OV 0.85 PT 0.88 ! SC 0.98 TS 1.00* WF 1.17 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.21 HH 1.00* HN 1.09 ! KP 0.82 LA 0.78 MC 1.00* MK 1.05 ! NL 0.86 NY 1.00* OV 0.85 PT 0.88 ! SC 1.06 TS 1.00* WF 1.00* WZ 1.00* ! ! ! Fit before self-cal, rms=0.220398Jy sigma=0.001736 ! Fit after self-cal, rms=0.165198Jy sigma=0.001243 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.245 mas, bmaj=8.948 mas, bpa=-2.894 degrees ! Estimated noise=911.862 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 (-9, -0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000914262 Jy ! Component: 100 - total flux cleaned = -0.00241143 Jy ! Component: 150 - total flux cleaned = -0.0020508 Jy ! Component: 200 - total flux cleaned = -0.00177599 Jy ! Total flux subtracted in 200 components = -0.00177599 Jy ! Clean residual min=-0.011251 max=0.006107 Jy/beam ! Clean residual mean=-0.000682 rms=0.002195 Jy/beam ! Combined flux in latest and established models = 0.992741 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 65 components and 0.992741 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162460Jy sigma=0.001220 ! Fit after self-cal, rms=0.162360Jy sigma=0.001218 ! 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.162360Jy sigma=0.001218 ! Fit after self-cal, rms=0.162404Jy sigma=0.001198 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.258 mas, bmaj=9.003 mas, bpa=-2.937 degrees ! Estimated noise=912.696 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00342305 Jy ! Component: 100 - total flux cleaned = -0.00372393 Jy ! Component: 150 - total flux cleaned = -0.00373298 Jy ! Component: 200 - total flux cleaned = -0.00353603 Jy ! Total flux subtracted in 200 components = -0.00353603 Jy ! Clean residual min=-0.006663 max=0.004990 Jy/beam ! Clean residual mean=-0.000462 rms=0.001725 Jy/beam ! Combined flux in latest and established models = 0.989205 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 83 components and 0.989205 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161441Jy sigma=0.001191 ! Fit after self-cal, rms=0.161429Jy sigma=0.001191 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.00089806 Jy ! Component: 100 - total flux cleaned = 0.00119344 Jy ! Component: 150 - total flux cleaned = 0.00119408 Jy ! Component: 200 - total flux cleaned = 0.00131651 Jy ! Total flux subtracted in 200 components = 0.00131651 Jy ! Clean residual min=-0.006201 max=0.004501 Jy/beam ! Clean residual mean=-0.000476 rms=0.001635 Jy/beam ! Combined flux in latest and established models = 0.990521 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 91 components and 0.990521 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.161129Jy sigma=0.001188 ! Fit after self-cal, rms=0.161094Jy sigma=0.001188 ! 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.161094Jy sigma=0.001188 ! Fit after self-cal, rms=0.158930Jy sigma=0.001166 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.261 mas, bmaj=9.004 mas, bpa=-2.897 degrees ! Estimated noise=913.072 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00332223 Jy ! Component: 100 - total flux cleaned = -0.00400177 Jy ! Component: 150 - total flux cleaned = -0.00414483 Jy ! Component: 200 - total flux cleaned = -0.00414463 Jy ! Total flux subtracted in 200 components = -0.00414463 Jy ! Clean residual min=-0.005790 max=0.004200 Jy/beam ! Clean residual mean=-0.000520 rms=0.001348 Jy/beam ! Combined flux in latest and established models = 0.986376 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 100 components and 0.986376 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158522Jy sigma=0.001163 ! Fit after self-cal, rms=0.158533Jy sigma=0.001163 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.000305311 Jy ! Component: 100 - total flux cleaned = 0.000408111 Jy ! Component: 150 - total flux cleaned = 0.000501378 Jy ! Component: 200 - total flux cleaned = 0.000588235 Jy ! Total flux subtracted in 200 components = 0.000588235 Jy ! Clean residual min=-0.005532 max=0.003876 Jy/beam ! Clean residual mean=-0.000527 rms=0.001306 Jy/beam ! Combined flux in latest and established models = 0.986965 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 105 components and 0.986965 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158394Jy sigma=0.001162 ! Fit after self-cal, rms=0.158379Jy sigma=0.001162 ! 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.158379Jy sigma=0.001162 ! Fit after self-cal, rms=0.158987Jy sigma=0.001161 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.673 mas, bmaj=22.1 mas, bpa=-3.535 degrees ! Estimated noise=1075.59 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00305426 Jy ! Component: 100 - total flux cleaned = -0.00322267 Jy ! Component: 150 - total flux cleaned = -0.0032255 Jy ! Component: 200 - total flux cleaned = -0.00322626 Jy ! Total flux subtracted in 200 components = -0.00322626 Jy ! Clean residual min=-0.006324 max=0.004483 Jy/beam ! Clean residual mean=-0.000871 rms=0.001605 Jy/beam ! Combined flux in latest and established models = 0.983738 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 110 components and 0.983739 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158920Jy sigma=0.001160 ! Fit after self-cal, rms=0.158916Jy sigma=0.001160 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.000270615 Jy ! Component: 100 - total flux cleaned = 0.000435744 Jy ! Component: 150 - total flux cleaned = 0.000484201 Jy ! Component: 200 - total flux cleaned = 0.000531974 Jy ! Total flux subtracted in 200 components = 0.000531974 Jy ! Clean residual min=-0.006179 max=0.004284 Jy/beam ! Clean residual mean=-0.000881 rms=0.001574 Jy/beam ! Combined flux in latest and established models = 0.984271 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 114 components and 0.984271 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.158893Jy sigma=0.001160 ! Fit after self-cal, rms=0.158886Jy sigma=0.001160 ! 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.264 mas, bmaj=9.002 mas, bpa=-2.857 degrees ! Estimated noise=913.383 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.264 x 9.002 at -2.857 degrees (North through East) ! Clean map min=-0.004965 max=0.75675 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.158886Jy sigma=0.001160 ! Fit after self-cal, rms=0.151431Jy sigma=0.001075 ! Inverting map and beam ! Estimated beam: bmin=3.268 mas, bmaj=9.006 mas, bpa=-2.784 degrees ! Estimated noise=913.799 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00279868 Jy ! Component: 100 - total flux cleaned = -0.00309781 Jy ! Component: 150 - total flux cleaned = -0.00310072 Jy ! Component: 200 - total flux cleaned = -0.0030999 Jy ! Total flux subtracted in 200 components = -0.0030999 Jy ! Clean residual min=-0.004961 max=0.003746 Jy/beam ! Clean residual mean=-0.000375 rms=0.001200 Jy/beam ! Combined flux in latest and established models = 0.981171 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 126 components and 0.981171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151265Jy sigma=0.001074 ! Fit after self-cal, rms=0.151254Jy sigma=0.001074 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.30789e-05 Jy ! Component: 100 - total flux cleaned = 0.000146256 Jy ! Component: 150 - total flux cleaned = 0.000215834 Jy ! Component: 200 - total flux cleaned = 0.000216191 Jy ! Total flux subtracted in 200 components = 0.000216191 Jy ! Clean residual min=-0.004870 max=0.003644 Jy/beam ! Clean residual mean=-0.000377 rms=0.001184 Jy/beam ! Combined flux in latest and established models = 0.981387 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 131 components and 0.981387 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151174Jy sigma=0.001074 ! Fit after self-cal, rms=0.151165Jy sigma=0.001074 ! 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.151165Jy sigma=0.001074 ! Fit after self-cal, rms=0.151616Jy sigma=0.001073 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.348 mas, bmaj=5.594 mas, bpa=0.7327 degrees ! Estimated noise=1472.74 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.60175 Jy ! Component: 100 - total flux cleaned = 0.740548 Jy ! Component: 150 - total flux cleaned = 0.797659 Jy ! Component: 200 - total flux cleaned = 0.83429 Jy ! Total flux subtracted in 200 components = 0.83429 Jy ! Clean residual min=-0.013675 max=0.024232 Jy/beam ! Clean residual mean=0.000051 rms=0.003352 Jy/beam ! Combined flux in latest and established models = 0.83429 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 23 model components to the UV plane model. ! The established model now contains 23 components and 0.83429 Jy ! Inverting map and beam ! Estimated beam: bmin=3.269 mas, bmaj=9.002 mas, bpa=-2.774 degrees ! Estimated noise=913.913 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-4, 4.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0550962 Jy ! Component: 100 - total flux cleaned = 0.0862238 Jy ! Component: 150 - total flux cleaned = 0.105781 Jy ! Component: 200 - total flux cleaned = 0.118578 Jy ! Total flux subtracted in 200 components = 0.118578 Jy ! Clean residual min=-0.005325 max=0.009073 Jy/beam ! Clean residual mean=-0.000167 rms=0.001943 Jy/beam ! Combined flux in latest and established models = 0.952868 Jy ! Adding 56 model components to the UV plane model. ! The established model now contains 70 components and 0.952867 Jy ! Inverting map ! Added new window around map position (3, 6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00938882 Jy ! Component: 100 - total flux cleaned = 0.0156127 Jy ! Component: 150 - total flux cleaned = 0.0192015 Jy ! Component: 200 - total flux cleaned = 0.0203559 Jy ! Total flux subtracted in 200 components = 0.0203559 Jy ! Clean residual min=-0.004422 max=0.004152 Jy/beam ! Clean residual mean=-0.000284 rms=0.001255 Jy/beam ! Combined flux in latest and established models = 0.973223 Jy ! Adding 34 model components to the UV plane model. ! The established model now contains 94 components and 0.973223 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.699 mas, bmaj=22.17 mas, bpa=-3.555 degrees ! Estimated noise=1073.09 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00377448 Jy ! Component: 100 - total flux cleaned = 0.0054335 Jy ! Component: 150 - total flux cleaned = 0.00590245 Jy ! Component: 200 - total flux cleaned = 0.00612857 Jy ! Total flux subtracted in 200 components = 0.00612857 Jy ! Clean residual min=-0.005390 max=0.004298 Jy/beam ! Clean residual mean=-0.000651 rms=0.001497 Jy/beam ! Combined flux in latest and established models = 0.979352 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.000149707 Jy ! Component: 100 - total flux cleaned = 0.000239484 Jy ! Component: 150 - total flux cleaned = 0.000239945 Jy ! Component: 200 - total flux cleaned = 0.00023987 Jy ! Total flux subtracted in 200 components = 0.00023987 Jy ! Clean residual min=-0.005367 max=0.004214 Jy/beam ! Clean residual mean=-0.000654 rms=0.001480 Jy/beam ! Combined flux in latest and established models = 0.979592 Jy ! Adding 31 model components to the UV plane model. ! The established model now contains 121 components and 0.979592 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.269 mas, bmaj=9.002 mas, bpa=-2.774 degrees ! Estimated noise=913.913 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000662921 Jy ! Component: 100 - total flux cleaned = -0.000577496 Jy ! Component: 150 - total flux cleaned = -0.000498331 Jy ! Component: 200 - total flux cleaned = -0.000495899 Jy ! Total flux subtracted in 200 components = -0.000495899 Jy ! Clean residual min=-0.004151 max=0.003368 Jy/beam ! Clean residual mean=-0.000324 rms=0.001139 Jy/beam ! Combined flux in latest and established models = 0.979096 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 134 components and 0.979096 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151285Jy sigma=0.001072 ! Fit after self-cal, rms=0.151176Jy sigma=0.001071 wmodel J1316-3338_S_2013_01_09_pus_map.mod ! Writing 134 model components to file: J1316-3338_S_2013_01_09_pus_map.mod wobs J1316-3338_S_2013_01_09_pus_uvs.fits ! Writing UV FITS file: J1316-3338_S_2013_01_09_pus_uvs.fits wwins J1316-3338_S_2013_01_09_pus_map.win ! wwins: Wrote 6 windows to J1316-3338_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.0019041 Jy ! Component: 100 - total flux cleaned = -0.00365937 Jy ! Component: 150 - total flux cleaned = -0.005314 Jy ! Component: 200 - total flux cleaned = -0.00688592 Jy ! Component: 250 - total flux cleaned = -0.00821371 Jy ! Component: 300 - total flux cleaned = -0.00926368 Jy ! Component: 350 - total flux cleaned = -0.0101139 Jy ! Component: 400 - total flux cleaned = -0.0109969 Jy ! Component: 450 - total flux cleaned = -0.0118035 Jy ! Component: 500 - total flux cleaned = -0.0125397 Jy ! Component: 550 - total flux cleaned = -0.0133624 Jy ! Component: 600 - total flux cleaned = -0.0141178 Jy ! Component: 650 - total flux cleaned = -0.0147116 Jy ! Component: 700 - total flux cleaned = -0.0153429 Jy ! Component: 750 - total flux cleaned = -0.0160115 Jy ! Component: 800 - total flux cleaned = -0.016575 Jy ! Component: 850 - total flux cleaned = -0.0171308 Jy ! Component: 900 - total flux cleaned = -0.01777 Jy ! Component: 950 - total flux cleaned = -0.018399 Jy ! Component: 1000 - total flux cleaned = -0.0187972 Jy ! Component: 1050 - total flux cleaned = -0.0193221 Jy ! Component: 1100 - total flux cleaned = -0.0198828 Jy ! Component: 1150 - total flux cleaned = -0.0203933 Jy ! Component: 1200 - total flux cleaned = -0.0207713 Jy ! Component: 1250 - total flux cleaned = -0.0211871 Jy ! Component: 1300 - total flux cleaned = -0.0215981 Jy ! Component: 1350 - total flux cleaned = -0.0220045 Jy ! Component: 1400 - total flux cleaned = -0.0223659 Jy ! Component: 1450 - total flux cleaned = -0.0226445 Jy ! Component: 1500 - total flux cleaned = -0.0230383 Jy ! Component: 1550 - total flux cleaned = -0.0233112 Jy ! Component: 1600 - total flux cleaned = -0.0236588 Jy ! Component: 1650 - total flux cleaned = -0.0239645 Jy ! Component: 1700 - total flux cleaned = -0.02423 Jy ! Component: 1750 - total flux cleaned = -0.0244929 Jy ! Component: 1800 - total flux cleaned = -0.0247912 Jy ! Component: 1850 - total flux cleaned = -0.0251236 Jy ! Component: 1900 - total flux cleaned = -0.0253066 Jy ! Component: 1950 - total flux cleaned = -0.0255972 Jy ! Component: 2000 - total flux cleaned = -0.0259569 Jy ! Total flux subtracted in 2048 components = -0.026135 Jy ! Clean residual min=-0.001902 max=0.002120 Jy/beam ! Clean residual mean=-0.000121 rms=0.000747 Jy/beam ! Combined flux in latest and established models = 0.952961 Jy keep ! Adding 133 model components to the UV plane model. ! The established model now contains 267 components and 0.952961 Jy ![Exited script file: muppet] wmap "J1316-3338_S_2013_01_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.269 x 9.002 at -2.774 degrees (North through East) ! Clean map min=-0.0039456 max=0.75711 Jy/beam ! Writing clean map to FITS file: J1316-3338_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.757106 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000745752 device "J1316-3338_S_2013_01_09_pus\_map.eps/VPS" ! Attempting to open device: 'J1316-3338_S_2013_01_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.295501 0.295501 0.591002 1.182 2.36401 4.72801 9.45603 18.9121 37.8241 75.6482 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0022 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.757 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.269 x 9.002 at -2.774 degrees (North through East) ! Clean map min=-0.0039456 max=0.75711 Jy/beam quit ! Quitting program ! Log file J1316-3338_S_2013_01_09_pus_uvs.log closed on Tue Nov 4 19:01:44 2014