! Started logfile: /scr/pima/br149ag_uvs/J1740+5211_X_uvs.log on Fri Feb 19 09:47:29 2016 obs /scr/pima/br149ag_uvs/J1740+5211_X_uva.fits ! Reading UV FITS file: /scr/pima/br149ag_uvs/J1740+5211_X_uva.fits ! AN table 1: 13 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.965812 visibilities/baseline/integration-bin. ! Found source: J1740+5211 ! ! 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 8.10949e+09 1.28e+08 1 1.28e+08 ! 02 2 8.12349e+09 1.28e+08 1 1.28e+08 ! 03 3 8.16549e+09 1.28e+08 1 1.28e+08 ! 04 4 8.23549e+09 1.28e+08 1 1.28e+08 ! 05 5 8.31949e+09 1.28e+08 1 1.28e+08 ! 06 6 8.41749e+09 1.28e+08 1 1.28e+08 ! 07 7 8.55649e+09 1.28e+08 1 1.28e+08 ! 08 8 8.58349e+09 1.28e+08 1 1.28e+08 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3616 visibilities. select rr, 1,2,3,4,6,7,8 ! Selecting polarization: RR, channels: 1..4 6..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: (none) ! 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.2 print "taper_size =", taper_size ! taper_size = 70 ![@/home/lpetrov/bin/01_muppet /scr/pima/br149ag_uvs/J1740+5211_X] 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.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 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. ! ! Not correcting unselected IF 5. ! ! 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.900525Jy sigma=44.937189 ! Fit after self-cal, rms=0.637999Jy sigma=15.338879 ! 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.6319 mas, bmaj=2.384 mas, bpa=-28.75 degrees ! Estimated noise=0.939962 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.19593 Jy ! Component: 100 - total flux cleaned = 1.47033 Jy ! Total flux subtracted in 100 components = 1.47033 Jy ! Clean residual min=-0.059628 max=0.104027 Jy/beam ! Clean residual mean=0.000013 rms=0.016273 Jy/beam ! Combined flux in latest and established models = 1.47033 Jy ! Performing phase self-cal ! Adding 5 model components to the UV plane model. ! The established model now contains 5 components and 1.47033 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.284244Jy sigma=6.839615 ! Fit after self-cal, rms=0.284233Jy sigma=6.839331 ! Inverting map ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.113069 Jy ! Component: 100 - total flux cleaned = 0.171973 Jy ! Total flux subtracted in 100 components = 0.171973 Jy ! Clean residual min=-0.043921 max=0.059974 Jy/beam ! Clean residual mean=-0.000007 rms=0.011934 Jy/beam ! Combined flux in latest and established models = 1.6423 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 10 components and 1.6423 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.230878Jy sigma=5.495231 ! Fit after self-cal, rms=0.230869Jy sigma=5.494865 ! 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.7937 mas, bmaj=2.45 mas, bpa=-29.47 degrees ! Estimated noise=0.541177 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.0426722 Jy ! Component: 100 - total flux cleaned = 0.0617615 Jy ! Component: 150 - total flux cleaned = 0.0738974 Jy ! Component: 200 - total flux cleaned = 0.0822802 Jy ! Total flux subtracted in 200 components = 0.0822802 Jy ! Clean residual min=-0.065033 max=0.084934 Jy/beam ! Clean residual mean=-0.000039 rms=0.016347 Jy/beam ! Combined flux in latest and established models = 1.72458 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 29 components and 1.72458 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.220797Jy sigma=5.236107 ! Fit after self-cal, rms=0.220753Jy sigma=5.234966 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! restore ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7937 x 2.45 at -29.47 degrees (North through East) ! Clean map min=-0.060622 max=1.5832 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* ! 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 1.00* FD 1.00* HN 1.00* KP 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! 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* ! 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 1.00 FD 1.02 HN 0.96 KP 0.91 ! MK 1.02 NL 0.97 OV 1.06 PT 1.00 ! SC 1.00 ! ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.00 HN 1.05 KP 0.93 ! MK 0.98 NL 0.96 OV 0.98 PT 0.98 ! SC 1.01 ! ! ! Correcting IF 7. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00* FD 1.00* HN 1.00* KP 1.00* ! MK 1.00* NL 1.00* OV 1.00* PT 1.00* ! SC 1.00* ! ! ! Correcting IF 8. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 0.99 HN 0.96 KP 0.86 ! MK 0.99 NL 0.97 OV 1.02 PT 0.98 ! SC 1.01 ! ! ! Fit before self-cal, rms=0.220753Jy sigma=5.234966 ! Fit after self-cal, rms=0.207143Jy sigma=4.868542 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.796 mas, bmaj=2.449 mas, bpa=-29.54 degrees ! Estimated noise=0.53178 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 (5.2, -23.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0618294 Jy ! Component: 100 - total flux cleaned = 0.069458 Jy ! Component: 150 - total flux cleaned = 0.0802989 Jy ! Component: 200 - total flux cleaned = 0.0914003 Jy ! Total flux subtracted in 200 components = 0.0914003 Jy ! Clean residual min=-0.039460 max=0.052349 Jy/beam ! Clean residual mean=-0.000013 rms=0.010759 Jy/beam ! Combined flux in latest and established models = 1.81598 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 62 components and 1.81598 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.180144Jy sigma=4.258123 ! Fit after self-cal, rms=0.178948Jy sigma=4.227463 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.178948Jy sigma=4.227463 ! Fit after self-cal, rms=0.176346Jy sigma=4.166428 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7963 mas, bmaj=2.451 mas, bpa=-29.55 degrees ! Estimated noise=0.530806 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0200348 Jy ! Component: 100 - total flux cleaned = 0.0304404 Jy ! Component: 150 - total flux cleaned = 0.0368737 Jy ! Component: 200 - total flux cleaned = 0.0410716 Jy ! Total flux subtracted in 200 components = 0.0410716 Jy ! Clean residual min=-0.037227 max=0.046782 Jy/beam ! Clean residual mean=-0.000011 rms=0.010198 Jy/beam ! Combined flux in latest and established models = 1.85705 Jy selfcal ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 87 components and 1.85705 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.167419Jy sigma=3.961963 ! Fit after self-cal, rms=0.167049Jy sigma=3.951027 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.00972942 Jy ! Component: 100 - total flux cleaned = 0.011563 Jy ! Component: 150 - total flux cleaned = 0.0126499 Jy ! Component: 200 - total flux cleaned = 0.013678 Jy ! Total flux subtracted in 200 components = 0.013678 Jy ! Clean residual min=-0.034602 max=0.040992 Jy/beam ! Clean residual mean=-0.000008 rms=0.009668 Jy/beam ! Combined flux in latest and established models = 1.87073 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 100 components and 1.87073 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.163395Jy sigma=3.866255 ! Fit after self-cal, rms=0.163318Jy sigma=3.863068 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.163318Jy sigma=3.863068 ! Fit after self-cal, rms=0.161653Jy sigma=3.827115 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.796 mas, bmaj=2.452 mas, bpa=-29.54 degrees ! Estimated noise=0.529468 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00886301 Jy ! Component: 100 - total flux cleaned = 0.0139137 Jy ! Component: 150 - total flux cleaned = 0.0154451 Jy ! Component: 200 - total flux cleaned = 0.0159111 Jy ! Total flux subtracted in 200 components = 0.0159111 Jy ! Clean residual min=-0.034894 max=0.042232 Jy/beam ! Clean residual mean=-0.000010 rms=0.009670 Jy/beam ! Combined flux in latest and established models = 1.88664 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 115 components and 1.88664 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.158298Jy sigma=3.748543 ! Fit after self-cal, rms=0.158212Jy sigma=3.745664 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.00308213 Jy ! Component: 100 - total flux cleaned = 0.00264063 Jy ! Component: 150 - total flux cleaned = 0.00346794 Jy ! Component: 200 - total flux cleaned = 0.00422791 Jy ! Total flux subtracted in 200 components = 0.00422791 Jy ! Clean residual min=-0.034274 max=0.039754 Jy/beam ! Clean residual mean=-0.000008 rms=0.009412 Jy/beam ! Combined flux in latest and established models = 1.89087 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 125 components and 1.89087 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.156279Jy sigma=3.700189 ! Fit after self-cal, rms=0.156261Jy sigma=3.699093 ! Inverting map ! Added new window around map position (-31.6, 39.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0361651 Jy ! Component: 100 - total flux cleaned = 0.0463243 Jy ! Component: 150 - total flux cleaned = 0.0515284 Jy ! Component: 200 - total flux cleaned = 0.0565885 Jy ! Total flux subtracted in 200 components = 0.0565885 Jy ! Clean residual min=-0.035640 max=0.032988 Jy/beam ! Clean residual mean=-0.000007 rms=0.008205 Jy/beam ! Combined flux in latest and established models = 1.94746 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 162 components and 1.94746 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.145029Jy sigma=3.444316 ! Fit after self-cal, rms=0.144497Jy sigma=3.427396 ! 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.144497Jy sigma=3.427396 ! Fit after self-cal, rms=0.138695Jy sigma=3.275597 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.846 mas, bmaj=2.596 mas, bpa=-34.18 degrees ! Estimated noise=0.706806 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0271706 Jy ! Component: 100 - total flux cleaned = 0.0309125 Jy ! Component: 150 - total flux cleaned = 0.0296786 Jy ! Component: 200 - total flux cleaned = 0.0285488 Jy ! Total flux subtracted in 200 components = 0.0285488 Jy ! Clean residual min=-0.038590 max=0.036946 Jy/beam ! Clean residual mean=-0.000015 rms=0.010835 Jy/beam ! Combined flux in latest and established models = 1.97601 Jy selfcal ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 190 components and 1.97601 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.135529Jy sigma=3.204263 ! Fit after self-cal, rms=0.135413Jy sigma=3.198112 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.00323573 Jy ! Component: 100 - total flux cleaned = 0.000694385 Jy ! Component: 150 - total flux cleaned = -0.000246439 Jy ! Component: 200 - total flux cleaned = -0.000693719 Jy ! Total flux subtracted in 200 components = -0.000693719 Jy ! Clean residual min=-0.036822 max=0.034044 Jy/beam ! Clean residual mean=-0.000017 rms=0.010187 Jy/beam ! Combined flux in latest and established models = 1.97531 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 207 components and 1.97531 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.133922Jy sigma=3.163094 ! Fit after self-cal, rms=0.133928Jy sigma=3.161770 ! 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.7949 mas, bmaj=2.45 mas, bpa=-29.55 degrees ! Estimated noise=0.529967 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7949 x 2.45 at -29.55 degrees (North through East) ! Clean map min=-0.029559 max=1.5154 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. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.133928Jy sigma=3.161770 ! Fit after self-cal, rms=0.132391Jy sigma=3.125458 ! Inverting map and beam ! Estimated beam: bmin=0.7941 mas, bmaj=2.451 mas, bpa=-29.56 degrees ! Estimated noise=0.529095 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000393823 Jy ! Component: 100 - total flux cleaned = -0.00120666 Jy ! Component: 150 - total flux cleaned = -0.000854453 Jy ! Component: 200 - total flux cleaned = -0.000519194 Jy ! Total flux subtracted in 200 components = -0.000519194 Jy ! Clean residual min=-0.030058 max=0.029221 Jy/beam ! Clean residual mean=-0.000010 rms=0.006982 Jy/beam ! Combined flux in latest and established models = 1.97479 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 233 components and 1.97479 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.130404Jy sigma=3.076651 ! Fit after self-cal, rms=0.130385Jy sigma=3.075584 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00100898 Jy ! Component: 100 - total flux cleaned = 0.000685593 Jy ! Component: 150 - total flux cleaned = 0.000691189 Jy ! Component: 200 - total flux cleaned = 0.000416323 Jy ! Total flux subtracted in 200 components = 0.000416323 Jy ! Clean residual min=-0.028891 max=0.028311 Jy/beam ! Clean residual mean=-0.000010 rms=0.006833 Jy/beam ! Combined flux in latest and established models = 1.97521 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 242 components and 1.97521 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.129197Jy sigma=3.046016 ! Fit after self-cal, rms=0.129195Jy sigma=3.045396 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.129195Jy sigma=3.045396 ! Fit after self-cal, rms=0.129027Jy sigma=3.041784 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.6323 mas, bmaj=2.392 mas, bpa=-28.79 degrees ! Estimated noise=0.911222 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.15859 Jy ! Component: 100 - total flux cleaned = 1.42466 Jy ! Component: 150 - total flux cleaned = 1.5448 Jy ! Component: 200 - total flux cleaned = 1.6066 Jy ! Total flux subtracted in 200 components = 1.6066 Jy ! Clean residual min=-0.045076 max=0.061338 Jy/beam ! Clean residual mean=0.000007 rms=0.011376 Jy/beam ! Combined flux in latest and established models = 1.6066 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 15 model components to the UV plane model. ! The established model now contains 15 components and 1.6066 Jy ! Inverting map and beam ! Estimated beam: bmin=0.794 mas, bmaj=2.451 mas, bpa=-29.57 degrees ! Estimated noise=0.528918 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 (-4.2, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0832536 Jy ! Component: 100 - total flux cleaned = 0.137824 Jy ! Component: 150 - total flux cleaned = 0.176969 Jy ! Component: 200 - total flux cleaned = 0.210478 Jy ! Total flux subtracted in 200 components = 0.210478 Jy ! Clean residual min=-0.039055 max=0.050355 Jy/beam ! Clean residual mean=0.000007 rms=0.010933 Jy/beam ! Combined flux in latest and established models = 1.81708 Jy ! Adding 57 model components to the UV plane model. ! The established model now contains 65 components and 1.81708 Jy ! Inverting map ! Added new window around map position (4.2, 17.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.031778 Jy ! Component: 100 - total flux cleaned = 0.0499152 Jy ! Component: 150 - total flux cleaned = 0.0636047 Jy ! Component: 200 - total flux cleaned = 0.0768475 Jy ! Total flux subtracted in 200 components = 0.0768475 Jy ! Clean residual min=-0.032182 max=0.033010 Jy/beam ! Clean residual mean=0.000003 rms=0.008322 Jy/beam ! Combined flux in latest and established models = 1.89393 Jy ! Adding 42 model components to the UV plane model. ! The established model now contains 101 components and 1.89393 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.845 mas, bmaj=2.595 mas, bpa=-34.15 degrees ! Estimated noise=0.706483 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.038205 Jy ! Component: 100 - total flux cleaned = 0.0543764 Jy ! Component: 150 - total flux cleaned = 0.0646795 Jy ! Component: 200 - total flux cleaned = 0.0750375 Jy ! Total flux subtracted in 200 components = 0.0750375 Jy ! Clean residual min=-0.033718 max=0.038677 Jy/beam ! Clean residual mean=0.000002 rms=0.009310 Jy/beam ! Combined flux in latest and established models = 1.96897 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.00809066 Jy ! Component: 100 - total flux cleaned = 0.0113622 Jy ! Component: 150 - total flux cleaned = 0.014337 Jy ! Component: 200 - total flux cleaned = 0.0154361 Jy ! Total flux subtracted in 200 components = 0.0154361 Jy ! Clean residual min=-0.030042 max=0.032841 Jy/beam ! Clean residual mean=-0.000002 rms=0.008240 Jy/beam ! Combined flux in latest and established models = 1.9844 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 168 components and 1.9844 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.794 mas, bmaj=2.451 mas, bpa=-29.57 degrees ! Estimated noise=0.528918 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00440276 Jy ! Component: 100 - total flux cleaned = -0.00491896 Jy ! Component: 150 - total flux cleaned = -0.00595213 Jy ! Component: 200 - total flux cleaned = -0.00691464 Jy ! Total flux subtracted in 200 components = -0.00691464 Jy ! Clean residual min=-0.024034 max=0.024851 Jy/beam ! Clean residual mean=-0.000003 rms=0.006197 Jy/beam ! Combined flux in latest and established models = 1.97749 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 194 components and 1.97749 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Not correcting unselected IF 5. ! ! Correcting IF 6. ! ! Correcting IF 7. ! ! Correcting IF 8. ! ! Fit before self-cal, rms=0.125346Jy sigma=2.930424 ! Fit after self-cal, rms=0.124761Jy sigma=2.917695 wmodel /scr/pima/br149ag_uvs/J1740+5211_X_map.mod ! Writing 194 model components to file: /scr/pima/br149ag_uvs/J1740+5211_X_map.mod wobs /scr/pima/br149ag_uvs/J1740+5211_X_uvs.fits ! Writing UV FITS file: /scr/pima/br149ag_uvs/J1740+5211_X_uvs.fits wwins /scr/pima/br149ag_uvs/J1740+5211_X_map.win ! wwins: Wrote 5 windows to /scr/pima/br149ag_uvs/J1740+5211_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.00175476 Jy ! Component: 100 - total flux cleaned = 0.00460084 Jy ! Component: 150 - total flux cleaned = 0.00805938 Jy ! Component: 200 - total flux cleaned = 0.0121182 Jy ! Component: 250 - total flux cleaned = 0.0149478 Jy ! Component: 300 - total flux cleaned = 0.0176774 Jy ! Component: 350 - total flux cleaned = 0.0196536 Jy ! Component: 400 - total flux cleaned = 0.0212595 Jy ! Component: 450 - total flux cleaned = 0.0234234 Jy ! Component: 500 - total flux cleaned = 0.0243171 Jy ! Component: 550 - total flux cleaned = 0.0248963 Jy ! Component: 600 - total flux cleaned = 0.0257481 Jy ! Component: 650 - total flux cleaned = 0.0271251 Jy ! Component: 700 - total flux cleaned = 0.0276582 Jy ! Component: 750 - total flux cleaned = 0.028969 Jy ! Component: 800 - total flux cleaned = 0.0299887 Jy ! Component: 850 - total flux cleaned = 0.0312321 Jy ! Component: 900 - total flux cleaned = 0.0314775 Jy ! Component: 950 - total flux cleaned = 0.0326662 Jy ! Component: 1000 - total flux cleaned = 0.0331294 Jy ! Component: 1050 - total flux cleaned = 0.0340421 Jy ! Component: 1100 - total flux cleaned = 0.0351504 Jy ! Component: 1150 - total flux cleaned = 0.0358031 Jy ! Component: 1200 - total flux cleaned = 0.0360201 Jy ! Component: 1250 - total flux cleaned = 0.0366444 Jy ! Component: 1300 - total flux cleaned = 0.0374661 Jy ! Component: 1350 - total flux cleaned = 0.0380671 Jy ! Component: 1400 - total flux cleaned = 0.0390547 Jy ! Component: 1450 - total flux cleaned = 0.0394438 Jy ! Component: 1500 - total flux cleaned = 0.0405867 Jy ! Component: 1550 - total flux cleaned = 0.0411442 Jy ! Component: 1600 - total flux cleaned = 0.0416946 Jy ! Component: 1650 - total flux cleaned = 0.0425951 Jy ! Component: 1700 - total flux cleaned = 0.0424169 Jy ! Component: 1750 - total flux cleaned = 0.0429385 Jy ! Component: 1800 - total flux cleaned = 0.0431102 Jy ! Component: 1850 - total flux cleaned = 0.044284 Jy ! Component: 1900 - total flux cleaned = 0.0451099 Jy ! Component: 1950 - total flux cleaned = 0.0455971 Jy ! Component: 2000 - total flux cleaned = 0.0468791 Jy ! Component: 2050 - total flux cleaned = 0.0473504 Jy ! Component: 2100 - total flux cleaned = 0.0482828 Jy ! Component: 2150 - total flux cleaned = 0.0492019 Jy ! Component: 2200 - total flux cleaned = 0.0498065 Jy ! Component: 2250 - total flux cleaned = 0.0505523 Jy ! Component: 2300 - total flux cleaned = 0.0509934 Jy ! Component: 2350 - total flux cleaned = 0.0518612 Jy ! Component: 2400 - total flux cleaned = 0.0528626 Jy ! Component: 2450 - total flux cleaned = 0.0535653 Jy ! Component: 2500 - total flux cleaned = 0.0539799 Jy ! Component: 2550 - total flux cleaned = 0.0550761 Jy ! Component: 2600 - total flux cleaned = 0.0556174 Jy ! Component: 2650 - total flux cleaned = 0.0565529 Jy ! Component: 2700 - total flux cleaned = 0.0573437 Jy ! Component: 2750 - total flux cleaned = 0.0579932 Jy ! Component: 2800 - total flux cleaned = 0.0586352 Jy ! Component: 2850 - total flux cleaned = 0.0590154 Jy ! Component: 2900 - total flux cleaned = 0.059516 Jy ! Component: 2950 - total flux cleaned = 0.0601345 Jy ! Component: 3000 - total flux cleaned = 0.0607446 Jy ! Component: 3050 - total flux cleaned = 0.0611054 Jy ! Component: 3100 - total flux cleaned = 0.0611051 Jy ! Component: 3150 - total flux cleaned = 0.061577 Jy ! Component: 3200 - total flux cleaned = 0.0622748 Jy ! Component: 3250 - total flux cleaned = 0.0626183 Jy ! Component: 3300 - total flux cleaned = 0.0630719 Jy ! Component: 3350 - total flux cleaned = 0.0631832 Jy ! Component: 3400 - total flux cleaned = 0.064186 Jy ! Component: 3450 - total flux cleaned = 0.0642957 Jy ! Component: 3500 - total flux cleaned = 0.0645143 Jy ! Component: 3550 - total flux cleaned = 0.0653797 Jy ! Component: 3600 - total flux cleaned = 0.0655926 Jy ! Component: 3650 - total flux cleaned = 0.0656986 Jy ! Component: 3700 - total flux cleaned = 0.0659096 Jy ! Component: 3750 - total flux cleaned = 0.0665335 Jy ! Component: 3800 - total flux cleaned = 0.0668431 Jy ! Component: 3850 - total flux cleaned = 0.067148 Jy ! Component: 3900 - total flux cleaned = 0.0674517 Jy ! Component: 3950 - total flux cleaned = 0.0676532 Jy ! Component: 4000 - total flux cleaned = 0.0679529 Jy ! Component: 4050 - total flux cleaned = 0.0676564 Jy ! Total flux subtracted in 4096 components = 0.0681504 Jy ! Clean residual min=-0.005346 max=0.005906 Jy/beam ! Clean residual mean=0.000002 rms=0.001760 Jy/beam ! Combined flux in latest and established models = 2.04564 Jy keep ! Adding 319 model components to the UV plane model. ! The established model now contains 513 components and 2.04564 Jy ![Exited script file: /home/lpetrov/bin/01_muppet] wmap "/scr/pima/br149ag_uvs/J1740+5211_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.794 x 2.451 at -29.57 degrees (North through East) ! Clean map min=-0.02321 max=1.4919 Jy/beam ! Writing clean map to FITS file: /scr/pima/br149ag_uvs/J1740+5211_X_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 1.49187 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00171206 loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.344279 0.344279 0.688559 1.37712 2.75423 5.50847 11.0169 22.0339 44.0678 88.1355 print "clev= ", int(3*image_rms*10000)/10000 ! clev= 0.0051 print "peak= ", int(1000*peak_flux)/1000 ! peak= 1.491 print "dyn_range= ", int(1000*peak_flux/image_rms)/1000 ! dyn_range= 871.385 quit ! Quitting program ! Log file /scr/pima/br149ag_uvs/J1740+5211_X_uvs.log closed on Fri Feb 19 09:47:33 2016