! Started logfile: /g0/scr/pima/npcsa_uvs/J0909+0121_X_uvs.log on Sun Oct 6 10:27:23 2013 obs /g0/scr/pima/npcsa_uvs/J0909+0121_X_uva.fits ! Reading UV FITS file: /g0/scr/pima/npcsa_uvs/J0909+0121_X_uva.fits ! AN table 1: 9 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.861111 visibilities/baseline/integration-bin. ! Found source: J0909+0121 ! ! 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 8.40599e+09 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 1116 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.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/01_muppet /g0/scr/pima/npcsa_uvs/J0909+0121_X] 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.200x0.200 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 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=2.018605Jy sigma=29.520673 ! Fit after self-cal, rms=0.814569Jy sigma=12.826316 ! 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=0.7273 mas, bmaj=1.827 mas, bpa=0.9177 degrees ! Estimated noise=1.97867 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.33612 Jy ! Component: 100 - total flux cleaned = 1.70527 Jy ! Total flux subtracted in 100 components = 1.70527 Jy ! Clean residual min=-0.078677 max=0.153179 Jy/beam ! Clean residual mean=-0.000055 rms=0.027746 Jy/beam ! Combined flux in latest and established models = 1.70527 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 1.70527 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.326694Jy sigma=5.110750 ! Fit after self-cal, rms=0.325997Jy sigma=5.106227 ! 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=0.8433 mas, bmaj=2.031 mas, bpa=0.455 degrees ! Estimated noise=1.66572 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.224564 Jy ! Component: 100 - total flux cleaned = 0.330944 Jy ! Component: 150 - total flux cleaned = 0.392794 Jy ! Component: 200 - total flux cleaned = 0.430428 Jy ! Total flux subtracted in 200 components = 0.430428 Jy ! Clean residual min=-0.048116 max=0.072635 Jy/beam ! Clean residual mean=0.000010 rms=0.013165 Jy/beam ! Combined flux in latest and established models = 2.1357 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 32 components and 2.1357 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.189683Jy sigma=2.429351 ! Fit after self-cal, rms=0.187836Jy sigma=2.418407 ! Inverting map ! Added new window around map position (1.2, 1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0636422 Jy ! Component: 100 - total flux cleaned = 0.0868665 Jy ! Component: 150 - total flux cleaned = 0.0894214 Jy ! Component: 200 - total flux cleaned = 0.0900591 Jy ! Total flux subtracted in 200 components = 0.0900591 Jy ! Clean residual min=-0.047361 max=0.039359 Jy/beam ! Clean residual mean=0.000027 rms=0.010155 Jy/beam ! Combined flux in latest and established models = 2.22576 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 56 components and 2.22576 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.162695Jy sigma=2.137672 ! Fit after self-cal, rms=0.159040Jy sigma=2.102307 ! 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 0.99 HN 1.08 KP 1.04 LA 1.02 ! MK 0.99 NL 1.02 OV 1.02 PT 1.02 ! SC 0.98 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 HN 1.05 KP 0.95 LA 0.96 ! MK 0.94 NL 1.00 OV 0.98 PT 1.05 ! SC 0.95 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.09 HN 1.04 KP 0.87 LA 0.98 ! MK 0.96 NL 1.00 OV 1.01 PT 1.01 ! SC 0.98 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 HN 1.03 KP 0.92 LA 0.95 ! MK 0.98 NL 0.99 OV 1.03 PT 1.09 ! SC 1.02 ! ! ! Fit before self-cal, rms=0.159040Jy sigma=2.102307 ! Fit after self-cal, rms=0.111653Jy sigma=1.355199 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.839 mas, bmaj=2.026 mas, bpa=0.4241 degrees ! Estimated noise=1.65864 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 (-3.6, -4.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00937705 Jy ! Component: 100 - total flux cleaned = 0.0116812 Jy ! Component: 150 - total flux cleaned = 0.0107795 Jy ! Component: 200 - total flux cleaned = 0.00971544 Jy ! Total flux subtracted in 200 components = 0.00971544 Jy ! Clean residual min=-0.014131 max=0.013683 Jy/beam ! Clean residual mean=0.000006 rms=0.003139 Jy/beam ! Combined flux in latest and established models = 2.23547 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 78 components and 2.23547 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104973Jy sigma=1.291799 ! Fit after self-cal, rms=0.104764Jy sigma=1.290313 ! 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.104764Jy sigma=1.290313 ! Fit after self-cal, rms=0.104654Jy sigma=1.288097 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8388 mas, bmaj=2.025 mas, bpa=0.4283 degrees ! Estimated noise=1.6599 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00315698 Jy ! Component: 100 - total flux cleaned = 0.00287328 Jy ! Component: 150 - total flux cleaned = 0.00238217 Jy ! Component: 200 - total flux cleaned = 0.0023841 Jy ! Total flux subtracted in 200 components = 0.0023841 Jy ! Clean residual min=-0.012691 max=0.011881 Jy/beam ! Clean residual mean=0.000003 rms=0.002796 Jy/beam ! Combined flux in latest and established models = 2.23786 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 94 components and 2.23786 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103153Jy sigma=1.273652 ! Fit after self-cal, rms=0.103102Jy sigma=1.273409 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.000998789 Jy ! Component: 100 - total flux cleaned = -0.00120787 Jy ! Component: 150 - total flux cleaned = -0.00160348 Jy ! Component: 200 - total flux cleaned = -0.00178621 Jy ! Total flux subtracted in 200 components = -0.00178621 Jy ! Clean residual min=-0.011822 max=0.011126 Jy/beam ! Clean residual mean=0.000003 rms=0.002657 Jy/beam ! Combined flux in latest and established models = 2.23607 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 99 components and 2.23607 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102268Jy sigma=1.265408 ! Fit after self-cal, rms=0.102234Jy sigma=1.265268 ! 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.102234Jy sigma=1.265268 ! Fit after self-cal, rms=0.102166Jy sigma=1.264887 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8388 mas, bmaj=2.025 mas, bpa=0.4368 degrees ! Estimated noise=1.65938 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000122544 Jy ! Component: 100 - total flux cleaned = 0.000126303 Jy ! Component: 150 - total flux cleaned = 0.000128172 Jy ! Component: 200 - total flux cleaned = -2.73628e-05 Jy ! Total flux subtracted in 200 components = -2.73628e-05 Jy ! Clean residual min=-0.010860 max=0.010933 Jy/beam ! Clean residual mean=0.000002 rms=0.002531 Jy/beam ! Combined flux in latest and established models = 2.23604 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 107 components and 2.23604 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101579Jy sigma=1.259203 ! Fit after self-cal, rms=0.101544Jy sigma=1.259075 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 (11.8, -5.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00892213 Jy ! Component: 100 - total flux cleaned = 0.00993746 Jy ! Component: 150 - total flux cleaned = 0.0106268 Jy ! Component: 200 - total flux cleaned = 0.0107762 Jy ! Total flux subtracted in 200 components = 0.0107762 Jy ! Clean residual min=-0.009528 max=0.010240 Jy/beam ! Clean residual mean=0.000003 rms=0.002246 Jy/beam ! Combined flux in latest and established models = 2.24682 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 127 components and 2.24682 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100255Jy sigma=1.246572 ! Fit after self-cal, rms=0.100253Jy sigma=1.246189 ! Inverting map ! Added new window around map position (-8, -3.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00542739 Jy ! Component: 100 - total flux cleaned = 0.00684734 Jy ! Component: 150 - total flux cleaned = 0.00630132 Jy ! Component: 200 - total flux cleaned = 0.00628344 Jy ! Total flux subtracted in 200 components = 0.00628344 Jy ! Clean residual min=-0.008256 max=0.007872 Jy/beam ! Clean residual mean=0.000002 rms=0.001860 Jy/beam ! Combined flux in latest and established models = 2.2531 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 143 components and 2.2531 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098986Jy sigma=1.231723 ! Fit after self-cal, rms=0.098945Jy sigma=1.231320 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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.098945Jy sigma=1.231320 ! Fit after self-cal, rms=0.098927Jy sigma=1.230580 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.291 mas, bmaj=3.399 mas, bpa=4.755 degrees ! Estimated noise=1.9183 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00230024 Jy ! Component: 100 - total flux cleaned = 0.00452893 Jy ! Component: 150 - total flux cleaned = 0.005026 Jy ! Component: 200 - total flux cleaned = 0.00527981 Jy ! Total flux subtracted in 200 components = 0.00527981 Jy ! Clean residual min=-0.005466 max=0.006223 Jy/beam ! Clean residual mean=0.000004 rms=0.001481 Jy/beam ! Combined flux in latest and established models = 2.25838 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 160 components and 2.25838 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098824Jy sigma=1.228678 ! Fit after self-cal, rms=0.098786Jy sigma=1.228450 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.00019954 Jy ! Component: 100 - total flux cleaned = 0.000278595 Jy ! Component: 150 - total flux cleaned = 0.000348675 Jy ! Component: 200 - total flux cleaned = 0.000587639 Jy ! Total flux subtracted in 200 components = 0.000587639 Jy ! Clean residual min=-0.005282 max=0.006191 Jy/beam ! Clean residual mean=0.000003 rms=0.001432 Jy/beam ! Combined flux in latest and established models = 2.25897 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 169 components and 2.25897 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098798Jy sigma=1.227870 ! Fit after self-cal, rms=0.098781Jy sigma=1.227770 ! Inverting map ! Added new window around map position (49.8, 19.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00282015 Jy ! Component: 100 - total flux cleaned = 0.00303638 Jy ! Component: 150 - total flux cleaned = 0.00313541 Jy ! Component: 200 - total flux cleaned = 0.00338564 Jy ! Total flux subtracted in 200 components = 0.00338564 Jy ! Clean residual min=-0.004798 max=0.004779 Jy/beam ! Clean residual mean=0.000002 rms=0.001231 Jy/beam ! Combined flux in latest and established models = 2.26236 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 188 components and 2.26235 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098616Jy sigma=1.225503 ! Fit after self-cal, rms=0.098596Jy sigma=1.225279 ! 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=0.839 mas, bmaj=2.025 mas, bpa=0.4265 degrees ! Estimated noise=1.65819 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.839 x 2.025 at 0.4265 degrees (North through East) ! Clean map min=-0.013303 max=1.7134 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.098596Jy sigma=1.225279 ! Fit after self-cal, rms=0.092605Jy sigma=1.110155 ! Inverting map and beam ! Estimated beam: bmin=0.8392 mas, bmaj=2.025 mas, bpa=0.4773 degrees ! Estimated noise=1.65718 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00215768 Jy ! Component: 100 - total flux cleaned = 0.00143889 Jy ! Component: 150 - total flux cleaned = 0.000967726 Jy ! Component: 200 - total flux cleaned = 0.000235165 Jy ! Total flux subtracted in 200 components = 0.000235165 Jy ! Clean residual min=-0.007058 max=0.007243 Jy/beam ! Clean residual mean=0.000001 rms=0.001531 Jy/beam ! Combined flux in latest and established models = 2.26259 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 215 components and 2.26259 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.091886Jy sigma=1.103047 ! Fit after self-cal, rms=0.091879Jy sigma=1.102914 ! Inverting map ! Added new window around map position (5.2, -29). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00262485 Jy ! Component: 100 - total flux cleaned = 0.0025977 Jy ! Component: 150 - total flux cleaned = 0.00366795 Jy ! Component: 200 - total flux cleaned = 0.00437477 Jy ! Total flux subtracted in 200 components = 0.00437477 Jy ! Clean residual min=-0.006609 max=0.005920 Jy/beam ! Clean residual mean=0.000000 rms=0.001372 Jy/beam ! Combined flux in latest and established models = 2.26696 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 234 components and 2.26696 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.090937Jy sigma=1.095221 ! Fit after self-cal, rms=0.090924Jy sigma=1.095065 ! 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.090924Jy sigma=1.095065 ! Fit after self-cal, rms=0.090916Jy sigma=1.094813 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=0.7231 mas, bmaj=1.817 mas, bpa=0.9811 degrees ! Estimated noise=1.96512 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.32739 Jy ! Component: 100 - total flux cleaned = 1.70294 Jy ! Component: 150 - total flux cleaned = 1.88055 Jy ! Component: 200 - total flux cleaned = 1.98043 Jy ! Total flux subtracted in 200 components = 1.98043 Jy ! Clean residual min=-0.033141 max=0.060868 Jy/beam ! Clean residual mean=-0.000035 rms=0.011072 Jy/beam ! Combined flux in latest and established models = 1.98043 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 20 model components to the UV plane model. ! The established model now contains 20 components and 1.98043 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8393 mas, bmaj=2.025 mas, bpa=0.4737 degrees ! Estimated noise=1.65656 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-5.8, 9). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0893948 Jy ! Component: 100 - total flux cleaned = 0.138276 Jy ! Component: 150 - total flux cleaned = 0.17073 Jy ! Component: 200 - total flux cleaned = 0.188612 Jy ! Total flux subtracted in 200 components = 0.188612 Jy ! Clean residual min=-0.015308 max=0.021149 Jy/beam ! Clean residual mean=-0.000009 rms=0.004552 Jy/beam ! Combined flux in latest and established models = 2.16904 Jy ! Adding 56 model components to the UV plane model. ! The established model now contains 71 components and 2.16904 Jy ! Inverting map ! Added new window around map position (13, -16.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0141122 Jy ! Component: 100 - total flux cleaned = 0.0213882 Jy ! Component: 150 - total flux cleaned = 0.0244215 Jy ! Component: 200 - total flux cleaned = 0.0292179 Jy ! Total flux subtracted in 200 components = 0.0292179 Jy ! Clean residual min=-0.010449 max=0.011614 Jy/beam ! Clean residual mean=-0.000003 rms=0.002842 Jy/beam ! Combined flux in latest and established models = 2.19826 Jy ! Adding 44 model components to the UV plane model. ! The established model now contains 109 components and 2.19826 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.292 mas, bmaj=3.4 mas, bpa=4.781 degrees ! Estimated noise=1.91559 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00971856 Jy ! Component: 100 - total flux cleaned = 0.017852 Jy ! Component: 150 - total flux cleaned = 0.0238694 Jy ! Component: 200 - total flux cleaned = 0.027803 Jy ! Total flux subtracted in 200 components = 0.027803 Jy ! Clean residual min=-0.009419 max=0.007644 Jy/beam ! Clean residual mean=-0.000001 rms=0.002245 Jy/beam ! Combined flux in latest and established models = 2.22606 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.00324638 Jy ! Component: 100 - total flux cleaned = 0.00572747 Jy ! Component: 150 - total flux cleaned = 0.00743669 Jy ! Component: 200 - total flux cleaned = 0.00867947 Jy ! Total flux subtracted in 200 components = 0.00867947 Jy ! Clean residual min=-0.007757 max=0.006258 Jy/beam ! Clean residual mean=-0.000000 rms=0.001766 Jy/beam ! Combined flux in latest and established models = 2.23474 Jy ! Adding 57 model components to the UV plane model. ! The established model now contains 162 components and 2.23474 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.8393 mas, bmaj=2.025 mas, bpa=0.4737 degrees ! Estimated noise=1.65656 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000639031 Jy ! Component: 100 - total flux cleaned = -5.75684e-05 Jy ! Component: 150 - total flux cleaned = 0.00199127 Jy ! Component: 200 - total flux cleaned = 0.00386014 Jy ! Total flux subtracted in 200 components = 0.00386014 Jy ! Clean residual min=-0.008292 max=0.006799 Jy/beam ! Clean residual mean=-0.000000 rms=0.001816 Jy/beam ! Combined flux in latest and established models = 2.2386 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 193 components and 2.2386 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092571Jy sigma=1.114710 ! Fit after self-cal, rms=0.092208Jy sigma=1.111849 wmodel /g0/scr/pima/npcsa_uvs/J0909+0121_X_map.mod ! Writing 193 model components to file: /g0/scr/pima/npcsa_uvs/J0909+0121_X_map.mod wobs /g0/scr/pima/npcsa_uvs/J0909+0121_X_uvs.fits ! Writing UV FITS file: /g0/scr/pima/npcsa_uvs/J0909+0121_X_uvs.fits wwins /g0/scr/pima/npcsa_uvs/J0909+0121_X_map.win ! wwins: Wrote 9 windows to /g0/scr/pima/npcsa_uvs/J0909+0121_X_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.00233508 Jy ! Component: 100 - total flux cleaned = -0.00443964 Jy ! Component: 150 - total flux cleaned = -0.0058201 Jy ! Component: 200 - total flux cleaned = -0.00668238 Jy ! Component: 250 - total flux cleaned = -0.0076038 Jy ! Component: 300 - total flux cleaned = -0.00819523 Jy ! Component: 350 - total flux cleaned = -0.00876768 Jy ! Component: 400 - total flux cleaned = -0.00904445 Jy ! Component: 450 - total flux cleaned = -0.00931412 Jy ! Component: 500 - total flux cleaned = -0.00948801 Jy ! Component: 550 - total flux cleaned = -0.00965703 Jy ! Component: 600 - total flux cleaned = -0.00965818 Jy ! Component: 650 - total flux cleaned = -0.00974051 Jy ! Component: 700 - total flux cleaned = -0.00982068 Jy ! Component: 750 - total flux cleaned = -0.00958833 Jy ! Component: 800 - total flux cleaned = -0.00981764 Jy ! Component: 850 - total flux cleaned = -0.00989276 Jy ! Component: 900 - total flux cleaned = -0.00996521 Jy ! Component: 950 - total flux cleaned = -0.00989447 Jy ! Component: 1000 - total flux cleaned = -0.00975429 Jy ! Component: 1050 - total flux cleaned = -0.00968562 Jy ! Component: 1100 - total flux cleaned = -0.00961881 Jy ! Component: 1150 - total flux cleaned = -0.00975309 Jy ! Component: 1200 - total flux cleaned = -0.00988583 Jy ! Component: 1250 - total flux cleaned = -0.00988721 Jy ! Component: 1300 - total flux cleaned = -0.00995157 Jy ! Component: 1350 - total flux cleaned = -0.00988717 Jy ! Component: 1400 - total flux cleaned = -0.0098245 Jy ! Component: 1450 - total flux cleaned = -0.00982362 Jy ! Component: 1500 - total flux cleaned = -0.00970297 Jy ! Component: 1550 - total flux cleaned = -0.00982135 Jy ! Component: 1600 - total flux cleaned = -0.00987895 Jy ! Component: 1650 - total flux cleaned = -0.00987852 Jy ! Component: 1700 - total flux cleaned = -0.00982075 Jy ! Component: 1750 - total flux cleaned = -0.00965153 Jy ! Component: 1800 - total flux cleaned = -0.00965119 Jy ! Component: 1850 - total flux cleaned = -0.0095966 Jy ! Component: 1900 - total flux cleaned = -0.00954253 Jy ! Component: 1950 - total flux cleaned = -0.0094892 Jy ! Component: 2000 - total flux cleaned = -0.00964865 Jy ! Component: 2050 - total flux cleaned = -0.00949206 Jy ! Component: 2100 - total flux cleaned = -0.00959597 Jy ! Component: 2150 - total flux cleaned = -0.00969857 Jy ! Component: 2200 - total flux cleaned = -0.00969907 Jy ! Component: 2250 - total flux cleaned = -0.00964928 Jy ! Component: 2300 - total flux cleaned = -0.00979812 Jy ! Component: 2350 - total flux cleaned = -0.00984746 Jy ! Component: 2400 - total flux cleaned = -0.00994419 Jy ! Component: 2450 - total flux cleaned = -0.00994436 Jy ! Component: 2500 - total flux cleaned = -0.0099918 Jy ! Component: 2550 - total flux cleaned = -0.0099922 Jy ! Component: 2600 - total flux cleaned = -0.00976078 Jy ! Component: 2650 - total flux cleaned = -0.00980673 Jy ! Component: 2700 - total flux cleaned = -0.00962398 Jy ! Component: 2750 - total flux cleaned = -0.00957853 Jy ! Component: 2800 - total flux cleaned = -0.00953349 Jy ! Component: 2850 - total flux cleaned = -0.00944415 Jy ! Component: 2900 - total flux cleaned = -0.00944414 Jy ! Component: 2950 - total flux cleaned = -0.00935584 Jy ! Component: 3000 - total flux cleaned = -0.00931206 Jy ! Component: 3050 - total flux cleaned = -0.00918182 Jy ! Component: 3100 - total flux cleaned = -0.00926812 Jy ! Component: 3150 - total flux cleaned = -0.00913977 Jy ! Component: 3200 - total flux cleaned = -0.00922486 Jy ! Component: 3250 - total flux cleaned = -0.00939425 Jy ! Component: 3300 - total flux cleaned = -0.0093526 Jy ! Component: 3350 - total flux cleaned = -0.0093944 Jy ! Component: 3400 - total flux cleaned = -0.00943595 Jy ! Component: 3450 - total flux cleaned = -0.00947754 Jy ! Component: 3500 - total flux cleaned = -0.0093535 Jy ! Component: 3550 - total flux cleaned = -0.00927121 Jy ! Component: 3600 - total flux cleaned = -0.00931206 Jy ! Component: 3650 - total flux cleaned = -0.00918956 Jy ! Component: 3700 - total flux cleaned = -0.00910798 Jy ! Component: 3750 - total flux cleaned = -0.00898638 Jy ! Component: 3800 - total flux cleaned = -0.00878369 Jy ! Component: 3850 - total flux cleaned = -0.0086219 Jy ! Component: 3900 - total flux cleaned = -0.00854117 Jy ! Component: 3950 - total flux cleaned = -0.00837977 Jy ! Component: 4000 - total flux cleaned = -0.00809777 Jy ! Component: 4050 - total flux cleaned = -0.00797701 Jy ! Total flux subtracted in 4096 components = -0.00777594 Jy ! Clean residual min=-0.002243 max=0.002475 Jy/beam ! Clean residual mean=-0.000005 rms=0.000631 Jy/beam ! Combined flux in latest and established models = 2.23083 Jy keep ! Adding 252 model components to the UV plane model. ! The established model now contains 445 components and 2.23083 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/g0/scr/pima/npcsa_uvs/J0909+0121_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8393 x 2.025 at 0.4737 degrees (North through East) ! Clean map min=-0.011348 max=1.7091 Jy/beam ! Writing clean map to FITS file: /g0/scr/pima/npcsa_uvs/J0909+0121_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.7091 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000609365 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.106963 0.106963 0.213925 0.42785 0.855701 1.7114 3.4228 6.84561 13.6912 27.3824 54.7649 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0018 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.709 quit ! Quitting program ! Log file /g0/scr/pima/npcsa_uvs/J0909+0121_X_uvs.log closed on Sun Oct 6 10:27:27 2013