! Started logfile: J1130-1449_S_2000_03_13_pus_uvs.log on Tue Apr 2 23:55:16 2013 obs J1130-1449_S_2000_03_13_pus_uva.fits ! Reading UV FITS file: J1130-1449_S_2000_03_13_pus_uva.fits ! AN table 1: 274 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.0246831 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1130-1449 ! ! 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.22497e+09 8e+06 1 8e+06 ! 02 2 2.24497e+09 8e+06 1 8e+06 ! 03 3 2.33497e+09 8e+06 1 8e+06 ! 04 4 2.36497e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 2251 lines of history. ! ! Reading 5140 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 J1130-1449_S_2000_03_13_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 4692 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 4779 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 4716 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 4776 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=2.626431Jy sigma=0.008590 ! Fit after self-cal, rms=2.502380Jy sigma=0.008249 ! 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=1.587 mas, bmaj=2.469 mas, bpa=7.98 degrees ! Estimated noise=4377.78 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.4654 Jy ! Component: 100 - total flux cleaned = 1.9227 Jy ! Total flux subtracted in 100 components = 1.9227 Jy ! Clean residual min=-0.189335 max=0.265152 Jy/beam ! Clean residual mean=-0.004771 rms=0.076851 Jy/beam ! Combined flux in latest and established models = 1.9227 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 1.9227 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.895985Jy sigma=0.006118 ! Fit after self-cal, rms=1.895769Jy sigma=0.006118 ! 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=2.065 mas, bmaj=3.314 mas, bpa=17.04 degrees ! Estimated noise=3207.26 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 = 1.0847 Jy ! Component: 100 - total flux cleaned = 1.61091 Jy ! Component: 150 - total flux cleaned = 1.87367 Jy ! Component: 200 - total flux cleaned = 2.00493 Jy ! Total flux subtracted in 200 components = 2.00493 Jy ! Clean residual min=-0.253420 max=0.306531 Jy/beam ! Clean residual mean=0.001151 rms=0.111900 Jy/beam ! Combined flux in latest and established models = 3.92763 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 16 components and 3.92763 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.103305Jy sigma=0.003711 ! Fit after self-cal, rms=1.102400Jy sigma=0.003706 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 447 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* GC 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.01 WF 1.22 HH 1.03 KK 1.00* ! MA 1.00* MC 1.06 MK 1.00* NY 1.00* ! ON 0.85 TS 1.00* WZ 0.89 HN 1.30 ! ! ! Correcting IF 2. ! A total of 372 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* GC 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.02 WF 1.27 HH 0.98 KK 1.00* ! MA 1.00* MC 0.85 MK 1.00* NY 1.00* ! ON 0.87 TS 1.00* WZ 0.90 HN 1.28 ! ! ! Correcting IF 3. ! A total of 435 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* GC 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.04 WF 1.61 HH 1.01 KK 1.00* ! MA 1.00* MC 1.54 MK 1.00* NY 1.00* ! ON 0.88 TS 1.00* WZ 0.88 HN 1.29 ! ! ! Correcting IF 4. ! A total of 363 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* GC 1.00* KP 1.00* ! LA 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.02 WF 1.23 HH 1.00 KK 1.00* ! MA 1.00* MC 1.29 MK 1.00* NY 1.00* ! ON 0.86 TS 1.00* WZ 0.90 HN 1.30 ! ! ! Fit before self-cal, rms=1.102400Jy sigma=0.003706 ! Fit after self-cal, rms=0.917348Jy sigma=0.001925 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.834 mas, bmaj=2.868 mas, bpa=9.757 degrees ! Estimated noise=4202.56 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.0150971 Jy ! Component: 100 - total flux cleaned = -0.0238203 Jy ! Component: 150 - total flux cleaned = -0.0323386 Jy ! Component: 200 - total flux cleaned = -0.038727 Jy ! Total flux subtracted in 200 components = -0.038727 Jy ! Clean residual min=-0.075803 max=0.069427 Jy/beam ! Clean residual mean=0.000323 rms=0.020591 Jy/beam ! Combined flux in latest and established models = 3.8889 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 17 components and 3.8889 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.913728Jy sigma=0.001912 ! Fit after self-cal, rms=0.919305Jy sigma=0.001900 ! 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.919305Jy sigma=0.001900 ! Fit after self-cal, rms=0.922241Jy sigma=0.001898 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.839 mas, bmaj=2.88 mas, bpa=9.379 degrees ! Estimated noise=4216.31 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.017043 Jy ! Component: 100 - total flux cleaned = -0.0210631 Jy ! Component: 150 - total flux cleaned = -0.0251187 Jy ! Component: 200 - total flux cleaned = -0.0286819 Jy ! Total flux subtracted in 200 components = -0.0286819 Jy ! Clean residual min=-0.067591 max=0.063072 Jy/beam ! Clean residual mean=0.000248 rms=0.018631 Jy/beam ! Combined flux in latest and established models = 3.86022 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 18 components and 3.86022 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.919815Jy sigma=0.001894 ! Fit after self-cal, rms=0.920866Jy sigma=0.001892 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.00756337 Jy ! Component: 100 - total flux cleaned = -0.00592415 Jy ! Component: 150 - total flux cleaned = -0.00545654 Jy ! Component: 200 - total flux cleaned = -0.00584811 Jy ! Total flux subtracted in 200 components = -0.00584811 Jy ! Clean residual min=-0.064779 max=0.061303 Jy/beam ! Clean residual mean=0.000209 rms=0.018418 Jy/beam ! Combined flux in latest and established models = 3.85437 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 19 components and 3.85437 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.919363Jy sigma=0.001892 ! Fit after self-cal, rms=0.920373Jy sigma=0.001891 ! 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.920373Jy sigma=0.001891 ! Fit after self-cal, rms=0.910541Jy sigma=0.001889 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.843 mas, bmaj=2.889 mas, bpa=9.222 degrees ! Estimated noise=4223.94 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0137029 Jy ! Component: 100 - total flux cleaned = -0.0147095 Jy ! Component: 150 - total flux cleaned = -0.0151551 Jy ! Component: 200 - total flux cleaned = -0.0159448 Jy ! Total flux subtracted in 200 components = -0.0159448 Jy ! Clean residual min=-0.063682 max=0.062729 Jy/beam ! Clean residual mean=0.000246 rms=0.018009 Jy/beam ! Combined flux in latest and established models = 3.83843 Jy selfcal ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 21 components and 3.83843 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.908945Jy sigma=0.001889 ! Fit after self-cal, rms=0.909642Jy sigma=0.001888 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.00141849 Jy ! Component: 100 - total flux cleaned = -0.000587568 Jy ! Component: 150 - total flux cleaned = -0.000247052 Jy ! Component: 200 - total flux cleaned = -0.000555962 Jy ! Total flux subtracted in 200 components = -0.000555962 Jy ! Clean residual min=-0.062690 max=0.061925 Jy/beam ! Clean residual mean=0.000242 rms=0.017905 Jy/beam ! Combined flux in latest and established models = 3.83787 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 22 components and 3.83787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.908356Jy sigma=0.001889 ! Fit after self-cal, rms=0.909344Jy sigma=0.001887 ! 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.909344Jy sigma=0.001887 ! Fit after self-cal, rms=0.905280Jy sigma=0.001887 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.3 mas, bmaj=18.66 mas, bpa=-82.18 degrees ! Estimated noise=7103.93 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019259 Jy ! Component: 100 - total flux cleaned = 0.00187485 Jy ! Component: 150 - total flux cleaned = 0.00185037 Jy ! Component: 200 - total flux cleaned = 0.0018363 Jy ! Total flux subtracted in 200 components = 0.0018363 Jy ! Clean residual min=-0.058232 max=0.049896 Jy/beam ! Clean residual mean=0.001202 rms=0.021016 Jy/beam ! Combined flux in latest and established models = 3.83971 Jy selfcal ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 22 components and 3.83971 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.904677Jy sigma=0.001889 ! Fit after self-cal, rms=0.905169Jy sigma=0.001887 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 = -6.00601e-05 Jy ! Component: 100 - total flux cleaned = -6.52895e-05 Jy ! Component: 150 - total flux cleaned = -5.85847e-05 Jy ! Component: 200 - total flux cleaned = -5.04098e-05 Jy ! Total flux subtracted in 200 components = -5.04098e-05 Jy ! Clean residual min=-0.057798 max=0.049464 Jy/beam ! Clean residual mean=0.001201 rms=0.021001 Jy/beam ! Combined flux in latest and established models = 3.83966 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 22 components and 3.83966 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.904861Jy sigma=0.001889 ! Fit after self-cal, rms=0.905411Jy sigma=0.001888 ! 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=1.846 mas, bmaj=2.898 mas, bpa=8.986 degrees ! Estimated noise=4228.34 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.846 x 2.898 at 8.986 degrees (North through East) ! Clean map min=-0.057261 max=2.4044 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.905411Jy sigma=0.001888 ! Fit after self-cal, rms=0.931815Jy sigma=0.001532 ! Inverting map and beam ! Estimated beam: bmin=1.842 mas, bmaj=2.898 mas, bpa=8.736 degrees ! Estimated noise=4204.33 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00512053 Jy ! Component: 100 - total flux cleaned = -0.0044119 Jy ! Component: 150 - total flux cleaned = -0.00324736 Jy ! Component: 200 - total flux cleaned = -0.00172248 Jy ! Total flux subtracted in 200 components = -0.00172248 Jy ! Clean residual min=-0.053087 max=0.052652 Jy/beam ! Clean residual mean=0.000148 rms=0.014659 Jy/beam ! Combined flux in latest and established models = 3.83794 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 24 components and 3.83794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.929706Jy sigma=0.001531 ! Fit after self-cal, rms=0.930074Jy sigma=0.001531 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0024984 Jy ! Component: 100 - total flux cleaned = 0.00391715 Jy ! Component: 150 - total flux cleaned = 0.00432885 Jy ! Component: 200 - total flux cleaned = 0.00472134 Jy ! Total flux subtracted in 200 components = 0.00472134 Jy ! Clean residual min=-0.052043 max=0.051155 Jy/beam ! Clean residual mean=0.000178 rms=0.014375 Jy/beam ! Combined flux in latest and established models = 3.84266 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 25 components and 3.84266 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.928898Jy sigma=0.001531 ! Fit after self-cal, rms=0.929190Jy sigma=0.001531 ! 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.929190Jy sigma=0.001531 ! Fit after self-cal, rms=0.930513Jy sigma=0.001530 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=1.617 mas, bmaj=2.436 mas, bpa=5.8 degrees ! Estimated noise=4901.14 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.57757 Jy ! Component: 100 - total flux cleaned = 2.0513 Jy ! Component: 150 - total flux cleaned = 2.3586 Jy ! Component: 200 - total flux cleaned = 2.59727 Jy ! Total flux subtracted in 200 components = 2.59727 Jy ! Clean residual min=-0.129159 max=0.156389 Jy/beam ! Clean residual mean=-0.003186 rms=0.051484 Jy/beam ! Combined flux in latest and established models = 2.59727 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 17 model components to the UV plane model. ! The established model now contains 17 components and 2.59727 Jy ! Inverting map and beam ! Estimated beam: bmin=1.844 mas, bmaj=2.909 mas, bpa=8.345 degrees ! Estimated noise=4209.65 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.447101 Jy ! Component: 100 - total flux cleaned = 0.726597 Jy ! Component: 150 - total flux cleaned = 0.902879 Jy ! Component: 200 - total flux cleaned = 1.01443 Jy ! Total flux subtracted in 200 components = 1.01443 Jy ! Clean residual min=-0.073872 max=0.104324 Jy/beam ! Clean residual mean=-0.001298 rms=0.026918 Jy/beam ! Combined flux in latest and established models = 3.6117 Jy ! Adding 15 model components to the UV plane model. ! The established model now contains 21 components and 3.6117 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=10.3 mas, bmaj=18.6 mas, bpa=-82.17 degrees ! Estimated noise=7142.34 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.174464 Jy ! Component: 100 - total flux cleaned = 0.218765 Jy ! Component: 150 - total flux cleaned = 0.230605 Jy ! Component: 200 - total flux cleaned = 0.230867 Jy ! Total flux subtracted in 200 components = 0.230867 Jy ! Clean residual min=-0.053559 max=0.056932 Jy/beam ! Clean residual mean=0.000758 rms=0.025061 Jy/beam ! Combined flux in latest and established models = 3.84256 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 = 5.8316e-06 Jy ! Component: 100 - total flux cleaned = 0.00025791 Jy ! Component: 150 - total flux cleaned = 0.000507477 Jy ! Component: 200 - total flux cleaned = 0.000745267 Jy ! Total flux subtracted in 200 components = 0.000745267 Jy ! Clean residual min=-0.052189 max=0.056067 Jy/beam ! Clean residual mean=0.000775 rms=0.024734 Jy/beam ! Combined flux in latest and established models = 3.84331 Jy ! Adding 8 model components to the UV plane model. ! The established model now contains 24 components and 3.84331 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.844 mas, bmaj=2.909 mas, bpa=8.345 degrees ! Estimated noise=4209.65 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0441874 Jy ! Component: 100 - total flux cleaned = -0.04541 Jy ! Component: 150 - total flux cleaned = -0.046528 Jy ! Component: 200 - total flux cleaned = -0.0465457 Jy ! Total flux subtracted in 200 components = -0.0465457 Jy ! Clean residual min=-0.066071 max=0.079790 Jy/beam ! Clean residual mean=-0.000088 rms=0.017643 Jy/beam ! Combined flux in latest and established models = 3.79676 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 25 components and 3.79676 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.934176Jy sigma=0.001553 ! Fit after self-cal, rms=0.934985Jy sigma=0.001547 wmodel J1130-1449_S_2000_03_13_pus_map.mod ! Writing 25 model components to file: J1130-1449_S_2000_03_13_pus_map.mod wobs J1130-1449_S_2000_03_13_pus_uvs.fits ! Writing UV FITS file: J1130-1449_S_2000_03_13_pus_uvs.fits wwins J1130-1449_S_2000_03_13_pus_map.win ! wwins: Wrote 1 windows to J1130-1449_S_2000_03_13_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.010237 Jy ! Component: 100 - total flux cleaned = 0.00810139 Jy ! Component: 150 - total flux cleaned = 0.0088157 Jy ! Component: 200 - total flux cleaned = 0.00729678 Jy ! Component: 250 - total flux cleaned = 0.00726133 Jy ! Component: 300 - total flux cleaned = 0.00725141 Jy ! Component: 350 - total flux cleaned = 0.00783661 Jy ! Component: 400 - total flux cleaned = 0.00886704 Jy ! Component: 450 - total flux cleaned = 0.0102733 Jy ! Component: 500 - total flux cleaned = 0.0115991 Jy ! Component: 550 - total flux cleaned = 0.0115958 Jy ! Component: 600 - total flux cleaned = 0.0119637 Jy ! Component: 650 - total flux cleaned = 0.0123233 Jy ! Component: 700 - total flux cleaned = 0.0144134 Jy ! Component: 750 - total flux cleaned = 0.0167197 Jy ! Component: 800 - total flux cleaned = 0.0195382 Jy ! Component: 850 - total flux cleaned = 0.022255 Jy ! Component: 900 - total flux cleaned = 0.0242718 Jy ! Component: 950 - total flux cleaned = 0.025938 Jy ! Component: 1000 - total flux cleaned = 0.0280887 Jy ! Component: 1050 - total flux cleaned = 0.0309545 Jy ! Component: 1100 - total flux cleaned = 0.0339788 Jy ! Component: 1150 - total flux cleaned = 0.0361837 Jy ! Component: 1200 - total flux cleaned = 0.0385658 Jy ! Component: 1250 - total flux cleaned = 0.0399626 Jy ! Component: 1300 - total flux cleaned = 0.0420158 Jy ! Component: 1350 - total flux cleaned = 0.0429061 Jy ! Component: 1400 - total flux cleaned = 0.0435634 Jy ! Component: 1450 - total flux cleaned = 0.0444255 Jy ! Component: 1500 - total flux cleaned = 0.0450604 Jy ! Component: 1550 - total flux cleaned = 0.0456858 Jy ! Component: 1600 - total flux cleaned = 0.046302 Jy ! Component: 1650 - total flux cleaned = 0.0467071 Jy ! Component: 1700 - total flux cleaned = 0.0475095 Jy ! Component: 1750 - total flux cleaned = 0.0482999 Jy ! Component: 1800 - total flux cleaned = 0.0484956 Jy ! Component: 1850 - total flux cleaned = 0.0492679 Jy ! Component: 1900 - total flux cleaned = 0.0498409 Jy ! Component: 1950 - total flux cleaned = 0.0504089 Jy ! Component: 2000 - total flux cleaned = 0.0505984 Jy ! Total flux subtracted in 2048 components = 0.0509732 Jy ! Clean residual min=-0.010879 max=0.010172 Jy/beam ! Clean residual mean=0.000161 rms=0.003283 Jy/beam ! Combined flux in latest and established models = 3.84774 Jy keep ! Adding 117 model components to the UV plane model. ! The established model now contains 142 components and 3.84774 Jy ![Exited script file: muppet] wmap "J1130-1449_S_2000_03_13_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.844 x 2.909 at 8.345 degrees (North through East) ! Clean map min=-0.024376 max=2.453 Jy/beam ! Writing clean map to FITS file: J1130-1449_S_2000_03_13_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 2.45298 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00322824 device "J1130-1449_S_2000_03_13_pus\_map.eps/VPS" ! Attempting to open device: 'J1130-1449_S_2000_03_13_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.394814 0.394814 0.789629 1.57926 3.15852 6.31703 12.6341 25.2681 50.5363 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0096 print "peak=",int(1000*peak_flux)/1000 ! peak= 2.452 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.844 x 2.909 at 8.345 degrees (North through East) ! Clean map min=-0.024376 max=2.453 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J1130-1449_S_2000_03_13_pus_uvs.log closed on Tue Apr 2 23:55:22 2013