! Started logfile: /d2/scr/pima/rv101_uvs/J0455-4615_S_uvs.log on Sat Nov 29 13:50:04 2014 obs /d2/scr/pima/rv101_uvs/J0455-4615_S_uva.fits ! Reading UV FITS file: /d2/scr/pima/rv101_uvs/J0455-4615_S_uva.fits ! AN table 1: 22 integrations on 21 of 21 possible baselines. ! AN table 2: 22 integrations on 15 of 15 possible baselines. ! Apparent sampling: 0.992424 visibilities/baseline/integration-bin. ! Found source: J0455-4615 ! ! 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.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3144 visibilities. select RR ! 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 ![@/home/lpetrov/bin/01_muppet /d2/scr/pima/rv101_uvs/J0455-4615_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 1 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=2.228510Jy sigma=16.808479 ! Fit after self-cal, rms=1.029213Jy sigma=7.078262 ! 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.04 mas, bmaj=7.757 mas, bpa=-4.377 degrees ! Estimated noise=3.23867 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.35482 Jy ! Component: 100 - total flux cleaned = 1.70932 Jy ! Total flux subtracted in 100 components = 1.70932 Jy ! Clean residual min=-0.113703 max=0.170315 Jy/beam ! Clean residual mean=-0.000585 rms=0.027656 Jy/beam ! Combined flux in latest and established models = 1.70931 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 1.70931 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.483654Jy sigma=3.052183 ! Fit after self-cal, rms=0.483593Jy sigma=3.052084 ! Inverting map ! Added new window around map position (-4, -11.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.177636 Jy ! Component: 100 - total flux cleaned = 0.271847 Jy ! Total flux subtracted in 100 components = 0.271847 Jy ! Clean residual min=-0.074485 max=0.115159 Jy/beam ! Clean residual mean=-0.000119 rms=0.018524 Jy/beam ! Combined flux in latest and established models = 1.98116 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 14 components and 1.98116 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.373940Jy sigma=2.554471 ! Fit after self-cal, rms=0.374009Jy sigma=2.551707 ! Inverting map ! Added new window around map position (24.5, 28.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0955268 Jy ! Component: 100 - total flux cleaned = 0.130382 Jy ! Total flux subtracted in 100 components = 0.130382 Jy ! Clean residual min=-0.073837 max=0.054748 Jy/beam ! Clean residual mean=-0.000095 rms=0.014528 Jy/beam ! Combined flux in latest and established models = 2.11154 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 32 components and 2.11154 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.338972Jy sigma=2.316102 ! Fit after self-cal, rms=0.337621Jy sigma=2.303631 ! 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.834 mas, bmaj=9.679 mas, bpa=-5.311 degrees ! Estimated noise=1.70488 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 (22, 4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0934034 Jy ! Component: 100 - total flux cleaned = 0.136631 Jy ! Component: 150 - total flux cleaned = 0.162691 Jy ! Component: 200 - total flux cleaned = 0.172878 Jy ! Total flux subtracted in 200 components = 0.172878 Jy ! Clean residual min=-0.046286 max=0.046565 Jy/beam ! Clean residual mean=-0.000693 rms=0.012188 Jy/beam ! Combined flux in latest and established models = 2.28442 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 76 components and 2.28442 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.303344Jy sigma=2.120875 ! Fit after self-cal, rms=0.299069Jy sigma=2.087518 ! 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.834 x 9.679 at -5.311 degrees (North through East) ! Clean map min=-0.042726 max=1.8144 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: ! FD 1.00* KP 1.00* LA 1.00* MK 1.00* ! OV 1.00* PT 1.00* SC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.86 KP 1.17 LA 1.17 MK 0.97 ! OV 1.00 PT 0.99 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 0.96 KP 1.14 LA 1.04 MK 0.99 ! OV 1.04 PT 1.01 SC 1.06 ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.01 KP 1.16 LA 1.06 MK 0.97 ! OV 1.03 PT 1.00 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 1.02 KP 0.98 LA 0.96 MK 0.99 ! OV 1.04 PT 1.01 SC 0.98 ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.06 KP 1.00 LA 0.94 MK 1.00 ! OV 1.02 PT 1.01 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 0.95 KP 0.89 LA 0.83 MK 1.03 ! OV 0.86 PT 0.91 SC 1.01 ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00* KP 1.00* LA 1.00* MK 1.00* ! OV 1.00* PT 1.00* ! ! ! Fit before self-cal, rms=0.299069Jy sigma=2.087518 ! Fit after self-cal, rms=0.216068Jy sigma=1.434436 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.828 mas, bmaj=9.603 mas, bpa=-5.389 degrees ! Estimated noise=1.71129 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.0321337 Jy ! Component: 100 - total flux cleaned = 0.0525135 Jy ! Component: 150 - total flux cleaned = 0.0668304 Jy ! Component: 200 - total flux cleaned = 0.0730141 Jy ! Total flux subtracted in 200 components = 0.0730141 Jy ! Clean residual min=-0.033634 max=0.036215 Jy/beam ! Clean residual mean=-0.000064 rms=0.009580 Jy/beam ! Combined flux in latest and established models = 2.35744 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 115 components and 2.35744 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206932Jy sigma=1.405534 ! Fit after self-cal, rms=0.206540Jy sigma=1.401164 ! 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.206540Jy sigma=1.401164 ! Fit after self-cal, rms=0.118668Jy sigma=0.883629 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.83 mas, bmaj=9.645 mas, bpa=-5.377 degrees ! Estimated noise=1.71268 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00468115 Jy ! Component: 100 - total flux cleaned = -0.00687375 Jy ! Component: 150 - total flux cleaned = -0.00576642 Jy ! Component: 200 - total flux cleaned = -0.00447337 Jy ! Total flux subtracted in 200 components = -0.00447337 Jy ! Clean residual min=-0.018318 max=0.021333 Jy/beam ! Clean residual mean=-0.000306 rms=0.005421 Jy/beam ! Combined flux in latest and established models = 2.35296 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 134 components and 2.35296 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.116004Jy sigma=0.858618 ! Fit after self-cal, rms=0.115694Jy sigma=0.857293 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.000605681 Jy ! Component: 100 - total flux cleaned = 0.00237278 Jy ! Component: 150 - total flux cleaned = 0.00343413 Jy ! Component: 200 - total flux cleaned = 0.00442579 Jy ! Total flux subtracted in 200 components = 0.00442579 Jy ! Clean residual min=-0.016201 max=0.020526 Jy/beam ! Clean residual mean=-0.000275 rms=0.005192 Jy/beam ! Combined flux in latest and established models = 2.35739 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 146 components and 2.35739 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114530Jy sigma=0.847182 ! Fit after self-cal, rms=0.114406Jy sigma=0.846780 ! 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.114406Jy sigma=0.846780 ! Fit after self-cal, rms=0.114004Jy sigma=0.843145 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.833 mas, bmaj=9.663 mas, bpa=-5.377 degrees ! Estimated noise=1.71308 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00448976 Jy ! Component: 100 - total flux cleaned = -0.00452771 Jy ! Component: 150 - total flux cleaned = -0.00355702 Jy ! Component: 200 - total flux cleaned = -0.00265676 Jy ! Total flux subtracted in 200 components = -0.00265676 Jy ! Clean residual min=-0.015336 max=0.019834 Jy/beam ! Clean residual mean=-0.000288 rms=0.004947 Jy/beam ! Combined flux in latest and established models = 2.35473 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 154 components and 2.35473 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112941Jy sigma=0.833763 ! Fit after self-cal, rms=0.112828Jy sigma=0.833420 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.000417647 Jy ! Component: 100 - total flux cleaned = 0.0014268 Jy ! Component: 150 - total flux cleaned = 0.00217741 Jy ! Component: 200 - total flux cleaned = 0.00287815 Jy ! Total flux subtracted in 200 components = 0.00287815 Jy ! Clean residual min=-0.015194 max=0.019512 Jy/beam ! Clean residual mean=-0.000271 rms=0.004841 Jy/beam ! Combined flux in latest and established models = 2.35761 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 161 components and 2.35761 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112235Jy sigma=0.828752 ! Fit after self-cal, rms=0.112175Jy sigma=0.828605 ! 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.112175Jy sigma=0.828605 ! Fit after self-cal, rms=0.112059Jy sigma=0.827062 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.949 mas, bmaj=33.48 mas, bpa=-0.1637 degrees ! Estimated noise=2.38629 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00341331 Jy ! Component: 100 - total flux cleaned = 0.00413553 Jy ! Component: 150 - total flux cleaned = 0.00480105 Jy ! Component: 200 - total flux cleaned = 0.00540366 Jy ! Total flux subtracted in 200 components = 0.00540366 Jy ! Clean residual min=-0.019259 max=0.029385 Jy/beam ! Clean residual mean=-0.000479 rms=0.008341 Jy/beam ! Combined flux in latest and established models = 2.36301 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 178 components and 2.36301 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.111514Jy sigma=0.824196 ! Fit after self-cal, rms=0.111440Jy sigma=0.823755 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.000839158 Jy ! Component: 100 - total flux cleaned = 0.00139732 Jy ! Component: 150 - total flux cleaned = 0.00174309 Jy ! Component: 200 - total flux cleaned = 0.00207178 Jy ! Total flux subtracted in 200 components = 0.00207178 Jy ! Clean residual min=-0.019064 max=0.028773 Jy/beam ! Clean residual mean=-0.000459 rms=0.008141 Jy/beam ! Combined flux in latest and established models = 2.36509 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 187 components and 2.36509 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.111136Jy sigma=0.822421 ! Fit after self-cal, rms=0.111113Jy sigma=0.822122 ! 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.836 mas, bmaj=9.674 mas, bpa=-5.379 degrees ! Estimated noise=1.71362 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.836 x 9.674 at -5.379 degrees (North through East) ! Clean map min=-0.015047 max=1.8184 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.111113Jy sigma=0.822122 ! Fit after self-cal, rms=0.091309Jy sigma=0.629428 ! Inverting map and beam ! Estimated beam: bmin=2.836 mas, bmaj=9.675 mas, bpa=-5.377 degrees ! Estimated noise=1.71278 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00508545 Jy ! Component: 100 - total flux cleaned = -0.00510663 Jy ! Component: 150 - total flux cleaned = -0.00493451 Jy ! Component: 200 - total flux cleaned = -0.00476865 Jy ! Total flux subtracted in 200 components = -0.00476865 Jy ! Clean residual min=-0.014424 max=0.017744 Jy/beam ! Clean residual mean=-0.000217 rms=0.004457 Jy/beam ! Combined flux in latest and established models = 2.36032 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 199 components and 2.36032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090540Jy sigma=0.622854 ! Fit after self-cal, rms=0.090508Jy sigma=0.622707 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000168499 Jy ! Component: 100 - total flux cleaned = 0.000310311 Jy ! Component: 150 - total flux cleaned = 0.000435973 Jy ! Component: 200 - total flux cleaned = 0.000682332 Jy ! Total flux subtracted in 200 components = 0.000682332 Jy ! Clean residual min=-0.014099 max=0.017806 Jy/beam ! Clean residual mean=-0.000209 rms=0.004414 Jy/beam ! Combined flux in latest and established models = 2.361 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 204 components and 2.361 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090124Jy sigma=0.619935 ! Fit after self-cal, rms=0.090113Jy sigma=0.619843 ! 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.090113Jy sigma=0.619843 ! Fit after self-cal, rms=0.090099Jy sigma=0.619024 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.046 mas, bmaj=7.764 mas, bpa=-4.483 degrees ! Estimated noise=3.25224 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.35874 Jy ! Component: 100 - total flux cleaned = 1.71553 Jy ! Component: 150 - total flux cleaned = 1.89904 Jy ! Component: 200 - total flux cleaned = 2.0127 Jy ! Total flux subtracted in 200 components = 2.0127 Jy ! Clean residual min=-0.070535 max=0.059166 Jy/beam ! Clean residual mean=-0.000367 rms=0.012439 Jy/beam ! Combined flux in latest and established models = 2.0127 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 18 model components to the UV plane model. ! The established model now contains 18 components and 2.0127 Jy ! Inverting map and beam ! Estimated beam: bmin=2.838 mas, bmaj=9.683 mas, bpa=-5.379 degrees ! Estimated noise=1.71351 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.103337 Jy ! Component: 100 - total flux cleaned = 0.170534 Jy ! Component: 150 - total flux cleaned = 0.217458 Jy ! Component: 200 - total flux cleaned = 0.251308 Jy ! Total flux subtracted in 200 components = 0.251308 Jy ! Clean residual min=-0.029213 max=0.028144 Jy/beam ! Clean residual mean=-0.000528 rms=0.007271 Jy/beam ! Combined flux in latest and established models = 2.264 Jy ! Adding 50 model components to the UV plane model. ! The established model now contains 62 components and 2.264 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=9.947 mas, bmaj=33.49 mas, bpa=-0.1588 degrees ! Estimated noise=2.38416 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0404875 Jy ! Component: 100 - total flux cleaned = 0.0645942 Jy ! Component: 150 - total flux cleaned = 0.0793965 Jy ! Component: 200 - total flux cleaned = 0.0867504 Jy ! Total flux subtracted in 200 components = 0.0867504 Jy ! Clean residual min=-0.018290 max=0.032025 Jy/beam ! Clean residual mean=-0.000437 rms=0.008533 Jy/beam ! Combined flux in latest and established models = 2.35075 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.00196917 Jy ! Component: 100 - total flux cleaned = 0.00336507 Jy ! Component: 150 - total flux cleaned = 0.0042105 Jy ! Component: 200 - total flux cleaned = 0.00503126 Jy ! Total flux subtracted in 200 components = 0.00503126 Jy ! Clean residual min=-0.018117 max=0.031918 Jy/beam ! Clean residual mean=-0.000402 rms=0.008266 Jy/beam ! Combined flux in latest and established models = 2.35578 Jy ! Adding 22 model components to the UV plane model. ! The established model now contains 80 components and 2.35579 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.838 mas, bmaj=9.683 mas, bpa=-5.379 degrees ! Estimated noise=1.71351 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00531279 Jy ! Component: 100 - total flux cleaned = -0.00781876 Jy ! Component: 150 - total flux cleaned = -0.00830576 Jy ! Component: 200 - total flux cleaned = -0.00833416 Jy ! Total flux subtracted in 200 components = -0.00833416 Jy ! Clean residual min=-0.013447 max=0.024319 Jy/beam ! Clean residual mean=-0.000250 rms=0.004911 Jy/beam ! Combined flux in latest and established models = 2.34745 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 94 components and 2.34745 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.094293Jy sigma=0.662561 ! Fit after self-cal, rms=0.093857Jy sigma=0.654689 wmodel /d2/scr/pima/rv101_uvs/J0455-4615_S_map.mod ! Writing 94 model components to file: /d2/scr/pima/rv101_uvs/J0455-4615_S_map.mod wobs /d2/scr/pima/rv101_uvs/J0455-4615_S_uvs.fits ! Writing UV FITS file: /d2/scr/pima/rv101_uvs/J0455-4615_S_uvs.fits wwins /d2/scr/pima/rv101_uvs/J0455-4615_S_map.win ! wwins: Wrote 4 windows to /d2/scr/pima/rv101_uvs/J0455-4615_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.00990261 Jy ! Component: 100 - total flux cleaned = 0.0174116 Jy ! Component: 150 - total flux cleaned = 0.0230061 Jy ! Component: 200 - total flux cleaned = 0.0264476 Jy ! Component: 250 - total flux cleaned = 0.0285918 Jy ! Component: 300 - total flux cleaned = 0.0301877 Jy ! Component: 350 - total flux cleaned = 0.0316988 Jy ! Component: 400 - total flux cleaned = 0.0333068 Jy ! Component: 450 - total flux cleaned = 0.0348346 Jy ! Component: 500 - total flux cleaned = 0.0359665 Jy ! Component: 550 - total flux cleaned = 0.0373516 Jy ! Component: 600 - total flux cleaned = 0.0380856 Jy ! Component: 650 - total flux cleaned = 0.0389283 Jy ! Component: 700 - total flux cleaned = 0.0394697 Jy ! Component: 750 - total flux cleaned = 0.0398611 Jy ! Component: 800 - total flux cleaned = 0.0404852 Jy ! Component: 850 - total flux cleaned = 0.0410807 Jy ! Component: 900 - total flux cleaned = 0.0416542 Jy ! Component: 950 - total flux cleaned = 0.0422121 Jy ! Component: 1000 - total flux cleaned = 0.0427498 Jy ! Component: 1050 - total flux cleaned = 0.0428559 Jy ! Component: 1100 - total flux cleaned = 0.0434622 Jy ! Component: 1150 - total flux cleaned = 0.043756 Jy ! Component: 1200 - total flux cleaned = 0.0441353 Jy ! Component: 1250 - total flux cleaned = 0.0442282 Jy ! Component: 1300 - total flux cleaned = 0.0448559 Jy ! Component: 1350 - total flux cleaned = 0.0452894 Jy ! Component: 1400 - total flux cleaned = 0.0457122 Jy ! Component: 1450 - total flux cleaned = 0.0460434 Jy ! Component: 1500 - total flux cleaned = 0.046205 Jy ! Component: 1550 - total flux cleaned = 0.0464416 Jy ! Component: 1600 - total flux cleaned = 0.0465941 Jy ! Component: 1650 - total flux cleaned = 0.0466675 Jy ! Component: 1700 - total flux cleaned = 0.0468867 Jy ! Component: 1750 - total flux cleaned = 0.0469581 Jy ! Component: 1800 - total flux cleaned = 0.0471676 Jy ! Component: 1850 - total flux cleaned = 0.046964 Jy ! Component: 1900 - total flux cleaned = 0.0472343 Jy ! Component: 1950 - total flux cleaned = 0.0472341 Jy ! Component: 2000 - total flux cleaned = 0.0471699 Jy ! Component: 2050 - total flux cleaned = 0.0469804 Jy ! Component: 2100 - total flux cleaned = 0.0466088 Jy ! Component: 2150 - total flux cleaned = 0.0464881 Jy ! Component: 2200 - total flux cleaned = 0.0461278 Jy ! Component: 2250 - total flux cleaned = 0.0458922 Jy ! Component: 2300 - total flux cleaned = 0.0455451 Jy ! Component: 2350 - total flux cleaned = 0.0453183 Jy ! Component: 2400 - total flux cleaned = 0.045319 Jy ! Component: 2450 - total flux cleaned = 0.0450434 Jy ! Component: 2500 - total flux cleaned = 0.0448813 Jy ! Component: 2550 - total flux cleaned = 0.0446692 Jy ! Component: 2600 - total flux cleaned = 0.0445122 Jy ! Component: 2650 - total flux cleaned = 0.0442537 Jy ! Component: 2700 - total flux cleaned = 0.0437964 Jy ! Component: 2750 - total flux cleaned = 0.0436468 Jy ! Component: 2800 - total flux cleaned = 0.0433988 Jy ! Component: 2850 - total flux cleaned = 0.0431071 Jy ! Component: 2900 - total flux cleaned = 0.0427232 Jy ! Component: 2950 - total flux cleaned = 0.042535 Jy ! Component: 3000 - total flux cleaned = 0.0422555 Jy ! Component: 3050 - total flux cleaned = 0.0420717 Jy ! Component: 3100 - total flux cleaned = 0.0417998 Jy ! Component: 3150 - total flux cleaned = 0.0415769 Jy ! Component: 3200 - total flux cleaned = 0.0412676 Jy ! Component: 3250 - total flux cleaned = 0.0411363 Jy ! Component: 3300 - total flux cleaned = 0.0407915 Jy ! Component: 3350 - total flux cleaned = 0.0404097 Jy ! Component: 3400 - total flux cleaned = 0.040326 Jy ! Component: 3450 - total flux cleaned = 0.0399118 Jy ! Component: 3500 - total flux cleaned = 0.0397477 Jy ! Component: 3550 - total flux cleaned = 0.039586 Jy ! Component: 3600 - total flux cleaned = 0.039267 Jy ! Component: 3650 - total flux cleaned = 0.0389525 Jy ! Component: 3700 - total flux cleaned = 0.0386809 Jy ! Component: 3750 - total flux cleaned = 0.0384112 Jy ! Component: 3800 - total flux cleaned = 0.038221 Jy ! Component: 3850 - total flux cleaned = 0.0379199 Jy ! Component: 3900 - total flux cleaned = 0.0376593 Jy ! Component: 3950 - total flux cleaned = 0.0373282 Jy ! Component: 4000 - total flux cleaned = 0.037037 Jy ! Component: 4050 - total flux cleaned = 0.036893 Jy ! Total flux subtracted in 4096 components = 0.0366433 Jy ! Clean residual min=-0.002100 max=0.001926 Jy/beam ! Clean residual mean=-0.000005 rms=0.000740 Jy/beam ! Combined flux in latest and established models = 2.38409 Jy keep ! Adding 332 model components to the UV plane model. ! The established model now contains 426 components and 2.3841 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/d2/scr/pima/rv101_uvs/J0455-4615_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.838 x 9.683 at -5.379 degrees (North through East) ! Clean map min=-0.0073822 max=1.818 Jy/beam ! Writing clean map to FITS file: /d2/scr/pima/rv101_uvs/J0455-4615_S_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.81803 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000733473 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.121033 0.121033 0.242066 0.484133 0.968266 1.93653 3.87306 7.74613 15.4923 30.9845 61.969 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0022 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.818 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2478.66 quit ! Quitting program ! Log file /d2/scr/pima/rv101_uvs/J0455-4615_S_uvs.log closed on Sat Nov 29 13:50:09 2014