! Started logfile: /g0/scr/pima//bp171af_uvs/J2037-0729_X_uvs.log on Sun Dec 8 11:59:40 2013 obs /g0/scr/pima//bp171af_uvs/J2037-0729_X_uva.fits ! Reading UV FITS file: /g0/scr/pima//bp171af_uvs/J2037-0729_X_uva.fits ! AN table 1: 4 integrations on 28 of 28 possible baselines. ! Apparent sampling: 0.526786 visibilities/baseline/integration-bin. ! Found source: J2037-0729 ! ! 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 472 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/J2037-0729_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. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 2 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.002441Jy sigma=1.424198 ! Fit after self-cal, rms=0.975831Jy sigma=1.340582 ! 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.8214 mas, bmaj=1.751 mas, bpa=-4.174 degrees ! Estimated noise=37.3114 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 = 0.0295367 Jy ! Component: 100 - total flux cleaned = 0.036524 Jy ! Total flux subtracted in 100 components = 0.036524 Jy ! Clean residual min=-0.012874 max=0.013662 Jy/beam ! Clean residual mean=-0.000001 rms=0.002606 Jy/beam ! Combined flux in latest and established models = 0.036524 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 11 components and 0.036524 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.035471Jy sigma=0.052644 ! Fit after self-cal, rms=0.033167Jy sigma=0.050655 ! 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.9447 mas, bmaj=1.918 mas, bpa=-1.767 degrees ! Estimated noise=28.2662 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.00377485 Jy ! Component: 100 - total flux cleaned = 0.00305866 Jy ! Component: 150 - total flux cleaned = 0.00264959 Jy ! Component: 200 - total flux cleaned = 0.00229671 Jy ! Total flux subtracted in 200 components = 0.00229671 Jy ! Clean residual min=-0.011405 max=0.011113 Jy/beam ! Clean residual mean=-0.000001 rms=0.002352 Jy/beam ! Combined flux in latest and established models = 0.0388207 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 20 components and 0.0388207 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.033450Jy sigma=0.041815 ! Fit after self-cal, rms=0.033350Jy sigma=0.041757 ! 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.61 HN 1.00 KP 0.67 MK 1.36 ! NL 1.87 OV 1.31 PT 0.97 SC 0.61 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.74 HN 0.95 KP 0.78 MK 0.85 ! NL 1.45 OV 1.37 PT 0.71 SC 0.61 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.65 HN 1.08 KP 0.86 MK 0.82 ! NL 0.90 OV 1.00 PT 0.82 SC 1.06 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.69 HN 2.34 KP 0.70 MK 0.97 ! NL 0.73 OV 0.86 PT 1.77 SC 2.54 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 0.82 HN 0.92 KP 0.79 MK 0.79 ! NL 1.31 OV 0.77 PT 0.82 SC 0.80 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 0.62 HN 1.08 KP 0.61 MK 1.14 ! NL 1.41 OV 1.07 PT 0.60 SC 0.80 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.66 HN 0.87 KP 2.04 MK 1.14 ! NL 1.11 OV 1.76 PT 0.77 SC 0.80 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 0.69 HN 0.65 KP 3.84 MK 1.38 ! NL 0.88 OV 1.26 PT 1.43 SC 1.08 ! ! ! Fit before self-cal, rms=0.033350Jy sigma=0.041757 ! Fit after self-cal, rms=0.037935Jy sigma=0.026822 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9388 mas, bmaj=1.876 mas, bpa=-0.5105 degrees ! Estimated noise=22.963 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.00488466 Jy ! Component: 100 - total flux cleaned = -0.00641498 Jy ! Component: 150 - total flux cleaned = -0.00687599 Jy ! Component: 200 - total flux cleaned = -0.00729745 Jy ! Total flux subtracted in 200 components = -0.00729745 Jy ! Clean residual min=-0.005253 max=0.005019 Jy/beam ! Clean residual mean=-0.000000 rms=0.000969 Jy/beam ! Combined flux in latest and established models = 0.0315233 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 26 components and 0.0315233 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.035724Jy sigma=0.027363 ! Fit after self-cal, rms=0.036822Jy sigma=0.027171 ! 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.036822Jy sigma=0.027171 ! Fit after self-cal, rms=0.034796Jy sigma=0.025493 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9297 mas, bmaj=1.862 mas, bpa=-0.1278 degrees ! Estimated noise=20.4223 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00374327 Jy ! Component: 100 - total flux cleaned = -0.00460271 Jy ! Component: 150 - total flux cleaned = -0.00479881 Jy ! Component: 200 - total flux cleaned = -0.00502759 Jy ! Total flux subtracted in 200 components = -0.00502759 Jy ! Clean residual min=-0.004458 max=0.004068 Jy/beam ! Clean residual mean=-0.000000 rms=0.000798 Jy/beam ! Combined flux in latest and established models = 0.0264957 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 30 components and 0.0264957 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.033182Jy sigma=0.026390 ! Fit after self-cal, rms=0.033169Jy sigma=0.026322 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.000473974 Jy ! Component: 100 - total flux cleaned = -0.000633482 Jy ! Component: 150 - total flux cleaned = -0.000740378 Jy ! Component: 200 - total flux cleaned = -0.00090896 Jy ! Total flux subtracted in 200 components = -0.00090896 Jy ! Clean residual min=-0.004471 max=0.003718 Jy/beam ! Clean residual mean=-0.000000 rms=0.000788 Jy/beam ! Combined flux in latest and established models = 0.0255867 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 33 components and 0.0255867 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.032988Jy sigma=0.026166 ! Fit after self-cal, rms=0.032970Jy sigma=0.026111 ! 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.032970Jy sigma=0.026111 ! Fit after self-cal, rms=0.033367Jy sigma=0.024977 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.912 mas, bmaj=1.802 mas, bpa=1.52 degrees ! Estimated noise=18.8321 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00130503 Jy ! Component: 100 - total flux cleaned = -0.00117452 Jy ! Component: 150 - total flux cleaned = -0.00124565 Jy ! Component: 200 - total flux cleaned = -0.00134507 Jy ! Total flux subtracted in 200 components = -0.00134507 Jy ! Clean residual min=-0.003676 max=0.003190 Jy/beam ! Clean residual mean=-0.000000 rms=0.000643 Jy/beam ! Combined flux in latest and established models = 0.0242417 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 37 components and 0.0242417 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.032762Jy sigma=0.025346 ! Fit after self-cal, rms=0.032624Jy sigma=0.025287 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.000282529 Jy ! Component: 100 - total flux cleaned = -0.00037254 Jy ! Component: 150 - total flux cleaned = -0.000428751 Jy ! Component: 200 - total flux cleaned = -0.000504889 Jy ! Total flux subtracted in 200 components = -0.000504889 Jy ! Clean residual min=-0.003657 max=0.003079 Jy/beam ! Clean residual mean=-0.000000 rms=0.000638 Jy/beam ! Combined flux in latest and established models = 0.0237368 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 39 components and 0.0237368 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.032538Jy sigma=0.025160 ! Fit after self-cal, rms=0.032523Jy sigma=0.025130 ! 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.032523Jy sigma=0.025130 ! Fit after self-cal, rms=0.032387Jy sigma=0.024363 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.719 mas, bmaj=2.952 mas, bpa=8.77 degrees ! Estimated noise=24.9552 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00168658 Jy ! Component: 100 - total flux cleaned = -0.00242963 Jy ! Component: 150 - total flux cleaned = -0.00250756 Jy ! Component: 200 - total flux cleaned = -0.00248671 Jy ! Total flux subtracted in 200 components = -0.00248671 Jy ! Clean residual min=-0.002477 max=0.002621 Jy/beam ! Clean residual mean=-0.000000 rms=0.000640 Jy/beam ! Combined flux in latest and established models = 0.0212501 Jy selfcal ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 40 components and 0.0212501 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.032108Jy sigma=0.025377 ! Fit after self-cal, rms=0.031885Jy sigma=0.025308 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 = 1.45714e-05 Jy ! Component: 100 - total flux cleaned = -2.27038e-05 Jy ! Component: 150 - total flux cleaned = -6.70419e-06 Jy ! Component: 200 - total flux cleaned = 2.03408e-05 Jy ! Total flux subtracted in 200 components = 2.03408e-05 Jy ! Clean residual min=-0.002364 max=0.002594 Jy/beam ! Clean residual mean=-0.000000 rms=0.000628 Jy/beam ! Combined flux in latest and established models = 0.0212704 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 41 components and 0.0212704 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.031782Jy sigma=0.026042 ! Fit after self-cal, rms=0.031765Jy sigma=0.026003 ! 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.9025 mas, bmaj=1.766 mas, bpa=2.37 degrees ! Estimated noise=17.6196 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9025 x 1.766 at 2.37 degrees (North through East) ! Clean map min=-0.010135 max=0.035105 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.7939 mas, bmaj=1.602 mas, bpa=0.7012 degrees ! Estimated noise=23.054 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.019317 Jy ! Component: 100 - total flux cleaned = 0.0247319 Jy ! Component: 150 - total flux cleaned = 0.0259349 Jy ! Component: 200 - total flux cleaned = 0.0266451 Jy ! Total flux subtracted in 200 components = 0.0266451 Jy ! Clean residual min=-0.006194 max=0.006194 Jy/beam ! Clean residual mean=-0.000000 rms=0.001123 Jy/beam ! Combined flux in latest and established models = 0.0266451 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 19 model components to the UV plane model. ! The established model now contains 19 components and 0.0266451 Jy ! Inverting map and beam ! Estimated beam: bmin=0.9025 mas, bmaj=1.766 mas, bpa=2.37 degrees ! Estimated noise=17.6196 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.00145088 Jy ! Component: 100 - total flux cleaned = 0.000945345 Jy ! Component: 150 - total flux cleaned = 0.000735601 Jy ! Component: 200 - total flux cleaned = 0.000355535 Jy ! Total flux subtracted in 200 components = 0.000355535 Jy ! Clean residual min=-0.004725 max=0.005280 Jy/beam ! Clean residual mean=-0.000000 rms=0.000935 Jy/beam ! Combined flux in latest and established models = 0.0270006 Jy ! Adding 21 model components to the UV plane model. ! The established model now contains 29 components and 0.0270006 Jy ! Inverting map ! Added new window around map position (-0.9, -3.75). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00202065 Jy ! Component: 100 - total flux cleaned = 0.00249463 Jy ! Component: 150 - total flux cleaned = 0.00287157 Jy ! Component: 200 - total flux cleaned = 0.00299708 Jy ! Total flux subtracted in 200 components = 0.00299708 Jy ! Clean residual min=-0.002488 max=0.002995 Jy/beam ! Clean residual mean=-0.000000 rms=0.000554 Jy/beam ! Combined flux in latest and established models = 0.0299977 Jy ! Adding 18 model components to the UV plane model. ! The established model now contains 45 components and 0.0299977 Jy ! Inverting map ! Added new window around map position (3.9, 12.9). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00151357 Jy ! Component: 100 - total flux cleaned = 0.00125078 Jy ! Component: 150 - total flux cleaned = 0.00102078 Jy ! Component: 200 - total flux cleaned = 0.000864239 Jy ! Total flux subtracted in 200 components = 0.000864239 Jy ! Clean residual min=-0.001840 max=0.001929 Jy/beam ! Clean residual mean=-0.000000 rms=0.000419 Jy/beam ! Combined flux in latest and established models = 0.030862 Jy ! Adding 19 model components to the UV plane model. ! The established model now contains 60 components and 0.030862 Jy ! Inverting map ! Added new window around map position (6.15, -73.05). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000476152 Jy ! Component: 100 - total flux cleaned = 0.00103002 Jy ! Component: 150 - total flux cleaned = 0.00136841 Jy ! Component: 200 - total flux cleaned = 0.0015772 Jy ! Total flux subtracted in 200 components = 0.0015772 Jy ! Clean residual min=-0.001747 max=0.001963 Jy/beam ! Clean residual mean=-0.000000 rms=0.000378 Jy/beam ! Combined flux in latest and established models = 0.0324391 Jy ! Adding 25 model components to the UV plane model. ! The established model now contains 82 components and 0.0324391 Jy ! Inverting map ! Added new window around map position (104.85, -64.35). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0016279 Jy ! Component: 100 - total flux cleaned = 0.00229961 Jy ! Component: 150 - total flux cleaned = 0.00279902 Jy ! Component: 200 - total flux cleaned = 0.0033149 Jy ! Total flux subtracted in 200 components = 0.0033149 Jy ! Clean residual min=-0.001525 max=0.001587 Jy/beam ! Clean residual mean=-0.000000 rms=0.000323 Jy/beam ! Combined flux in latest and established models = 0.0357541 Jy ! Adding 27 model components to the UV plane model. ! The established model now contains 106 components and 0.0357541 Jy ! Inverting map ! Added new window around map position (-70.2, -22.95). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00145292 Jy ! Component: 100 - total flux cleaned = 0.00207243 Jy ! Component: 150 - total flux cleaned = 0.00239262 Jy ! Component: 200 - total flux cleaned = 0.00278949 Jy ! Total flux subtracted in 200 components = 0.00278949 Jy ! Clean residual min=-0.001492 max=0.001335 Jy/beam ! Clean residual mean=-0.000000 rms=0.000305 Jy/beam ! Combined flux in latest and established models = 0.0385435 Jy ! Adding 34 model components to the UV plane model. ! The established model now contains 134 components and 0.0385435 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=1.719 mas, bmaj=2.952 mas, bpa=8.77 degrees ! Estimated noise=24.9552 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 2.8731e-05 Jy ! Component: 100 - total flux cleaned = 0.000297814 Jy ! Component: 150 - total flux cleaned = 0.000439449 Jy ! Component: 200 - total flux cleaned = 0.000484309 Jy ! Total flux subtracted in 200 components = 0.000484309 Jy ! Clean residual min=-0.001538 max=0.001396 Jy/beam ! Clean residual mean=-0.000000 rms=0.000331 Jy/beam ! Combined flux in latest and established models = 0.0390278 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 = -7.67252e-05 Jy ! Component: 100 - total flux cleaned = -0.00016822 Jy ! Component: 150 - total flux cleaned = -0.000221158 Jy ! Component: 200 - total flux cleaned = -0.000303856 Jy ! Total flux subtracted in 200 components = -0.000303856 Jy ! Clean residual min=-0.001561 max=0.001382 Jy/beam ! Clean residual mean=-0.000000 rms=0.000318 Jy/beam ! Combined flux in latest and established models = 0.038724 Jy ! Adding 37 model components to the UV plane model. ! The established model now contains 163 components and 0.038724 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.9025 mas, bmaj=1.766 mas, bpa=2.37 degrees ! Estimated noise=17.6196 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00032516 Jy ! Component: 100 - total flux cleaned = 0.0006552 Jy ! Component: 150 - total flux cleaned = 0.000884809 Jy ! Component: 200 - total flux cleaned = 0.000980089 Jy ! Total flux subtracted in 200 components = 0.000980089 Jy ! Clean residual min=-0.001382 max=0.001249 Jy/beam ! Clean residual mean=-0.000000 rms=0.000286 Jy/beam ! Combined flux in latest and established models = 0.0397041 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 181 components and 0.0397041 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.029864Jy sigma=0.022210 ! Fit after self-cal, rms=0.029587Jy sigma=0.021738 wmodel /g0/scr/pima//bp171af_uvs/J2037-0729_X_map.mod ! Writing 181 model components to file: /g0/scr/pima//bp171af_uvs/J2037-0729_X_map.mod wobs /g0/scr/pima//bp171af_uvs/J2037-0729_X_uvs.fits ! Writing UV FITS file: /g0/scr/pima//bp171af_uvs/J2037-0729_X_uvs.fits wwins /g0/scr/pima//bp171af_uvs/J2037-0729_X_map.win ! wwins: Wrote 6 windows to /g0/scr/pima//bp171af_uvs/J2037-0729_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 = -3.5963e-05 Jy ! Component: 100 - total flux cleaned = 6.27004e-06 Jy ! Component: 150 - total flux cleaned = 0.000190489 Jy ! Component: 200 - total flux cleaned = 0.000287792 Jy ! Component: 250 - total flux cleaned = 0.00036284 Jy ! Component: 300 - total flux cleaned = 0.000362915 Jy ! Component: 350 - total flux cleaned = 0.000415547 Jy ! Component: 400 - total flux cleaned = 0.000466392 Jy ! Component: 450 - total flux cleaned = 0.000400207 Jy ! Component: 500 - total flux cleaned = 0.000400647 Jy ! Component: 550 - total flux cleaned = 0.000416439 Jy ! Component: 600 - total flux cleaned = 0.000416695 Jy ! Component: 650 - total flux cleaned = 0.000355267 Jy ! Component: 700 - total flux cleaned = 0.000310014 Jy ! Component: 750 - total flux cleaned = 0.00022086 Jy ! Component: 800 - total flux cleaned = 0.000191543 Jy ! Component: 850 - total flux cleaned = 0.000119333 Jy ! Component: 900 - total flux cleaned = 9.10388e-05 Jy ! Component: 950 - total flux cleaned = 6.29006e-05 Jy ! Component: 1000 - total flux cleaned = -6.63081e-06 Jy ! Component: 1050 - total flux cleaned = -2.00954e-05 Jy ! Component: 1100 - total flux cleaned = -3.37716e-05 Jy ! Component: 1150 - total flux cleaned = -2.02381e-05 Jy ! Component: 1200 - total flux cleaned = -2.02614e-05 Jy ! Component: 1250 - total flux cleaned = -2.03816e-05 Jy ! Component: 1300 - total flux cleaned = -2.03285e-05 Jy ! Component: 1350 - total flux cleaned = 5.92427e-06 Jy ! Component: 1400 - total flux cleaned = 4.4823e-05 Jy ! Component: 1450 - total flux cleaned = 7.06247e-05 Jy ! Component: 1500 - total flux cleaned = 0.000109018 Jy ! Component: 1550 - total flux cleaned = 0.000109142 Jy ! Component: 1600 - total flux cleaned = 0.000147068 Jy ! Component: 1650 - total flux cleaned = 0.000184645 Jy ! Component: 1700 - total flux cleaned = 0.000197157 Jy ! Component: 1750 - total flux cleaned = 0.000160005 Jy ! Component: 1800 - total flux cleaned = 0.000135288 Jy ! Component: 1850 - total flux cleaned = 0.000122957 Jy ! Component: 1900 - total flux cleaned = 8.63268e-05 Jy ! Component: 1950 - total flux cleaned = 6.20452e-05 Jy ! Component: 2000 - total flux cleaned = 7.42169e-05 Jy ! Component: 2050 - total flux cleaned = 2.57888e-05 Jy ! Component: 2100 - total flux cleaned = 2.57248e-05 Jy ! Component: 2150 - total flux cleaned = 1.37227e-05 Jy ! Component: 2200 - total flux cleaned = -1.03188e-05 Jy ! Component: 2250 - total flux cleaned = -1.03118e-05 Jy ! Component: 2300 - total flux cleaned = 1.73517e-06 Jy ! Component: 2350 - total flux cleaned = -1.0275e-05 Jy ! Component: 2400 - total flux cleaned = -3.42992e-05 Jy ! Component: 2450 - total flux cleaned = -9.43649e-05 Jy ! Component: 2500 - total flux cleaned = -0.000118446 Jy ! Component: 2550 - total flux cleaned = -0.000142468 Jy ! Component: 2600 - total flux cleaned = -0.0001786 Jy ! Component: 2650 - total flux cleaned = -0.000214636 Jy ! Component: 2700 - total flux cleaned = -0.000250818 Jy ! Component: 2750 - total flux cleaned = -0.000299099 Jy ! Component: 2800 - total flux cleaned = -0.000420309 Jy ! Component: 2850 - total flux cleaned = -0.000541842 Jy ! Component: 2900 - total flux cleaned = -0.000639485 Jy ! Component: 2950 - total flux cleaned = -0.000700751 Jy ! Component: 3000 - total flux cleaned = -0.000860714 Jy ! Component: 3050 - total flux cleaned = -0.000922604 Jy ! Component: 3100 - total flux cleaned = -0.00102186 Jy ! Component: 3150 - total flux cleaned = -0.00110924 Jy ! Component: 3200 - total flux cleaned = -0.00113434 Jy ! Component: 3250 - total flux cleaned = -0.00114699 Jy ! Component: 3300 - total flux cleaned = -0.00117246 Jy ! Component: 3350 - total flux cleaned = -0.0011723 Jy ! Component: 3400 - total flux cleaned = -0.00118519 Jy ! Component: 3450 - total flux cleaned = -0.00121115 Jy ! Component: 3500 - total flux cleaned = -0.00121112 Jy ! Component: 3550 - total flux cleaned = -0.00123755 Jy ! Component: 3600 - total flux cleaned = -0.00127774 Jy ! Component: 3650 - total flux cleaned = -0.00130477 Jy ! Component: 3700 - total flux cleaned = -0.00130491 Jy ! Component: 3750 - total flux cleaned = -0.00134598 Jy ! Component: 3800 - total flux cleaned = -0.0014014 Jy ! Component: 3850 - total flux cleaned = -0.00140159 Jy ! Component: 3900 - total flux cleaned = -0.00141574 Jy ! Component: 3950 - total flux cleaned = -0.00141594 Jy ! Component: 4000 - total flux cleaned = -0.00138734 Jy ! Component: 4050 - total flux cleaned = -0.00138734 Jy ! Component: 4100 - total flux cleaned = -0.00137263 Jy ! Component: 4150 - total flux cleaned = -0.00137265 Jy ! Component: 4200 - total flux cleaned = -0.00135753 Jy ! Component: 4250 - total flux cleaned = -0.00132686 Jy ! Component: 4300 - total flux cleaned = -0.0013267 Jy ! Component: 4350 - total flux cleaned = -0.00127946 Jy ! Component: 4400 - total flux cleaned = -0.00126362 Jy ! Component: 4450 - total flux cleaned = -0.00119893 Jy ! Component: 4500 - total flux cleaned = -0.00118235 Jy ! Component: 4550 - total flux cleaned = -0.00113256 Jy ! Component: 4600 - total flux cleaned = -0.00113258 Jy ! Component: 4650 - total flux cleaned = -0.00108163 Jy ! Component: 4700 - total flux cleaned = -0.00102984 Jy ! Component: 4750 - total flux cleaned = -0.000994729 Jy ! Component: 4800 - total flux cleaned = -0.000941602 Jy ! Component: 4850 - total flux cleaned = -0.000887658 Jy ! Component: 4900 - total flux cleaned = -0.000887769 Jy ! Component: 4950 - total flux cleaned = -0.00086936 Jy ! Component: 5000 - total flux cleaned = -0.000907168 Jy ! Component: 5050 - total flux cleaned = -0.00088822 Jy ! Component: 5100 - total flux cleaned = -0.000907332 Jy ! Component: 5150 - total flux cleaned = -0.000966048 Jy ! Component: 5200 - total flux cleaned = -0.000966274 Jy ! Component: 5250 - total flux cleaned = -0.000965921 Jy ! Component: 5300 - total flux cleaned = -0.00102724 Jy ! Component: 5350 - total flux cleaned = -0.00111053 Jy ! Component: 5400 - total flux cleaned = -0.00113184 Jy ! Component: 5450 - total flux cleaned = -0.00117485 Jy ! Component: 5500 - total flux cleaned = -0.00117502 Jy ! Component: 5550 - total flux cleaned = -0.00117523 Jy ! Component: 5600 - total flux cleaned = -0.00124188 Jy ! Component: 5650 - total flux cleaned = -0.00126442 Jy ! Component: 5700 - total flux cleaned = -0.00135684 Jy ! Component: 5750 - total flux cleaned = -0.00140437 Jy ! Component: 5800 - total flux cleaned = -0.00145236 Jy ! Component: 5850 - total flux cleaned = -0.00152527 Jy ! Component: 5900 - total flux cleaned = -0.00162364 Jy ! Component: 5950 - total flux cleaned = -0.00164856 Jy ! Component: 6000 - total flux cleaned = -0.00172471 Jy ! Component: 6050 - total flux cleaned = -0.00169882 Jy ! Component: 6100 - total flux cleaned = -0.00172462 Jy ! Component: 6150 - total flux cleaned = -0.0017781 Jy ! Component: 6200 - total flux cleaned = -0.00177739 Jy ! Component: 6250 - total flux cleaned = -0.001832 Jy ! Component: 6300 - total flux cleaned = -0.00194438 Jy ! Component: 6350 - total flux cleaned = -0.00200165 Jy ! Component: 6400 - total flux cleaned = -0.00200128 Jy ! Component: 6450 - total flux cleaned = -0.00208955 Jy ! Component: 6500 - total flux cleaned = -0.00211918 Jy ! Component: 6550 - total flux cleaned = -0.00218031 Jy ! Component: 6600 - total flux cleaned = -0.00230402 Jy ! Component: 6650 - total flux cleaned = -0.00227234 Jy ! Component: 6700 - total flux cleaned = -0.00246565 Jy ! Component: 6750 - total flux cleaned = -0.00243294 Jy ! Component: 6800 - total flux cleaned = -0.00263365 Jy ! Component: 6850 - total flux cleaned = -0.00263509 Jy ! Component: 6900 - total flux cleaned = -0.00260018 Jy ! Component: 6950 - total flux cleaned = -0.00270568 Jy ! Component: 7000 - total flux cleaned = -0.00274166 Jy ! Component: 7050 - total flux cleaned = -0.00277742 Jy ! Component: 7100 - total flux cleaned = -0.00299987 Jy ! Component: 7150 - total flux cleaned = -0.00303835 Jy ! Component: 7200 - total flux cleaned = -0.00319204 Jy ! Component: 7250 - total flux cleaned = -0.00319252 Jy ! Component: 7300 - total flux cleaned = -0.00327176 Jy ! Component: 7350 - total flux cleaned = -0.00339266 Jy ! Component: 7400 - total flux cleaned = -0.00355596 Jy ! Component: 7450 - total flux cleaned = -0.00364044 Jy ! Component: 7500 - total flux cleaned = -0.00372597 Jy ! Component: 7550 - total flux cleaned = -0.00389914 Jy ! Component: 7600 - total flux cleaned = -0.00398807 Jy ! Component: 7650 - total flux cleaned = -0.00407759 Jy ! Component: 7700 - total flux cleaned = -0.00421413 Jy ! Component: 7750 - total flux cleaned = -0.00426068 Jy ! Component: 7800 - total flux cleaned = -0.00440178 Jy ! Component: 7850 - total flux cleaned = -0.00459355 Jy ! Component: 7900 - total flux cleaned = -0.00474022 Jy ! Component: 7950 - total flux cleaned = -0.00484053 Jy ! Component: 8000 - total flux cleaned = -0.00488876 Jy ! Component: 8050 - total flux cleaned = -0.00519926 Jy ! Component: 8100 - total flux cleaned = -0.0054097 Jy ! Component: 8150 - total flux cleaned = -0.00557106 Jy ! Component: 8200 - total flux cleaned = -0.00584446 Jy ! Component: 8250 - total flux cleaned = -0.00601282 Jy ! Component: 8300 - total flux cleaned = -0.00618258 Jy ! Component: 8350 - total flux cleaned = -0.00629747 Jy ! Component: 8400 - total flux cleaned = -0.00664936 Jy ! Component: 8450 - total flux cleaned = -0.00658926 Jy ! Component: 8500 - total flux cleaned = -0.00689231 Jy ! Component: 8550 - total flux cleaned = -0.00695315 Jy ! Component: 8600 - total flux cleaned = -0.00695291 Jy ! Component: 8650 - total flux cleaned = -0.00727022 Jy ! Component: 8700 - total flux cleaned = -0.00740073 Jy ! Component: 8750 - total flux cleaned = -0.00760088 Jy ! Component: 8800 - total flux cleaned = -0.00746377 Jy ! Component: 8850 - total flux cleaned = -0.00773949 Jy ! Component: 8900 - total flux cleaned = -0.00766936 Jy ! Component: 8950 - total flux cleaned = -0.00752656 Jy ! Component: 9000 - total flux cleaned = -0.00767201 Jy ! Component: 9050 - total flux cleaned = -0.00797095 Jy ! Component: 9100 - total flux cleaned = -0.00759194 Jy ! Component: 9150 - total flux cleaned = -0.00766862 Jy ! Component: 9200 - total flux cleaned = -0.00766853 Jy ! Component: 9250 - total flux cleaned = -0.00726651 Jy ! Component: 9300 - total flux cleaned = -0.00743046 Jy ! Component: 9350 - total flux cleaned = -0.00709857 Jy ! Component: 9400 - total flux cleaned = -0.00718147 Jy ! Component: 9450 - total flux cleaned = -0.00683419 Jy ! Component: 9500 - total flux cleaned = -0.00683624 Jy ! Component: 9550 - total flux cleaned = -0.00692444 Jy ! Component: 9600 - total flux cleaned = -0.00683481 Jy ! Component: 9650 - total flux cleaned = -0.00673987 Jy ! Component: 9700 - total flux cleaned = -0.00692635 Jy ! Component: 9750 - total flux cleaned = -0.00682648 Jy ! Component: 9800 - total flux cleaned = -0.00692301 Jy ! Component: 9850 - total flux cleaned = -0.00732317 Jy ! Component: 9900 - total flux cleaned = -0.00691665 Jy ! Component: 9950 - total flux cleaned = -0.00754099 Jy ! Component: 10000 - total flux cleaned = -0.00701034 Jy ! Component: 10050 - total flux cleaned = -0.00722137 Jy ! Component: 10100 - total flux cleaned = -0.00721959 Jy ! Component: 10150 - total flux cleaned = -0.00777685 Jy ! Component: 10200 - total flux cleaned = -0.00766361 Jy ! Component: 10250 - total flux cleaned = -0.00812462 Jy ! Component: 10300 - total flux cleaned = -0.00836473 Jy ! Component: 10350 - total flux cleaned = -0.00812541 Jy ! Component: 10400 - total flux cleaned = -0.00787812 Jy ! Component: 10450 - total flux cleaned = -0.00800341 Jy ! Component: 10500 - total flux cleaned = -0.00826085 Jy ! Component: 10550 - total flux cleaned = -0.00839322 Jy ! Component: 10600 - total flux cleaned = -0.00839458 Jy ! Component: 10650 - total flux cleaned = -0.00852489 Jy ! Component: 10700 - total flux cleaned = -0.00852654 Jy ! Component: 10750 - total flux cleaned = -0.00881025 Jy ! Component: 10800 - total flux cleaned = -0.00881088 Jy ! Component: 10850 - total flux cleaned = -0.00910673 Jy ! Component: 10900 - total flux cleaned = -0.00940558 Jy ! Component: 10950 - total flux cleaned = -0.00910225 Jy ! Component: 11000 - total flux cleaned = -0.00925167 Jy ! Component: 11050 - total flux cleaned = -0.00925101 Jy ! Component: 11100 - total flux cleaned = -0.00973587 Jy ! Component: 11150 - total flux cleaned = -0.0102279 Jy ! Component: 11200 - total flux cleaned = -0.00939526 Jy ! Component: 11250 - total flux cleaned = -0.0100696 Jy ! Component: 11300 - total flux cleaned = -0.0104171 Jy ! Component: 11350 - total flux cleaned = -0.0104122 Jy ! Component: 11400 - total flux cleaned = -0.0109505 Jy ! Component: 11450 - total flux cleaned = -0.0111317 Jy ! Component: 11500 - total flux cleaned = -0.0113135 Jy ! Component: 11550 - total flux cleaned = -0.0111188 Jy ! Component: 11600 - total flux cleaned = -0.0115008 Jy ! Component: 11650 - total flux cleaned = -0.0118936 Jy ! Component: 11700 - total flux cleaned = -0.0116939 Jy ! Component: 11750 - total flux cleaned = -0.0116862 Jy ! Component: 11800 - total flux cleaned = -0.01148 Jy ! Component: 11850 - total flux cleaned = -0.0119079 Jy ! Component: 11900 - total flux cleaned = -0.0123354 Jy ! Component: 11950 - total flux cleaned = -0.0132197 Jy ! Component: 12000 - total flux cleaned = -0.0123244 Jy ! Component: 12050 - total flux cleaned = -0.0125612 Jy ! Component: 12100 - total flux cleaned = -0.013023 Jy ! Component: 12150 - total flux cleaned = -0.0139789 Jy ! Component: 12200 - total flux cleaned = -0.013737 Jy ! Component: 12250 - total flux cleaned = -0.0129997 Jy ! Component: 12300 - total flux cleaned = -0.0137552 Jy ! Component: 12350 - total flux cleaned = -0.0134957 Jy ! Component: 12400 - total flux cleaned = -0.0137508 Jy ! Component: 12450 - total flux cleaned = -0.013744 Jy ! Component: 12500 - total flux cleaned = -0.0148363 Jy ! Component: 12550 - total flux cleaned = -0.0153973 Jy ! Component: 12600 - total flux cleaned = -0.0142696 Jy ! Component: 12650 - total flux cleaned = -0.0142636 Jy ! Component: 12700 - total flux cleaned = -0.0148399 Jy ! Component: 12750 - total flux cleaned = -0.0160428 Jy ! Component: 12800 - total flux cleaned = -0.0151377 Jy ! Component: 12850 - total flux cleaned = -0.0148283 Jy ! Component: 12900 - total flux cleaned = -0.0154537 Jy ! Component: 12950 - total flux cleaned = -0.0144788 Jy ! Component: 13000 - total flux cleaned = -0.0154602 Jy ! Component: 13050 - total flux cleaned = -0.0154649 Jy ! Component: 13100 - total flux cleaned = -0.0144483 Jy ! Component: 13150 - total flux cleaned = -0.0151352 Jy ! Component: 13200 - total flux cleaned = -0.0158333 Jy ! Component: 13250 - total flux cleaned = -0.0154569 Jy ! Component: 13300 - total flux cleaned = -0.0169235 Jy ! Component: 13350 - total flux cleaned = -0.016922 Jy ! Component: 13400 - total flux cleaned = -0.0173105 Jy ! Component: 13450 - total flux cleaned = -0.0161296 Jy ! Component: 13500 - total flux cleaned = -0.0181108 Jy ! Component: 13550 - total flux cleaned = -0.0177155 Jy ! Component: 13600 - total flux cleaned = -0.0185371 Jy ! Component: 13650 - total flux cleaned = -0.0172737 Jy ! Component: 13700 - total flux cleaned = -0.0189913 Jy ! Component: 13750 - total flux cleaned = -0.0181133 Jy ! Component: 13800 - total flux cleaned = -0.0185653 Jy ! Component: 13850 - total flux cleaned = -0.0176502 Jy ! Component: 13900 - total flux cleaned = -0.0171903 Jy ! Component: 13950 - total flux cleaned = -0.0186029 Jy ! Component: 14000 - total flux cleaned = -0.0186041 Jy ! Component: 14050 - total flux cleaned = -0.0190999 Jy ! Component: 14100 - total flux cleaned = -0.0185887 Jy ! Component: 14150 - total flux cleaned = -0.0180812 Jy ! Component: 14200 - total flux cleaned = -0.0191246 Jy ! Component: 14250 - total flux cleaned = -0.017022 Jy ! Component: 14300 - total flux cleaned = -0.0191654 Jy ! Component: 14350 - total flux cleaned = -0.0169568 Jy ! Component: 14400 - total flux cleaned = -0.0186289 Jy ! Component: 14450 - total flux cleaned = -0.019774 Jy ! Component: 14500 - total flux cleaned = -0.0209389 Jy ! Component: 14550 - total flux cleaned = -0.0203476 Jy ! Component: 14600 - total flux cleaned = -0.0197396 Jy ! Component: 14650 - total flux cleaned = -0.0203638 Jy ! Component: 14700 - total flux cleaned = -0.0197397 Jy ! Component: 14750 - total flux cleaned = -0.0190965 Jy ! Component: 14800 - total flux cleaned = -0.017784 Jy ! Component: 14850 - total flux cleaned = -0.0191167 Jy ! Component: 14900 - total flux cleaned = -0.0198127 Jy ! Component: 14950 - total flux cleaned = -0.01707 Jy ! Component: 15000 - total flux cleaned = -0.017772 Jy ! Component: 15050 - total flux cleaned = -0.0156064 Jy ! Component: 15100 - total flux cleaned = -0.0163166 Jy ! Component: 15150 - total flux cleaned = -0.0185462 Jy ! Component: 15200 - total flux cleaned = -0.017008 Jy ! Component: 15250 - total flux cleaned = -0.0169932 Jy ! Component: 15300 - total flux cleaned = -0.0162226 Jy ! Component: 15350 - total flux cleaned = -0.0162036 Jy ! Component: 15400 - total flux cleaned = -0.0161981 Jy ! Component: 15450 - total flux cleaned = -0.0170296 Jy ! Component: 15500 - total flux cleaned = -0.0170186 Jy ! Component: 15550 - total flux cleaned = -0.0170121 Jy ! Component: 15600 - total flux cleaned = -0.0161268 Jy ! Component: 15650 - total flux cleaned = -0.0179206 Jy ! Component: 15700 - total flux cleaned = -0.0142809 Jy ! Component: 15750 - total flux cleaned = -0.0142966 Jy ! Component: 15800 - total flux cleaned = -0.0161822 Jy ! Component: 15850 - total flux cleaned = -0.0171572 Jy ! Component: 15900 - total flux cleaned = -0.0112263 Jy ! Component: 15950 - total flux cleaned = -0.0132338 Jy ! Component: 16000 - total flux cleaned = -0.00911416 Jy ! Component: 16050 - total flux cleaned = -0.012236 Jy ! Component: 16100 - total flux cleaned = -0.0111737 Jy ! Component: 16150 - total flux cleaned = -0.0122569 Jy ! Component: 16200 - total flux cleaned = -0.012252 Jy ! Component: 16250 - total flux cleaned = -0.0133719 Jy ! Component: 16300 - total flux cleaned = -0.0099057 Jy ! Component: 16350 - total flux cleaned = -0.00987811 Jy ! Total flux subtracted in 16384 components = -0.0110512 Jy ! Clean residual min=-0.064730 max=0.063920 Jy/beam ! Clean residual mean=0.000000 rms=0.011381 Jy/beam ! Combined flux in latest and established models = 0.0286528 Jy keep ! Adding 474 model components to the UV plane model. ! The established model now contains 655 components and 0.0286528 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/g0/scr/pima//bp171af_uvs/J2037-0729_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9025 x 1.766 at 2.37 degrees (North through East) ! Clean map min=-0.068419 max=0.065835 Jy/beam ! Writing clean map to FITS file: /g0/scr/pima//bp171af_uvs/J2037-0729_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.0658351 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00954661 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -43.5024 43.5024 87.0048 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0286 print "peak= ", int(1000*peak_flux)/1000 ! peak= 0.065 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 6.896 quit ! Quitting program ! Log file /g0/scr/pima//bp171af_uvs/J2037-0729_X_uvs.log closed on Sun Dec 8 12:05:07 2013