! Started logfile: /scr/pima/rdv94_uvs/J1321+8316_S_uvs.log on Tue Oct 8 12:20:11 2013 obs /scr/pima/rdv94_uvs/J1321+8316_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J1321+8316_S_uva.fits ! AN table 1: 36 integrations on 28 of 28 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J1321+8316 ! ! 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 2.23299e+09 1.6e+07 1 1.6e+07 ! 02 2 2.26299e+09 1.6e+07 1 1.6e+07 ! 03 3 2.35299e+09 1.6e+07 1 1.6e+07 ! 04 4 2.37299e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 4032 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.5 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv94_uvs/J1321+8316_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 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. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.013720Jy sigma=4.218877 ! Fit after self-cal, rms=0.845341Jy sigma=3.476340 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.748 mas, bmaj=9.875 mas, bpa=-36.25 degrees ! Estimated noise=5.5287 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.120888 Jy ! Component: 100 - total flux cleaned = 0.14936 Jy ! Total flux subtracted in 100 components = 0.14936 Jy ! Clean residual min=-0.008802 max=0.010577 Jy/beam ! Clean residual mean=-0.000005 rms=0.002566 Jy/beam ! Combined flux in latest and established models = 0.14936 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 0.14936 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.055685Jy sigma=0.247388 ! Fit after self-cal, rms=0.055667Jy sigma=0.247004 ! 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=3.569 mas, bmaj=10.39 mas, bpa=-35.93 degrees ! Estimated noise=2.74301 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.0166783 Jy ! Component: 100 - total flux cleaned = 0.0257711 Jy ! Component: 150 - total flux cleaned = 0.0314376 Jy ! Component: 200 - total flux cleaned = 0.0338201 Jy ! Total flux subtracted in 200 components = 0.0338201 Jy ! Clean residual min=-0.008212 max=0.010540 Jy/beam ! Clean residual mean=0.000026 rms=0.002443 Jy/beam ! Combined flux in latest and established models = 0.183181 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 38 components and 0.183181 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.051668Jy sigma=0.220301 ! Fit after self-cal, rms=0.051450Jy sigma=0.219569 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.88 FD 1.03 HN 1.08 KP 1.18 ! LA 1.23 MK 0.94 NL 0.90 OV 0.94 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.92 FD 1.10 HN 1.12 KP 1.20 ! LA 1.13 MK 1.00 NL 0.94 OV 0.97 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.10 HN 1.09 KP 1.01 ! LA 0.89 MK 1.02 NL 1.00 OV 1.00 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.87 FD 1.03 HN 0.97 KP 0.95 ! LA 0.74 MK 1.13 NL 0.92 OV 0.87 ! ! ! Fit before self-cal, rms=0.051450Jy sigma=0.219569 ! Fit after self-cal, rms=0.043539Jy sigma=0.183178 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.623 mas, bmaj=10.54 mas, bpa=-35.66 degrees ! Estimated noise=2.68689 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 (-9.5, -2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00941751 Jy ! Component: 100 - total flux cleaned = 0.0147358 Jy ! Component: 150 - total flux cleaned = 0.0180951 Jy ! Component: 200 - total flux cleaned = 0.0195966 Jy ! Total flux subtracted in 200 components = 0.0195966 Jy ! Clean residual min=-0.005300 max=0.007098 Jy/beam ! Clean residual mean=0.000038 rms=0.001317 Jy/beam ! Combined flux in latest and established models = 0.202777 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 84 components and 0.202777 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.041632Jy sigma=0.171523 ! Fit after self-cal, rms=0.041576Jy sigma=0.170923 ! Inverting map ! Added new window around map position (-25, -23). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0062192 Jy ! Component: 100 - total flux cleaned = 0.00988306 Jy ! Component: 150 - total flux cleaned = 0.0129124 Jy ! Component: 200 - total flux cleaned = 0.0155069 Jy ! Total flux subtracted in 200 components = 0.0155069 Jy ! Clean residual min=-0.003822 max=0.003212 Jy/beam ! Clean residual mean=0.000036 rms=0.000981 Jy/beam ! Combined flux in latest and established models = 0.218284 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 136 components and 0.218284 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.040556Jy sigma=0.165222 ! Fit after self-cal, rms=0.040484Jy sigma=0.164924 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.040484Jy sigma=0.164924 ! Fit after self-cal, rms=0.040262Jy sigma=0.164724 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.622 mas, bmaj=10.56 mas, bpa=-35.59 degrees ! Estimated noise=2.67766 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00267203 Jy ! Component: 100 - total flux cleaned = 0.00493431 Jy ! Component: 150 - total flux cleaned = 0.00687245 Jy ! Component: 200 - total flux cleaned = 0.00832479 Jy ! Total flux subtracted in 200 components = 0.00832479 Jy ! Clean residual min=-0.003267 max=0.002960 Jy/beam ! Clean residual mean=0.000034 rms=0.000832 Jy/beam ! Combined flux in latest and established models = 0.226609 Jy selfcal ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 177 components and 0.226609 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039973Jy sigma=0.163050 ! Fit after self-cal, rms=0.039954Jy sigma=0.162981 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.00159678 Jy ! Component: 100 - total flux cleaned = 0.00276501 Jy ! Component: 150 - total flux cleaned = 0.00376846 Jy ! Component: 200 - total flux cleaned = 0.00468339 Jy ! Total flux subtracted in 200 components = 0.00468339 Jy ! Clean residual min=-0.002991 max=0.002759 Jy/beam ! Clean residual mean=0.000034 rms=0.000765 Jy/beam ! Combined flux in latest and established models = 0.231292 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 207 components and 0.231292 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039805Jy sigma=0.162147 ! Fit after self-cal, rms=0.039792Jy sigma=0.162109 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039792Jy sigma=0.162109 ! Fit after self-cal, rms=0.039687Jy sigma=0.160836 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.62 mas, bmaj=10.56 mas, bpa=-35.6 degrees ! Estimated noise=2.67868 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000720497 Jy ! Component: 100 - total flux cleaned = 0.00158685 Jy ! Component: 150 - total flux cleaned = 0.00225324 Jy ! Component: 200 - total flux cleaned = 0.00281172 Jy ! Total flux subtracted in 200 components = 0.00281172 Jy ! Clean residual min=-0.002437 max=0.002453 Jy/beam ! Clean residual mean=0.000017 rms=0.000592 Jy/beam ! Combined flux in latest and established models = 0.234104 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 233 components and 0.234104 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039577Jy sigma=0.160291 ! Fit after self-cal, rms=0.039567Jy sigma=0.160228 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.000264975 Jy ! Component: 100 - total flux cleaned = 0.000679105 Jy ! Component: 150 - total flux cleaned = 0.00108663 Jy ! Component: 200 - total flux cleaned = 0.00152384 Jy ! Total flux subtracted in 200 components = 0.00152384 Jy ! Clean residual min=-0.002257 max=0.002209 Jy/beam ! Clean residual mean=0.000017 rms=0.000560 Jy/beam ! Combined flux in latest and established models = 0.235628 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 256 components and 0.235628 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039506Jy sigma=0.159948 ! Fit after self-cal, rms=0.039498Jy sigma=0.159915 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039498Jy sigma=0.159915 ! Fit after self-cal, rms=0.039429Jy sigma=0.159870 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.22 mas, bmaj=12.68 mas, bpa=-41.64 degrees ! Estimated noise=2.99432 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000802007 Jy ! Component: 100 - total flux cleaned = 0.00141695 Jy ! Component: 150 - total flux cleaned = 0.00199638 Jy ! Component: 200 - total flux cleaned = 0.00245183 Jy ! Total flux subtracted in 200 components = 0.00245183 Jy ! Clean residual min=-0.002493 max=0.002719 Jy/beam ! Clean residual mean=0.000030 rms=0.000768 Jy/beam ! Combined flux in latest and established models = 0.23808 Jy selfcal ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 304 components and 0.23808 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039385Jy sigma=0.159621 ! Fit after self-cal, rms=0.039374Jy sigma=0.159585 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.000127085 Jy ! Component: 100 - total flux cleaned = 0.000551056 Jy ! Component: 150 - total flux cleaned = 0.0010157 Jy ! Component: 200 - total flux cleaned = 0.00144508 Jy ! Total flux subtracted in 200 components = 0.00144508 Jy ! Clean residual min=-0.002374 max=0.002543 Jy/beam ! Clean residual mean=0.000029 rms=0.000742 Jy/beam ! Combined flux in latest and established models = 0.239525 Jy ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 332 components and 0.239525 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039355Jy sigma=0.159474 ! Fit after self-cal, rms=0.039351Jy sigma=0.159455 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039351Jy sigma=0.159455 ! Fit after self-cal, rms=0.039351Jy sigma=0.159455 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039351Jy sigma=0.159455 ! Fit after self-cal, rms=0.039288Jy sigma=0.159420 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.789 mas, bmaj=10.05 mas, bpa=-36.19 degrees ! Estimated noise=5.37979 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 (50.5, 93). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019289 Jy ! Component: 100 - total flux cleaned = 0.00193208 Jy ! Component: 150 - total flux cleaned = 0.00176796 Jy ! Component: 200 - total flux cleaned = 0.00169086 Jy ! Total flux subtracted in 200 components = 0.00169086 Jy ! Clean residual min=-0.003000 max=0.002668 Jy/beam ! Clean residual mean=0.000006 rms=0.000697 Jy/beam ! Combined flux in latest and established models = 0.241216 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 359 components and 0.241216 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039238Jy sigma=0.159319 ! Fit after self-cal, rms=0.039231Jy sigma=0.159228 ! Inverting map 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=3.619 mas, bmaj=10.57 mas, bpa=-35.61 degrees ! Estimated noise=2.67256 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00105638 Jy ! Component: 100 - total flux cleaned = 0.00167906 Jy ! Component: 150 - total flux cleaned = 0.00217969 Jy ! Component: 200 - total flux cleaned = 0.00255797 Jy ! Component: 250 - total flux cleaned = 0.00279466 Jy ! Component: 300 - total flux cleaned = 0.0030862 Jy ! Component: 350 - total flux cleaned = 0.003246 Jy ! Component: 400 - total flux cleaned = 0.0034181 Jy ! Total flux subtracted in 400 components = 0.0034181 Jy ! Clean residual min=-0.001701 max=0.001709 Jy/beam ! Clean residual mean=0.000007 rms=0.000434 Jy/beam ! Combined flux in latest and established models = 0.244634 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 414 components and 0.244634 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039130Jy sigma=0.158871 ! Fit after self-cal, rms=0.039124Jy sigma=0.158837 ! Inverting map dynam = dynam - 0.5 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) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000259778 Jy ! Component: 100 - total flux cleaned = 0.00053559 Jy ! Component: 150 - total flux cleaned = 0.000746712 Jy ! Component: 200 - total flux cleaned = 0.000845426 Jy ! Component: 250 - total flux cleaned = 0.000995891 Jy ! Component: 300 - total flux cleaned = 0.00112077 Jy ! Component: 350 - total flux cleaned = 0.00118963 Jy ! Component: 400 - total flux cleaned = 0.001289 Jy ! Total flux subtracted in 400 components = 0.001289 Jy ! Clean residual min=-0.001644 max=0.001683 Jy/beam ! Clean residual mean=0.000006 rms=0.000420 Jy/beam ! Combined flux in latest and established models = 0.245923 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 452 components and 0.245923 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.039082Jy sigma=0.158704 ! Fit after self-cal, rms=0.039079Jy sigma=0.158688 ! Inverting map ! Added new window around map position (-107, 114). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000951105 Jy ! Component: 100 - total flux cleaned = 0.00119658 Jy ! Component: 150 - total flux cleaned = 0.00123736 Jy ! Component: 200 - total flux cleaned = 0.00133419 Jy ! Component: 250 - total flux cleaned = 0.00145063 Jy ! Component: 300 - total flux cleaned = 0.00153286 Jy ! Component: 350 - total flux cleaned = 0.0016895 Jy ! Component: 400 - total flux cleaned = 0.00181238 Jy ! Total flux subtracted in 400 components = 0.00181238 Jy ! Clean residual min=-0.001320 max=0.001231 Jy/beam ! Clean residual mean=0.000005 rms=0.000343 Jy/beam ! Combined flux in latest and established models = 0.247735 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 503 components and 0.247735 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038943Jy sigma=0.158038 ! Fit after self-cal, rms=0.038940Jy sigma=0.158017 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038940Jy sigma=0.158017 ! Fit after self-cal, rms=0.038886Jy sigma=0.157982 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.615 mas, bmaj=10.57 mas, bpa=-35.63 degrees ! Estimated noise=2.66936 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000142424 Jy ! Component: 100 - total flux cleaned = 0.00021674 Jy ! Component: 150 - total flux cleaned = 0.000286891 Jy ! Component: 200 - total flux cleaned = 0.000354365 Jy ! Component: 250 - total flux cleaned = 0.000440921 Jy ! Component: 300 - total flux cleaned = 0.000544397 Jy ! Component: 350 - total flux cleaned = 0.000684288 Jy ! Component: 400 - total flux cleaned = 0.000800738 Jy ! Total flux subtracted in 400 components = 0.000800738 Jy ! Clean residual min=-0.001159 max=0.001204 Jy/beam ! Clean residual mean=0.000003 rms=0.000313 Jy/beam ! Combined flux in latest and established models = 0.248536 Jy selfcal ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 541 components and 0.248536 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038839Jy sigma=0.157792 ! Fit after self-cal, rms=0.038839Jy sigma=0.157786 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 (-87, -113). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000670664 Jy ! Component: 100 - total flux cleaned = 0.00123813 Jy ! Component: 150 - total flux cleaned = 0.00151863 Jy ! Component: 200 - total flux cleaned = 0.00182263 Jy ! Component: 250 - total flux cleaned = 0.00212136 Jy ! Component: 300 - total flux cleaned = 0.00242533 Jy ! Component: 350 - total flux cleaned = 0.00273373 Jy ! Component: 400 - total flux cleaned = 0.00292399 Jy ! Total flux subtracted in 400 components = 0.00292399 Jy ! Clean residual min=-0.001098 max=0.001083 Jy/beam ! Clean residual mean=0.000004 rms=0.000270 Jy/beam ! Combined flux in latest and established models = 0.25146 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 612 components and 0.25146 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038745Jy sigma=0.157412 ! Fit after self-cal, rms=0.038743Jy sigma=0.157403 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038743Jy sigma=0.157403 ! Fit after self-cal, rms=0.038699Jy sigma=0.157376 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.616 mas, bmaj=10.57 mas, bpa=-35.65 degrees ! Estimated noise=2.66694 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000420148 Jy ! Component: 100 - total flux cleaned = 0.000728829 Jy ! Component: 150 - total flux cleaned = 0.00106498 Jy ! Component: 200 - total flux cleaned = 0.00140636 Jy ! Component: 250 - total flux cleaned = 0.00167665 Jy ! Component: 300 - total flux cleaned = 0.00188232 Jy ! Component: 350 - total flux cleaned = 0.00206367 Jy ! Component: 400 - total flux cleaned = 0.00231108 Jy ! Total flux subtracted in 400 components = 0.00231108 Jy ! Clean residual min=-0.001084 max=0.001040 Jy/beam ! Clean residual mean=0.000004 rms=0.000251 Jy/beam ! Combined flux in latest and established models = 0.253771 Jy selfcal ! Performing phase self-cal ! Adding 84 model components to the UV plane model. ! The established model now contains 656 components and 0.253771 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038660Jy sigma=0.157222 ! Fit after self-cal, rms=0.038660Jy sigma=0.157218 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000203799 Jy ! Component: 100 - total flux cleaned = 0.000343785 Jy ! Component: 150 - total flux cleaned = 0.000514498 Jy ! Component: 200 - total flux cleaned = 0.000665253 Jy ! Component: 250 - total flux cleaned = 0.000795691 Jy ! Component: 300 - total flux cleaned = 0.000972121 Jy ! Component: 350 - total flux cleaned = 0.00112899 Jy ! Component: 400 - total flux cleaned = 0.00126758 Jy ! Total flux subtracted in 400 components = 0.00126758 Jy ! Clean residual min=-0.001056 max=0.001016 Jy/beam ! Clean residual mean=0.000004 rms=0.000242 Jy/beam ! Combined flux in latest and established models = 0.255039 Jy ! Performing phase self-cal ! Adding 86 model components to the UV plane model. ! The established model now contains 684 components and 0.255039 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038632Jy sigma=0.157120 ! Fit after self-cal, rms=0.038632Jy sigma=0.157118 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038632Jy sigma=0.157118 ! Fit after self-cal, rms=0.038313Jy sigma=0.154416 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.217 mas, bmaj=12.69 mas, bpa=-41.74 degrees ! Estimated noise=2.98875 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000539873 Jy ! Component: 100 - total flux cleaned = 0.000900541 Jy ! Component: 150 - total flux cleaned = 0.00124352 Jy ! Component: 200 - total flux cleaned = 0.00156115 Jy ! Component: 250 - total flux cleaned = 0.00188857 Jy ! Component: 300 - total flux cleaned = 0.00209998 Jy ! Component: 350 - total flux cleaned = 0.00232939 Jy ! Component: 400 - total flux cleaned = 0.00253712 Jy ! Total flux subtracted in 400 components = 0.00253712 Jy ! Clean residual min=-0.000811 max=0.000967 Jy/beam ! Clean residual mean=0.000004 rms=0.000220 Jy/beam ! Combined flux in latest and established models = 0.257576 Jy selfcal ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 750 components and 0.257576 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038301Jy sigma=0.154349 ! Fit after self-cal, rms=0.038305Jy sigma=0.154341 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 (88, -73). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -9.83056e-05 Jy ! Component: 100 - total flux cleaned = -0.00012876 Jy ! Component: 150 - total flux cleaned = -0.000152457 Jy ! Component: 200 - total flux cleaned = -9.47411e-05 Jy ! Component: 250 - total flux cleaned = -9.26605e-06 Jy ! Component: 300 - total flux cleaned = 2.20694e-05 Jy ! Component: 350 - total flux cleaned = 9.69479e-05 Jy ! Component: 400 - total flux cleaned = 0.000182467 Jy ! Total flux subtracted in 400 components = 0.000182467 Jy ! Clean residual min=-0.000662 max=0.000730 Jy/beam ! Clean residual mean=0.000003 rms=0.000173 Jy/beam ! Combined flux in latest and established models = 0.257758 Jy ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 820 components and 0.257758 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038272Jy sigma=0.154164 ! Fit after self-cal, rms=0.038270Jy sigma=0.154158 ! Inverting map ! Added new window around map position (-119.5, -46.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000544481 Jy ! Component: 100 - total flux cleaned = 0.000565307 Jy ! Component: 150 - total flux cleaned = 0.000581813 Jy ! Component: 200 - total flux cleaned = 0.000565905 Jy ! Component: 250 - total flux cleaned = 0.000579493 Jy ! Component: 300 - total flux cleaned = 0.00060729 Jy ! Component: 350 - total flux cleaned = 0.000620808 Jy ! Component: 400 - total flux cleaned = 0.000647544 Jy ! Total flux subtracted in 400 components = 0.000647544 Jy ! Clean residual min=-0.000658 max=0.000553 Jy/beam ! Clean residual mean=0.000001 rms=0.000145 Jy/beam ! Combined flux in latest and established models = 0.258406 Jy ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 889 components and 0.258406 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038251Jy sigma=0.154051 ! Fit after self-cal, rms=0.038250Jy sigma=0.154045 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=3.618 mas, bmaj=10.57 mas, bpa=-35.67 degrees ! Estimated noise=2.66906 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.618 x 10.57 at -35.67 degrees (North through East) ! Clean map min=-0.0018548 max=0.15712 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.782 mas, bmaj=10.05 mas, bpa=-36.22 degrees ! Estimated noise=5.36932 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.117252 Jy ! Component: 100 - total flux cleaned = 0.145322 Jy ! Component: 150 - total flux cleaned = 0.157146 Jy ! Component: 200 - total flux cleaned = 0.165103 Jy ! Component: 250 - total flux cleaned = 0.1712 Jy ! Component: 300 - total flux cleaned = 0.175806 Jy ! Component: 350 - total flux cleaned = 0.179731 Jy ! Component: 400 - total flux cleaned = 0.182771 Jy ! Total flux subtracted in 400 components = 0.182771 Jy ! Clean residual min=-0.005078 max=0.004287 Jy/beam ! Clean residual mean=-0.000005 rms=0.001196 Jy/beam ! Combined flux in latest and established models = 0.182771 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 70 model components to the UV plane model. ! The established model now contains 70 components and 0.182771 Jy ! Inverting map and beam ! Estimated beam: bmin=3.618 mas, bmaj=10.57 mas, bpa=-35.67 degrees ! Estimated noise=2.66906 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00772261 Jy ! Component: 100 - total flux cleaned = 0.0133196 Jy ! Component: 150 - total flux cleaned = 0.0180286 Jy ! Component: 200 - total flux cleaned = 0.0217566 Jy ! Component: 250 - total flux cleaned = 0.0250316 Jy ! Component: 300 - total flux cleaned = 0.0278221 Jy ! Component: 350 - total flux cleaned = 0.0298793 Jy ! Component: 400 - total flux cleaned = 0.0316641 Jy ! Total flux subtracted in 400 components = 0.0316641 Jy ! Clean residual min=-0.002982 max=0.002319 Jy/beam ! Clean residual mean=0.000002 rms=0.000791 Jy/beam ! Combined flux in latest and established models = 0.214435 Jy ! Adding 166 model components to the UV plane model. ! The established model now contains 220 components and 0.214435 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.217 mas, bmaj=12.69 mas, bpa=-41.74 degrees ! Estimated noise=2.98875 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00395047 Jy ! Component: 100 - total flux cleaned = 0.00653563 Jy ! Component: 150 - total flux cleaned = 0.00892321 Jy ! Component: 200 - total flux cleaned = 0.0111513 Jy ! Component: 250 - total flux cleaned = 0.013058 Jy ! Component: 300 - total flux cleaned = 0.0148613 Jy ! Component: 350 - total flux cleaned = 0.0162638 Jy ! Component: 400 - total flux cleaned = 0.0175454 Jy ! Total flux subtracted in 400 components = 0.0175454 Jy ! Clean residual min=-0.002127 max=0.001868 Jy/beam ! Clean residual mean=0.000000 rms=0.000648 Jy/beam ! Combined flux in latest and established models = 0.231981 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00111785 Jy ! Component: 100 - total flux cleaned = 0.00215028 Jy ! Component: 150 - total flux cleaned = 0.00310605 Jy ! Component: 200 - total flux cleaned = 0.00394537 Jy ! Component: 250 - total flux cleaned = 0.00477382 Jy ! Component: 300 - total flux cleaned = 0.00541544 Jy ! Component: 350 - total flux cleaned = 0.00610839 Jy ! Component: 400 - total flux cleaned = 0.00684051 Jy ! Total flux subtracted in 400 components = 0.00684051 Jy ! Clean residual min=-0.001391 max=0.001301 Jy/beam ! Clean residual mean=-0.000001 rms=0.000412 Jy/beam ! Combined flux in latest and established models = 0.238821 Jy ! Adding 300 model components to the UV plane model. ! The established model now contains 492 components and 0.238821 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.618 mas, bmaj=10.57 mas, bpa=-35.67 degrees ! Estimated noise=2.66906 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000223349 Jy ! Component: 100 - total flux cleaned = 0.000793367 Jy ! Component: 150 - total flux cleaned = 0.00128337 Jy ! Component: 200 - total flux cleaned = 0.00178237 Jy ! Component: 250 - total flux cleaned = 0.00211448 Jy ! Component: 300 - total flux cleaned = 0.00246776 Jy ! Component: 350 - total flux cleaned = 0.00280757 Jy ! Component: 400 - total flux cleaned = 0.00316896 Jy ! Total flux subtracted in 400 components = 0.00316896 Jy ! Clean residual min=-0.001197 max=0.001160 Jy/beam ! Clean residual mean=-0.000002 rms=0.000289 Jy/beam ! Combined flux in latest and established models = 0.24199 Jy selfcal ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 562 components and 0.24199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.038489Jy sigma=0.154964 ! Fit after self-cal, rms=0.038471Jy sigma=0.154889 wmodel /scr/pima/rdv94_uvs/J1321+8316_S_map.mod ! Writing 562 model components to file: /scr/pima/rdv94_uvs/J1321+8316_S_map.mod wobs /scr/pima/rdv94_uvs/J1321+8316_S_uvs.fits ! Applying 2016 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J1321+8316_S_uvs.fits wwins /scr/pima/rdv94_uvs/J1321+8316_S_map.win ! wwins: Wrote 8 windows to /scr/pima/rdv94_uvs/J1321+8316_S_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.000116324 Jy ! Component: 100 - total flux cleaned = 1.62904e-05 Jy ! Component: 150 - total flux cleaned = 0.000105259 Jy ! Component: 200 - total flux cleaned = 0.000105185 Jy ! Component: 250 - total flux cleaned = 0.000153275 Jy ! Component: 300 - total flux cleaned = 0.00015346 Jy ! Component: 350 - total flux cleaned = 0.000168614 Jy ! Component: 400 - total flux cleaned = 0.000154116 Jy ! Component: 450 - total flux cleaned = 0.000168727 Jy ! Component: 500 - total flux cleaned = 0.00018274 Jy ! Component: 550 - total flux cleaned = 0.000196341 Jy ! Component: 600 - total flux cleaned = 0.000209902 Jy ! Component: 650 - total flux cleaned = 0.000170114 Jy ! Component: 700 - total flux cleaned = 0.000143951 Jy ! Component: 750 - total flux cleaned = 5.39788e-05 Jy ! Component: 800 - total flux cleaned = 3.25186e-06 Jy ! Component: 850 - total flux cleaned = -0.000146588 Jy ! Component: 900 - total flux cleaned = -0.000232726 Jy ! Component: 950 - total flux cleaned = -0.000318152 Jy ! Component: 1000 - total flux cleaned = -0.000426763 Jy ! Component: 1050 - total flux cleaned = -0.000546323 Jy ! Component: 1100 - total flux cleaned = -0.000664641 Jy ! Component: 1150 - total flux cleaned = -0.000805081 Jy ! Component: 1200 - total flux cleaned = -0.000920945 Jy ! Component: 1250 - total flux cleaned = -0.00104711 Jy ! Component: 1300 - total flux cleaned = -0.00111528 Jy ! Component: 1350 - total flux cleaned = -0.00125035 Jy ! Component: 1400 - total flux cleaned = -0.00137297 Jy ! Component: 1450 - total flux cleaned = -0.00142819 Jy ! Component: 1500 - total flux cleaned = -0.00150512 Jy ! Component: 1550 - total flux cleaned = -0.00163567 Jy ! Component: 1600 - total flux cleaned = -0.00171118 Jy ! Component: 1650 - total flux cleaned = -0.00179668 Jy ! Component: 1700 - total flux cleaned = -0.00184974 Jy ! Component: 1750 - total flux cleaned = -0.00190243 Jy ! Component: 1800 - total flux cleaned = -0.00195476 Jy ! Component: 1850 - total flux cleaned = -0.00204814 Jy ! Component: 1900 - total flux cleaned = -0.002079 Jy ! Component: 1950 - total flux cleaned = -0.00210966 Jy ! Component: 2000 - total flux cleaned = -0.00215032 Jy ! Component: 2050 - total flux cleaned = -0.00220068 Jy ! Component: 2100 - total flux cleaned = -0.00224068 Jy ! Component: 2150 - total flux cleaned = -0.00224063 Jy ! Component: 2200 - total flux cleaned = -0.00230994 Jy ! Component: 2250 - total flux cleaned = -0.00234926 Jy ! Component: 2300 - total flux cleaned = -0.00236894 Jy ! Component: 2350 - total flux cleaned = -0.00247576 Jy ! Component: 2400 - total flux cleaned = -0.00249502 Jy ! Component: 2450 - total flux cleaned = -0.00253343 Jy ! Component: 2500 - total flux cleaned = -0.00259059 Jy ! Component: 2550 - total flux cleaned = -0.00260949 Jy ! Component: 2600 - total flux cleaned = -0.00265662 Jy ! Component: 2650 - total flux cleaned = -0.00266597 Jy ! Component: 2700 - total flux cleaned = -0.00270329 Jy ! Component: 2750 - total flux cleaned = -0.00271262 Jy ! Component: 2800 - total flux cleaned = -0.00282311 Jy ! Component: 2850 - total flux cleaned = -0.00282304 Jy ! Component: 2900 - total flux cleaned = -0.0028595 Jy ! Component: 2950 - total flux cleaned = -0.00288665 Jy ! Component: 3000 - total flux cleaned = -0.00295873 Jy ! Component: 3050 - total flux cleaned = -0.00297663 Jy ! Component: 3100 - total flux cleaned = -0.00302119 Jy ! Component: 3150 - total flux cleaned = -0.00308323 Jy ! Component: 3200 - total flux cleaned = -0.00311845 Jy ! Component: 3250 - total flux cleaned = -0.00314477 Jy ! Component: 3300 - total flux cleaned = -0.00318845 Jy ! Component: 3350 - total flux cleaned = -0.00324925 Jy ! Component: 3400 - total flux cleaned = -0.00324924 Jy ! Component: 3450 - total flux cleaned = -0.003318 Jy ! Component: 3500 - total flux cleaned = -0.00334363 Jy ! Component: 3550 - total flux cleaned = -0.00341173 Jy ! Component: 3600 - total flux cleaned = -0.00342022 Jy ! Component: 3650 - total flux cleaned = -0.00347091 Jy ! Component: 3700 - total flux cleaned = -0.00351285 Jy ! Component: 3750 - total flux cleaned = -0.00352957 Jy ! Component: 3800 - total flux cleaned = -0.00356285 Jy ! Component: 3850 - total flux cleaned = -0.00358777 Jy ! Component: 3900 - total flux cleaned = -0.00360424 Jy ! Component: 3950 - total flux cleaned = -0.00362063 Jy ! Component: 4000 - total flux cleaned = -0.00363706 Jy ! Component: 4050 - total flux cleaned = -0.00366162 Jy ! Total flux subtracted in 4096 components = -0.00369414 Jy ! Clean residual min=-0.000521 max=0.000485 Jy/beam ! Clean residual mean=-0.000000 rms=0.000149 Jy/beam ! Combined flux in latest and established models = 0.238296 Jy keep ! Adding 262 model components to the UV plane model. ! The established model now contains 820 components and 0.238296 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J1321+8316_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.618 x 10.57 at -35.67 degrees (North through East) ! Clean map min=-0.0015649 max=0.15533 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J1321+8316_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.1553 rms= 0.00014 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J1321+8316_S_uvs.log closed on Tue Oct 8 12:20:18 2013