! Started logfile: /d2/scr/pima/bp177i_uvs/J0854+2006_C_uvs.log on Tue Apr 15 00:20:21 2014 obs /d2/scr/pima/bp177i_uvs/J0854+2006_C_uva.fits ! Reading UV FITS file: /d2/scr/pima/bp177i_uvs/J0854+2006_C_uva.fits ! AN table 1: 4 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.755556 visibilities/baseline/integration-bin. ! Found source: J0854+2006 ! ! 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 4.128e+09 3.2e+07 1 3.2e+07 ! 02 2 4.16e+09 3.2e+07 1 3.2e+07 ! 03 3 4.192e+09 3.2e+07 1 3.2e+07 ! 04 4 4.224e+09 3.2e+07 1 3.2e+07 ! 05 5 4.416e+09 3.2e+07 1 3.2e+07 ! 06 6 4.512e+09 3.2e+07 1 3.2e+07 ! 07 7 4.544e+09 3.2e+07 1 3.2e+07 ! 08 8 4.576e+09 3.2e+07 1 3.2e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 1088 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 = 1024 print "field_cell =",field_cell ! field_cell = 0.3 print "taper_size =",taper_size ! taper_size = 40 ![@/home/lpetrov/bin/01_muppet /d2/scr/pima/bp177i_uvs/J0854+2006_C] float lim_self lim_self = 120.0 float lim_val lim_va = 1.2 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.300x0.300 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 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 5. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 6. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 7. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 8. ! A total of 10 telescope corrections were flagged in sub-array 1. ! ! Fit before self-cal, rms=2.440628Jy sigma=114.148508 ! Fit after self-cal, rms=1.335518Jy sigma=62.327087 ! 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.334 mas, bmaj=3.813 mas, bpa=1.423 degrees ! Estimated noise=0.647106 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 = 1.52731 Jy ! Component: 100 - total flux cleaned = 2.0572 Jy ! Total flux subtracted in 100 components = 2.0572 Jy ! Clean residual min=-0.088767 max=0.242102 Jy/beam ! Clean residual mean=-0.000096 rms=0.030048 Jy/beam ! Combined flux in latest and established models = 2.0572 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 2.0572 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.658711Jy sigma=26.793708 ! Fit after self-cal, rms=0.655782Jy sigma=26.745120 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.267845 Jy ! Component: 100 - total flux cleaned = 0.411781 Jy ! Total flux subtracted in 100 components = 0.411781 Jy ! Clean residual min=-0.070336 max=0.109166 Jy/beam ! Clean residual mean=-0.000065 rms=0.017644 Jy/beam ! Combined flux in latest and established models = 2.46898 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 17 components and 2.46898 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.520926Jy sigma=17.687711 ! Fit after self-cal, rms=0.519904Jy sigma=17.669426 ! 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=1.593 mas, bmaj=4.031 mas, bpa=2.846 degrees ! Estimated noise=0.509053 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 (-3.6, -2.1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.143046 Jy ! Component: 100 - total flux cleaned = 0.227065 Jy ! Component: 150 - total flux cleaned = 0.28426 Jy ! Component: 200 - total flux cleaned = 0.31304 Jy ! Total flux subtracted in 200 components = 0.31304 Jy ! Clean residual min=-0.073496 max=0.089376 Jy/beam ! Clean residual mean=-0.000108 rms=0.017305 Jy/beam ! Combined flux in latest and established models = 2.78203 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 48 components and 2.78203 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.484508Jy sigma=14.523575 ! Fit after self-cal, rms=0.481576Jy sigma=14.353764 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.593 x 4.031 at 2.846 degrees (North through East) ! Clean map min=-0.068028 max=2.0647 Jy/beam if ( peak(flux,max)/imstat(rms) < lim_self ) selflim lim_val ! Selfcal amplitude corrections will be limited to 0.833333 -> 1.2. ! Selfcal phase corrections will not be limited. gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* MK 1.00* NL 1.00* OV 1.00* ! PT 1.00* SC 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.19 HN 0.96 KP 0.91 ! LA 1.10 MK 1.20 NL 1.15 OV 1.00 ! PT 0.97 SC 1.00 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! LA 1.00* MK 1.00* NL 1.00* OV 1.00* ! PT 1.00* SC 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.08 HN 0.95 KP 0.96 ! LA 1.01 MK 1.19 NL 1.06 OV 0.94 ! PT 0.99 SC 1.01 ! ! ! Correcting IF 5. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 0.99 HN 0.96 KP 0.97 ! LA 1.05 MK 1.15 NL 1.08 OV 0.92 ! PT 0.91 SC 1.03 ! ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 0.99 HN 0.98 KP 0.97 ! LA 1.10 MK 1.09 NL 1.05 OV 0.93 ! PT 1.01 SC 1.07 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 0.92 HN 0.92 KP 0.89 ! LA 1.01 MK 1.06 NL 1.07 OV 0.92 ! PT 0.93 SC 0.97 ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 0.97 HN 0.96 KP 0.96 ! LA 1.02 MK 1.07 NL 1.07 OV 0.96 ! PT 0.92 SC 1.02 ! ! ! Fit before self-cal, rms=0.481576Jy sigma=14.353764 ! Fit after self-cal, rms=0.448498Jy sigma=11.482273 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.616 mas, bmaj=4.05 mas, bpa=3.184 degrees ! Estimated noise=0.512634 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.0174781 Jy ! Component: 100 - total flux cleaned = 0.0259966 Jy ! Component: 150 - total flux cleaned = 0.0281701 Jy ! Component: 200 - total flux cleaned = 0.0286335 Jy ! Total flux subtracted in 200 components = 0.0286335 Jy ! Clean residual min=-0.038054 max=0.044595 Jy/beam ! Clean residual mean=-0.000055 rms=0.010167 Jy/beam ! Combined flux in latest and established models = 2.81066 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 68 components and 2.81066 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.446953Jy sigma=11.353378 ! Fit after self-cal, rms=0.446781Jy sigma=11.347345 ! 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.446781Jy sigma=11.347345 ! Fit after self-cal, rms=0.446265Jy sigma=11.328350 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.617 mas, bmaj=4.049 mas, bpa=3.191 degrees ! Estimated noise=0.515395 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0176141 Jy ! Component: 100 - total flux cleaned = 0.0205898 Jy ! Component: 150 - total flux cleaned = 0.0214429 Jy ! Component: 200 - total flux cleaned = 0.0222366 Jy ! Total flux subtracted in 200 components = 0.0222366 Jy ! Clean residual min=-0.038190 max=0.044953 Jy/beam ! Clean residual mean=-0.000057 rms=0.010140 Jy/beam ! Combined flux in latest and established models = 2.8329 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 80 components and 2.8329 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.445518Jy sigma=11.243565 ! Fit after self-cal, rms=0.445384Jy sigma=11.240089 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) ! Added new window around map position (-27.6, -22.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0440348 Jy ! Component: 100 - total flux cleaned = 0.068122 Jy ! Component: 150 - total flux cleaned = 0.0784817 Jy ! Component: 200 - total flux cleaned = 0.083462 Jy ! Total flux subtracted in 200 components = 0.083462 Jy ! Clean residual min=-0.034884 max=0.032151 Jy/beam ! Clean residual mean=-0.000036 rms=0.008325 Jy/beam ! Combined flux in latest and established models = 2.91636 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 117 components and 2.91636 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.433286Jy sigma=11.220960 ! Fit after self-cal, rms=0.432944Jy sigma=11.192347 ! 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.432944Jy sigma=11.192347 ! Fit after self-cal, rms=0.432960Jy sigma=11.064388 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.615 mas, bmaj=4.046 mas, bpa=3.206 degrees ! Estimated noise=0.518254 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0301163 Jy ! Component: 100 - total flux cleaned = 0.0473485 Jy ! Component: 150 - total flux cleaned = 0.0605631 Jy ! Component: 200 - total flux cleaned = 0.0677549 Jy ! Total flux subtracted in 200 components = 0.0677549 Jy ! Clean residual min=-0.037494 max=0.032487 Jy/beam ! Clean residual mean=-0.000036 rms=0.008349 Jy/beam ! Combined flux in latest and established models = 2.98411 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 132 components and 2.98411 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.426055Jy sigma=11.098002 ! Fit after self-cal, rms=0.425970Jy sigma=11.084144 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.00821436 Jy ! Component: 100 - total flux cleaned = 0.0102106 Jy ! Component: 150 - total flux cleaned = 0.0116813 Jy ! Component: 200 - total flux cleaned = 0.0127154 Jy ! Total flux subtracted in 200 components = 0.0127154 Jy ! Clean residual min=-0.036827 max=0.031170 Jy/beam ! Clean residual mean=-0.000033 rms=0.008082 Jy/beam ! Combined flux in latest and established models = 2.99683 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 143 components and 2.99683 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.424025Jy sigma=11.098471 ! Fit after self-cal, rms=0.423987Jy sigma=11.096187 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 40 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.423987Jy sigma=11.096187 ! Fit after self-cal, rms=0.424495Jy sigma=11.022833 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.645 mas, bmaj=5.327 mas, bpa=0.8927 degrees ! Estimated noise=0.644041 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0344585 Jy ! Component: 100 - total flux cleaned = 0.0518111 Jy ! Component: 150 - total flux cleaned = 0.0586218 Jy ! Component: 200 - total flux cleaned = 0.0617509 Jy ! Total flux subtracted in 200 components = 0.0617509 Jy ! Clean residual min=-0.052254 max=0.047989 Jy/beam ! Clean residual mean=-0.000070 rms=0.012976 Jy/beam ! Combined flux in latest and established models = 3.05858 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 158 components and 3.05858 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.419627Jy sigma=11.139105 ! Fit after self-cal, rms=0.419559Jy sigma=11.130342 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.00583113 Jy ! Component: 100 - total flux cleaned = 0.00740445 Jy ! Component: 150 - total flux cleaned = 0.0084492 Jy ! Component: 200 - total flux cleaned = 0.00974278 Jy ! Total flux subtracted in 200 components = 0.00974278 Jy ! Clean residual min=-0.051251 max=0.045924 Jy/beam ! Clean residual mean=-0.000066 rms=0.012626 Jy/beam ! Combined flux in latest and established models = 3.06832 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 163 components and 3.06832 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.418606Jy sigma=11.152827 ! Fit after self-cal, rms=0.418542Jy sigma=11.150686 ! 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=1.613 mas, bmaj=4.043 mas, bpa=3.195 degrees ! Estimated noise=0.519635 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.613 x 4.043 at 3.195 degrees (North through East) ! Clean map min=-0.03623 max=2.0772 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 ! Performing amp+phase self-cal ! ! 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.418542Jy sigma=11.150686 ! Fit after self-cal, rms=0.416792Jy sigma=10.908739 ! Inverting map and beam ! Estimated beam: bmin=1.614 mas, bmaj=4.047 mas, bpa=3.254 degrees ! Estimated noise=0.517243 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00938639 Jy ! Component: 100 - total flux cleaned = 0.0111813 Jy ! Component: 150 - total flux cleaned = 0.0130884 Jy ! Component: 200 - total flux cleaned = 0.0126591 Jy ! Total flux subtracted in 200 components = 0.0126591 Jy ! Clean residual min=-0.040641 max=0.037031 Jy/beam ! Clean residual mean=-0.000030 rms=0.008987 Jy/beam ! Combined flux in latest and established models = 3.08098 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 182 components and 3.08098 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.412384Jy sigma=10.935771 ! Fit after self-cal, rms=0.412380Jy sigma=10.929567 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00299305 Jy ! Component: 100 - total flux cleaned = 0.00337637 Jy ! Component: 150 - total flux cleaned = 0.00337603 Jy ! Component: 200 - total flux cleaned = 0.00338918 Jy ! Total flux subtracted in 200 components = 0.00338918 Jy ! Clean residual min=-0.040429 max=0.035282 Jy/beam ! Clean residual mean=-0.000029 rms=0.008803 Jy/beam ! Combined flux in latest and established models = 3.08437 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 189 components and 3.08437 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.410621Jy sigma=10.941656 ! Fit after self-cal, rms=0.410631Jy sigma=10.940286 ! Inverting map ! Performing amp+phase self-cal ! ! 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.410631Jy sigma=10.940286 ! Fit after self-cal, rms=0.410488Jy sigma=10.911414 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.354 mas, bmaj=3.82 mas, bpa=1.839 degrees ! Estimated noise=0.666166 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.5673 Jy ! Component: 100 - total flux cleaned = 2.09489 Jy ! Component: 150 - total flux cleaned = 2.35709 Jy ! Component: 200 - total flux cleaned = 2.50574 Jy ! Total flux subtracted in 200 components = 2.50574 Jy ! Clean residual min=-0.046188 max=0.081285 Jy/beam ! Clean residual mean=-0.000047 rms=0.015473 Jy/beam ! Combined flux in latest and established models = 2.50574 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 11 model components to the UV plane model. ! The established model now contains 11 components and 2.50574 Jy ! Inverting map and beam ! Estimated beam: bmin=1.616 mas, bmaj=4.049 mas, bpa=3.29 degrees ! Estimated noise=0.516777 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.138159 Jy ! Component: 100 - total flux cleaned = 0.232253 Jy ! Component: 150 - total flux cleaned = 0.298567 Jy ! Component: 200 - total flux cleaned = 0.349947 Jy ! Total flux subtracted in 200 components = 0.349947 Jy ! Clean residual min=-0.044905 max=0.052351 Jy/beam ! Clean residual mean=-0.000054 rms=0.012225 Jy/beam ! Combined flux in latest and established models = 2.85568 Jy ! Adding 47 model components to the UV plane model. ! The established model now contains 54 components and 2.85568 Jy ! Inverting map ! Added new window around map position (-5.4, 66.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0496496 Jy ! Component: 100 - total flux cleaned = 0.0859644 Jy ! Component: 150 - total flux cleaned = 0.1165 Jy ! Component: 200 - total flux cleaned = 0.142335 Jy ! Total flux subtracted in 200 components = 0.142335 Jy ! Clean residual min=-0.040926 max=0.038419 Jy/beam ! Clean residual mean=-0.000034 rms=0.009753 Jy/beam ! Combined flux in latest and established models = 2.99802 Jy ! Adding 40 model components to the UV plane model. ! The established model now contains 85 components and 2.99802 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 40 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.652 mas, bmaj=5.339 mas, bpa=0.9498 degrees ! Estimated noise=0.638582 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0389946 Jy ! Component: 100 - total flux cleaned = 0.058169 Jy ! Component: 150 - total flux cleaned = 0.0692594 Jy ! Component: 200 - total flux cleaned = 0.0792789 Jy ! Total flux subtracted in 200 components = 0.0792789 Jy ! Clean residual min=-0.049062 max=0.048814 Jy/beam ! Clean residual mean=-0.000053 rms=0.013589 Jy/beam ! Combined flux in latest and established models = 3.0773 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.00752446 Jy ! Component: 100 - total flux cleaned = 0.0136934 Jy ! Component: 150 - total flux cleaned = 0.0192098 Jy ! Component: 200 - total flux cleaned = 0.0225402 Jy ! Total flux subtracted in 200 components = 0.0225402 Jy ! Clean residual min=-0.050065 max=0.048260 Jy/beam ! Clean residual mean=-0.000047 rms=0.013159 Jy/beam ! Combined flux in latest and established models = 3.09984 Jy ! Adding 60 model components to the UV plane model. ! The established model now contains 137 components and 3.09984 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.616 mas, bmaj=4.049 mas, bpa=3.29 degrees ! Estimated noise=0.516777 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0040829 Jy ! Component: 100 - total flux cleaned = -0.00112881 Jy ! Component: 150 - total flux cleaned = -7.14579e-05 Jy ! Component: 200 - total flux cleaned = 0.000860288 Jy ! Total flux subtracted in 200 components = 0.000860288 Jy ! Clean residual min=-0.038546 max=0.032969 Jy/beam ! Clean residual mean=-0.000020 rms=0.008502 Jy/beam ! Combined flux in latest and established models = 3.1007 Jy selfcal ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 157 components and 3.1007 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.405287Jy sigma=10.947178 ! Fit after self-cal, rms=0.404794Jy sigma=10.928936 wmodel /d2/scr/pima/bp177i_uvs/J0854+2006_C_map.mod ! Writing 157 model components to file: /d2/scr/pima/bp177i_uvs/J0854+2006_C_map.mod wobs /d2/scr/pima/bp177i_uvs/J0854+2006_C_uvs.fits ! Writing UV FITS file: /d2/scr/pima/bp177i_uvs/J0854+2006_C_uvs.fits wwins /d2/scr/pima/bp177i_uvs/J0854+2006_C_map.win ! wwins: Wrote 4 windows to /d2/scr/pima/bp177i_uvs/J0854+2006_C_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.0129174 Jy ! Component: 100 - total flux cleaned = -0.0181128 Jy ! Component: 150 - total flux cleaned = -0.0202424 Jy ! Component: 200 - total flux cleaned = -0.0192568 Jy ! Component: 250 - total flux cleaned = -0.0178164 Jy ! Component: 300 - total flux cleaned = -0.0155049 Jy ! Component: 350 - total flux cleaned = -0.015064 Jy ! Component: 400 - total flux cleaned = -0.0129112 Jy ! Component: 450 - total flux cleaned = -0.0108263 Jy ! Component: 500 - total flux cleaned = -0.00799944 Jy ! Component: 550 - total flux cleaned = -0.00603884 Jy ! Component: 600 - total flux cleaned = -0.00336176 Jy ! Component: 650 - total flux cleaned = -0.00112387 Jy ! Component: 700 - total flux cleaned = 0.00141038 Jy ! Component: 750 - total flux cleaned = 0.00247329 Jy ! Component: 800 - total flux cleaned = 0.00420577 Jy ! Component: 850 - total flux cleaned = 0.00522125 Jy ! Component: 900 - total flux cleaned = 0.00554472 Jy ! Component: 950 - total flux cleaned = 0.00619658 Jy ! Component: 1000 - total flux cleaned = 0.0061977 Jy ! Component: 1050 - total flux cleaned = 0.00745527 Jy ! Component: 1100 - total flux cleaned = 0.00683609 Jy ! Component: 1150 - total flux cleaned = 0.0068448 Jy ! Component: 1200 - total flux cleaned = 0.00804611 Jy ! Component: 1250 - total flux cleaned = 0.00892895 Jy ! Component: 1300 - total flux cleaned = 0.00950673 Jy ! Component: 1350 - total flux cleaned = 0.00978869 Jy ! Component: 1400 - total flux cleaned = 0.00978808 Jy ! Component: 1450 - total flux cleaned = 0.0100658 Jy ! Component: 1500 - total flux cleaned = 0.00979405 Jy ! Component: 1550 - total flux cleaned = 0.0100558 Jy ! Component: 1600 - total flux cleaned = 0.00952856 Jy ! Component: 1650 - total flux cleaned = 0.0100489 Jy ! Component: 1700 - total flux cleaned = 0.00876014 Jy ! Component: 1750 - total flux cleaned = 0.00901506 Jy ! Component: 1800 - total flux cleaned = 0.00851402 Jy ! Component: 1850 - total flux cleaned = 0.00727961 Jy ! Component: 1900 - total flux cleaned = 0.00728319 Jy ! Component: 1950 - total flux cleaned = 0.00776094 Jy ! Component: 2000 - total flux cleaned = 0.00704066 Jy ! Component: 2050 - total flux cleaned = 0.00656783 Jy ! Component: 2100 - total flux cleaned = 0.00610336 Jy ! Component: 2150 - total flux cleaned = 0.00633008 Jy ! Component: 2200 - total flux cleaned = 0.00519642 Jy ! Component: 2250 - total flux cleaned = 0.00542392 Jy ! Component: 2300 - total flux cleaned = 0.00587013 Jy ! Component: 2350 - total flux cleaned = 0.0056462 Jy ! Component: 2400 - total flux cleaned = 0.00478159 Jy ! Component: 2450 - total flux cleaned = 0.00542639 Jy ! Component: 2500 - total flux cleaned = 0.00457883 Jy ! Component: 2550 - total flux cleaned = 0.00479091 Jy ! Component: 2600 - total flux cleaned = 0.00499767 Jy ! Component: 2650 - total flux cleaned = 0.00478531 Jy ! Component: 2700 - total flux cleaned = 0.003768 Jy ! Component: 2750 - total flux cleaned = 0.00296586 Jy ! Component: 2800 - total flux cleaned = 0.00316443 Jy ! Component: 2850 - total flux cleaned = 0.00316274 Jy ! Component: 2900 - total flux cleaned = 0.00257777 Jy ! Component: 2950 - total flux cleaned = 0.00161367 Jy ! Component: 3000 - total flux cleaned = 0.00180686 Jy ! Component: 3050 - total flux cleaned = 0.00142247 Jy ! Component: 3100 - total flux cleaned = 0.000670625 Jy ! Component: 3150 - total flux cleaned = 0.000299632 Jy ! Component: 3200 - total flux cleaned = -0.000252547 Jy ! Component: 3250 - total flux cleaned = -0.00116634 Jy ! Component: 3300 - total flux cleaned = -0.00134555 Jy ! Component: 3350 - total flux cleaned = -0.00224296 Jy ! Component: 3400 - total flux cleaned = -0.00241923 Jy ! Component: 3450 - total flux cleaned = -0.00259374 Jy ! Component: 3500 - total flux cleaned = -0.00398985 Jy ! Component: 3550 - total flux cleaned = -0.00433642 Jy ! Component: 3600 - total flux cleaned = -0.00467836 Jy ! Component: 3650 - total flux cleaned = -0.00485074 Jy ! Component: 3700 - total flux cleaned = -0.00535882 Jy ! Component: 3750 - total flux cleaned = -0.00636246 Jy ! Component: 3800 - total flux cleaned = -0.00669279 Jy ! Component: 3850 - total flux cleaned = -0.0076814 Jy ! Component: 3900 - total flux cleaned = -0.00849652 Jy ! Component: 3950 - total flux cleaned = -0.00898064 Jy ! Component: 4000 - total flux cleaned = -0.00962058 Jy ! Component: 4050 - total flux cleaned = -0.0102578 Jy ! Total flux subtracted in 4096 components = -0.0107295 Jy ! Clean residual min=-0.009616 max=0.008496 Jy/beam ! Clean residual mean=0.000016 rms=0.002976 Jy/beam ! Combined flux in latest and established models = 3.08997 Jy keep ! Adding 321 model components to the UV plane model. ! The established model now contains 478 components and 3.08997 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/d2/scr/pima/bp177i_uvs/J0854+2006_C\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.616 x 4.049 at 3.29 degrees (North through East) ! Clean map min=-0.021154 max=2.0639 Jy/beam ! Writing clean map to FITS file: /d2/scr/pima/bp177i_uvs/J0854+2006_C_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 2.06385 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00294172 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.427607 0.427607 0.855213 1.71043 3.42085 6.84171 13.6834 27.3668 54.7337 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0088 print "peak= ", int(1000*peak_flux)/1000 ! peak= 2.063 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 701.579 quit ! Quitting program ! Log file /d2/scr/pima/bp177i_uvs/J0854+2006_C_uvs.log closed on Tue Apr 15 00:20:27 2014