! Started logfile: /g0/scr/pima/npcsb_uvs/J1904+7648_X_uvs.log on Sat Oct 5 18:54:45 2013 obs /g0/scr/pima/npcsb_uvs/J1904+7648_X_uva.fits ! Reading UV FITS file: /g0/scr/pima/npcsb_uvs/J1904+7648_X_uva.fits ! AN table 1: 47 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.614995 visibilities/baseline/integration-bin. ! Found source: J1904+7648 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2428 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/J1904+7648_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 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 48 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=1.001717Jy sigma=0.708179 ! Fit after self-cal, rms=0.954538Jy sigma=0.648834 ! 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=1.905 mas, bmaj=5.948 mas, bpa=45.96 degrees ! Estimated noise=50.319 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.0557434 Jy ! Component: 100 - total flux cleaned = 0.0655354 Jy ! Total flux subtracted in 100 components = 0.0655354 Jy ! Clean residual min=-0.008110 max=0.008112 Jy/beam ! Clean residual mean=0.000011 rms=0.002816 Jy/beam ! Combined flux in latest and established models = 0.0655354 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 18 components and 0.0655354 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068933Jy sigma=0.048096 ! Fit after self-cal, rms=0.067165Jy sigma=0.047101 ! 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.414 mas, bmaj=5.629 mas, bpa=42.56 degrees ! Estimated noise=23.7587 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.00137852 Jy ! Component: 100 - total flux cleaned = 0.00274239 Jy ! Component: 150 - total flux cleaned = 0.00328927 Jy ! Component: 200 - total flux cleaned = 0.00376708 Jy ! Total flux subtracted in 200 components = 0.00376708 Jy ! Clean residual min=-0.001890 max=0.001905 Jy/beam ! Clean residual mean=0.000005 rms=0.000684 Jy/beam ! Combined flux in latest and established models = 0.0693025 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 64 components and 0.0693025 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.067798Jy sigma=0.045499 ! Fit after self-cal, rms=0.067665Jy sigma=0.045424 ! 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.85 FD 0.87 KP 0.91 LA 0.93 ! OV 0.87 PT 1.08 SC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.82 FD 0.87 KP 0.98 LA 0.88 ! OV 0.90 PT 0.93 SC 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.09 KP 0.76 LA 0.90 ! OV 0.84 PT 0.88 SC 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.84 FD 0.81 KP 0.85 LA 0.90 ! OV 0.92 PT 0.82 SC 1.00* ! ! ! Fit before self-cal, rms=0.067665Jy sigma=0.045424 ! Fit after self-cal, rms=0.058772Jy sigma=0.042755 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.393 mas, bmaj=5.608 mas, bpa=43.08 degrees ! Estimated noise=18.6535 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.0105137 Jy ! Component: 100 - total flux cleaned = -0.0108393 Jy ! Component: 150 - total flux cleaned = -0.0106133 Jy ! Component: 200 - total flux cleaned = -0.0106728 Jy ! Total flux subtracted in 200 components = -0.0106728 Jy ! Clean residual min=-0.001002 max=0.001013 Jy/beam ! Clean residual mean=-0.000001 rms=0.000337 Jy/beam ! Combined flux in latest and established models = 0.0586297 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 96 components and 0.0586297 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.053638Jy sigma=0.044132 ! Fit after self-cal, rms=0.053631Jy sigma=0.044065 ! 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.053631Jy sigma=0.044065 ! Fit after self-cal, rms=0.047834Jy sigma=0.042509 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.387 mas, bmaj=5.603 mas, bpa=43.04 degrees ! Estimated noise=15.1911 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00850394 Jy ! Component: 100 - total flux cleaned = -0.0103032 Jy ! Component: 150 - total flux cleaned = -0.0102144 Jy ! Component: 200 - total flux cleaned = -0.00995104 Jy ! Total flux subtracted in 200 components = -0.00995104 Jy ! Clean residual min=-0.000670 max=0.000674 Jy/beam ! Clean residual mean=-0.000002 rms=0.000231 Jy/beam ! Combined flux in latest and established models = 0.0486787 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 140 components and 0.0486787 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.043703Jy sigma=0.043966 ! Fit after self-cal, rms=0.043680Jy sigma=0.043938 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.000916441 Jy ! Component: 100 - total flux cleaned = 0.00139298 Jy ! Component: 150 - total flux cleaned = 0.00158824 Jy ! Component: 200 - total flux cleaned = 0.00187665 Jy ! Total flux subtracted in 200 components = 0.00187665 Jy ! Clean residual min=-0.000584 max=0.000585 Jy/beam ! Clean residual mean=-0.000001 rms=0.000193 Jy/beam ! Combined flux in latest and established models = 0.0505553 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 184 components and 0.0505553 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.043546Jy sigma=0.043941 ! Fit after self-cal, rms=0.043533Jy sigma=0.043924 ! 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.043533Jy sigma=0.043924 ! Fit after self-cal, rms=0.038802Jy sigma=0.042273 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.381 mas, bmaj=5.588 mas, bpa=42.96 degrees ! Estimated noise=12.4227 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00645487 Jy ! Component: 100 - total flux cleaned = -0.00787429 Jy ! Component: 150 - total flux cleaned = -0.00807074 Jy ! Component: 200 - total flux cleaned = -0.00809786 Jy ! Total flux subtracted in 200 components = -0.00809786 Jy ! Clean residual min=-0.000462 max=0.000458 Jy/beam ! Clean residual mean=-0.000002 rms=0.000160 Jy/beam ! Combined flux in latest and established models = 0.0424575 Jy selfcal ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 230 components and 0.0424575 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.035731Jy sigma=0.043743 ! Fit after self-cal, rms=0.035718Jy sigma=0.043709 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.0005162 Jy ! Component: 100 - total flux cleaned = 0.000636635 Jy ! Component: 150 - total flux cleaned = 0.000582857 Jy ! Component: 200 - total flux cleaned = 0.000506711 Jy ! Total flux subtracted in 200 components = 0.000506711 Jy ! Clean residual min=-0.000403 max=0.000401 Jy/beam ! Clean residual mean=-0.000001 rms=0.000136 Jy/beam ! Combined flux in latest and established models = 0.0429642 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 269 components and 0.0429642 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.035619Jy sigma=0.043745 ! Fit after self-cal, rms=0.035579Jy sigma=0.043727 ! Inverting map ! Added new window around map position (35.4, -44.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000385497 Jy ! Component: 100 - total flux cleaned = 0.000186884 Jy ! Component: 150 - total flux cleaned = 4.74745e-05 Jy ! Component: 200 - total flux cleaned = -1.88313e-05 Jy ! Total flux subtracted in 200 components = -1.88313e-05 Jy ! Clean residual min=-0.000597 max=0.000630 Jy/beam ! Clean residual mean=-0.000001 rms=0.000175 Jy/beam ! Combined flux in latest and established models = 0.0429453 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 303 components and 0.0429453 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.035559Jy sigma=0.043720 ! Fit after self-cal, rms=0.035562Jy sigma=0.043713 ! 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.035562Jy sigma=0.043713 ! Fit after self-cal, rms=0.031681Jy sigma=0.042038 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=5.638 mas, bpa=40.62 degrees ! Estimated noise=10.3946 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00236264 Jy ! Component: 100 - total flux cleaned = -0.0028283 Jy ! Component: 150 - total flux cleaned = -0.00300225 Jy ! Component: 200 - total flux cleaned = -0.0030974 Jy ! Total flux subtracted in 200 components = -0.0030974 Jy ! Clean residual min=-0.004803 max=0.002774 Jy/beam ! Clean residual mean=-0.000009 rms=0.000922 Jy/beam ! Combined flux in latest and established models = 0.0398479 Jy selfcal ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 344 components and 0.0398479 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.030992Jy sigma=0.042360 ! Fit after self-cal, rms=0.031026Jy sigma=0.042204 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.00145174 Jy ! Component: 100 - total flux cleaned = -0.00163337 Jy ! Component: 150 - total flux cleaned = -0.00173941 Jy ! Component: 200 - total flux cleaned = -0.00194761 Jy ! Total flux subtracted in 200 components = -0.00194761 Jy ! Clean residual min=-0.003929 max=0.002129 Jy/beam ! Clean residual mean=-0.000007 rms=0.000734 Jy/beam ! Combined flux in latest and established models = 0.0379003 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 361 components and 0.0379003 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.030719Jy sigma=0.042444 ! Fit after self-cal, rms=0.030723Jy sigma=0.042399 ! 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.376 mas, bmaj=5.57 mas, bpa=42.91 degrees ! Estimated noise=10.1276 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.376 x 5.57 at 42.91 degrees (North through East) ! Clean map min=-0.0043778 max=0.029915 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=1.882 mas, bmaj=5.951 mas, bpa=46.04 degrees ! Estimated noise=20.9004 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00698665 Jy ! Component: 100 - total flux cleaned = 0.00926609 Jy ! Component: 150 - total flux cleaned = 0.010399 Jy ! Component: 200 - total flux cleaned = 0.0115094 Jy ! Total flux subtracted in 200 components = 0.0115094 Jy ! Clean residual min=-0.011981 max=0.022207 Jy/beam ! Clean residual mean=0.000030 rms=0.004172 Jy/beam ! Combined flux in latest and established models = 0.0115093 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 40 model components to the UV plane model. ! The established model now contains 40 components and 0.0115093 Jy ! Inverting map and beam ! Estimated beam: bmin=2.376 mas, bmaj=5.57 mas, bpa=42.91 degrees ! Estimated noise=10.1276 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0175431 Jy ! Component: 100 - total flux cleaned = 0.0170737 Jy ! Component: 150 - total flux cleaned = 0.017057 Jy ! Component: 200 - total flux cleaned = 0.0178677 Jy ! Total flux subtracted in 200 components = 0.0178677 Jy ! Clean residual min=-0.006731 max=0.004874 Jy/beam ! Clean residual mean=0.000012 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 0.029377 Jy ! Adding 18 model components to the UV plane model. ! The established model now contains 57 components and 0.029377 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.684 mas, bmaj=5.638 mas, bpa=40.62 degrees ! Estimated noise=10.3946 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000552968 Jy ! Component: 100 - total flux cleaned = 0.000791337 Jy ! Component: 150 - total flux cleaned = 0.00143286 Jy ! Component: 200 - total flux cleaned = 0.00211441 Jy ! Total flux subtracted in 200 components = 0.00211441 Jy ! Clean residual min=-0.005651 max=0.004097 Jy/beam ! Clean residual mean=0.000011 rms=0.001245 Jy/beam ! Combined flux in latest and established models = 0.0314914 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.000442315 Jy ! Component: 100 - total flux cleaned = 0.000847692 Jy ! Component: 150 - total flux cleaned = 0.0010737 Jy ! Component: 200 - total flux cleaned = 0.00163308 Jy ! Total flux subtracted in 200 components = 0.00163308 Jy ! Clean residual min=-0.004682 max=0.003324 Jy/beam ! Clean residual mean=0.000011 rms=0.001016 Jy/beam ! Combined flux in latest and established models = 0.0331245 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 103 components and 0.0331245 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.376 mas, bmaj=5.57 mas, bpa=42.91 degrees ! Estimated noise=10.1276 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000248652 Jy ! Component: 100 - total flux cleaned = 0.000489406 Jy ! Component: 150 - total flux cleaned = 0.00088901 Jy ! Component: 200 - total flux cleaned = 0.00132118 Jy ! Total flux subtracted in 200 components = 0.00132118 Jy ! Clean residual min=-0.004339 max=0.003122 Jy/beam ! Clean residual mean=0.000008 rms=0.000901 Jy/beam ! Combined flux in latest and established models = 0.0344457 Jy selfcal ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 133 components and 0.0344457 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.029725Jy sigma=0.046936 ! Fit after self-cal, rms=0.029282Jy sigma=0.045569 wmodel /g0/scr/pima/npcsb_uvs/J1904+7648_X_map.mod ! Writing 133 model components to file: /g0/scr/pima/npcsb_uvs/J1904+7648_X_map.mod wobs /g0/scr/pima/npcsb_uvs/J1904+7648_X_uvs.fits ! Writing UV FITS file: /g0/scr/pima/npcsb_uvs/J1904+7648_X_uvs.fits wwins /g0/scr/pima/npcsb_uvs/J1904+7648_X_map.win ! wwins: Wrote 2 windows to /g0/scr/pima/npcsb_uvs/J1904+7648_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00063731 Jy ! Component: 100 - total flux cleaned = -0.000665248 Jy ! Component: 150 - total flux cleaned = -0.000641674 Jy ! Component: 200 - total flux cleaned = -0.000686763 Jy ! Component: 250 - total flux cleaned = -0.000730577 Jy ! Component: 300 - total flux cleaned = -0.000580764 Jy ! Component: 350 - total flux cleaned = -0.000280219 Jy ! Component: 400 - total flux cleaned = -1.27239e-05 Jy ! Component: 450 - total flux cleaned = 0.000258265 Jy ! Component: 500 - total flux cleaned = 0.000531827 Jy ! Component: 550 - total flux cleaned = 0.000778329 Jy ! Component: 600 - total flux cleaned = 0.00101394 Jy ! Component: 650 - total flux cleaned = 0.00125417 Jy ! Component: 700 - total flux cleaned = 0.00152525 Jy ! Component: 750 - total flux cleaned = 0.00175062 Jy ! Component: 800 - total flux cleaned = 0.00195738 Jy ! Component: 850 - total flux cleaned = 0.00213476 Jy ! Component: 900 - total flux cleaned = 0.00234218 Jy ! Component: 950 - total flux cleaned = 0.00254479 Jy ! Component: 1000 - total flux cleaned = 0.00261101 Jy ! Component: 1050 - total flux cleaned = 0.00268748 Jy ! Component: 1100 - total flux cleaned = 0.00276291 Jy ! Component: 1150 - total flux cleaned = 0.00283733 Jy ! Component: 1200 - total flux cleaned = 0.00288984 Jy ! Component: 1250 - total flux cleaned = 0.00293143 Jy ! Component: 1300 - total flux cleaned = 0.00298284 Jy ! Component: 1350 - total flux cleaned = 0.00304383 Jy ! Component: 1400 - total flux cleaned = 0.00312433 Jy ! Component: 1450 - total flux cleaned = 0.00319415 Jy ! Component: 1500 - total flux cleaned = 0.00325346 Jy ! Component: 1550 - total flux cleaned = 0.00330258 Jy ! Component: 1600 - total flux cleaned = 0.00336114 Jy ! Component: 1650 - total flux cleaned = 0.00341935 Jy ! Component: 1700 - total flux cleaned = 0.00346758 Jy ! Component: 1750 - total flux cleaned = 0.00352523 Jy ! Component: 1800 - total flux cleaned = 0.00357307 Jy ! Component: 1850 - total flux cleaned = 0.00361123 Jy ! Component: 1900 - total flux cleaned = 0.0036875 Jy ! Component: 1950 - total flux cleaned = 0.00374458 Jy ! Component: 2000 - total flux cleaned = 0.0038016 Jy ! Component: 2050 - total flux cleaned = 0.00387745 Jy ! Component: 2100 - total flux cleaned = 0.00395331 Jy ! Component: 2150 - total flux cleaned = 0.00400078 Jy ! Component: 2200 - total flux cleaned = 0.00407672 Jy ! Component: 2250 - total flux cleaned = 0.00411467 Jy ! Component: 2300 - total flux cleaned = 0.0041622 Jy ! Component: 2350 - total flux cleaned = 0.00420028 Jy ! Component: 2400 - total flux cleaned = 0.00424803 Jy ! Component: 2450 - total flux cleaned = 0.00427673 Jy ! Component: 2500 - total flux cleaned = 0.00431509 Jy ! Component: 2550 - total flux cleaned = 0.00436311 Jy ! Component: 2600 - total flux cleaned = 0.00440162 Jy ! Component: 2650 - total flux cleaned = 0.00444022 Jy ! Component: 2700 - total flux cleaned = 0.00444993 Jy ! Component: 2750 - total flux cleaned = 0.00448891 Jy ! Component: 2800 - total flux cleaned = 0.00450857 Jy ! Component: 2850 - total flux cleaned = 0.00455773 Jy ! Component: 2900 - total flux cleaned = 0.00458739 Jy ! Component: 2950 - total flux cleaned = 0.00462731 Jy ! Component: 3000 - total flux cleaned = 0.00466748 Jy ! Component: 3050 - total flux cleaned = 0.00469776 Jy ! Component: 3100 - total flux cleaned = 0.00473846 Jy ! Component: 3150 - total flux cleaned = 0.00477947 Jy ! Component: 3200 - total flux cleaned = 0.00480011 Jy ! Component: 3250 - total flux cleaned = 0.00485238 Jy ! Component: 3300 - total flux cleaned = 0.00486296 Jy ! Component: 3350 - total flux cleaned = 0.00489487 Jy ! Component: 3400 - total flux cleaned = 0.00491635 Jy ! Component: 3450 - total flux cleaned = 0.00494883 Jy ! Component: 3500 - total flux cleaned = 0.00497064 Jy ! Component: 3550 - total flux cleaned = 0.00500368 Jy ! Component: 3600 - total flux cleaned = 0.00504802 Jy ! Component: 3650 - total flux cleaned = 0.00507031 Jy ! Component: 3700 - total flux cleaned = 0.00510419 Jy ! Component: 3750 - total flux cleaned = 0.0051384 Jy ! Component: 3800 - total flux cleaned = 0.00517285 Jy ! Component: 3850 - total flux cleaned = 0.00520771 Jy ! Component: 3900 - total flux cleaned = 0.00524287 Jy ! Component: 3950 - total flux cleaned = 0.00526654 Jy ! Component: 4000 - total flux cleaned = 0.00530252 Jy ! Component: 4050 - total flux cleaned = 0.00535094 Jy ! Total flux subtracted in 4096 components = 0.00537548 Jy ! Clean residual min=-0.000718 max=0.000652 Jy/beam ! Clean residual mean=0.000005 rms=0.000224 Jy/beam ! Combined flux in latest and established models = 0.0398212 Jy keep ! Adding 213 model components to the UV plane model. ! The established model now contains 346 components and 0.0398212 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/g0/scr/pima/npcsb_uvs/J1904+7648_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.376 x 5.57 at 42.91 degrees (North through East) ! Clean map min=-0.0043228 max=0.026934 Jy/beam ! Writing clean map to FITS file: /g0/scr/pima/npcsb_uvs/J1904+7648_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.0269339 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.000258558 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -2.87992 2.87992 5.75984 11.5197 23.0394 46.0787 92.1575 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0007 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.026 quit ! Quitting program ! Log file /g0/scr/pima/npcsb_uvs/J1904+7648_X_uvs.log closed on Sat Oct 5 18:54:49 2013