! Started logfile: /scr/pima/rdv93_uvs/J0914+0245_S_uvs.log on Tue Oct 8 15:27:20 2013 obs /scr/pima/rdv93_uvs/J0914+0245_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv93_uvs/J0914+0245_S_uva.fits ! AN table 1: 21 integrations on 55 of 55 possible baselines. ! AN table 2: 2 integrations on 3 of 3 possible baselines. ! AN table 3: 3 integrations on 6 of 6 possible baselines. ! AN table 4: 2 integrations on 1 of 1 possible baselines. ! AN table 5: 3 integrations on 21 of 21 possible baselines. ! AN table 6: 2 integrations on 6 of 6 possible baselines. ! AN table 7: 8 integrations on 1 of 1 possible baselines. ! AN table 8: 26 integrations on 45 of 45 possible baselines. ! AN table 9: 26 integrations on 55 of 55 possible baselines. ! AN table 10: 27 integrations on 78 of 78 possible baselines. ! Apparent sampling: 0.71206 visibilities/baseline/integration-bin. ! Found source: J0914+0245 ! ! 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 17004 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.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv93_uvs/J0914+0245_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:FT 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:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:FT self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:MA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:FT self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:MA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:FT self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MA self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:SC self-cal status=correctable weight=1 ! 10:WZ self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 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 53 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 78 telescope corrections were flagged in sub-array 8. ! A total of 98 telescope corrections were flagged in sub-array 9. ! A total of 126 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 2. ! A total of 53 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 101 telescope corrections were flagged in sub-array 8. ! A total of 118 telescope corrections were flagged in sub-array 9. ! A total of 147 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 3. ! A total of 53 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 78 telescope corrections were flagged in sub-array 8. ! A total of 98 telescope corrections were flagged in sub-array 9. ! A total of 126 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 4. ! A total of 53 telescope corrections were flagged in sub-array 1. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 101 telescope corrections were flagged in sub-array 8. ! A total of 118 telescope corrections were flagged in sub-array 9. ! A total of 147 telescope corrections were flagged in sub-array 10. ! ! Fit before self-cal, rms=1.217446Jy sigma=10.971722 ! Fit after self-cal, rms=0.438180Jy sigma=3.543709 ! 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 \] flag 1:FT ! Flagging baselines of antenna 1:FT in the currently selected channels. flag 1:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 7:NY ! Flagging baselines of antenna 7:NY in the currently selected channels. flag 7:WZ ! Flagging baselines of antenna 7:WZ in the currently selected channels. flag 8:FT ! Flagging baselines of antenna 8:FT in the currently selected channels. flag 8:MA ! Flagging baselines of antenna 8:MA in the currently selected channels. flag 8:NY ! Flagging baselines of antenna 8:NY in the currently selected channels. flag 8:WZ ! Flagging baselines of antenna 8:WZ in the currently selected channels. flag 9:FT ! Flagging baselines of antenna 9:FT in the currently selected channels. flag 9:MA ! Flagging baselines of antenna 9:MA in the currently selected channels. flag 9:NY ! Flagging baselines of antenna 9:NY in the currently selected channels. flag 9:WZ ! Flagging baselines of antenna 9:WZ in the currently selected channels. flag 10:FT ! Flagging baselines of antenna 10:FT in the currently selected channels. flag 10:MA ! Flagging baselines of antenna 10:MA in the currently selected channels. flag 10:NY ! Flagging baselines of antenna 10:NY in the currently selected channels. flag 10:WZ ! Flagging baselines of antenna 10:WZ in the currently selected channels. ![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=4.326 mas, bmaj=6.607 mas, bpa=0.4689 degrees ! Estimated noise=1.44637 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.48295 Jy ! Component: 100 - total flux cleaned = 0.590223 Jy ! Total flux subtracted in 100 components = 0.590223 Jy ! Clean residual min=-0.018592 max=0.035801 Jy/beam ! Clean residual mean=0.000022 rms=0.004563 Jy/beam ! Combined flux in latest and established models = 0.590223 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.590223 Jy ! ! Correcting IF 1. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 63 telescope corrections were flagged in sub-array 8. ! A total of 60 telescope corrections were flagged in sub-array 9. ! A total of 66 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 40 telescope corrections were flagged in sub-array 9. ! A total of 45 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 63 telescope corrections were flagged in sub-array 8. ! A total of 60 telescope corrections were flagged in sub-array 9. ! A total of 66 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 40 telescope corrections were flagged in sub-array 9. ! A total of 45 telescope corrections were flagged in sub-array 10. ! ! Fit before self-cal, rms=0.124556Jy sigma=1.368578 ! Fit after self-cal, rms=0.124305Jy sigma=1.366710 ! Inverting map and beam ! Estimated beam: bmin=4.326 mas, bmaj=6.607 mas, bpa=0.4689 degrees ! Estimated noise=1.4464 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.037779 Jy ! Component: 100 - total flux cleaned = 0.057032 Jy ! Total flux subtracted in 100 components = 0.057032 Jy ! Clean residual min=-0.016218 max=0.017121 Jy/beam ! Clean residual mean=0.000017 rms=0.003611 Jy/beam ! Combined flux in latest and established models = 0.647255 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 22 components and 0.647255 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113408Jy sigma=1.202885 ! Fit after self-cal, rms=0.113344Jy sigma=1.202182 ! 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=5.824 mas, bmaj=8.768 mas, bpa=5.2 degrees ! Estimated noise=0.75873 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.0124091 Jy ! Component: 100 - total flux cleaned = 0.0148932 Jy ! Component: 150 - total flux cleaned = 0.0147306 Jy ! Component: 200 - total flux cleaned = 0.0144295 Jy ! Total flux subtracted in 200 components = 0.0144295 Jy ! Clean residual min=-0.012656 max=0.015019 Jy/beam ! Clean residual mean=-0.000015 rms=0.003277 Jy/beam ! Combined flux in latest and established models = 0.661684 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 37 components and 0.661684 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112546Jy sigma=1.181728 ! Fit after self-cal, rms=0.112502Jy sigma=1.181308 ! 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. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 0.99 FT 1.00* HN 0.92 ! KP 1.01 LA 1.12 NL 0.99 OV 0.98 ! PT 1.03 SC 0.97 WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.07 MK 0.87 NL 1.05 OV 1.03 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.00 LA 1.24 MK 0.89 ! NL 1.41 OV 0.87 PT 1.02 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.08 LA 1.20 OV 1.00 PT 1.06 ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* WZ 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.04 FT 1.00* HN 0.94 LA 1.13 ! MA 1.00* NL 1.03 NY 1.00* PT 1.08 ! SC 0.97 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! FD 1.02 FT 1.00* HN 0.93 KP 1.04 ! LA 1.14 MA 1.00* NL 0.99 NY 1.00* ! PT 1.04 SC 0.96 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.99 FD 1.02 FT 1.00* HN 0.92 ! KP 1.05 LA 1.12 MA 1.00* NL 1.00 ! NY 1.00* OV 0.99 PT 1.03 SC 0.98 ! WZ 1.00* ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.08 FT 1.00* HN 1.04 ! KP 1.00 LA 1.02 NL 1.02 OV 1.00 ! PT 1.03 SC 1.05 WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.10 MK 1.04 NL 1.07 OV 1.03 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.09 FD 1.14 LA 1.05 MK 1.05 ! NL 1.32 OV 1.01 PT 1.05 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.07 LA 1.10 OV 1.04 PT 1.14 ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* WZ 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.13 FT 1.00* HN 1.06 LA 1.06 ! MA 1.00* NL 1.07 NY 1.00* PT 1.09 ! SC 1.04 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! FD 1.10 FT 1.00* HN 1.07 KP 1.05 ! LA 1.04 MA 1.00* NL 1.03 NY 1.00* ! PT 1.06 SC 1.04 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.05 FD 1.09 FT 1.00* HN 1.07 ! KP 1.03 LA 1.02 MA 1.00* NL 1.03 ! NY 1.00* OV 1.03 PT 1.06 SC 1.05 ! WZ 1.00* ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.09 FD 1.08 FT 1.00* HN 0.95 ! KP 0.88 LA 0.96 NL 1.09 OV 1.02 ! PT 1.06 SC 0.95 WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.15 MK 0.98 NL 1.18 OV 1.02 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.11 FD 1.16 LA 1.02 MK 1.01 ! NL 1.44 OV 1.02 PT 1.09 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.10 LA 1.04 OV 1.01 PT 1.12 ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* WZ 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.14 FT 1.00* HN 0.96 LA 1.00 ! MA 1.00* NL 1.12 NY 1.00* PT 1.13 ! SC 0.95 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! FD 1.10 FT 1.00* HN 0.96 KP 0.92 ! LA 0.99 MA 1.00* NL 1.10 NY 1.00* ! PT 1.09 SC 0.94 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.00 FD 1.09 FT 1.00* HN 0.97 ! KP 0.91 LA 0.97 MA 1.00* NL 1.12 ! NY 1.00* OV 1.01 PT 1.07 SC 0.95 ! WZ 1.00* ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.90 FD 1.05 FT 1.00* HN 0.93 ! KP 0.82 LA 0.80 NL 0.99 OV 0.86 ! PT 0.95 SC 0.99 WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! FD 1.08 MK 1.11 NL 1.07 OV 0.90 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.12 FD 1.11 LA 0.84 MK 1.11 ! NL 1.26 OV 0.89 PT 0.92 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 LA 0.88 OV 0.92 PT 1.01 ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* WZ 0.00* ! ! A total of 12 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.13 FT 1.00* HN 0.94 LA 0.85 ! MA 1.00* NL 1.04 NY 1.00* PT 0.98 ! SC 0.98 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! FD 1.09 FT 1.00* HN 0.93 KP 0.84 ! LA 0.83 MA 1.00* NL 1.00 NY 1.00* ! PT 0.97 SC 1.00 WZ 1.00* ! ! A total of 3 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.92 FD 1.07 FT 1.00* HN 0.94 ! KP 0.83 LA 0.82 MA 1.00* NL 1.01 ! NY 1.00* OV 0.90 PT 0.96 SC 1.01 ! WZ 1.00* ! ! ! Fit before self-cal, rms=0.112502Jy sigma=1.181308 ! Fit after self-cal, rms=0.080408Jy sigma=0.723310 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.75 mas, bmaj=8.766 mas, bpa=5.796 degrees ! Estimated noise=0.76039 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 (2.4, 10.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0113699 Jy ! Component: 100 - total flux cleaned = 0.013182 Jy ! Component: 150 - total flux cleaned = 0.0130979 Jy ! Component: 200 - total flux cleaned = 0.0128886 Jy ! Total flux subtracted in 200 components = 0.0128886 Jy ! Clean residual min=-0.002018 max=0.002539 Jy/beam ! Clean residual mean=-0.000010 rms=0.000603 Jy/beam ! Combined flux in latest and established models = 0.674573 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 50 components and 0.674573 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079028Jy sigma=0.704032 ! Fit after self-cal, rms=0.079040Jy sigma=0.702994 ! 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.079040Jy sigma=0.702994 ! Fit after self-cal, rms=0.079200Jy sigma=0.701912 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.745 mas, bmaj=8.757 mas, bpa=5.829 degrees ! Estimated noise=0.759241 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0023268 Jy ! Component: 100 - total flux cleaned = 0.002858 Jy ! Component: 150 - total flux cleaned = 0.00293224 Jy ! Component: 200 - total flux cleaned = 0.00300162 Jy ! Total flux subtracted in 200 components = 0.00300162 Jy ! Clean residual min=-0.001914 max=0.002512 Jy/beam ! Clean residual mean=-0.000011 rms=0.000535 Jy/beam ! Combined flux in latest and established models = 0.677574 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 72 components and 0.677574 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079134Jy sigma=0.700879 ! Fit after self-cal, rms=0.079122Jy sigma=0.700819 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 (-36.8, -51.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00217256 Jy ! Component: 100 - total flux cleaned = 0.00286438 Jy ! Component: 150 - total flux cleaned = 0.00310331 Jy ! Component: 200 - total flux cleaned = 0.00334947 Jy ! Total flux subtracted in 200 components = 0.00334947 Jy ! Clean residual min=-0.001729 max=0.001825 Jy/beam ! Clean residual mean=-0.000008 rms=0.000482 Jy/beam ! Combined flux in latest and established models = 0.680924 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 100 components and 0.680924 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079052Jy sigma=0.699864 ! Fit after self-cal, rms=0.079049Jy sigma=0.699809 ! 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.079049Jy sigma=0.699809 ! Fit after self-cal, rms=0.079129Jy sigma=0.699591 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.744 mas, bmaj=8.754 mas, bpa=5.817 degrees ! Estimated noise=0.759309 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00107655 Jy ! Component: 100 - total flux cleaned = 0.00148192 Jy ! Component: 150 - total flux cleaned = 0.00169013 Jy ! Component: 200 - total flux cleaned = 0.00185445 Jy ! Total flux subtracted in 200 components = 0.00185445 Jy ! Clean residual min=-0.001671 max=0.001763 Jy/beam ! Clean residual mean=-0.000008 rms=0.000466 Jy/beam ! Combined flux in latest and established models = 0.682778 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 120 components and 0.682778 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079102Jy sigma=0.699259 ! Fit after self-cal, rms=0.079102Jy sigma=0.699246 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.000315329 Jy ! Component: 100 - total flux cleaned = 0.000391529 Jy ! Component: 150 - total flux cleaned = 0.000416299 Jy ! Component: 200 - total flux cleaned = 0.000462464 Jy ! Total flux subtracted in 200 components = 0.000462464 Jy ! Clean residual min=-0.001649 max=0.001765 Jy/beam ! Clean residual mean=-0.000008 rms=0.000461 Jy/beam ! Combined flux in latest and established models = 0.683241 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 136 components and 0.683241 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079089Jy sigma=0.699072 ! Fit after self-cal, rms=0.079089Jy sigma=0.699069 ! 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.079089Jy sigma=0.699069 ! Fit after self-cal, rms=0.079134Jy sigma=0.698912 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.705 mas, bmaj=11.98 mas, bpa=-14.55 degrees ! Estimated noise=0.84244 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000577347 Jy ! Component: 100 - total flux cleaned = 0.000902458 Jy ! Component: 150 - total flux cleaned = 0.00113822 Jy ! Component: 200 - total flux cleaned = 0.00120294 Jy ! Total flux subtracted in 200 components = 0.00120294 Jy ! Clean residual min=-0.001532 max=0.001711 Jy/beam ! Clean residual mean=-0.000013 rms=0.000471 Jy/beam ! Combined flux in latest and established models = 0.684444 Jy selfcal ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 163 components and 0.684444 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079125Jy sigma=0.698798 ! Fit after self-cal, rms=0.079127Jy sigma=0.698793 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.000215957 Jy ! Component: 100 - total flux cleaned = 0.000258638 Jy ! Component: 150 - total flux cleaned = 0.00031218 Jy ! Component: 200 - total flux cleaned = 0.000350203 Jy ! Total flux subtracted in 200 components = 0.000350203 Jy ! Clean residual min=-0.001497 max=0.001708 Jy/beam ! Clean residual mean=-0.000013 rms=0.000468 Jy/beam ! Combined flux in latest and established models = 0.684794 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 176 components and 0.684794 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079121Jy sigma=0.698743 ! Fit after self-cal, rms=0.079122Jy sigma=0.698742 ! 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. ! Applying 2560 buffered edits. ! Applying 2560 buffered edits. 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:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:OV,true ! 6:OV self-cal status=fixed weight=1 selfant 6:PT,true ! 6:PT self-cal status=fixed weight=1 selfant 8:FD,true ! 8:FD self-cal status=fixed weight=1 selfant 8:HN,true ! 8:HN self-cal status=fixed weight=1 selfant 8:LA,true ! 8:LA self-cal status=fixed weight=1 selfant 8:NL,true ! 8:NL self-cal status=fixed weight=1 selfant 8:PT,true ! 8:PT self-cal status=fixed weight=1 selfant 8:SC,true ! 8:SC self-cal status=fixed weight=1 selfant 9:FD,true ! 9:FD self-cal status=fixed weight=1 selfant 9:HN,true ! 9:HN self-cal status=fixed weight=1 selfant 9:KP,true ! 9:KP self-cal status=fixed weight=1 selfant 9:LA,true ! 9:LA self-cal status=fixed weight=1 selfant 9:NL,true ! 9:NL self-cal status=fixed weight=1 selfant 9:PT,true ! 9:PT self-cal status=fixed weight=1 selfant 9:SC,true ! 9:SC self-cal status=fixed weight=1 selfant 10:BR,true ! 10:BR self-cal status=fixed weight=1 selfant 10:FD,true ! 10:FD self-cal status=fixed weight=1 selfant 10:HN,true ! 10:HN self-cal status=fixed weight=1 selfant 10:KP,true ! 10:KP self-cal status=fixed weight=1 selfant 10:LA,true ! 10:LA self-cal status=fixed weight=1 selfant 10:NL,true ! 10:NL self-cal status=fixed weight=1 selfant 10:OV,true ! 10:OV self-cal status=fixed weight=1 selfant 10:PT,true ! 10:PT self-cal status=fixed weight=1 selfant 10:SC,true ! 10:SC 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. ! A total of 26 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 41 telescope corrections were flagged in sub-array 8. ! A total of 44 telescope corrections were flagged in sub-array 9. ! A total of 46 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 2. ! A total of 26 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 64 telescope corrections were flagged in sub-array 8. ! A total of 64 telescope corrections were flagged in sub-array 9. ! A total of 68 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 3. ! A total of 26 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 41 telescope corrections were flagged in sub-array 8. ! A total of 44 telescope corrections were flagged in sub-array 9. ! A total of 46 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 4. ! A total of 26 telescope corrections were flagged in sub-array 1. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 64 telescope corrections were flagged in sub-array 8. ! A total of 64 telescope corrections were flagged in sub-array 9. ! A total of 68 telescope corrections were flagged in sub-array 10. ! ! Fit before self-cal, rms=0.235715Jy sigma=1.062031 ! Fit after self-cal, rms=0.126115Jy sigma=0.770155 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:FT 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:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:FT self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:MA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:FT self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:MA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:FT self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MA self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:SC self-cal status=correctable weight=1 ! 10:WZ self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 37 telescope corrections were flagged in sub-array 8. ! A total of 54 telescope corrections were flagged in sub-array 9. ! A total of 80 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 2. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 46 telescope corrections were flagged in sub-array 8. ! A total of 54 telescope corrections were flagged in sub-array 9. ! A total of 82 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 3. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 37 telescope corrections were flagged in sub-array 8. ! A total of 54 telescope corrections were flagged in sub-array 9. ! A total of 80 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 4. ! A total of 30 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 46 telescope corrections were flagged in sub-array 8. ! A total of 54 telescope corrections were flagged in sub-array 9. ! A total of 82 telescope corrections were flagged in sub-array 10. ! ! Fit before self-cal, rms=0.126115Jy sigma=0.770155 ! Fit after self-cal, rms=0.124704Jy sigma=0.740164 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.107 mas, bmaj=6.714 mas, bpa=-8.351 degrees ! Estimated noise=1.23124 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.00115118 Jy ! Component: 100 - total flux cleaned = -0.00087971 Jy ! Component: 150 - total flux cleaned = -0.000900371 Jy ! Component: 200 - total flux cleaned = -0.000490011 Jy ! Total flux subtracted in 200 components = -0.000490011 Jy ! Clean residual min=-0.004007 max=0.004033 Jy/beam ! Clean residual mean=-0.000002 rms=0.001009 Jy/beam ! Combined flux in latest and established models = 0.684304 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 196 components and 0.684304 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123155Jy sigma=0.729255 ! Fit after self-cal, rms=0.122708Jy sigma=0.724173 ! 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=2.42 mas, bmaj=8.224 mas, bpa=-9.04 degrees ! Estimated noise=0.710898 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.0021267 Jy ! Component: 100 - total flux cleaned = -0.00212481 Jy ! Component: 150 - total flux cleaned = -0.0020438 Jy ! Component: 200 - total flux cleaned = -0.00196935 Jy ! Component: 250 - total flux cleaned = -0.0020392 Jy ! Component: 300 - total flux cleaned = -0.00216938 Jy ! Component: 350 - total flux cleaned = -0.00223008 Jy ! Component: 400 - total flux cleaned = -0.00234691 Jy ! Total flux subtracted in 400 components = -0.00234691 Jy ! Clean residual min=-0.002030 max=0.002273 Jy/beam ! Clean residual mean=-0.000006 rms=0.000519 Jy/beam ! Combined flux in latest and established models = 0.681957 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 234 components and 0.681957 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122408Jy sigma=0.722614 ! Fit after self-cal, rms=0.122400Jy sigma=0.722554 ! Inverting map ! Added new window around map position (60.8, 55.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000909734 Jy ! Component: 100 - total flux cleaned = 0.00092136 Jy ! Component: 150 - total flux cleaned = 0.000742334 Jy ! Component: 200 - total flux cleaned = 0.000631155 Jy ! Component: 250 - total flux cleaned = 0.000524202 Jy ! Component: 300 - total flux cleaned = 0.000524432 Jy ! Component: 350 - total flux cleaned = 0.000473886 Jy ! Component: 400 - total flux cleaned = 0.000475518 Jy ! Total flux subtracted in 400 components = 0.000475518 Jy ! Clean residual min=-0.001898 max=0.002020 Jy/beam ! Clean residual mean=-0.000005 rms=0.000493 Jy/beam ! Combined flux in latest and established models = 0.682433 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 268 components and 0.682432 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122225Jy sigma=0.721545 ! Fit after self-cal, rms=0.122219Jy sigma=0.721505 ! Inverting map ! Added new window around map position (60, -49.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00072451 Jy ! Component: 100 - total flux cleaned = 0.000601247 Jy ! Component: 150 - total flux cleaned = 0.00080926 Jy ! Component: 200 - total flux cleaned = 0.000809429 Jy ! Component: 250 - total flux cleaned = 0.000808676 Jy ! Component: 300 - total flux cleaned = 0.000808956 Jy ! Component: 350 - total flux cleaned = 0.000766693 Jy ! Component: 400 - total flux cleaned = 0.000766362 Jy ! Total flux subtracted in 400 components = 0.000766362 Jy ! Clean residual min=-0.001914 max=0.001798 Jy/beam ! Clean residual mean=-0.000005 rms=0.000470 Jy/beam ! Combined flux in latest and established models = 0.683199 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 304 components and 0.683199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122078Jy sigma=0.720716 ! Fit after self-cal, rms=0.122071Jy sigma=0.720684 ! 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.000128549 Jy ! Component: 100 - total flux cleaned = 4.36633e-05 Jy ! Component: 150 - total flux cleaned = 4.35946e-05 Jy ! Component: 200 - total flux cleaned = 3.05292e-06 Jy ! Component: 250 - total flux cleaned = -3.60007e-05 Jy ! Component: 300 - total flux cleaned = -7.37667e-05 Jy ! Component: 350 - total flux cleaned = -0.000111296 Jy ! Component: 400 - total flux cleaned = -0.000110746 Jy ! Total flux subtracted in 400 components = -0.000110746 Jy ! Clean residual min=-0.001860 max=0.001757 Jy/beam ! Clean residual mean=-0.000005 rms=0.000459 Jy/beam ! Combined flux in latest and established models = 0.683088 Jy ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 318 components and 0.683088 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121989Jy sigma=0.720252 ! Fit after self-cal, rms=0.121986Jy sigma=0.720245 ! 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. ! Applying 2560 buffered edits. ! Applying 2560 buffered edits. 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. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 78 telescope corrections were flagged in sub-array 8. ! A total of 98 telescope corrections were flagged in sub-array 9. ! A total of 126 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 2. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 110 telescope corrections were flagged in sub-array 8. ! A total of 118 telescope corrections were flagged in sub-array 9. ! A total of 150 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 3. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 78 telescope corrections were flagged in sub-array 8. ! A total of 98 telescope corrections were flagged in sub-array 9. ! A total of 126 telescope corrections were flagged in sub-array 10. ! ! Correcting IF 4. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 7. ! A total of 110 telescope corrections were flagged in sub-array 8. ! A total of 118 telescope corrections were flagged in sub-array 9. ! A total of 150 telescope corrections were flagged in sub-array 10. ! ! Fit before self-cal, rms=0.132253Jy sigma=0.777067 ! Fit after self-cal, rms=0.122312Jy sigma=0.719617 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.426 mas, bmaj=8.224 mas, bpa=-9.01 degrees ! Estimated noise=0.711317 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000720921 Jy ! Component: 100 - total flux cleaned = -0.000621886 Jy ! Component: 150 - total flux cleaned = -0.000578319 Jy ! Component: 200 - total flux cleaned = -0.000372864 Jy ! Component: 250 - total flux cleaned = -0.000295101 Jy ! Component: 300 - total flux cleaned = -0.000219626 Jy ! Component: 350 - total flux cleaned = -0.000147602 Jy ! Component: 400 - total flux cleaned = -0.000147703 Jy ! Total flux subtracted in 400 components = -0.000147703 Jy ! Clean residual min=-0.001787 max=0.001716 Jy/beam ! Clean residual mean=-0.000006 rms=0.000436 Jy/beam ! Combined flux in latest and established models = 0.68294 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 349 components and 0.68294 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122212Jy sigma=0.719132 ! Fit after self-cal, rms=0.122212Jy sigma=0.719120 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 = 7.00008e-05 Jy ! Component: 100 - total flux cleaned = 0.000136806 Jy ! Component: 150 - total flux cleaned = 0.00010327 Jy ! Component: 200 - total flux cleaned = 0.000136548 Jy ! Component: 250 - total flux cleaned = 0.000136207 Jy ! Component: 300 - total flux cleaned = 0.000105043 Jy ! Component: 350 - total flux cleaned = 0.000167593 Jy ! Component: 400 - total flux cleaned = 0.000137256 Jy ! Total flux subtracted in 400 components = 0.000137256 Jy ! Clean residual min=-0.001792 max=0.001712 Jy/beam ! Clean residual mean=-0.000006 rms=0.000430 Jy/beam ! Combined flux in latest and established models = 0.683078 Jy ! Performing phase self-cal ! Adding 89 model components to the UV plane model. ! The established model now contains 381 components and 0.683078 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122154Jy sigma=0.718850 ! Fit after self-cal, rms=0.122153Jy sigma=0.718847 ! 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.122153Jy sigma=0.718847 ! Fit after self-cal, rms=0.122487Jy sigma=0.718478 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.429 mas, bmaj=8.223 mas, bpa=-9.006 degrees ! Estimated noise=0.711756 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00020393 Jy ! Component: 100 - total flux cleaned = -4.84818e-05 Jy ! Component: 150 - total flux cleaned = 0.000158309 Jy ! Component: 200 - total flux cleaned = 0.000222868 Jy ! Component: 250 - total flux cleaned = 0.000222318 Jy ! Component: 300 - total flux cleaned = 0.00025219 Jy ! Component: 350 - total flux cleaned = 0.000252276 Jy ! Component: 400 - total flux cleaned = 0.000280998 Jy ! Total flux subtracted in 400 components = 0.000280998 Jy ! Clean residual min=-0.001755 max=0.001663 Jy/beam ! Clean residual mean=-0.000007 rms=0.000419 Jy/beam ! Combined flux in latest and established models = 0.683358 Jy selfcal ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 405 components and 0.683359 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122426Jy sigma=0.718172 ! Fit after self-cal, rms=0.122424Jy sigma=0.718164 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 = 2.59682e-05 Jy ! Component: 100 - total flux cleaned = 8.28271e-05 Jy ! Component: 150 - total flux cleaned = 8.31051e-05 Jy ! Component: 200 - total flux cleaned = 0.000138219 Jy ! Component: 250 - total flux cleaned = 0.000111023 Jy ! Component: 300 - total flux cleaned = 0.000165128 Jy ! Component: 350 - total flux cleaned = 0.000191409 Jy ! Component: 400 - total flux cleaned = 0.000244011 Jy ! Total flux subtracted in 400 components = 0.000244011 Jy ! Clean residual min=-0.001748 max=0.001648 Jy/beam ! Clean residual mean=-0.000007 rms=0.000414 Jy/beam ! Combined flux in latest and established models = 0.683603 Jy ! Performing phase self-cal ! Adding 91 model components to the UV plane model. ! The established model now contains 421 components and 0.683603 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122381Jy sigma=0.717972 ! Fit after self-cal, rms=0.122381Jy sigma=0.717970 ! 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.122381Jy sigma=0.717970 ! Fit after self-cal, rms=0.122208Jy sigma=0.703260 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.634 mas, bmaj=11.76 mas, bpa=-12.15 degrees ! Estimated noise=0.833875 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -2.60176e-05 Jy ! Component: 100 - total flux cleaned = 9.18948e-05 Jy ! Component: 150 - total flux cleaned = 0.000145663 Jy ! Component: 200 - total flux cleaned = 0.000228501 Jy ! Component: 250 - total flux cleaned = 0.000290088 Jy ! Component: 300 - total flux cleaned = 0.000349521 Jy ! Component: 350 - total flux cleaned = 0.00039243 Jy ! Component: 400 - total flux cleaned = 0.000433562 Jy ! Total flux subtracted in 400 components = 0.000433562 Jy ! Clean residual min=-0.001786 max=0.001491 Jy/beam ! Clean residual mean=-0.000013 rms=0.000440 Jy/beam ! Combined flux in latest and established models = 0.684036 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 460 components and 0.684036 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122207Jy sigma=0.703208 ! Fit after self-cal, rms=0.122219Jy sigma=0.703185 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 = 4.5751e-05 Jy ! Component: 100 - total flux cleaned = 8.66105e-05 Jy ! Component: 150 - total flux cleaned = 0.000112545 Jy ! Component: 200 - total flux cleaned = 0.000137691 Jy ! Component: 250 - total flux cleaned = 0.00014988 Jy ! Component: 300 - total flux cleaned = 0.000162123 Jy ! Component: 350 - total flux cleaned = 0.00017442 Jy ! Component: 400 - total flux cleaned = 0.000209488 Jy ! Total flux subtracted in 400 components = 0.000209488 Jy ! Clean residual min=-0.001800 max=0.001516 Jy/beam ! Clean residual mean=-0.000014 rms=0.000438 Jy/beam ! Combined flux in latest and established models = 0.684246 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 484 components and 0.684246 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122216Jy sigma=0.703174 ! Fit after self-cal, rms=0.122218Jy sigma=0.703171 ! 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=2.431 mas, bmaj=8.224 mas, bpa=-9.002 degrees ! Estimated noise=0.712173 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.431 x 8.224 at -9.002 degrees (North through East) ! Clean map min=-0.003275 max=0.60728 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. ! ! Fit before self-cal, rms=0.122218Jy sigma=0.703171 ! Fit after self-cal, rms=0.119156Jy sigma=0.656919 ! Inverting map and beam ! Estimated beam: bmin=2.431 mas, bmaj=8.226 mas, bpa=-8.994 degrees ! Estimated noise=0.711918 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 8.0907e-05 Jy ! Component: 100 - total flux cleaned = 0.000275601 Jy ! Component: 150 - total flux cleaned = 0.000366093 Jy ! Component: 200 - total flux cleaned = 0.000424169 Jy ! Component: 250 - total flux cleaned = 0.000451653 Jy ! Component: 300 - total flux cleaned = 0.000452072 Jy ! Component: 350 - total flux cleaned = 0.000505381 Jy ! Component: 400 - total flux cleaned = 0.000531318 Jy ! Total flux subtracted in 400 components = 0.000531318 Jy ! Clean residual min=-0.001662 max=0.001632 Jy/beam ! Clean residual mean=-0.000007 rms=0.000400 Jy/beam ! Combined flux in latest and established models = 0.684777 Jy ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 524 components and 0.684777 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119108Jy sigma=0.656673 ! Fit after self-cal, rms=0.119106Jy sigma=0.656670 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 7.92822e-05 Jy ! Component: 100 - total flux cleaned = 0.000129884 Jy ! Component: 150 - total flux cleaned = 0.000130313 Jy ! Component: 200 - total flux cleaned = 0.000154736 Jy ! Component: 250 - total flux cleaned = 0.000179056 Jy ! Component: 300 - total flux cleaned = 0.000178979 Jy ! Component: 350 - total flux cleaned = 0.000179414 Jy ! Component: 400 - total flux cleaned = 0.000225901 Jy ! Total flux subtracted in 400 components = 0.000225901 Jy ! Clean residual min=-0.001649 max=0.001622 Jy/beam ! Clean residual mean=-0.000007 rms=0.000397 Jy/beam ! Combined flux in latest and established models = 0.685003 Jy ! Performing phase self-cal ! Adding 90 model components to the UV plane model. ! The established model now contains 544 components and 0.685003 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119071Jy sigma=0.656513 ! Fit after self-cal, rms=0.119069Jy sigma=0.656513 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119069Jy sigma=0.656513 ! Fit after self-cal, rms=0.119142Jy sigma=0.656482 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.113 mas, bmaj=6.717 mas, bpa=-8.278 degrees ! Estimated noise=1.23556 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.465326 Jy ! Component: 100 - total flux cleaned = 0.569476 Jy ! Component: 150 - total flux cleaned = 0.606868 Jy ! Component: 200 - total flux cleaned = 0.628453 Jy ! Component: 250 - total flux cleaned = 0.642565 Jy ! Component: 300 - total flux cleaned = 0.652257 Jy ! Component: 350 - total flux cleaned = 0.659224 Jy ! Component: 400 - total flux cleaned = 0.664299 Jy ! Total flux subtracted in 400 components = 0.664299 Jy ! Clean residual min=-0.004108 max=0.004142 Jy/beam ! Clean residual mean=-0.000003 rms=0.001014 Jy/beam ! Combined flux in latest and established models = 0.664299 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 49 model components to the UV plane model. ! The established model now contains 49 components and 0.664299 Jy ! Inverting map and beam ! Estimated beam: bmin=2.432 mas, bmaj=8.226 mas, bpa=-8.988 degrees ! Estimated noise=0.712143 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.00489312 Jy ! Component: 100 - total flux cleaned = 0.00828625 Jy ! Component: 150 - total flux cleaned = 0.0106617 Jy ! Component: 200 - total flux cleaned = 0.0122252 Jy ! Component: 250 - total flux cleaned = 0.0133521 Jy ! Component: 300 - total flux cleaned = 0.0143856 Jy ! Component: 350 - total flux cleaned = 0.0151199 Jy ! Component: 400 - total flux cleaned = 0.0155342 Jy ! Total flux subtracted in 400 components = 0.0155342 Jy ! Clean residual min=-0.001977 max=0.002131 Jy/beam ! Clean residual mean=-0.000009 rms=0.000494 Jy/beam ! Combined flux in latest and established models = 0.679833 Jy ! Adding 62 model components to the UV plane model. ! The established model now contains 105 components and 0.679833 Jy ! Inverting map ! Added new window around map position (-47.2, -58). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00132839 Jy ! Component: 100 - total flux cleaned = 0.00164676 Jy ! Component: 150 - total flux cleaned = 0.00195111 Jy ! Component: 200 - total flux cleaned = 0.00224458 Jy ! Component: 250 - total flux cleaned = 0.00247265 Jy ! Component: 300 - total flux cleaned = 0.00269329 Jy ! Component: 350 - total flux cleaned = 0.00285438 Jy ! Component: 400 - total flux cleaned = 0.00300914 Jy ! Total flux subtracted in 400 components = 0.00300914 Jy ! Clean residual min=-0.001886 max=0.001801 Jy/beam ! Clean residual mean=-0.000009 rms=0.000454 Jy/beam ! Combined flux in latest and established models = 0.682842 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 138 components and 0.682842 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=8.633 mas, bmaj=11.76 mas, bpa=-12.14 degrees ! Estimated noise=0.833694 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00149407 Jy ! Component: 100 - total flux cleaned = 0.0020278 Jy ! Component: 150 - total flux cleaned = 0.00210449 Jy ! Component: 200 - total flux cleaned = 0.00221189 Jy ! Component: 250 - total flux cleaned = 0.00221217 Jy ! Component: 300 - total flux cleaned = 0.00217885 Jy ! Component: 350 - total flux cleaned = 0.00208766 Jy ! Component: 400 - total flux cleaned = 0.00202948 Jy ! Total flux subtracted in 400 components = 0.00202948 Jy ! Clean residual min=-0.001852 max=0.001558 Jy/beam ! Clean residual mean=-0.000019 rms=0.000474 Jy/beam ! Combined flux in latest and established models = 0.684871 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 = -5.47117e-05 Jy ! Component: 100 - total flux cleaned = -8.10821e-05 Jy ! Component: 150 - total flux cleaned = -0.000158481 Jy ! Component: 200 - total flux cleaned = -0.000208752 Jy ! Component: 250 - total flux cleaned = -0.000232699 Jy ! Component: 300 - total flux cleaned = -0.000256944 Jy ! Component: 350 - total flux cleaned = -0.00028005 Jy ! Component: 400 - total flux cleaned = -0.00025824 Jy ! Total flux subtracted in 400 components = -0.00025824 Jy ! Clean residual min=-0.001846 max=0.001486 Jy/beam ! Clean residual mean=-0.000018 rms=0.000459 Jy/beam ! Combined flux in latest and established models = 0.684613 Jy ! Adding 102 model components to the UV plane model. ! The established model now contains 236 components and 0.684613 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.432 mas, bmaj=8.226 mas, bpa=-8.988 degrees ! Estimated noise=0.712143 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000699296 Jy ! Component: 100 - total flux cleaned = -0.000698841 Jy ! Component: 150 - total flux cleaned = -0.000475957 Jy ! Component: 200 - total flux cleaned = -0.000433908 Jy ! Component: 250 - total flux cleaned = -0.000391185 Jy ! Component: 300 - total flux cleaned = -0.000309906 Jy ! Component: 350 - total flux cleaned = -0.000270669 Jy ! Component: 400 - total flux cleaned = -0.000308503 Jy ! Total flux subtracted in 400 components = -0.000308503 Jy ! Clean residual min=-0.001792 max=0.001712 Jy/beam ! Clean residual mean=-0.000009 rms=0.000415 Jy/beam ! Combined flux in latest and established models = 0.684304 Jy selfcal ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 272 components and 0.684304 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119159Jy sigma=0.657167 ! Fit after self-cal, rms=0.119157Jy sigma=0.657062 wmodel /scr/pima/rdv93_uvs/J0914+0245_S_map.mod ! Writing 272 model components to file: /scr/pima/rdv93_uvs/J0914+0245_S_map.mod wobs /scr/pima/rdv93_uvs/J0914+0245_S_uvs.fits ! Applying 2074 buffered edits. ! Writing UV FITS file: /scr/pima/rdv93_uvs/J0914+0245_S_uvs.fits wwins /scr/pima/rdv93_uvs/J0914+0245_S_map.win ! wwins: Wrote 6 windows to /scr/pima/rdv93_uvs/J0914+0245_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.000388624 Jy ! Component: 100 - total flux cleaned = -0.000532669 Jy ! Component: 150 - total flux cleaned = -0.000754308 Jy ! Component: 200 - total flux cleaned = -0.000888146 Jy ! Component: 250 - total flux cleaned = -0.00104428 Jy ! Component: 300 - total flux cleaned = -0.00124766 Jy ! Component: 350 - total flux cleaned = -0.00142143 Jy ! Component: 400 - total flux cleaned = -0.00154347 Jy ! Component: 450 - total flux cleaned = -0.00175975 Jy ! Component: 500 - total flux cleaned = -0.0019015 Jy ! Component: 550 - total flux cleaned = -0.00204145 Jy ! Component: 600 - total flux cleaned = -0.00217941 Jy ! Component: 650 - total flux cleaned = -0.00231491 Jy ! Component: 700 - total flux cleaned = -0.00242655 Jy ! Component: 750 - total flux cleaned = -0.00253705 Jy ! Component: 800 - total flux cleaned = -0.00266848 Jy ! Component: 850 - total flux cleaned = -0.00277686 Jy ! Component: 900 - total flux cleaned = -0.00292712 Jy ! Component: 950 - total flux cleaned = -0.00299103 Jy ! Component: 1000 - total flux cleaned = -0.00309665 Jy ! Component: 1050 - total flux cleaned = -0.00324288 Jy ! Component: 1100 - total flux cleaned = -0.00330489 Jy ! Component: 1150 - total flux cleaned = -0.0034074 Jy ! Component: 1200 - total flux cleaned = -0.00346866 Jy ! Component: 1250 - total flux cleaned = -0.0035093 Jy ! Component: 1300 - total flux cleaned = -0.0036698 Jy ! Component: 1350 - total flux cleaned = -0.00370949 Jy ! Component: 1400 - total flux cleaned = -0.00376869 Jy ! Component: 1450 - total flux cleaned = -0.00382755 Jy ! Component: 1500 - total flux cleaned = -0.00392477 Jy ! Component: 1550 - total flux cleaned = -0.0038862 Jy ! Component: 1600 - total flux cleaned = -0.00400122 Jy ! Component: 1650 - total flux cleaned = -0.00398234 Jy ! Component: 1700 - total flux cleaned = -0.00403904 Jy ! Component: 1750 - total flux cleaned = -0.00409546 Jy ! Component: 1800 - total flux cleaned = -0.0040767 Jy ! Component: 1850 - total flux cleaned = -0.00411396 Jy ! Component: 1900 - total flux cleaned = -0.00415069 Jy ! Component: 1950 - total flux cleaned = -0.00418726 Jy ! Component: 2000 - total flux cleaned = -0.00415093 Jy ! Component: 2050 - total flux cleaned = -0.00424164 Jy ! Component: 2100 - total flux cleaned = -0.00425977 Jy ! Component: 2150 - total flux cleaned = -0.00429552 Jy ! Component: 2200 - total flux cleaned = -0.0043135 Jy ! Component: 2250 - total flux cleaned = -0.00438429 Jy ! Component: 2300 - total flux cleaned = -0.00443725 Jy ! Component: 2350 - total flux cleaned = -0.00438465 Jy ! Component: 2400 - total flux cleaned = -0.00450689 Jy ! Component: 2450 - total flux cleaned = -0.00448938 Jy ! Component: 2500 - total flux cleaned = -0.00450668 Jy ! Component: 2550 - total flux cleaned = -0.00457561 Jy ! Component: 2600 - total flux cleaned = -0.00462698 Jy ! Component: 2650 - total flux cleaned = -0.00462713 Jy ! Component: 2700 - total flux cleaned = -0.00464413 Jy ! Component: 2750 - total flux cleaned = -0.00472842 Jy ! Component: 2800 - total flux cleaned = -0.00469482 Jy ! Component: 2850 - total flux cleaned = -0.00477835 Jy ! Component: 2900 - total flux cleaned = -0.00476166 Jy ! Component: 2950 - total flux cleaned = -0.00487752 Jy ! Component: 3000 - total flux cleaned = -0.0049104 Jy ! Component: 3050 - total flux cleaned = -0.00494321 Jy ! Component: 3100 - total flux cleaned = -0.00495955 Jy ! Component: 3150 - total flux cleaned = -0.00500846 Jy ! Component: 3200 - total flux cleaned = -0.00502457 Jy ! Component: 3250 - total flux cleaned = -0.00504072 Jy ! Component: 3300 - total flux cleaned = -0.00502471 Jy ! Component: 3350 - total flux cleaned = -0.00510464 Jy ! Component: 3400 - total flux cleaned = -0.00512064 Jy ! Component: 3450 - total flux cleaned = -0.00518407 Jy ! Component: 3500 - total flux cleaned = -0.00519982 Jy ! Component: 3550 - total flux cleaned = -0.00516835 Jy ! Component: 3600 - total flux cleaned = -0.00532463 Jy ! Component: 3650 - total flux cleaned = -0.00530894 Jy ! Component: 3700 - total flux cleaned = -0.00526243 Jy ! Component: 3750 - total flux cleaned = -0.00535504 Jy ! Component: 3800 - total flux cleaned = -0.0053704 Jy ! Component: 3850 - total flux cleaned = -0.00543159 Jy ! Component: 3900 - total flux cleaned = -0.00544672 Jy ! Component: 3950 - total flux cleaned = -0.0053556 Jy ! Component: 4000 - total flux cleaned = -0.0054465 Jy ! Component: 4050 - total flux cleaned = -0.00552189 Jy ! Total flux subtracted in 4096 components = -0.00547677 Jy ! Clean residual min=-0.001111 max=0.000805 Jy/beam ! Clean residual mean=-0.000007 rms=0.000289 Jy/beam ! Combined flux in latest and established models = 0.678828 Jy keep ! Adding 290 model components to the UV plane model. ! The established model now contains 561 components and 0.678827 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv93_uvs/J0914+0245_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.432 x 8.226 at -8.988 degrees (North through East) ! Clean map min=-0.0019714 max=0.60542 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J0914+0245_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.6054 rms= 0.00028 Jy quit ! Quitting program ! Log file /scr/pima/rdv93_uvs/J0914+0245_S_uvs.log closed on Tue Oct 8 15:27:28 2013