! Started logfile: /g0/scr/pima/npcsb_uvs/J2325+7917_X_uvs.log on Wed Sep 25 00:08:25 2013 obs /g0/scr/pima/npcsb_uvs/J2325+7917_X_uva.fits ! Reading UV FITS file: /g0/scr/pima/npcsb_uvs/J2325+7917_X_uva.fits ! AN table 1: 46 integrations on 15 of 15 possible baselines. ! AN table 2: 46 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.741848 visibilities/baseline/integration-bin. ! Found source: J2325+7917 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40993e+09 8e+06 1 8e+06 ! 02 2 8.47993e+09 8e+06 1 8e+06 ! 03 3 8.79493e+09 8e+06 1 8e+06 ! 04 4 8.89993e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2184 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/01_muppet /g0/scr/pima/npcsb_uvs/J2325+7917_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 14 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 14 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 14 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 14 telescope corrections were flagged in sub-array 1. ! A total of 92 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=1.038518Jy sigma=2.514624 ! Fit after self-cal, rms=0.892073Jy sigma=2.212474 ! 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.392 mas, bmaj=7.333 mas, bpa=39.47 degrees ! Estimated noise=11.357 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.0816325 Jy ! Component: 100 - total flux cleaned = 0.114893 Jy ! Total flux subtracted in 100 components = 0.114893 Jy ! Clean residual min=-0.013992 max=0.016503 Jy/beam ! Clean residual mean=0.000041 rms=0.005721 Jy/beam ! Combined flux in latest and established models = 0.114893 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 38 components and 0.114893 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084743Jy sigma=0.190876 ! Fit after self-cal, rms=0.084168Jy sigma=0.189418 ! Inverting map print "************** FINISHED UNIFORM WEIGHTING CLEAN **************" ! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** ! uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.77 mas, bmaj=7.497 mas, bpa=34.16 degrees ! Estimated noise=6.4729 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.0232537 Jy ! Component: 100 - total flux cleaned = 0.0317482 Jy ! Component: 150 - total flux cleaned = 0.0361276 Jy ! Component: 200 - total flux cleaned = 0.0370799 Jy ! Total flux subtracted in 200 components = 0.0370799 Jy ! Clean residual min=-0.004751 max=0.004742 Jy/beam ! Clean residual mean=0.000005 rms=0.001986 Jy/beam ! Combined flux in latest and established models = 0.151973 Jy ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 110 components and 0.151973 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.076421Jy sigma=0.169950 ! Fit after self-cal, rms=0.075850Jy sigma=0.168728 ! 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 1.06 FD 0.92 KP 1.02 LA 0.97 ! OV 1.00 PT 0.97 ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* SC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.89 KP 0.95 LA 0.96 ! OV 1.00 PT 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* SC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.03 KP 0.89 LA 0.96 ! OV 0.94 PT 1.01 ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* SC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 1.00 KP 0.94 LA 0.94 ! OV 1.00 PT 1.11 ! ! Telescope amplitude corrections in sub-array 2: ! HN 0.00* SC 0.00* ! ! ! Fit before self-cal, rms=0.075850Jy sigma=0.168728 ! Fit after self-cal, rms=0.071342Jy sigma=0.166729 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.789 mas, bmaj=7.51 mas, bpa=33.69 degrees ! Estimated noise=6.21552 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.00604231 Jy ! Component: 100 - total flux cleaned = 0.00820457 Jy ! Component: 150 - total flux cleaned = 0.00948379 Jy ! Component: 200 - total flux cleaned = 0.0105428 Jy ! Total flux subtracted in 200 components = 0.0105428 Jy ! Clean residual min=-0.001461 max=0.001467 Jy/beam ! Clean residual mean=0.000003 rms=0.000544 Jy/beam ! Combined flux in latest and established models = 0.162516 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 176 components and 0.162516 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070807Jy sigma=0.165027 ! Fit after self-cal, rms=0.070848Jy sigma=0.164796 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070848Jy sigma=0.164796 ! Fit after self-cal, rms=0.070152Jy sigma=0.164749 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.78 mas, bmaj=7.499 mas, bpa=33.68 degrees ! Estimated noise=6.17119 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00173815 Jy ! Component: 100 - total flux cleaned = 0.00173453 Jy ! Component: 150 - total flux cleaned = 0.00195184 Jy ! Component: 200 - total flux cleaned = 0.00175532 Jy ! Total flux subtracted in 200 components = 0.00175532 Jy ! Clean residual min=-0.001067 max=0.001069 Jy/beam ! Clean residual mean=0.000002 rms=0.000397 Jy/beam ! Combined flux in latest and established models = 0.164271 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 232 components and 0.164271 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070000Jy sigma=0.164313 ! Fit after self-cal, rms=0.069901Jy sigma=0.164277 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 = 2.55431e-05 Jy ! Component: 100 - total flux cleaned = -3.96719e-05 Jy ! Component: 150 - total flux cleaned = -0.000513796 Jy ! Component: 200 - total flux cleaned = -0.000916049 Jy ! Total flux subtracted in 200 components = -0.000916049 Jy ! Clean residual min=-0.000942 max=0.000937 Jy/beam ! Clean residual mean=0.000001 rms=0.000336 Jy/beam ! Combined flux in latest and established models = 0.163355 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 281 components and 0.163355 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069802Jy sigma=0.164027 ! Fit after self-cal, rms=0.069736Jy sigma=0.164015 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069736Jy sigma=0.164015 ! Fit after self-cal, rms=0.069606Jy sigma=0.163966 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.773 mas, bmaj=7.497 mas, bpa=33.7 degrees ! Estimated noise=6.15179 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000847236 Jy ! Component: 100 - total flux cleaned = 0.00053227 Jy ! Component: 150 - total flux cleaned = 0.000140267 Jy ! Component: 200 - total flux cleaned = -5.18876e-05 Jy ! Total flux subtracted in 200 components = -5.18876e-05 Jy ! Clean residual min=-0.000780 max=0.000786 Jy/beam ! Clean residual mean=0.000000 rms=0.000283 Jy/beam ! Combined flux in latest and established models = 0.163303 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 321 components and 0.163303 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069515Jy sigma=0.163807 ! Fit after self-cal, rms=0.069498Jy sigma=0.163798 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.000165459 Jy ! Component: 100 - total flux cleaned = -9.44351e-06 Jy ! Component: 150 - total flux cleaned = -9.36885e-05 Jy ! Component: 200 - total flux cleaned = -0.000135808 Jy ! Total flux subtracted in 200 components = -0.000135808 Jy ! Clean residual min=-0.000699 max=0.000702 Jy/beam ! Clean residual mean=0.000000 rms=0.000248 Jy/beam ! Combined flux in latest and established models = 0.163168 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 361 components and 0.163167 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069439Jy sigma=0.163681 ! Fit after self-cal, rms=0.069445Jy sigma=0.163674 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069445Jy sigma=0.163674 ! Fit after self-cal, rms=0.069286Jy sigma=0.163635 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.026 mas, bmaj=7.815 mas, bpa=31.36 degrees ! Estimated noise=6.20198 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000545136 Jy ! Component: 100 - total flux cleaned = 0.000353825 Jy ! Component: 150 - total flux cleaned = 0.000132597 Jy ! Component: 200 - total flux cleaned = -1.36046e-05 Jy ! Total flux subtracted in 200 components = -1.36046e-05 Jy ! Clean residual min=-0.000609 max=0.000600 Jy/beam ! Clean residual mean=0.000000 rms=0.000221 Jy/beam ! Combined flux in latest and established models = 0.163154 Jy selfcal ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 400 components and 0.163154 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069226Jy sigma=0.163547 ! Fit after self-cal, rms=0.069205Jy sigma=0.163541 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.000146889 Jy ! Component: 100 - total flux cleaned = 0.000116029 Jy ! Component: 150 - total flux cleaned = 0.000183156 Jy ! Component: 200 - total flux cleaned = 0.000216671 Jy ! Total flux subtracted in 200 components = 0.000216671 Jy ! Clean residual min=-0.000551 max=0.000550 Jy/beam ! Clean residual mean=0.000000 rms=0.000203 Jy/beam ! Combined flux in latest and established models = 0.163371 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 429 components and 0.163371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069166Jy sigma=0.163470 ! Fit after self-cal, rms=0.069170Jy sigma=0.163465 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=2.769 mas, bmaj=7.497 mas, bpa=33.65 degrees ! Estimated noise=6.13069 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.769 x 7.497 at 33.65 degrees (North through East) ! Clean map min=-0.0029784 max=0.096546 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=2.403 mas, bmaj=7.323 mas, bpa=38.64 degrees ! Estimated noise=10.7692 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0774654 Jy ! Component: 100 - total flux cleaned = 0.110312 Jy ! Component: 150 - total flux cleaned = 0.129662 Jy ! Component: 200 - total flux cleaned = 0.13815 Jy ! Total flux subtracted in 200 components = 0.13815 Jy ! Clean residual min=-0.007415 max=0.007586 Jy/beam ! Clean residual mean=0.000014 rms=0.003050 Jy/beam ! Combined flux in latest and established models = 0.13815 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 69 model components to the UV plane model. ! The established model now contains 69 components and 0.13815 Jy ! Inverting map and beam ! Estimated beam: bmin=2.769 mas, bmaj=7.497 mas, bpa=33.65 degrees ! Estimated noise=6.13069 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.00894838 Jy ! Component: 100 - total flux cleaned = 0.0130371 Jy ! Component: 150 - total flux cleaned = 0.0132896 Jy ! Component: 200 - total flux cleaned = 0.0136614 Jy ! Total flux subtracted in 200 components = 0.0136614 Jy ! Clean residual min=-0.002881 max=0.002882 Jy/beam ! Clean residual mean=0.000001 rms=0.001179 Jy/beam ! Combined flux in latest and established models = 0.151811 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 129 components and 0.151811 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=3.026 mas, bmaj=7.815 mas, bpa=31.36 degrees ! Estimated noise=6.20198 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000365184 Jy ! Component: 100 - total flux cleaned = 0.000628984 Jy ! Component: 150 - total flux cleaned = 0.000180031 Jy ! Component: 200 - total flux cleaned = -2.50878e-05 Jy ! Total flux subtracted in 200 components = -2.50878e-05 Jy ! Clean residual min=-0.001624 max=0.001625 Jy/beam ! Clean residual mean=0.000002 rms=0.000622 Jy/beam ! Combined flux in latest and established models = 0.151786 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.000380164 Jy ! Component: 100 - total flux cleaned = -0.000471781 Jy ! Component: 150 - total flux cleaned = -0.000470791 Jy ! Component: 200 - total flux cleaned = -0.000665309 Jy ! Total flux subtracted in 200 components = -0.000665309 Jy ! Clean residual min=-0.001602 max=0.001610 Jy/beam ! Clean residual mean=0.000001 rms=0.000562 Jy/beam ! Combined flux in latest and established models = 0.151121 Jy ! Adding 108 model components to the UV plane model. ! The established model now contains 231 components and 0.151121 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.769 mas, bmaj=7.497 mas, bpa=33.65 degrees ! Estimated noise=6.13069 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000267719 Jy ! Component: 100 - total flux cleaned = -0.000472908 Jy ! Component: 150 - total flux cleaned = -0.000869476 Jy ! Component: 200 - total flux cleaned = -0.00117983 Jy ! Total flux subtracted in 200 components = -0.00117983 Jy ! Clean residual min=-0.000841 max=0.000834 Jy/beam ! Clean residual mean=-0.000001 rms=0.000295 Jy/beam ! Combined flux in latest and established models = 0.149941 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 283 components and 0.149941 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069187Jy sigma=0.163532 ! Fit after self-cal, rms=0.069127Jy sigma=0.163487 wmodel /g0/scr/pima/npcsb_uvs/J2325+7917_X_map.mod ! Writing 283 model components to file: /g0/scr/pima/npcsb_uvs/J2325+7917_X_map.mod wobs /g0/scr/pima/npcsb_uvs/J2325+7917_X_uvs.fits ! Writing UV FITS file: /g0/scr/pima/npcsb_uvs/J2325+7917_X_uvs.fits wwins /g0/scr/pima/npcsb_uvs/J2325+7917_X_map.win ! wwins: No CLEAN windows to write. ! Error occured in command: wwins ![Exited script file: /home/lpetrov/bin/01_muppet] ! wmap "/g0/scr/pima/npcsb_uvs/J2325+7917_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.769 x 7.497 at 33.65 degrees (North through East) ! Clean map min=-0.0041073 max=0.089318 Jy/beam ! Writing clean map to FITS file: /g0/scr/pima/npcsb_uvs/J2325+7917_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.0893181 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000292478 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.982369 0.982369 1.96474 3.92948 7.85896 15.7179 31.4358 62.8716 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0008 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.089 quit ! Quitting program ! Log file /g0/scr/pima/npcsb_uvs/J2325+7917_X_uvs.log closed on Wed Sep 25 00:08:28 2013