! Started logfile: J0828-3731_S_2004_05_08_yyk_uvs.log on Wed Dec 7 00:11:09 2005 obs J0828-3731_S_2004_05_08_yyk_uve.fits ! Reading UV FITS file: J0828-3731_S_2004_05_08_yyk_uve.fits ! AN table 1: 5 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.684444 visibilities/baseline/integration-bin. ! Found source: J0828-3731 ! ! 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.23697e+09 8e+06 1 8e+06 ! 02 2 2.26697e+09 8e+06 1 8e+06 ! 03 3 2.35697e+09 8e+06 1 8e+06 ! 04 4 2.37697e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1432 lines of history. ! ! Reading 616 visibilities. select I ! Polarization I is unavailable. ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 ![@muppet J0828-3731_S_2004_05_08_yyk] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.04 float dynam dynam = 6.0 float amp_soltime1 amp_soltime1 = 10 float amp_soltime2 amp_soltime2 = 1 float phase_soltime1 phase_soltime1=10 float phase_soltime2 phase_soltime2=0.5 float phase_soltime phase_soltime=phase_soltime1 float ampcor_selflim ampcor_selflim=1.20 float thresh thresh = 0.10 float win_mult win_mult = 1.8 float time_av time_av = 30 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 512x512 pixels with 0.500x0.500 milli-arcsec cellsize. phase_soltime=phase_soltime1 print "uvstat(rms)=",uvstat(rms),"Jy" ! uvstat(rms)= 2.12475 Jy startmod ,phase_soltime ! Applying default point source starting model. ! Performing phase self-cal over 10 minute time intervals ! 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 14 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 14 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 14 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 23 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.211571Jy sigma=16.709748 ! Fit after self-cal, rms=1.211507Jy sigma=16.708661 ! 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.364 mas, bmaj=6.892 mas, bpa=6.573 degrees ! Estimated noise=3.56702 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 1.70114 Jy ! Component: 100 - total flux cleaned = 2.17873 Jy ! Total flux subtracted in 100 components = 2.17873 Jy ! Clean residual min=-0.086527 max=0.161490 Jy/beam ! Clean residual mean=0.000693 rms=0.041926 Jy/beam ! Combined flux in latest and established models = 2.17873 Jy ! Performing phase self-cal over 10 minute time intervals ! Adding 27 model components to the UV plane model. ! The established model now contains 27 components and 2.17873 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.416510Jy sigma=5.848085 ! Fit after self-cal, rms=0.416376Jy sigma=5.847047 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 0,-2 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -2. ! 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.876 mas, bmaj=8.101 mas, bpa=2.197 degrees ! Estimated noise=2.70144 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.322172 Jy ! Component: 100 - total flux cleaned = 0.437279 Jy ! Component: 150 - total flux cleaned = 0.480986 Jy ! Component: 200 - total flux cleaned = 0.496537 Jy ! Total flux subtracted in 200 components = 0.496537 Jy ! Clean residual min=-0.073577 max=0.100829 Jy/beam ! Clean residual mean=0.001252 rms=0.031888 Jy/beam ! Combined flux in latest and established models = 2.67527 Jy ! Performing phase self-cal over 10 minute time intervals ! Adding 49 model components to the UV plane model. ! The established model now contains 69 components and 2.67527 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.298630Jy sigma=3.645816 ! Fit after self-cal, rms=0.298678Jy sigma=3.645477 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! selflims ampcor_selflim ! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2. ! Selfcal phase corrections will not be limited. gscale ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Amplitude normalization factor in sub-array 1: 0.99079 ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.06 HN 1.00* KP 0.98 ! LA 1.01 MK 1.01 NL 1.18 OV 0.96 ! PT 0.96 SC 0.93 ! ! ! Correcting IF 2. ! Amplitude normalization factor in sub-array 1: 1 ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* MK 1.00* NL 1.00* OV 1.00* ! PT 1.00* SC 1.00* ! ! ! Correcting IF 3. ! Amplitude normalization factor in sub-array 1: 1 ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* MK 1.00* NL 1.00* OV 1.00* ! PT 1.00* SC 1.00* ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Amplitude normalization factor in sub-array 1: 0.99728 ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 0.97 HN 1.00* KP 0.88 ! LA 1.00* MK 0.96 NL 1.18 OV 1.00* ! PT 0.96 SC 1.09 ! ! ! Fit before self-cal, rms=0.298678Jy sigma=3.645477 ! Fit after self-cal, rms=0.267390Jy sigma=3.230219 selfcal false, false, phase_soltime ! Performing phase self-cal over 10 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.267390Jy sigma=3.230219 ! Fit after self-cal, rms=0.267416Jy sigma=3.230017 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.877 mas, bmaj=8.074 mas, bpa=2.35 degrees ! Estimated noise=2.67309 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = -0.0169372 Jy ! Component: 100 - total flux cleaned = -0.016945 Jy ! Component: 150 - total flux cleaned = -0.0165898 Jy ! Component: 200 - total flux cleaned = -0.0162576 Jy ! Total flux subtracted in 200 components = -0.0162576 Jy ! Clean residual min=-0.059766 max=0.099118 Jy/beam ! Clean residual mean=0.001840 rms=0.025639 Jy/beam ! Combined flux in latest and established models = 2.65901 Jy ! Performing phase self-cal over 10 minute time intervals ! Adding 10 model components to the UV plane model. ! The established model now contains 77 components and 2.65901 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.264829Jy sigma=3.213856 ! Fit after self-cal, rms=0.264823Jy sigma=3.213846 ! Inverting map print "************** FINISHED AMPLITUDE GSCALE SELF-CAL **************" ! ************** FINISHED AMPLITUDE GSCALE SELF-CAL ************** ! phase_soltime=phase_soltime2 selfcal false, false, phase_soltime ! Performing phase self-cal over 0.5 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.264823Jy sigma=3.213846 ! Fit after self-cal, rms=0.262806Jy sigma=3.188991 dynam = dynam - 0.5 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 false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.000282445 Jy ! Component: 100 - total flux cleaned = 0.000290828 Jy ! Component: 150 - total flux cleaned = 0.000544369 Jy ! Component: 200 - total flux cleaned = 0.000566749 Jy ! Total flux subtracted in 200 components = 0.000566749 Jy ! Clean residual min=-0.058713 max=0.097848 Jy/beam ! Clean residual mean=0.001829 rms=0.025435 Jy/beam ! Combined flux in latest and established models = 2.65958 Jy ! Performing phase self-cal over 0.5 minute time intervals ! Adding 8 model components to the UV plane model. ! The established model now contains 80 components and 2.65958 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262448Jy sigma=3.185157 ! Fit after self-cal, rms=0.262361Jy sigma=3.183296 ! Inverting map print "*** FINISHED NATURAL WEIGHTING CLEAN with shorter phase_soltime ***" ! *** FINISHED NATURAL WEIGHTING CLEAN with shorter phase_soltime *** ! dynam = dynam - 0.25 restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.877 x 8.074 at 2.35 degrees (North through East) ! Clean map min=-0.057846 max=1.9737 Jy/beam if(peak(flux) > thresh) selfcal true, true, amp_soltime1 selfcal false, false, phase_soltime 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 false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true, amp_soltime2 selfcal false, false, phase_soltime 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 false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) print "************** FINISHED AMPLITUDE SELF-CAL **************" end if ! Performing amp+phase self-cal over 10 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262361Jy sigma=3.183296 ! Fit after self-cal, rms=0.262345Jy sigma=3.176249 ! Performing phase self-cal over 0.5 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262345Jy sigma=3.176249 ! Fit after self-cal, rms=0.262349Jy sigma=3.176016 ! Inverting map and beam ! Estimated beam: bmin=2.866 mas, bmaj=8.061 mas, bpa=2.228 degrees ! Estimated noise=2.67847 mJy/beam. ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.00915166 Jy ! Component: 100 - total flux cleaned = 0.01003 Jy ! Component: 150 - total flux cleaned = 0.0104277 Jy ! Component: 200 - total flux cleaned = 0.0108121 Jy ! Total flux subtracted in 200 components = 0.0108121 Jy ! Clean residual min=-0.058904 max=0.098232 Jy/beam ! Clean residual mean=0.001934 rms=0.025510 Jy/beam ! Combined flux in latest and established models = 2.67039 Jy ! Performing phase self-cal over 0.5 minute time intervals ! Adding 11 model components to the UV plane model. ! The established model now contains 83 components and 2.67039 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262746Jy sigma=3.171568 ! Fit after self-cal, rms=0.262742Jy sigma=3.171493 ! Inverting map ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262742Jy sigma=3.171493 ! Fit after self-cal, rms=0.262482Jy sigma=3.164595 ! Performing phase self-cal over 0.5 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262482Jy sigma=3.164595 ! Fit after self-cal, rms=0.262482Jy sigma=3.164582 ! Inverting map and beam ! Estimated beam: bmin=2.862 mas, bmaj=8.055 mas, bpa=2.182 degrees ! Estimated noise=2.68012 mJy/beam. ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.00428638 Jy ! Component: 100 - total flux cleaned = 0.00484594 Jy ! Component: 150 - total flux cleaned = 0.00520502 Jy ! Component: 200 - total flux cleaned = 0.00554855 Jy ! Total flux subtracted in 200 components = 0.00554855 Jy ! Clean residual min=-0.058988 max=0.097955 Jy/beam ! Clean residual mean=0.001945 rms=0.025481 Jy/beam ! Combined flux in latest and established models = 2.67594 Jy ! Performing phase self-cal over 0.5 minute time intervals ! Adding 13 model components to the UV plane model. ! The established model now contains 84 components and 2.67594 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.262618Jy sigma=3.162057 ! Fit after self-cal, rms=0.262618Jy sigma=3.162039 ! Inverting map ! ************** FINISHED AMPLITUDE SELF-CAL ************** uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.4 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=9.421 mas, bmaj=30.74 mas, bpa=5.76 degrees ! Estimated noise=3.56239 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal false,false,phase_soltime new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.0144502 Jy ! Component: 100 - total flux cleaned = 0.0145089 Jy ! Component: 150 - total flux cleaned = 0.0149075 Jy ! Component: 200 - total flux cleaned = 0.0149792 Jy ! Total flux subtracted in 200 components = 0.0149792 Jy ! Clean residual min=-0.088070 max=0.144148 Jy/beam ! Clean residual mean=0.003720 rms=0.044788 Jy/beam ! Combined flux in latest and established models = 2.69092 Jy ! Performing phase self-cal over 0.5 minute time intervals ! Adding 7 model components to the UV plane model. ! The established model now contains 87 components and 2.69092 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.264116Jy sigma=3.168922 ! Fit after self-cal, rms=0.263977Jy sigma=3.167396 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.4 print "************** FINISHED TAPERING SELF-CAL CLEAN **************" ! ************** FINISHED TAPERING SELF-CAL CLEAN ************** ! print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true,true,true ! clrmod: Cleared the established, tentative and continuum 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.351 mas, bmaj=6.883 mas, bpa=6.517 degrees ! Estimated noise=3.54895 mJy/beam. ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = 1.69636 Jy ! Component: 100 - total flux cleaned = 2.18032 Jy ! Component: 150 - total flux cleaned = 2.40438 Jy ! Component: 200 - total flux cleaned = 2.51507 Jy ! Total flux subtracted in 200 components = 2.51507 Jy ! Clean residual min=-0.058370 max=0.088641 Jy/beam ! Clean residual mean=0.001002 rms=0.023019 Jy/beam ! Combined flux in latest and established models = 2.51507 Jy uvw 0,-2 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -2. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Adding 54 model components to the UV plane model. ! The established model now contains 54 components and 2.51507 Jy ! Inverting map and beam ! Estimated beam: bmin=2.862 mas, bmaj=8.055 mas, bpa=2.182 degrees ! Estimated noise=2.68012 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.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.0992703 Jy ! Component: 100 - total flux cleaned = 0.137026 Jy ! Component: 150 - total flux cleaned = 0.15195 Jy ! Component: 200 - total flux cleaned = 0.155478 Jy ! Total flux subtracted in 200 components = 0.155478 Jy ! Clean residual min=-0.059042 max=0.097617 Jy/beam ! Clean residual mean=0.001980 rms=0.025629 Jy/beam ! Combined flux in latest and established models = 2.67055 Jy ! Adding 39 model components to the UV plane model. ! The established model now contains 88 components and 2.67055 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=9.421 mas, bmaj=30.74 mas, bpa=5.76 degrees ! Estimated noise=3.56239 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.04 cutoff=0 ! Component: 050 - total flux cleaned = 0.0162096 Jy ! Component: 100 - total flux cleaned = 0.0162344 Jy ! Component: 150 - total flux cleaned = 0.0163076 Jy ! Component: 200 - total flux cleaned = 0.0163783 Jy ! Total flux subtracted in 200 components = 0.0163783 Jy ! Clean residual min=-0.089098 max=0.143634 Jy/beam ! Clean residual mean=0.003779 rms=0.045009 Jy/beam ! Combined flux in latest and established models = 2.68693 Jy ! Adding 7 model components to the UV plane model. ! The established model now contains 91 components and 2.68693 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.862 mas, bmaj=8.055 mas, bpa=2.182 degrees ! Estimated noise=2.68012 mJy/beam. ! clean: niter=200 gain=0.04 cutoff=0 ! Component: 050 - total flux cleaned = -0.0123213 Jy ! Component: 100 - total flux cleaned = -0.0123225 Jy ! Component: 150 - total flux cleaned = -0.0123294 Jy ! Component: 200 - total flux cleaned = -0.0120327 Jy ! Total flux subtracted in 200 components = -0.0120327 Jy ! Clean residual min=-0.058677 max=0.097923 Jy/beam ! Clean residual mean=0.001959 rms=0.025447 Jy/beam ! Combined flux in latest and established models = 2.6749 Jy selfcal false, false, phase_soltime ! Performing phase self-cal over 0.5 minute time intervals ! Adding 14 model components to the UV plane model. ! The established model now contains 102 components and 2.6749 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.263848Jy sigma=3.168948 ! Fit after self-cal, rms=0.263851Jy sigma=3.168879 wmodel J0828-3731_S_2004_05_08_yyk\_map.mod ! Writing 102 model components to file: J0828-3731_S_2004_05_08_yyk_map.mod wobs J0828-3731_S_2004_05_08_yyk\_uvs.fits ! Writing UV FITS file: J0828-3731_S_2004_05_08_yyk_uvs.fits wwins J0828-3731_S_2004_05_08_yyk\_map.win ! wwins: Wrote 1 windows to J0828-3731_S_2004_05_08_yyk_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=2048 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.0394963 Jy ! Component: 100 - total flux cleaned = 0.0679655 Jy ! Component: 150 - total flux cleaned = 0.0881603 Jy ! Component: 200 - total flux cleaned = 0.104487 Jy ! Component: 250 - total flux cleaned = 0.116958 Jy ! Component: 300 - total flux cleaned = 0.125426 Jy ! Component: 350 - total flux cleaned = 0.132665 Jy ! Component: 400 - total flux cleaned = 0.142975 Jy ! Component: 450 - total flux cleaned = 0.150145 Jy ! Component: 500 - total flux cleaned = 0.15943 Jy ! Component: 550 - total flux cleaned = 0.166507 Jy ! Component: 600 - total flux cleaned = 0.174386 Jy ! Component: 650 - total flux cleaned = 0.18136 Jy ! Component: 700 - total flux cleaned = 0.187006 Jy ! Component: 750 - total flux cleaned = 0.190944 Jy ! Component: 800 - total flux cleaned = 0.195662 Jy ! Component: 850 - total flux cleaned = 0.200195 Jy ! Component: 900 - total flux cleaned = 0.204137 Jy ! Component: 950 - total flux cleaned = 0.207522 Jy ! Component: 1000 - total flux cleaned = 0.211622 Jy ! Component: 1050 - total flux cleaned = 0.215588 Jy ! Component: 1100 - total flux cleaned = 0.217901 Jy ! Component: 1150 - total flux cleaned = 0.221286 Jy ! Component: 1200 - total flux cleaned = 0.224938 Jy ! Component: 1250 - total flux cleaned = 0.228142 Jy ! Component: 1300 - total flux cleaned = 0.229537 Jy ! Component: 1350 - total flux cleaned = 0.23158 Jy ! Component: 1400 - total flux cleaned = 0.233244 Jy ! Component: 1450 - total flux cleaned = 0.234216 Jy ! Component: 1500 - total flux cleaned = 0.235176 Jy ! Component: 1550 - total flux cleaned = 0.236438 Jy ! Component: 1600 - total flux cleaned = 0.237979 Jy ! Component: 1650 - total flux cleaned = 0.238582 Jy ! Component: 1700 - total flux cleaned = 0.240074 Jy ! Component: 1750 - total flux cleaned = 0.240073 Jy ! Component: 1800 - total flux cleaned = 0.240364 Jy ! Component: 1850 - total flux cleaned = 0.240362 Jy ! Component: 1900 - total flux cleaned = 0.239797 Jy ! Component: 1950 - total flux cleaned = 0.23952 Jy ! Component: 2000 - total flux cleaned = 0.23897 Jy ! Total flux subtracted in 2048 components = 0.238154 Jy ! Clean residual min=-0.016719 max=0.016027 Jy/beam ! Clean residual mean=-0.000197 rms=0.006671 Jy/beam ! Combined flux in latest and established models = 2.91305 Jy keep ! Adding 206 model components to the UV plane model. ! The established model now contains 308 components and 2.91305 Jy ![Exited script file: muppet] wmap "J0828-3731_S_2004_05_08_yyk\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.862 x 8.055 at 2.182 degrees (North through East) ! Clean map min=-0.044678 max=1.9698 Jy/beam ! Writing clean map to FITS file: J0828-3731_S_2004_05_08_yyk_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.96981 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00665026 device "J0828-3731_S_2004_05_08_yyk\_map.eps/VPS" ! Attempting to open device: 'J0828-3731_S_2004_05_08_yyk_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.01283 1.01283 2.02566 4.05131 8.10262 16.2052 32.4105 64.821 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0199 print "peak=",int(peak_flux*1000)/1000 ! peak= 1.969 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.862 x 8.055 at 2.182 degrees (North through East) ! Clean map min=-0.044678 max=1.9698 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J0828-3731_S_2004_05_08_yyk\_map.gif/gif" ! Attempting to open device: 'J0828-3731_S_2004_05_08_yyk_map.gif/gif' mapplot cln quit ! Quitting program ! Log file J0828-3731_S_2004_05_08_yyk_uvs.log closed on Wed Dec 7 00:11:13 2005