! Started logfile: ./J1924-2914_S_uvs.log on Sat Jan 2 10:00:27 2016 obs J1924-2914_S_uva.fits ! Reading UV FITS file: J1924-2914_S_uva.fits ! AN table 1: 11 integrations on 28 of 28 possible baselines. ! AN table 2: 30 integrations on 36 of 36 possible baselines. ! AN table 3: 2 integrations on 21 of 21 possible baselines. ! AN table 4: 12 integrations on 6 of 6 possible baselines. ! Apparent sampling: 0.768309 visibilities/baseline/integration-bin. ! Found source: J1924-2914 ! ! There are 8 IFs, and a total of 8 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.22499e+09 8e+06 1 8e+06 ! 02 2 2.23299e+09 8e+06 1 8e+06 ! 03 3 2.25499e+09 8e+06 1 8e+06 ! 04 4 2.26299e+09 8e+06 1 8e+06 ! 05 5 2.34499e+09 8e+06 1 8e+06 ! 06 6 2.35299e+09 8e+06 1 8e+06 ! 07 7 2.36499e+09 8e+06 1 8e+06 ! 08 8 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 9232 visibilities. select rr, 1,2,3,4,6,7,8 ! Selecting polarization: RR, channels: 1..4 6..8 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ! Reading IF 5 channels: (none) ! Reading IF 6 channels: 6..6 ! Reading IF 7 channels: 7..7 ! Reading IF 8 channels: 8..8 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 ![@/home/lpetrov/bin/01_muppet ./J1924-2914_S] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff 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 73 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 73 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 95 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 95 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! A total of 11 telescope corrections were flagged in sub-array 1. ! A total of 95 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 15 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 7. ! A total of 73 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 48 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 8. ! A total of 73 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 48 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=9.414323Jy sigma=20.410812 ! Fit after self-cal, rms=8.423892Jy sigma=18.027404 ! 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.349 mas, bmaj=4.389 mas, bpa=3.913 degrees ! Estimated noise=6.41493 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 = 5.88262 Jy ! Component: 100 - total flux cleaned = 7.78993 Jy ! Total flux subtracted in 100 components = 7.78993 Jy ! Clean residual min=-0.450527 max=0.827024 Jy/beam ! Clean residual mean=-0.000978 rms=0.105925 Jy/beam ! Combined flux in latest and established models = 7.78993 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 15 components and 7.78993 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=2.682616Jy sigma=6.060783 ! Fit after self-cal, rms=2.681387Jy sigma=6.060259 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.9441 Jy ! Component: 100 - total flux cleaned = 1.4883 Jy ! Total flux subtracted in 100 components = 1.4883 Jy ! Clean residual min=-0.273740 max=0.408346 Jy/beam ! Clean residual mean=-0.000891 rms=0.069040 Jy/beam ! Combined flux in latest and established models = 9.27823 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 24 components and 9.27823 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.939300Jy sigma=5.003358 ! Fit after self-cal, rms=1.937845Jy sigma=5.002492 ! 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.771 mas, bmaj=5.066 mas, bpa=5.414 degrees ! Estimated noise=4.21982 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.479307 Jy ! Component: 100 - total flux cleaned = 0.703134 Jy ! Component: 150 - total flux cleaned = 0.814233 Jy ! Component: 200 - total flux cleaned = 0.828704 Jy ! Total flux subtracted in 200 components = 0.828704 Jy ! Clean residual min=-0.291712 max=0.372884 Jy/beam ! Clean residual mean=-0.003840 rms=0.116226 Jy/beam ! Combined flux in latest and established models = 10.1069 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 33 components and 10.1069 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.702663Jy sigma=5.008751 ! Fit after self-cal, rms=1.701572Jy sigma=5.006780 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.771 x 5.066 at 5.414 degrees (North through East) ! Clean map min=-0.29109 max=7.8773 Jy/beam if ( peak(flux,max)/imstat(rms) < lim_self ) selflim lim_val ! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2. ! Selfcal phase corrections will not be limited. gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 0.95 KP 1.00 LA 1.00 ! NL 1.06 OV 0.99 PT 0.98 SC 1.02 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.90 FD 0.95 KP 0.97 LA 0.97 ! MK 0.99 NL 1.08 OV 0.94 PT 0.98 ! SC 0.98 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.97 FD 0.96 KP 1.03 LA 1.01 ! MK 1.08 NL 1.08 PT 1.00* ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 0.93 KP 1.01 LA 1.03 ! NL 1.20 OV 1.01 PT 1.01 SC 1.08 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.90 FD 0.93 KP 1.00 LA 1.03 ! MK 0.95 NL 1.13 OV 0.93 PT 0.99 ! SC 1.01 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.96 FD 0.98 KP 1.03 LA 1.10 ! MK 1.04 NL 1.18 PT 1.00* ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.96 KP 1.00 LA 0.95 ! NL 1.15 OV 1.00* PT 0.95 SC 1.14 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 0.98 KP 0.98 LA 0.95 ! MK 1.08 NL 1.15 OV 1.00* PT 0.96 ! SC 1.07 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.06 FD 1.02 KP 1.04 LA 1.00 ! MK 1.16 NL 1.14 PT 1.00* ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 0.95 KP 1.00 LA 0.99 ! NL 1.19 OV 1.00* PT 0.96 SC 1.11 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* PT 1.00* ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* PT 1.00* ! ! A total of 33 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! NL 1.00* OV 1.00* PT 1.00* SC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.00* KP 1.00* LA 1.00* ! MK 1.00* NL 1.00* PT 1.00* ! ! Telescope amplitude corrections in sub-array 4: ! EF 0.00* HH 0.00* MC 0.00* ZC 0.00* ! ! ! Fit before self-cal, rms=1.701572Jy sigma=5.006780 ! Fit after self-cal, rms=1.429949Jy sigma=2.384066 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.71 mas, bmaj=7.473 mas, bpa=0.4879 degrees ! Estimated noise=4.49815 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 (7.5, 1.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.391599 Jy ! Component: 100 - total flux cleaned = 0.592829 Jy ! Component: 150 - total flux cleaned = 0.710354 Jy ! Component: 200 - total flux cleaned = 0.748513 Jy ! Total flux subtracted in 200 components = 0.748513 Jy ! Clean residual min=-0.133197 max=0.205892 Jy/beam ! Clean residual mean=-0.000343 rms=0.037532 Jy/beam ! Combined flux in latest and established models = 10.8554 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 89 components and 10.8554 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.249009Jy sigma=2.194654 ! Fit after self-cal, rms=1.225315Jy sigma=2.121902 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.225315Jy sigma=2.121902 ! Fit after self-cal, rms=1.138286Jy sigma=1.895303 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.693 mas, bmaj=7.44 mas, bpa=0.3452 degrees ! Estimated noise=4.53527 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.204974 Jy ! Component: 100 - total flux cleaned = 0.326728 Jy ! Component: 150 - total flux cleaned = 0.387296 Jy ! Component: 200 - total flux cleaned = 0.417789 Jy ! Total flux subtracted in 200 components = 0.417789 Jy ! Clean residual min=-0.116472 max=0.140502 Jy/beam ! Clean residual mean=-0.000398 rms=0.031010 Jy/beam ! Combined flux in latest and established models = 11.2732 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 140 components and 11.2732 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.062680Jy sigma=1.847120 ! Fit after self-cal, rms=1.055873Jy sigma=1.825789 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.0866199 Jy ! Component: 100 - total flux cleaned = 0.113787 Jy ! Component: 150 - total flux cleaned = 0.131256 Jy ! Component: 200 - total flux cleaned = 0.142787 Jy ! Total flux subtracted in 200 components = 0.142787 Jy ! Clean residual min=-0.107739 max=0.125771 Jy/beam ! Clean residual mean=-0.000213 rms=0.029634 Jy/beam ! Combined flux in latest and established models = 11.416 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 162 components and 11.416 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.037873Jy sigma=1.821263 ! Fit after self-cal, rms=1.036874Jy sigma=1.816759 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.036874Jy sigma=1.816759 ! Fit after self-cal, rms=0.932712Jy sigma=1.613051 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=7.426 mas, bpa=0.2799 degrees ! Estimated noise=4.58385 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.117436 Jy ! Component: 100 - total flux cleaned = 0.188994 Jy ! Component: 150 - total flux cleaned = 0.236157 Jy ! Component: 200 - total flux cleaned = 0.263754 Jy ! Total flux subtracted in 200 components = 0.263754 Jy ! Clean residual min=-0.080623 max=0.102706 Jy/beam ! Clean residual mean=-0.000183 rms=0.023007 Jy/beam ! Combined flux in latest and established models = 11.6798 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 183 components and 11.6798 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.897799Jy sigma=1.608810 ! Fit after self-cal, rms=0.896435Jy sigma=1.603321 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.0432753 Jy ! Component: 100 - total flux cleaned = 0.0610645 Jy ! Component: 150 - total flux cleaned = 0.0646066 Jy ! Component: 200 - total flux cleaned = 0.0668846 Jy ! Total flux subtracted in 200 components = 0.0668846 Jy ! Clean residual min=-0.079785 max=0.096405 Jy/beam ! Clean residual mean=-0.000097 rms=0.022556 Jy/beam ! Combined flux in latest and established models = 11.7467 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 193 components and 11.7467 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.890755Jy sigma=1.602676 ! Fit after self-cal, rms=0.890630Jy sigma=1.601420 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.890630Jy sigma=1.601420 ! Fit after self-cal, rms=0.838297Jy sigma=1.465409 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.11 mas, bmaj=19.63 mas, bpa=3.303 degrees ! Estimated noise=7.18109 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.142863 Jy ! Component: 100 - total flux cleaned = 0.205973 Jy ! Component: 150 - total flux cleaned = 0.238504 Jy ! Component: 200 - total flux cleaned = 0.243603 Jy ! Total flux subtracted in 200 components = 0.243603 Jy ! Clean residual min=-0.105351 max=0.113086 Jy/beam ! Clean residual mean=-0.000202 rms=0.035938 Jy/beam ! Combined flux in latest and established models = 11.9903 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 214 components and 11.9903 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.817178Jy sigma=1.465474 ! Fit after self-cal, rms=0.817195Jy sigma=1.462457 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.0319249 Jy ! Component: 100 - total flux cleaned = 0.0312335 Jy ! Component: 150 - total flux cleaned = 0.0312262 Jy ! Component: 200 - total flux cleaned = 0.0317775 Jy ! Total flux subtracted in 200 components = 0.0317775 Jy ! Clean residual min=-0.109259 max=0.112946 Jy/beam ! Clean residual mean=-0.000107 rms=0.035717 Jy/beam ! Combined flux in latest and established models = 12.022 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 217 components and 12.022 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.816740Jy sigma=1.466263 ! Fit after self-cal, rms=0.816814Jy sigma=1.465414 ! 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=2.677 mas, bmaj=7.415 mas, bpa=0.2208 degrees ! Estimated noise=4.62477 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.677 x 7.415 at 0.2208 degrees (North through East) ! Clean map min=-0.12263 max=8.3848 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.816814Jy sigma=1.465414 ! Fit after self-cal, rms=0.786679Jy sigma=1.407890 ! Inverting map and beam ! Estimated beam: bmin=2.667 mas, bmaj=7.396 mas, bpa=0.249 degrees ! Estimated noise=4.64805 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0793364 Jy ! Component: 100 - total flux cleaned = 0.11854 Jy ! Component: 150 - total flux cleaned = 0.136585 Jy ! Component: 200 - total flux cleaned = 0.141944 Jy ! Total flux subtracted in 200 components = 0.141944 Jy ! Clean residual min=-0.065048 max=0.077303 Jy/beam ! Clean residual mean=-0.000180 rms=0.019247 Jy/beam ! Combined flux in latest and established models = 12.164 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 235 components and 12.164 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.770321Jy sigma=1.406078 ! Fit after self-cal, rms=0.769907Jy sigma=1.404870 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0184894 Jy ! Component: 100 - total flux cleaned = 0.0225259 Jy ! Component: 150 - total flux cleaned = 0.023098 Jy ! Component: 200 - total flux cleaned = 0.0231444 Jy ! Total flux subtracted in 200 components = 0.0231444 Jy ! Clean residual min=-0.064968 max=0.074392 Jy/beam ! Clean residual mean=-0.000150 rms=0.019087 Jy/beam ! Combined flux in latest and established models = 12.1871 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 242 components and 12.1871 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.768250Jy sigma=1.402553 ! Fit after self-cal, rms=0.768059Jy sigma=1.402073 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.768059Jy sigma=1.402073 ! Fit after self-cal, rms=0.672964Jy sigma=1.278363 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.255 mas, bmaj=6.673 mas, bpa=-0.834 degrees ! Estimated noise=7.15843 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 6.37604 Jy ! Component: 100 - total flux cleaned = 8.24582 Jy ! Component: 150 - total flux cleaned = 9.18243 Jy ! Component: 200 - total flux cleaned = 9.77761 Jy ! Total flux subtracted in 200 components = 9.77761 Jy ! Clean residual min=-0.195390 max=0.349036 Jy/beam ! Clean residual mean=-0.000975 rms=0.052644 Jy/beam ! Combined flux in latest and established models = 9.7776 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 35 model components to the UV plane model. ! The established model now contains 35 components and 9.7776 Jy ! Inverting map and beam ! Estimated beam: bmin=2.663 mas, bmaj=7.391 mas, bpa=0.2112 degrees ! Estimated noise=4.66791 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.803205 Jy ! Component: 100 - total flux cleaned = 1.32493 Jy ! Component: 150 - total flux cleaned = 1.68191 Jy ! Component: 200 - total flux cleaned = 1.92952 Jy ! Total flux subtracted in 200 components = 1.92952 Jy ! Clean residual min=-0.077244 max=0.155813 Jy/beam ! Clean residual mean=-0.001021 rms=0.026761 Jy/beam ! Combined flux in latest and established models = 11.7071 Jy ! Adding 69 model components to the UV plane model. ! The established model now contains 94 components and 11.7071 Jy ! Inverting map ! Added new window around map position (2.5, -10). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.179339 Jy ! Component: 100 - total flux cleaned = 0.30452 Jy ! Component: 150 - total flux cleaned = 0.394032 Jy ! Component: 200 - total flux cleaned = 0.460703 Jy ! Total flux subtracted in 200 components = 0.460703 Jy ! Clean residual min=-0.060362 max=0.074511 Jy/beam ! Clean residual mean=-0.000430 rms=0.018013 Jy/beam ! Combined flux in latest and established models = 12.1678 Jy ! Adding 42 model components to the UV plane model. ! The established model now contains 127 components and 12.1678 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.078 mas, bmaj=19.55 mas, bpa=3.225 degrees ! Estimated noise=7.3093 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.101614 Jy ! Component: 100 - total flux cleaned = 0.14908 Jy ! Component: 150 - total flux cleaned = 0.162193 Jy ! Component: 200 - total flux cleaned = 0.172481 Jy ! Total flux subtracted in 200 components = 0.172481 Jy ! Clean residual min=-0.091431 max=0.111560 Jy/beam ! Clean residual mean=-0.000530 rms=0.030867 Jy/beam ! Combined flux in latest and established models = 12.3403 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.00715346 Jy ! Component: 100 - total flux cleaned = 0.0121573 Jy ! Component: 150 - total flux cleaned = 0.0143283 Jy ! Component: 200 - total flux cleaned = 0.0163494 Jy ! Total flux subtracted in 200 components = 0.0163494 Jy ! Clean residual min=-0.091954 max=0.110466 Jy/beam ! Clean residual mean=-0.000485 rms=0.030755 Jy/beam ! Combined flux in latest and established models = 12.3567 Jy ! Adding 24 model components to the UV plane model. ! The established model now contains 145 components and 12.3567 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.663 mas, bmaj=7.391 mas, bpa=0.2112 degrees ! Estimated noise=4.66791 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00590816 Jy ! Component: 100 - total flux cleaned = -0.0133221 Jy ! Component: 150 - total flux cleaned = -0.0172652 Jy ! Component: 200 - total flux cleaned = -0.0196875 Jy ! Total flux subtracted in 200 components = -0.0196875 Jy ! Clean residual min=-0.058344 max=0.072558 Jy/beam ! Clean residual mean=-0.000223 rms=0.016995 Jy/beam ! Combined flux in latest and established models = 12.337 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 167 components and 12.337 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.662794Jy sigma=1.292987 ! Fit after self-cal, rms=0.659631Jy sigma=1.287548 wmodel ./J1924-2914_S_map.mod ! Writing 167 model components to file: ./J1924-2914_S_map.mod wobs ./J1924-2914_S_uvs.fits ! Writing UV FITS file: ./J1924-2914_S_uvs.fits wwins ./J1924-2914_S_map.win ! wwins: Wrote 3 windows to ./J1924-2914_S_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.0324313 Jy ! Component: 100 - total flux cleaned = 0.0575424 Jy ! Component: 150 - total flux cleaned = 0.064931 Jy ! Component: 200 - total flux cleaned = 0.0699978 Jy ! Component: 250 - total flux cleaned = 0.0747604 Jy ! Component: 300 - total flux cleaned = 0.0793535 Jy ! Component: 350 - total flux cleaned = 0.082004 Jy ! Component: 400 - total flux cleaned = 0.0819909 Jy ! Component: 450 - total flux cleaned = 0.0836739 Jy ! Component: 500 - total flux cleaned = 0.0844972 Jy ! Component: 550 - total flux cleaned = 0.0828904 Jy ! Component: 600 - total flux cleaned = 0.08447 Jy ! Component: 650 - total flux cleaned = 0.0821583 Jy ! Component: 700 - total flux cleaned = 0.0814014 Jy ! Component: 750 - total flux cleaned = 0.0791676 Jy ! Component: 800 - total flux cleaned = 0.0777032 Jy ! Component: 850 - total flux cleaned = 0.0762573 Jy ! Component: 900 - total flux cleaned = 0.0748506 Jy ! Component: 950 - total flux cleaned = 0.0727684 Jy ! Component: 1000 - total flux cleaned = 0.0720887 Jy ! Component: 1050 - total flux cleaned = 0.0707302 Jy ! Component: 1100 - total flux cleaned = 0.0693987 Jy ! Component: 1150 - total flux cleaned = 0.0687464 Jy ! Component: 1200 - total flux cleaned = 0.0668161 Jy ! Component: 1250 - total flux cleaned = 0.0687314 Jy ! Component: 1300 - total flux cleaned = 0.0674669 Jy ! Component: 1350 - total flux cleaned = 0.0662219 Jy ! Component: 1400 - total flux cleaned = 0.0662252 Jy ! Component: 1450 - total flux cleaned = 0.0644106 Jy ! Component: 1500 - total flux cleaned = 0.0650124 Jy ! Component: 1550 - total flux cleaned = 0.0626383 Jy ! Component: 1600 - total flux cleaned = 0.0620495 Jy ! Component: 1650 - total flux cleaned = 0.0608835 Jy ! Component: 1700 - total flux cleaned = 0.0614598 Jy ! Component: 1750 - total flux cleaned = 0.0597432 Jy ! Component: 1800 - total flux cleaned = 0.0580575 Jy ! Component: 1850 - total flux cleaned = 0.057503 Jy ! Component: 1900 - total flux cleaned = 0.0569533 Jy ! Component: 1950 - total flux cleaned = 0.0585876 Jy ! Component: 2000 - total flux cleaned = 0.0564127 Jy ! Component: 2050 - total flux cleaned = 0.0569401 Jy ! Component: 2100 - total flux cleaned = 0.0548148 Jy ! Component: 2150 - total flux cleaned = 0.054289 Jy ! Component: 2200 - total flux cleaned = 0.0516877 Jy ! Component: 2250 - total flux cleaned = 0.0532357 Jy ! Component: 2300 - total flux cleaned = 0.0522043 Jy ! Component: 2350 - total flux cleaned = 0.0491581 Jy ! Component: 2400 - total flux cleaned = 0.0496681 Jy ! Component: 2450 - total flux cleaned = 0.0476715 Jy ! Component: 2500 - total flux cleaned = 0.0476713 Jy ! Component: 2550 - total flux cleaned = 0.0457129 Jy ! Component: 2600 - total flux cleaned = 0.046688 Jy ! Component: 2650 - total flux cleaned = 0.0452413 Jy ! Component: 2700 - total flux cleaned = 0.0438116 Jy ! Component: 2750 - total flux cleaned = 0.0438142 Jy ! Component: 2800 - total flux cleaned = 0.0428778 Jy ! Component: 2850 - total flux cleaned = 0.0428811 Jy ! Component: 2900 - total flux cleaned = 0.0410346 Jy ! Component: 2950 - total flux cleaned = 0.0410358 Jy ! Component: 3000 - total flux cleaned = 0.0396679 Jy ! Component: 3050 - total flux cleaned = 0.0383124 Jy ! Component: 3100 - total flux cleaned = 0.0374166 Jy ! Component: 3150 - total flux cleaned = 0.0369695 Jy ! Component: 3200 - total flux cleaned = 0.0352033 Jy ! Component: 3250 - total flux cleaned = 0.0330156 Jy ! Component: 3300 - total flux cleaned = 0.0325807 Jy ! Component: 3350 - total flux cleaned = 0.0312876 Jy ! Component: 3400 - total flux cleaned = 0.0282896 Jy ! Component: 3450 - total flux cleaned = 0.0287193 Jy ! Component: 3500 - total flux cleaned = 0.0266104 Jy ! Component: 3550 - total flux cleaned = 0.0253555 Jy ! Component: 3600 - total flux cleaned = 0.0249346 Jy ! Component: 3650 - total flux cleaned = 0.0232861 Jy ! Component: 3700 - total flux cleaned = 0.0228777 Jy ! Component: 3750 - total flux cleaned = 0.0200336 Jy ! Component: 3800 - total flux cleaned = 0.022053 Jy ! Component: 3850 - total flux cleaned = 0.0200437 Jy ! Component: 3900 - total flux cleaned = 0.0196403 Jy ! Component: 3950 - total flux cleaned = 0.0196425 Jy ! Component: 4000 - total flux cleaned = 0.0180684 Jy ! Component: 4050 - total flux cleaned = 0.017288 Jy ! Total flux subtracted in 4096 components = 0.0165099 Jy ! Clean residual min=-0.021815 max=0.024956 Jy/beam ! Clean residual mean=-0.000092 rms=0.007834 Jy/beam ! Combined flux in latest and established models = 12.3535 Jy keep ! Adding 395 model components to the UV plane model. ! The established model now contains 562 components and 12.3535 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "./J1924-2914_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.663 x 7.391 at 0.2112 degrees (North through East) ! Clean map min=-0.047526 max=8.5044 Jy/beam ! Writing clean map to FITS file: ./J1924-2914_S_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 8.5044 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00779468 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.274964 0.274964 0.549928 1.09986 2.19971 4.39942 8.79884 17.5977 35.1954 70.3907 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0233 print "peak= ", int(1000*peak_flux)/1000 ! peak= 8.504 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 1091.05 quit ! Quitting program ! Log file ./J1924-2914_S_uvs.log closed on Sat Jan 2 10:00:31 2016