! Started logfile: J1733-1304_S_2017_10_21_pus_uvs.log on Sat May 11 21:03:10 2019 obs J1733-1304_S_2017_10_21_pus_uva.fits ! Reading UV FITS file: J1733-1304_S_2017_10_21_pus_uva.fits ! AN table 1: 35 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.511111 visibilities/baseline/integration-bin. ! Found source: J1733-1304 ! ! 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.23588e+09 3.2e+07 1 3.2e+07 ! 02 2 2.26788e+09 3.2e+07 1 3.2e+07 ! 03 3 2.29988e+09 3.2e+07 1 3.2e+07 ! 04 4 2.36388e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 644 lines of history. ! ! Reading 2576 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 = 1024 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J1733-1304_S_2017_10_21_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 = 1024x1024 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 148 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 168 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 148 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 230 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=3.298741Jy sigma=0.027302 ! Fit after self-cal, rms=3.268122Jy sigma=0.027150 ! 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.791 mas, bmaj=14.19 mas, bpa=-17.19 degrees ! Estimated noise=2478.79 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 = 2.83979 Jy ! Component: 100 - total flux cleaned = 3.64141 Jy ! Total flux subtracted in 100 components = 3.64141 Jy ! Clean residual min=-0.203436 max=0.361078 Jy/beam ! Clean residual mean=0.001119 rms=0.065793 Jy/beam ! Combined flux in latest and established models = 3.64141 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 12 components and 3.64141 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.227080Jy sigma=0.009670 ! Fit after self-cal, rms=1.227238Jy sigma=0.009670 ! 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.499 mas, bmaj=16.09 mas, bpa=-17.09 degrees ! Estimated noise=1986.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) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.644227 Jy ! Component: 100 - total flux cleaned = 1.01721 Jy ! Component: 150 - total flux cleaned = 1.26589 Jy ! Component: 200 - total flux cleaned = 1.43828 Jy ! Total flux subtracted in 200 components = 1.43828 Jy ! Clean residual min=-0.145434 max=0.249521 Jy/beam ! Clean residual mean=0.000721 rms=0.046686 Jy/beam ! Combined flux in latest and established models = 5.07968 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 71 components and 5.07968 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.806866Jy sigma=0.005589 ! Fit after self-cal, rms=0.802140Jy sigma=0.005548 ! Inverting map ! Added new window around map position (-2, 25). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.237191 Jy ! Component: 100 - total flux cleaned = 0.367337 Jy ! Component: 150 - total flux cleaned = 0.467215 Jy ! Component: 200 - total flux cleaned = 0.547801 Jy ! Total flux subtracted in 200 components = 0.547801 Jy ! Clean residual min=-0.111810 max=0.128158 Jy/beam ! Clean residual mean=0.000294 rms=0.031338 Jy/beam ! Combined flux in latest and established models = 5.62748 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 113 components and 5.62749 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.715220Jy sigma=0.004755 ! Fit after self-cal, rms=0.707913Jy sigma=0.004688 ! 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.16 FD 0.94 HN 1.01 KP 1.07 ! LA 1.00 MK 0.99 NL 0.89 OV 0.98 ! PT 0.99 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 0.99 HN 1.06 KP 1.04 ! LA 0.97 MK 1.00 NL 1.00 OV 1.00 ! PT 0.97 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.05 FD 1.00 HN 1.04 KP 1.04 ! LA 0.76 MK 1.02 NL 1.21 OV 0.92 ! PT 0.96 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.11 HN 1.03 KP 1.00* ! LA 1.01 MK 1.19 NL 1.00* OV 1.00* ! PT 1.02 ! ! ! Fit before self-cal, rms=0.707913Jy sigma=0.004688 ! Fit after self-cal, rms=0.504456Jy sigma=0.003004 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.544 mas, bmaj=16.38 mas, bpa=-17.12 degrees ! Estimated noise=1998.77 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.0255566 Jy ! Component: 100 - total flux cleaned = 0.0249588 Jy ! Component: 150 - total flux cleaned = 0.0286451 Jy ! Component: 200 - total flux cleaned = 0.0383656 Jy ! Total flux subtracted in 200 components = 0.0383656 Jy ! Clean residual min=-0.070648 max=0.090905 Jy/beam ! Clean residual mean=0.000080 rms=0.019741 Jy/beam ! Combined flux in latest and established models = 5.66585 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 147 components and 5.66585 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.481619Jy sigma=0.002809 ! Fit after self-cal, rms=0.478438Jy sigma=0.002784 ! 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.478438Jy sigma=0.002784 ! Fit after self-cal, rms=0.225923Jy sigma=0.001476 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.55 mas, bmaj=16.25 mas, bpa=-17.08 degrees ! Estimated noise=1992.49 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0276668 Jy ! Component: 100 - total flux cleaned = 0.0355745 Jy ! Component: 150 - total flux cleaned = 0.0376178 Jy ! Component: 200 - total flux cleaned = 0.0385084 Jy ! Total flux subtracted in 200 components = 0.0385084 Jy ! Clean residual min=-0.036389 max=0.041114 Jy/beam ! Clean residual mean=-0.000038 rms=0.008225 Jy/beam ! Combined flux in latest and established models = 5.70436 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 173 components and 5.70436 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210173Jy sigma=0.001327 ! Fit after self-cal, rms=0.208168Jy sigma=0.001312 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) ! Added new window around map position (22, 0.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0141177 Jy ! Component: 100 - total flux cleaned = 0.0174689 Jy ! Component: 150 - total flux cleaned = 0.017503 Jy ! Component: 200 - total flux cleaned = 0.0183318 Jy ! Total flux subtracted in 200 components = 0.0183318 Jy ! Clean residual min=-0.034373 max=0.026727 Jy/beam ! Clean residual mean=-0.000026 rms=0.006763 Jy/beam ! Combined flux in latest and established models = 5.72269 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 195 components and 5.72269 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.194909Jy sigma=0.001205 ! Fit after self-cal, rms=0.194250Jy sigma=0.001200 ! 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.194250Jy sigma=0.001200 ! Fit after self-cal, rms=0.192370Jy sigma=0.001186 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.552 mas, bmaj=16.26 mas, bpa=-17.07 degrees ! Estimated noise=1988.92 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00292641 Jy ! Component: 100 - total flux cleaned = 0.00379607 Jy ! Component: 150 - total flux cleaned = 0.0044754 Jy ! Component: 200 - total flux cleaned = 0.00320161 Jy ! Total flux subtracted in 200 components = 0.00320161 Jy ! Clean residual min=-0.029060 max=0.021239 Jy/beam ! Clean residual mean=-0.000012 rms=0.005835 Jy/beam ! Combined flux in latest and established models = 5.72589 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 217 components and 5.72589 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185844Jy sigma=0.001135 ! Fit after self-cal, rms=0.185266Jy sigma=0.001132 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.00183657 Jy ! Component: 100 - total flux cleaned = 0.00306343 Jy ! Component: 150 - total flux cleaned = 0.00362536 Jy ! Component: 200 - total flux cleaned = 0.0030865 Jy ! Total flux subtracted in 200 components = 0.0030865 Jy ! Clean residual min=-0.026617 max=0.018517 Jy/beam ! Clean residual mean=-0.000009 rms=0.005380 Jy/beam ! Combined flux in latest and established models = 5.72898 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 236 components and 5.72898 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.181061Jy sigma=0.001100 ! Fit after self-cal, rms=0.180680Jy sigma=0.001098 ! 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.180680Jy sigma=0.001098 ! Fit after self-cal, rms=0.180392Jy sigma=0.001095 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.39 mas, bmaj=19.52 mas, bpa=-6.71 degrees ! Estimated noise=2347.97 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00347538 Jy ! Component: 100 - total flux cleaned = 0.00304642 Jy ! Component: 150 - total flux cleaned = 0.00162366 Jy ! Component: 200 - total flux cleaned = 0.000777876 Jy ! Total flux subtracted in 200 components = 0.000777876 Jy ! Clean residual min=-0.022847 max=0.016081 Jy/beam ! Clean residual mean=-0.000005 rms=0.005376 Jy/beam ! Combined flux in latest and established models = 5.72976 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 260 components and 5.72976 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.178444Jy sigma=0.001076 ! Fit after self-cal, rms=0.178072Jy sigma=0.001073 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.00103285 Jy ! Component: 100 - total flux cleaned = 0.000197548 Jy ! Component: 150 - total flux cleaned = 0.000559705 Jy ! Component: 200 - total flux cleaned = 0.000222234 Jy ! Total flux subtracted in 200 components = 0.000222234 Jy ! Clean residual min=-0.021118 max=0.015373 Jy/beam ! Clean residual mean=-0.000002 rms=0.005061 Jy/beam ! Combined flux in latest and established models = 5.72998 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 279 components and 5.72998 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176882Jy sigma=0.001063 ! Fit after self-cal, rms=0.176665Jy sigma=0.001061 ! 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.554 mas, bmaj=16.26 mas, bpa=-17.07 degrees ! Estimated noise=1987.59 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.554 x 16.26 at -17.07 degrees (North through East) ! Clean map min=-0.033528 max=3.7635 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.176665Jy sigma=0.001061 ! Fit after self-cal, rms=0.129962Jy sigma=0.000860 ! Inverting map and beam ! Estimated beam: bmin=3.554 mas, bmaj=16.28 mas, bpa=-17.07 degrees ! Estimated noise=1987.43 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00663424 Jy ! Component: 100 - total flux cleaned = -0.00814678 Jy ! Component: 150 - total flux cleaned = -0.00906695 Jy ! Component: 200 - total flux cleaned = -0.00950184 Jy ! Total flux subtracted in 200 components = -0.00950184 Jy ! Clean residual min=-0.022118 max=0.015920 Jy/beam ! Clean residual mean=0.000004 rms=0.004333 Jy/beam ! Combined flux in latest and established models = 5.72048 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 293 components and 5.72048 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125597Jy sigma=0.000834 ! Fit after self-cal, rms=0.125277Jy sigma=0.000832 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000330219 Jy ! Component: 100 - total flux cleaned = 0.000313047 Jy ! Component: 150 - total flux cleaned = -0.000821871 Jy ! Component: 200 - total flux cleaned = -0.00189798 Jy ! Total flux subtracted in 200 components = -0.00189798 Jy ! Clean residual min=-0.020994 max=0.015771 Jy/beam ! Clean residual mean=0.000006 rms=0.004135 Jy/beam ! Combined flux in latest and established models = 5.71858 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 309 components and 5.71858 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122393Jy sigma=0.000815 ! Fit after self-cal, rms=0.122210Jy sigma=0.000814 ! 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.122210Jy sigma=0.000814 ! Fit after self-cal, rms=0.121781Jy sigma=0.000812 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.821 mas, bmaj=14.35 mas, bpa=-17.19 degrees ! Estimated noise=2515.08 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 2.87876 Jy ! Component: 100 - total flux cleaned = 3.66487 Jy ! Component: 150 - total flux cleaned = 4.06413 Jy ! Component: 200 - total flux cleaned = 4.33272 Jy ! Total flux subtracted in 200 components = 4.33272 Jy ! Clean residual min=-0.078047 max=0.152635 Jy/beam ! Clean residual mean=0.000593 rms=0.030007 Jy/beam ! Combined flux in latest and established models = 4.33271 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 54 model components to the UV plane model. ! The established model now contains 54 components and 4.33271 Jy ! Inverting map and beam ! Estimated beam: bmin=3.554 mas, bmaj=16.29 mas, bpa=-17.07 degrees ! Estimated noise=1986.76 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.346165 Jy ! Component: 100 - total flux cleaned = 0.571739 Jy ! Component: 150 - total flux cleaned = 0.733442 Jy ! Component: 200 - total flux cleaned = 0.854018 Jy ! Total flux subtracted in 200 components = 0.854018 Jy ! Clean residual min=-0.059907 max=0.069541 Jy/beam ! Clean residual mean=0.000316 rms=0.020480 Jy/beam ! Combined flux in latest and established models = 5.18673 Jy ! Adding 143 model components to the UV plane model. ! The established model now contains 187 components and 5.18673 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=8.391 mas, bmaj=19.55 mas, bpa=-6.741 degrees ! Estimated noise=2346.74 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.136293 Jy ! Component: 100 - total flux cleaned = 0.226797 Jy ! Component: 150 - total flux cleaned = 0.291985 Jy ! Component: 200 - total flux cleaned = 0.340455 Jy ! Total flux subtracted in 200 components = 0.340455 Jy ! Clean residual min=-0.031861 max=0.034892 Jy/beam ! Clean residual mean=0.000133 rms=0.011302 Jy/beam ! Combined flux in latest and established models = 5.52719 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.0372842 Jy ! Component: 100 - total flux cleaned = 0.0668006 Jy ! Component: 150 - total flux cleaned = 0.0906095 Jy ! Component: 200 - total flux cleaned = 0.110071 Jy ! Total flux subtracted in 200 components = 0.110071 Jy ! Clean residual min=-0.021912 max=0.020422 Jy/beam ! Clean residual mean=0.000044 rms=0.006276 Jy/beam ! Combined flux in latest and established models = 5.63726 Jy ! Adding 179 model components to the UV plane model. ! The established model now contains 362 components and 5.63726 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.554 mas, bmaj=16.29 mas, bpa=-17.07 degrees ! Estimated noise=1986.76 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0127464 Jy ! Component: 100 - total flux cleaned = 0.0181802 Jy ! Component: 150 - total flux cleaned = 0.0236801 Jy ! Component: 200 - total flux cleaned = 0.0280964 Jy ! Total flux subtracted in 200 components = 0.0280964 Jy ! Clean residual min=-0.022826 max=0.021628 Jy/beam ! Clean residual mean=0.000009 rms=0.005243 Jy/beam ! Combined flux in latest and established models = 5.66536 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 394 components and 5.66536 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134620Jy sigma=0.000911 ! Fit after self-cal, rms=0.132760Jy sigma=0.000895 wmodel J1733-1304_S_2017_10_21_pus_map.mod ! Writing 394 model components to file: J1733-1304_S_2017_10_21_pus_map.mod wobs J1733-1304_S_2017_10_21_pus_uvs.fits ! Writing UV FITS file: J1733-1304_S_2017_10_21_pus_uvs.fits wwins J1733-1304_S_2017_10_21_pus_map.win ! wwins: Wrote 3 windows to J1733-1304_S_2017_10_21_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=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00636956 Jy ! Component: 100 - total flux cleaned = -0.0077689 Jy ! Component: 150 - total flux cleaned = -0.00710443 Jy ! Component: 200 - total flux cleaned = -0.0074237 Jy ! Component: 250 - total flux cleaned = -0.00652605 Jy ! Component: 300 - total flux cleaned = -0.00566301 Jy ! Component: 350 - total flux cleaned = -0.00429033 Jy ! Component: 400 - total flux cleaned = -0.00376077 Jy ! Component: 450 - total flux cleaned = -0.00299134 Jy ! Component: 500 - total flux cleaned = -0.00373139 Jy ! Component: 550 - total flux cleaned = -0.00469424 Jy ! Component: 600 - total flux cleaned = -0.00586829 Jy ! Component: 650 - total flux cleaned = -0.0067817 Jy ! Component: 700 - total flux cleaned = -0.00855497 Jy ! Component: 750 - total flux cleaned = -0.00985161 Jy ! Component: 800 - total flux cleaned = -0.0113249 Jy ! Component: 850 - total flux cleaned = -0.0127636 Jy ! Component: 900 - total flux cleaned = -0.014175 Jy ! Component: 950 - total flux cleaned = -0.0155461 Jy ! Component: 1000 - total flux cleaned = -0.0167038 Jy ! Component: 1050 - total flux cleaned = -0.0178404 Jy ! Component: 1100 - total flux cleaned = -0.0191346 Jy ! Component: 1150 - total flux cleaned = -0.019862 Jy ! Component: 1200 - total flux cleaned = -0.0202233 Jy ! Component: 1250 - total flux cleaned = -0.0209212 Jy ! Component: 1300 - total flux cleaned = -0.0216077 Jy ! Component: 1350 - total flux cleaned = -0.0217749 Jy ! Component: 1400 - total flux cleaned = -0.0222723 Jy ! Component: 1450 - total flux cleaned = -0.022434 Jy ! Component: 1500 - total flux cleaned = -0.0230713 Jy ! Component: 1550 - total flux cleaned = -0.0232263 Jy ! Component: 1600 - total flux cleaned = -0.0233781 Jy ! Component: 1650 - total flux cleaned = -0.0238346 Jy ! Component: 1700 - total flux cleaned = -0.024283 Jy ! Component: 1750 - total flux cleaned = -0.0248692 Jy ! Component: 1800 - total flux cleaned = -0.0251575 Jy ! Component: 1850 - total flux cleaned = -0.0254443 Jy ! Component: 1900 - total flux cleaned = -0.0258664 Jy ! Component: 1950 - total flux cleaned = -0.0265572 Jy ! Component: 2000 - total flux cleaned = -0.0265574 Jy ! Component: 2050 - total flux cleaned = -0.0268291 Jy ! Component: 2100 - total flux cleaned = -0.0272238 Jy ! Component: 2150 - total flux cleaned = -0.0274843 Jy ! Component: 2200 - total flux cleaned = -0.0276111 Jy ! Component: 2250 - total flux cleaned = -0.0282478 Jy ! Component: 2300 - total flux cleaned = -0.0284994 Jy ! Component: 2350 - total flux cleaned = -0.0288688 Jy ! Component: 2400 - total flux cleaned = -0.0294765 Jy ! Component: 2450 - total flux cleaned = -0.0297151 Jy ! Component: 2500 - total flux cleaned = -0.0301916 Jy ! Component: 2550 - total flux cleaned = -0.0304251 Jy ! Component: 2600 - total flux cleaned = -0.0305432 Jy ! Component: 2650 - total flux cleaned = -0.0313417 Jy ! Component: 2700 - total flux cleaned = -0.031567 Jy ! Component: 2750 - total flux cleaned = -0.0317897 Jy ! Component: 2800 - total flux cleaned = -0.0319015 Jy ! Component: 2850 - total flux cleaned = -0.0324435 Jy ! Component: 2900 - total flux cleaned = -0.0326558 Jy ! Component: 2950 - total flux cleaned = -0.0331824 Jy ! Component: 3000 - total flux cleaned = -0.0330767 Jy ! Component: 3050 - total flux cleaned = -0.0333888 Jy ! Component: 3100 - total flux cleaned = -0.0336967 Jy ! Component: 3150 - total flux cleaned = -0.0337981 Jy ! Component: 3200 - total flux cleaned = -0.0342978 Jy ! Component: 3250 - total flux cleaned = -0.0341994 Jy ! Component: 3300 - total flux cleaned = -0.0345916 Jy ! Component: 3350 - total flux cleaned = -0.0347856 Jy ! Component: 3400 - total flux cleaned = -0.0348821 Jy ! Component: 3450 - total flux cleaned = -0.0349765 Jy ! Component: 3500 - total flux cleaned = -0.0348835 Jy ! Component: 3550 - total flux cleaned = -0.0354416 Jy ! Component: 3600 - total flux cleaned = -0.0349817 Jy ! Component: 3650 - total flux cleaned = -0.0350744 Jy ! Component: 3700 - total flux cleaned = -0.0351651 Jy ! Component: 3750 - total flux cleaned = -0.0351648 Jy ! Component: 3800 - total flux cleaned = -0.0349007 Jy ! Component: 3850 - total flux cleaned = -0.0350762 Jy ! Component: 3900 - total flux cleaned = -0.0346408 Jy ! Component: 3950 - total flux cleaned = -0.0347278 Jy ! Component: 4000 - total flux cleaned = -0.0347272 Jy ! Component: 4050 - total flux cleaned = -0.034559 Jy ! Total flux subtracted in 4096 components = -0.0346421 Jy ! Clean residual min=-0.004175 max=0.004571 Jy/beam ! Clean residual mean=0.000022 rms=0.001633 Jy/beam ! Combined flux in latest and established models = 5.63072 Jy keep ! Adding 298 model components to the UV plane model. ! The established model now contains 692 components and 5.63072 Jy ![Exited script file: muppet] wmap "J1733-1304_S_2017_10_21_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.554 x 16.29 at -17.07 degrees (North through East) ! Clean map min=-0.019935 max=3.8478 Jy/beam ! Writing clean map to FITS file: J1733-1304_S_2017_10_21_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 3.84778 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00160926 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.125469 0.125469 0.250938 0.501876 1.00375 2.0075 4.01501 8.03001 16.06 32.1201 64.2401 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0048 print "peak=",int(1000*peak_flux)/1000 ! peak= 3.847 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.554 x 16.29 at -17.07 degrees (North through East) ! Clean map min=-0.019935 max=3.8478 Jy/beam ! Error occured in command: mapplot ! ! ! ! ! ! 4 consecutive EOF's on stdin - 10 will exit program ! ! 5 consecutive EOF's on stdin - 10 will exit program ! ! 6 consecutive EOF's on stdin - 10 will exit program ! ! 7 consecutive EOF's on stdin - 10 will exit program ! ! 8 consecutive EOF's on stdin - 10 will exit program ! ! 9 consecutive EOF's on stdin - 10 will exit program ! ! 10 consecutive EOF's on stdin - 10 will exit program ! Log file J1733-1304_S_2017_10_21_pus_uvs.log closed on Sat May 11 21:03:15 2019