! Started logfile: /g0/scr/pima//bp171af_uvs/J2219-1221_X_uvs.log on Sun Dec 8 13:13:51 2013 obs /g0/scr/pima//bp171af_uvs/J2219-1221_X_uva.fits ! Reading UV FITS file: /g0/scr/pima//bp171af_uvs/J2219-1221_X_uva.fits ! AN table 1: 3 integrations on 28 of 28 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J2219-1221 ! ! 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 7.392e+09 3.2e+07 1 3.2e+07 ! 02 2 7.424e+09 3.2e+07 1 3.2e+07 ! 03 3 7.456e+09 3.2e+07 1 3.2e+07 ! 04 4 7.552e+09 3.2e+07 1 3.2e+07 ! 05 5 7.744e+09 3.2e+07 1 3.2e+07 ! 06 6 7.776e+09 3.2e+07 1 3.2e+07 ! 07 7 7.808e+09 3.2e+07 1 3.2e+07 ! 08 8 7.84e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 672 visibilities. select RR ! Selecting polarization: RR, channels: 1..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: 5..5 ! 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 = 4096 print "field_cell =",field_cell ! field_cell = 0.15 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/01_muppet /g0/scr/pima//bp171af_uvs/J2219-1221_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 = 4096x4096 pixels with 0.150x0.150 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=1.001103Jy sigma=1.176586 ! Fit after self-cal, rms=0.988163Jy sigma=1.116443 ! 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.6915 mas, bmaj=1.733 mas, bpa=-3.815 degrees ! Estimated noise=37.0613 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=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.015239 Jy ! Component: 100 - total flux cleaned = 0.0195557 Jy ! Total flux subtracted in 100 components = 0.0195557 Jy ! Clean residual min=-0.007476 max=0.007493 Jy/beam ! Clean residual mean=0.000000 rms=0.002379 Jy/beam ! Combined flux in latest and established models = 0.0195557 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 41 components and 0.0195557 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.026300Jy sigma=0.056155 ! Fit after self-cal, rms=0.025276Jy sigma=0.055474 ! 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.8199 mas, bmaj=1.993 mas, bpa=-3.029 degrees ! Estimated noise=29.0958 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.00700376 Jy ! Component: 100 - total flux cleaned = 0.00712716 Jy ! Component: 150 - total flux cleaned = 0.00860721 Jy ! Component: 200 - total flux cleaned = 0.00782958 Jy ! Total flux subtracted in 200 components = 0.00782958 Jy ! Clean residual min=-0.003968 max=0.003965 Jy/beam ! Clean residual mean=0.000000 rms=0.001979 Jy/beam ! Combined flux in latest and established models = 0.0273853 Jy ! Performing phase self-cal ! Adding 128 model components to the UV plane model. ! The established model now contains 165 components and 0.0273853 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.020324Jy sigma=0.028537 ! Fit after self-cal, rms=0.020063Jy sigma=0.028348 ! 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.90 HN 0.77 KP 0.90 MK 1.10 ! NL 1.24 OV 0.71 PT 0.79 SC 0.71 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 HN 0.95 KP 0.81 MK 0.99 ! NL 0.81 OV 1.23 PT 0.83 SC 0.98 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.77 HN 0.63 KP 0.91 MK 1.26 ! NL 0.97 OV 1.21 PT 0.88 SC 0.89 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.59 HN 0.97 KP 0.87 MK 1.57 ! NL 0.88 OV 0.82 PT 0.80 SC 1.03 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 0.78 HN 0.83 KP 1.44 MK 1.20 ! NL 0.85 OV 0.83 PT 0.64 SC 1.05 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.34 HN 0.73 KP 0.77 MK 0.79 ! NL 1.09 OV 0.85 PT 1.12 SC 0.84 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.77 HN 0.89 KP 0.90 MK 1.11 ! NL 0.98 OV 1.63 PT 0.76 SC 0.87 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 HN 0.88 KP 0.90 MK 0.96 ! NL 0.96 OV 0.80 PT 0.89 SC 0.90 ! ! ! Fit before self-cal, rms=0.020063Jy sigma=0.028348 ! Fit after self-cal, rms=0.017771Jy sigma=0.019471 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8454 mas, bmaj=2.039 mas, bpa=-1.152 degrees ! Estimated noise=23.666 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.000949983 Jy ! Component: 100 - total flux cleaned = -0.00107125 Jy ! Component: 150 - total flux cleaned = -0.000832298 Jy ! Component: 200 - total flux cleaned = -0.000780416 Jy ! Total flux subtracted in 200 components = -0.000780416 Jy ! Clean residual min=-0.000867 max=0.000866 Jy/beam ! Clean residual mean=-0.000000 rms=0.000230 Jy/beam ! Combined flux in latest and established models = 0.0266049 Jy ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 231 components and 0.0266049 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.017122Jy sigma=0.019187 ! Fit after self-cal, rms=0.017093Jy sigma=0.019138 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.017093Jy sigma=0.019138 ! Fit after self-cal, rms=0.016503Jy sigma=0.018992 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8442 mas, bmaj=2.042 mas, bpa=-0.9163 degrees ! Estimated noise=22.4677 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000999194 Jy ! Component: 100 - total flux cleaned = -0.000898107 Jy ! Component: 150 - total flux cleaned = -0.000848166 Jy ! Component: 200 - total flux cleaned = -0.000625979 Jy ! Total flux subtracted in 200 components = -0.000625979 Jy ! Clean residual min=-0.000705 max=0.000708 Jy/beam ! Clean residual mean=0.000000 rms=0.000195 Jy/beam ! Combined flux in latest and established models = 0.0259789 Jy selfcal ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 274 components and 0.0259789 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.016060Jy sigma=0.018817 ! Fit after self-cal, rms=0.016031Jy sigma=0.018791 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.000314363 Jy ! Component: 100 - total flux cleaned = 0.000602732 Jy ! Component: 150 - total flux cleaned = 0.000799446 Jy ! Component: 200 - total flux cleaned = 0.000908781 Jy ! Total flux subtracted in 200 components = 0.000908781 Jy ! Clean residual min=-0.000601 max=0.000602 Jy/beam ! Clean residual mean=-0.000000 rms=0.000166 Jy/beam ! Combined flux in latest and established models = 0.0268877 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 324 components and 0.0268877 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015821Jy sigma=0.018636 ! Fit after self-cal, rms=0.015800Jy sigma=0.018620 ! Inverting map ! Added new window around map position (-28.65, 56.7). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000596713 Jy ! Component: 100 - total flux cleaned = 0.000674094 Jy ! Component: 150 - total flux cleaned = 0.000703359 Jy ! Component: 200 - total flux cleaned = 0.000728661 Jy ! Total flux subtracted in 200 components = 0.000728661 Jy ! Clean residual min=-0.000767 max=0.000745 Jy/beam ! Clean residual mean=0.000000 rms=0.000179 Jy/beam ! Combined flux in latest and established models = 0.0276164 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 351 components and 0.0276164 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015772Jy sigma=0.018609 ! Fit after self-cal, rms=0.015758Jy sigma=0.018592 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015758Jy sigma=0.018592 ! Fit after self-cal, rms=0.015067Jy sigma=0.018408 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8449 mas, bmaj=2.048 mas, bpa=-0.4703 degrees ! Estimated noise=21.1247 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000544924 Jy ! Component: 100 - total flux cleaned = 0.000882756 Jy ! Component: 150 - total flux cleaned = 0.00103925 Jy ! Component: 200 - total flux cleaned = 0.00110908 Jy ! Total flux subtracted in 200 components = 0.00110908 Jy ! Clean residual min=-0.001398 max=0.000994 Jy/beam ! Clean residual mean=0.000000 rms=0.000243 Jy/beam ! Combined flux in latest and established models = 0.0287255 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 364 components and 0.0287255 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015028Jy sigma=0.018432 ! Fit after self-cal, rms=0.015019Jy sigma=0.018410 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.000299505 Jy ! Component: 100 - total flux cleaned = 0.000455549 Jy ! Component: 150 - total flux cleaned = 0.000503406 Jy ! Component: 200 - total flux cleaned = 0.000538947 Jy ! Total flux subtracted in 200 components = 0.000538947 Jy ! Clean residual min=-0.001401 max=0.001034 Jy/beam ! Clean residual mean=0.000000 rms=0.000244 Jy/beam ! Combined flux in latest and established models = 0.0292644 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 370 components and 0.0292644 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015006Jy sigma=0.018427 ! Fit after self-cal, rms=0.015005Jy sigma=0.018420 ! 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. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.015005Jy sigma=0.018420 ! Fit after self-cal, rms=0.014932Jy sigma=0.018382 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.055 mas, bmaj=4.085 mas, bpa=2.399 degrees ! Estimated noise=24.6643 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000546482 Jy ! Component: 100 - total flux cleaned = 0.000806414 Jy ! Component: 150 - total flux cleaned = 0.000889792 Jy ! Component: 200 - total flux cleaned = 0.000898879 Jy ! Total flux subtracted in 200 components = 0.000898879 Jy ! Clean residual min=-0.001487 max=0.001373 Jy/beam ! Clean residual mean=0.000000 rms=0.000390 Jy/beam ! Combined flux in latest and established models = 0.0301633 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 389 components and 0.0301633 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.014917Jy sigma=0.018411 ! Fit after self-cal, rms=0.014915Jy sigma=0.018399 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.000281221 Jy ! Component: 100 - total flux cleaned = 0.000365281 Jy ! Component: 150 - total flux cleaned = 0.00038287 Jy ! Component: 200 - total flux cleaned = 0.000394165 Jy ! Total flux subtracted in 200 components = 0.000394165 Jy ! Clean residual min=-0.001486 max=0.001358 Jy/beam ! Clean residual mean=0.000000 rms=0.000391 Jy/beam ! Combined flux in latest and established models = 0.0305575 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 399 components and 0.0305575 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.014914Jy sigma=0.018418 ! Fit after self-cal, rms=0.014918Jy sigma=0.018413 ! 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.8452 mas, bmaj=2.049 mas, bpa=-0.39 degrees ! Estimated noise=20.9584 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8452 x 2.049 at -0.39 degrees (North through East) ! Clean map min=-0.0025721 max=0.0095879 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 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.7172 mas, bmaj=1.77 mas, bpa=-1.341 degrees ! Estimated noise=27.2414 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00766633 Jy ! Component: 100 - total flux cleaned = 0.0123274 Jy ! Component: 150 - total flux cleaned = 0.0147801 Jy ! Component: 200 - total flux cleaned = 0.0161645 Jy ! Total flux subtracted in 200 components = 0.0161645 Jy ! Clean residual min=-0.010960 max=0.016119 Jy/beam ! Clean residual mean=0.000000 rms=0.002527 Jy/beam ! Combined flux in latest and established models = 0.0161644 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 22 model components to the UV plane model. ! The established model now contains 22 components and 0.0161644 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8452 mas, bmaj=2.049 mas, bpa=-0.39 degrees ! Estimated noise=20.9584 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.00521187 Jy ! Component: 100 - total flux cleaned = 0.0070412 Jy ! Component: 150 - total flux cleaned = 0.00833949 Jy ! Component: 200 - total flux cleaned = 0.00920711 Jy ! Total flux subtracted in 200 components = 0.00920711 Jy ! Clean residual min=-0.013133 max=0.020421 Jy/beam ! Clean residual mean=0.000000 rms=0.004519 Jy/beam ! Combined flux in latest and established models = 0.0253716 Jy ! Adding 17 model components to the UV plane model. ! The established model now contains 36 components and 0.0253716 Jy ! Inverting map ! Added new window around map position (0, 0). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0177089 Jy ! Component: 100 - total flux cleaned = 0.0269231 Jy ! Component: 150 - total flux cleaned = 0.0330615 Jy ! Component: 200 - total flux cleaned = 0.0373311 Jy ! Total flux subtracted in 200 components = 0.0373311 Jy ! Clean residual min=-0.005837 max=0.006474 Jy/beam ! Clean residual mean=-0.000000 rms=0.001700 Jy/beam ! Combined flux in latest and established models = 0.0627027 Jy ! Adding 30 model components to the UV plane model. ! The established model now contains 66 components and 0.0627027 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.055 mas, bmaj=4.085 mas, bpa=2.399 degrees ! Estimated noise=24.6643 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00590691 Jy ! Component: 100 - total flux cleaned = 0.00885962 Jy ! Component: 150 - total flux cleaned = 0.00973314 Jy ! Component: 200 - total flux cleaned = 0.0104037 Jy ! Total flux subtracted in 200 components = 0.0104037 Jy ! Clean residual min=-0.006079 max=0.007482 Jy/beam ! Clean residual mean=-0.000001 rms=0.002107 Jy/beam ! Combined flux in latest and established models = 0.0731064 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.000497588 Jy ! Component: 100 - total flux cleaned = 0.000831915 Jy ! Component: 150 - total flux cleaned = 0.00100437 Jy ! Component: 200 - total flux cleaned = 0.00108623 Jy ! Total flux subtracted in 200 components = 0.00108623 Jy ! Clean residual min=-0.005848 max=0.007203 Jy/beam ! Clean residual mean=-0.000001 rms=0.001998 Jy/beam ! Combined flux in latest and established models = 0.0741927 Jy ! Adding 33 model components to the UV plane model. ! The established model now contains 92 components and 0.0741927 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.8452 mas, bmaj=2.049 mas, bpa=-0.39 degrees ! Estimated noise=20.9584 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0012629 Jy ! Component: 100 - total flux cleaned = -0.00135228 Jy ! Component: 150 - total flux cleaned = -0.00134953 Jy ! Component: 200 - total flux cleaned = -0.00127577 Jy ! Total flux subtracted in 200 components = -0.00127577 Jy ! Clean residual min=-0.004855 max=0.005466 Jy/beam ! Clean residual mean=-0.000000 rms=0.001185 Jy/beam ! Combined flux in latest and established models = 0.0729169 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 101 components and 0.0729169 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Correcting IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.019111Jy sigma=0.028604 ! Fit after self-cal, rms=0.018042Jy sigma=0.024514 wmodel /g0/scr/pima//bp171af_uvs/J2219-1221_X_map.mod ! Writing 101 model components to file: /g0/scr/pima//bp171af_uvs/J2219-1221_X_map.mod wobs /g0/scr/pima//bp171af_uvs/J2219-1221_X_uvs.fits ! Writing UV FITS file: /g0/scr/pima//bp171af_uvs/J2219-1221_X_uvs.fits wwins /g0/scr/pima//bp171af_uvs/J2219-1221_X_map.win ! wwins: Wrote 2 windows to /g0/scr/pima//bp171af_uvs/J2219-1221_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=16384 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -5.33301e-05 Jy ! Component: 100 - total flux cleaned = -5.5015e-05 Jy ! Component: 150 - total flux cleaned = 4.27367e-05 Jy ! Component: 200 - total flux cleaned = 0.000496805 Jy ! Component: 250 - total flux cleaned = 0.000665837 Jy ! Component: 300 - total flux cleaned = 0.000787165 Jy ! Component: 350 - total flux cleaned = 0.000787641 Jy ! Component: 400 - total flux cleaned = 0.000823121 Jy ! Component: 450 - total flux cleaned = 0.000856567 Jy ! Component: 500 - total flux cleaned = 0.000951936 Jy ! Component: 550 - total flux cleaned = 0.00092089 Jy ! Component: 600 - total flux cleaned = 0.000861838 Jy ! Component: 650 - total flux cleaned = 0.000861435 Jy ! Component: 700 - total flux cleaned = 0.000916134 Jy ! Component: 750 - total flux cleaned = 0.00088897 Jy ! Component: 800 - total flux cleaned = 0.000941496 Jy ! Component: 850 - total flux cleaned = 0.00101732 Jy ! Component: 900 - total flux cleaned = 0.00106648 Jy ! Component: 950 - total flux cleaned = 0.00116163 Jy ! Component: 1000 - total flux cleaned = 0.0011384 Jy ! Component: 1050 - total flux cleaned = 0.00116078 Jy ! Component: 1100 - total flux cleaned = 0.00118326 Jy ! Component: 1150 - total flux cleaned = 0.00116175 Jy ! Component: 1200 - total flux cleaned = 0.0012464 Jy ! Component: 1250 - total flux cleaned = 0.00132961 Jy ! Component: 1300 - total flux cleaned = 0.00147094 Jy ! Component: 1350 - total flux cleaned = 0.00153017 Jy ! Component: 1400 - total flux cleaned = 0.00162745 Jy ! Component: 1450 - total flux cleaned = 0.00174164 Jy ! Component: 1500 - total flux cleaned = 0.00181619 Jy ! Component: 1550 - total flux cleaned = 0.00187148 Jy ! Component: 1600 - total flux cleaned = 0.00194302 Jy ! Component: 1650 - total flux cleaned = 0.00197821 Jy ! Component: 1700 - total flux cleaned = 0.00199577 Jy ! Component: 1750 - total flux cleaned = 0.00201303 Jy ! Component: 1800 - total flux cleaned = 0.00197906 Jy ! Component: 1850 - total flux cleaned = 0.00194536 Jy ! Component: 1900 - total flux cleaned = 0.0018955 Jy ! Component: 1950 - total flux cleaned = 0.00189542 Jy ! Component: 2000 - total flux cleaned = 0.00179725 Jy ! Component: 2050 - total flux cleaned = 0.00176481 Jy ! Component: 2100 - total flux cleaned = 0.00171652 Jy ! Component: 2150 - total flux cleaned = 0.00162045 Jy ! Component: 2200 - total flux cleaned = 0.00150882 Jy ! Component: 2250 - total flux cleaned = 0.00141358 Jy ! Component: 2300 - total flux cleaned = 0.00133434 Jy ! Component: 2350 - total flux cleaned = 0.00120768 Jy ! Component: 2400 - total flux cleaned = 0.00111267 Jy ! Component: 2450 - total flux cleaned = 0.00103366 Jy ! Component: 2500 - total flux cleaned = 0.000922877 Jy ! Component: 2550 - total flux cleaned = 0.000827761 Jy ! Component: 2600 - total flux cleaned = 0.000700487 Jy ! Component: 2650 - total flux cleaned = 0.000668496 Jy ! Component: 2700 - total flux cleaned = 0.000540079 Jy ! Component: 2750 - total flux cleaned = 0.000394884 Jy ! Component: 2800 - total flux cleaned = 0.000264931 Jy ! Component: 2850 - total flux cleaned = 0.000183202 Jy ! Component: 2900 - total flux cleaned = 0.000100885 Jy ! Component: 2950 - total flux cleaned = 1.77771e-05 Jy ! Component: 3000 - total flux cleaned = 1.77363e-05 Jy ! Component: 3050 - total flux cleaned = -3.28035e-05 Jy ! Component: 3100 - total flux cleaned = -6.68222e-05 Jy ! Component: 3150 - total flux cleaned = -0.000135689 Jy ! Component: 3200 - total flux cleaned = -0.000205128 Jy ! Component: 3250 - total flux cleaned = -0.000275324 Jy ! Component: 3300 - total flux cleaned = -0.000293133 Jy ! Component: 3350 - total flux cleaned = -0.00034743 Jy ! Component: 3400 - total flux cleaned = -0.000401791 Jy ! Component: 3450 - total flux cleaned = -0.000420235 Jy ! Component: 3500 - total flux cleaned = -0.000438949 Jy ! Component: 3550 - total flux cleaned = -0.000438715 Jy ! Component: 3600 - total flux cleaned = -0.000514638 Jy ! Component: 3650 - total flux cleaned = -0.000572513 Jy ! Component: 3700 - total flux cleaned = -0.000572548 Jy ! Component: 3750 - total flux cleaned = -0.000532461 Jy ! Component: 3800 - total flux cleaned = -0.000612388 Jy ! Component: 3850 - total flux cleaned = -0.000673801 Jy ! Component: 3900 - total flux cleaned = -0.000632422 Jy ! Component: 3950 - total flux cleaned = -0.000611592 Jy ! Component: 4000 - total flux cleaned = -0.000590107 Jy ! Component: 4050 - total flux cleaned = -0.000633483 Jy ! Component: 4100 - total flux cleaned = -0.000677751 Jy ! Component: 4150 - total flux cleaned = -0.000744932 Jy ! Component: 4200 - total flux cleaned = -0.000745065 Jy ! Component: 4250 - total flux cleaned = -0.000814484 Jy ! Component: 4300 - total flux cleaned = -0.000838138 Jy ! Component: 4350 - total flux cleaned = -0.000862885 Jy ! Component: 4400 - total flux cleaned = -0.000887952 Jy ! Component: 4450 - total flux cleaned = -0.000837657 Jy ! Component: 4500 - total flux cleaned = -0.000837102 Jy ! Component: 4550 - total flux cleaned = -0.000836627 Jy ! Component: 4600 - total flux cleaned = -0.000809823 Jy ! Component: 4650 - total flux cleaned = -0.000754386 Jy ! Component: 4700 - total flux cleaned = -0.000697592 Jy ! Component: 4750 - total flux cleaned = -0.00063836 Jy ! Component: 4800 - total flux cleaned = -0.000577778 Jy ! Component: 4850 - total flux cleaned = -0.000545872 Jy ! Component: 4900 - total flux cleaned = -0.000513808 Jy ! Component: 4950 - total flux cleaned = -0.000480607 Jy ! Component: 5000 - total flux cleaned = -0.000479879 Jy ! Component: 5050 - total flux cleaned = -0.000515227 Jy ! Component: 5100 - total flux cleaned = -0.000553205 Jy ! Component: 5150 - total flux cleaned = -0.000479169 Jy ! Component: 5200 - total flux cleaned = -0.000440935 Jy ! Component: 5250 - total flux cleaned = -0.000401508 Jy ! Component: 5300 - total flux cleaned = -0.00027767 Jy ! Component: 5350 - total flux cleaned = -0.000277424 Jy ! Component: 5400 - total flux cleaned = -0.000278614 Jy ! Component: 5450 - total flux cleaned = -9.61228e-05 Jy ! Component: 5500 - total flux cleaned = -4.81355e-05 Jy ! Component: 5550 - total flux cleaned = -4.80714e-05 Jy ! Component: 5600 - total flux cleaned = 0.000104209 Jy ! Component: 5650 - total flux cleaned = 0.000209399 Jy ! Component: 5700 - total flux cleaned = 0.000264184 Jy ! Component: 5750 - total flux cleaned = 0.000320354 Jy ! Component: 5800 - total flux cleaned = 0.000435946 Jy ! Component: 5850 - total flux cleaned = 0.000436258 Jy ! Component: 5900 - total flux cleaned = 0.000436368 Jy ! Component: 5950 - total flux cleaned = 0.000562947 Jy ! Component: 6000 - total flux cleaned = 0.000562913 Jy ! Component: 6050 - total flux cleaned = 0.000630189 Jy ! Component: 6100 - total flux cleaned = 0.000701716 Jy ! Component: 6150 - total flux cleaned = 0.000701048 Jy ! Component: 6200 - total flux cleaned = 0.000698909 Jy ! Component: 6250 - total flux cleaned = 0.000855355 Jy ! Component: 6300 - total flux cleaned = 0.000937112 Jy ! Component: 6350 - total flux cleaned = 0.000854192 Jy ! Component: 6400 - total flux cleaned = 0.00093792 Jy ! Component: 6450 - total flux cleaned = 0.00102646 Jy ! Component: 6500 - total flux cleaned = 0.00111691 Jy ! Component: 6550 - total flux cleaned = 0.00112089 Jy ! Component: 6600 - total flux cleaned = 0.00131908 Jy ! Component: 6650 - total flux cleaned = 0.00142353 Jy ! Component: 6700 - total flux cleaned = 0.00132111 Jy ! Component: 6750 - total flux cleaned = 0.0013196 Jy ! Component: 6800 - total flux cleaned = 0.00143168 Jy ! Component: 6850 - total flux cleaned = 0.00119738 Jy ! Component: 6900 - total flux cleaned = 0.001313 Jy ! Component: 6950 - total flux cleaned = 0.00156047 Jy ! Component: 7000 - total flux cleaned = 0.00156433 Jy ! Component: 7050 - total flux cleaned = 0.00183542 Jy ! Component: 7100 - total flux cleaned = 0.00170302 Jy ! Component: 7150 - total flux cleaned = 0.00170554 Jy ! Component: 7200 - total flux cleaned = 0.00185626 Jy ! Component: 7250 - total flux cleaned = 0.0020177 Jy ! Component: 7300 - total flux cleaned = 0.00202136 Jy ! Component: 7350 - total flux cleaned = 0.00219052 Jy ! Component: 7400 - total flux cleaned = 0.00202466 Jy ! Component: 7450 - total flux cleaned = 0.00219539 Jy ! Component: 7500 - total flux cleaned = 0.00237818 Jy ! Component: 7550 - total flux cleaned = 0.00201069 Jy ! Component: 7600 - total flux cleaned = 0.00219777 Jy ! Component: 7650 - total flux cleaned = 0.00239723 Jy ! Component: 7700 - total flux cleaned = 0.00280817 Jy ! Component: 7750 - total flux cleaned = 0.00302422 Jy ! Component: 7800 - total flux cleaned = 0.00281266 Jy ! Component: 7850 - total flux cleaned = 0.00280715 Jy ! Component: 7900 - total flux cleaned = 0.0032784 Jy ! Component: 7950 - total flux cleaned = 0.00329067 Jy ! Component: 8000 - total flux cleaned = 0.00328953 Jy ! Component: 8050 - total flux cleaned = 0.00381409 Jy ! Component: 8100 - total flux cleaned = 0.00354312 Jy ! Component: 8150 - total flux cleaned = 0.00410722 Jy ! Component: 8200 - total flux cleaned = 0.00381226 Jy ! Component: 8250 - total flux cleaned = 0.00440581 Jy ! Component: 8300 - total flux cleaned = 0.00410305 Jy ! Component: 8350 - total flux cleaned = 0.00536441 Jy ! Component: 8400 - total flux cleaned = 0.005699 Jy ! Component: 8450 - total flux cleaned = 0.00501986 Jy ! Component: 8500 - total flux cleaned = 0.00501687 Jy ! Component: 8550 - total flux cleaned = 0.00538984 Jy ! Component: 8600 - total flux cleaned = 0.00614138 Jy ! Component: 8650 - total flux cleaned = 0.00614702 Jy ! Component: 8700 - total flux cleaned = 0.00655082 Jy ! Component: 8750 - total flux cleaned = 0.00612867 Jy ! Component: 8800 - total flux cleaned = 0.00656193 Jy ! Component: 8850 - total flux cleaned = 0.00789007 Jy ! Component: 8900 - total flux cleaned = 0.00694847 Jy ! Component: 8950 - total flux cleaned = 0.00789226 Jy ! Component: 9000 - total flux cleaned = 0.00740927 Jy ! Component: 9050 - total flux cleaned = 0.00844233 Jy ! Component: 9100 - total flux cleaned = 0.00791478 Jy ! Component: 9150 - total flux cleaned = 0.00898593 Jy ! Component: 9200 - total flux cleaned = 0.00955494 Jy ! Component: 9250 - total flux cleaned = 0.00897694 Jy ! Component: 9300 - total flux cleaned = 0.00956625 Jy ! Component: 9350 - total flux cleaned = 0.00892483 Jy ! Component: 9400 - total flux cleaned = 0.00956793 Jy ! Component: 9450 - total flux cleaned = 0.0102514 Jy ! Component: 9500 - total flux cleaned = 0.0116256 Jy ! Component: 9550 - total flux cleaned = 0.0109208 Jy ! Component: 9600 - total flux cleaned = 0.0124048 Jy ! Component: 9650 - total flux cleaned = 0.0116557 Jy ! Component: 9700 - total flux cleaned = 0.0132365 Jy ! Component: 9750 - total flux cleaned = 0.0132666 Jy ! Component: 9800 - total flux cleaned = 0.013244 Jy ! Component: 9850 - total flux cleaned = 0.0141471 Jy ! Component: 9900 - total flux cleaned = 0.0168647 Jy ! Component: 9950 - total flux cleaned = 0.0149936 Jy ! Component: 10000 - total flux cleaned = 0.0159239 Jy ! Component: 10050 - total flux cleaned = 0.0189258 Jy ! Component: 10100 - total flux cleaned = 0.0179013 Jy ! Component: 10150 - total flux cleaned = 0.0189547 Jy ! Component: 10200 - total flux cleaned = 0.0200585 Jy ! Component: 10250 - total flux cleaned = 0.0188812 Jy ! Component: 10300 - total flux cleaned = 0.0247961 Jy ! Component: 10350 - total flux cleaned = 0.0235977 Jy ! Component: 10400 - total flux cleaned = 0.0235983 Jy ! Component: 10450 - total flux cleaned = 0.0262445 Jy ! Component: 10500 - total flux cleaned = 0.0248908 Jy ! Component: 10550 - total flux cleaned = 0.0248214 Jy ! Component: 10600 - total flux cleaned = 0.0291469 Jy ! Component: 10650 - total flux cleaned = 0.0276926 Jy ! Component: 10700 - total flux cleaned = 0.0277042 Jy ! Component: 10750 - total flux cleaned = 0.0292976 Jy ! Component: 10800 - total flux cleaned = 0.0342245 Jy ! Component: 10850 - total flux cleaned = 0.032497 Jy ! Component: 10900 - total flux cleaned = 0.0342716 Jy ! Component: 10950 - total flux cleaned = 0.0342481 Jy ! Component: 11000 - total flux cleaned = 0.0323837 Jy ! Component: 11050 - total flux cleaned = 0.0362525 Jy ! Component: 11100 - total flux cleaned = 0.0341765 Jy ! Component: 11150 - total flux cleaned = 0.0383589 Jy ! Component: 11200 - total flux cleaned = 0.0383495 Jy ! Component: 11250 - total flux cleaned = 0.03833 Jy ! Component: 11300 - total flux cleaned = 0.0382844 Jy ! Component: 11350 - total flux cleaned = 0.0405915 Jy ! Component: 11400 - total flux cleaned = 0.0405557 Jy ! Component: 11450 - total flux cleaned = 0.0405879 Jy ! Component: 11500 - total flux cleaned = 0.0432465 Jy ! Component: 11550 - total flux cleaned = 0.0431935 Jy ! Component: 11600 - total flux cleaned = 0.0488483 Jy ! Component: 11650 - total flux cleaned = 0.0576886 Jy ! Component: 11700 - total flux cleaned = 0.0515553 Jy ! Component: 11750 - total flux cleaned = 0.0483752 Jy ! Component: 11800 - total flux cleaned = 0.0548794 Jy ! Component: 11850 - total flux cleaned = 0.0514786 Jy ! Component: 11900 - total flux cleaned = 0.0550484 Jy ! Component: 11950 - total flux cleaned = 0.0585644 Jy ! Component: 12000 - total flux cleaned = 0.0585209 Jy ! Component: 12050 - total flux cleaned = 0.0661296 Jy ! Component: 12100 - total flux cleaned = 0.066322 Jy ! Component: 12150 - total flux cleaned = 0.0620785 Jy ! Component: 12200 - total flux cleaned = 0.0704251 Jy ! Component: 12250 - total flux cleaned = 0.0748222 Jy ! Component: 12300 - total flux cleaned = 0.0749456 Jy ! Component: 12350 - total flux cleaned = 0.0703071 Jy ! Component: 12400 - total flux cleaned = 0.0606593 Jy ! Component: 12450 - total flux cleaned = 0.0708733 Jy ! Component: 12500 - total flux cleaned = 0.0654305 Jy ! Component: 12550 - total flux cleaned = 0.0653989 Jy ! Component: 12600 - total flux cleaned = 0.0766355 Jy ! Component: 12650 - total flux cleaned = 0.0825712 Jy ! Component: 12700 - total flux cleaned = 0.0827147 Jy ! Component: 12750 - total flux cleaned = 0.0703294 Jy ! Component: 12800 - total flux cleaned = 0.0765656 Jy ! Component: 12850 - total flux cleaned = 0.0765401 Jy ! Component: 12900 - total flux cleaned = 0.0832403 Jy ! Component: 12950 - total flux cleaned = 0.0759803 Jy ! Component: 13000 - total flux cleaned = 0.083353 Jy ! Component: 13050 - total flux cleaned = 0.0678734 Jy ! Component: 13100 - total flux cleaned = 0.0833837 Jy ! Component: 13150 - total flux cleaned = 0.0751678 Jy ! Component: 13200 - total flux cleaned = 0.0753424 Jy ! Component: 13250 - total flux cleaned = 0.0840711 Jy ! Component: 13300 - total flux cleaned = 0.0927907 Jy ! Component: 13350 - total flux cleaned = 0.0837596 Jy ! Component: 13400 - total flux cleaned = 0.0837908 Jy ! Component: 13450 - total flux cleaned = 0.0840978 Jy ! Component: 13500 - total flux cleaned = 0.0837913 Jy ! Component: 13550 - total flux cleaned = 0.105298 Jy ! Component: 13600 - total flux cleaned = 0.105365 Jy ! Component: 13650 - total flux cleaned = 0.105609 Jy ! Component: 13700 - total flux cleaned = 0.0820924 Jy ! Component: 13750 - total flux cleaned = 0.105761 Jy ! Component: 13800 - total flux cleaned = 0.0932042 Jy ! Component: 13850 - total flux cleaned = 0.0799407 Jy ! Component: 13900 - total flux cleaned = 0.0931418 Jy ! Component: 13950 - total flux cleaned = 0.107099 Jy ! Component: 14000 - total flux cleaned = 0.121809 Jy ! Component: 14050 - total flux cleaned = 0.151428 Jy ! Component: 14100 - total flux cleaned = 0.0897993 Jy ! Component: 14150 - total flux cleaned = 0.0895587 Jy ! Component: 14200 - total flux cleaned = 0.122261 Jy ! Component: 14250 - total flux cleaned = 0.105926 Jy ! Component: 14300 - total flux cleaned = 0.106025 Jy ! Component: 14350 - total flux cleaned = 0.106396 Jy ! Component: 14400 - total flux cleaned = 0.144386 Jy ! Component: 14450 - total flux cleaned = 0.0662178 Jy ! Component: 14500 - total flux cleaned = 0.0856981 Jy ! Component: 14550 - total flux cleaned = 0.128271 Jy ! Component: 14600 - total flux cleaned = 0.128899 Jy ! Component: 14650 - total flux cleaned = 0.106189 Jy ! Component: 14700 - total flux cleaned = 0.15295 Jy ! Component: 14750 - total flux cleaned = 0.201319 Jy ! Component: 14800 - total flux cleaned = 0.17689 Jy ! Component: 14850 - total flux cleaned = 0.151001 Jy ! Component: 14900 - total flux cleaned = 0.150306 Jy ! Component: 14950 - total flux cleaned = 0.177888 Jy ! Component: 15000 - total flux cleaned = 0.177944 Jy ! Component: 15050 - total flux cleaned = 0.120905 Jy ! Component: 15100 - total flux cleaned = 0.150267 Jy ! Component: 15150 - total flux cleaned = 0.211974 Jy ! Component: 15200 - total flux cleaned = 0.113175 Jy ! Component: 15250 - total flux cleaned = 0.214994 Jy ! Component: 15300 - total flux cleaned = 0.180315 Jy ! Component: 15350 - total flux cleaned = 0.1814 Jy ! Component: 15400 - total flux cleaned = 0.144481 Jy ! Component: 15450 - total flux cleaned = 0.182351 Jy ! Component: 15500 - total flux cleaned = 0.0231278 Jy ! Component: 15550 - total flux cleaned = 0.186583 Jy ! Component: 15600 - total flux cleaned = 0.271673 Jy ! Component: 15650 - total flux cleaned = 0.318043 Jy ! Component: 15700 - total flux cleaned = 0.273754 Jy ! Component: 15750 - total flux cleaned = 0.228267 Jy ! Component: 15800 - total flux cleaned = 0.278376 Jy ! Component: 15850 - total flux cleaned = 0.278925 Jy ! Component: 15900 - total flux cleaned = 0.226407 Jy ! Component: 15950 - total flux cleaned = 0.282231 Jy ! Component: 16000 - total flux cleaned = 0.116199 Jy ! Component: 16050 - total flux cleaned = 0.289689 Jy ! Component: 16100 - total flux cleaned = 0.110281 Jy ! Component: 16150 - total flux cleaned = 0.173704 Jy ! Component: 16200 - total flux cleaned = 0.110056 Jy ! Component: 16250 - total flux cleaned = 0.0441465 Jy ! Component: 16300 - total flux cleaned = 0.109682 Jy ! Component: 16350 - total flux cleaned = 0.178561 Jy ! Total flux subtracted in 16384 components = 0.106856 Jy ! Clean residual min=-3.965342 max=3.674196 Jy/beam ! Clean residual mean=-0.000208 rms=0.967099 Jy/beam ! Combined flux in latest and established models = 0.179776 Jy keep ! Adding 570 model components to the UV plane model. ! The established model now contains 671 components and 0.179776 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/g0/scr/pima//bp171af_uvs/J2219-1221_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8452 x 2.049 at -0.39 degrees (North through East) ! Clean map min=-3.1587 max=3.2658 Jy/beam ! Writing clean map to FITS file: /g0/scr/pima//bp171af_uvs/J2219-1221_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 3.26584 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 1.09413 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -100 100 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 3.2823 print "peak= ", int(1000*peak_flux)/1000 ! peak= 3.265 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 2.984 quit ! Quitting program ! Log file /g0/scr/pima//bp171af_uvs/J2219-1221_X_uvs.log closed on Sun Dec 8 13:19:17 2013