! Started logfile: /scr/pima/rdv73_uvs/J0607+6720_S_uvs.log on Sat Jan 9 17:45:57 2016 obs /scr/pima/rdv73_uvs/J0607+6720_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J0607+6720_S_uva.fits ! AN table 1: 28 integrations on 45 of 45 possible baselines. ! AN table 2: 25 integrations on 66 of 66 possible baselines. ! AN table 3: 7 integrations on 3 of 3 possible baselines. ! AN table 4: 31 integrations on 91 of 91 possible baselines. ! AN table 5: 2 integrations on 1 of 1 possible baselines. ! AN table 6: 13 integrations on 3 of 3 possible baselines. ! AN table 7: 60 integrations on 6 of 6 possible baselines. ! Apparent sampling: 0.781732 visibilities/baseline/integration-bin. ! Found source: J0607+6720 ! ! 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 19240 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/rdv73_uvs/J0607+6720_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:NY 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:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC 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 50 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 102 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 50 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 77 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 180 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 50 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 102 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 50 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 102 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=1.226907Jy sigma=9.251169 ! Fit after self-cal, rms=0.419245Jy sigma=3.446851 ! 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:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 3:KK ! Flagging baselines of antenna 3:KK in the currently selected channels. flag 4:KK ! Flagging baselines of antenna 4:KK in the currently selected channels. flag 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 4:ZC ! Flagging baselines of antenna 4:ZC in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 5:ZC ! Flagging baselines of antenna 5:ZC in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. flag 7:KK ! Flagging baselines of antenna 7:KK 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 7:ZC ! Flagging baselines of antenna 7:ZC 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=2.659 mas, bmaj=5.817 mas, bpa=-23.95 degrees ! Estimated noise=1.71687 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.557619 Jy ! Component: 100 - total flux cleaned = 0.693973 Jy ! Total flux subtracted in 100 components = 0.693973 Jy ! Clean residual min=-0.028250 max=0.059671 Jy/beam ! Clean residual mean=0.000270 rms=0.007818 Jy/beam ! Combined flux in latest and established models = 0.693973 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 8 components and 0.693973 Jy ! ! Correcting IF 1. ! A total of 46 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 138 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 46 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 71 telescope corrections were flagged in sub-array 4. ! A total of 60 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 46 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 138 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 46 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 138 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.222197Jy sigma=1.505067 ! Fit after self-cal, rms=0.222009Jy sigma=1.503081 ! Inverting map and beam ! Estimated beam: bmin=2.646 mas, bmaj=5.917 mas, bpa=-23.99 degrees ! Estimated noise=1.70521 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0695568 Jy ! Component: 100 - total flux cleaned = 0.112004 Jy ! Total flux subtracted in 100 components = 0.112004 Jy ! Clean residual min=-0.026570 max=0.033462 Jy/beam ! Clean residual mean=0.000207 rms=0.006143 Jy/beam ! Combined flux in latest and established models = 0.805977 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 25 components and 0.805977 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180554Jy sigma=1.211696 ! Fit after self-cal, rms=0.180546Jy sigma=1.211573 ! 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.5 mas, bmaj=7.497 mas, bpa=-24.82 degrees ! Estimated noise=0.941649 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 (-4, -3.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0484148 Jy ! Component: 100 - total flux cleaned = 0.0782333 Jy ! Component: 150 - total flux cleaned = 0.0942831 Jy ! Component: 200 - total flux cleaned = 0.101545 Jy ! Total flux subtracted in 200 components = 0.101545 Jy ! Clean residual min=-0.019018 max=0.020880 Jy/beam ! Clean residual mean=0.000254 rms=0.005114 Jy/beam ! Combined flux in latest and established models = 0.907522 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 62 components and 0.907522 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160860Jy sigma=1.113451 ! Fit after self-cal, rms=0.160628Jy sigma=1.111791 ! 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 1.01 FD 0.90 HN 1.17 KP 1.02 ! LA 0.98 NY 1.00* OV 0.98 PT 1.07 ! SC 0.87 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 0.87 HN 1.14 KP 1.01 ! LA 0.98 NL 0.94 NY 1.00* OV 0.98 ! PT 1.07 SC 0.97 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.94 FD 0.86 HN 1.22 KK 1.00* ! KP 0.97 LA 0.98 MK 0.93 NL 0.98 ! NY 1.00* OV 0.98 PT 1.05 SC 0.96 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 0.96 HN 1.21 KP 1.06 ! LA 0.92 NY 1.00* OV 1.03 PT 1.09 ! SC 0.95 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.94 HN 1.21 KP 1.03 ! LA 0.92 NL 0.97 NY 1.00* OV 1.03 ! PT 1.07 SC 1.05 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.98 FD 0.93 HN 1.30 KK 1.00* ! KP 1.00 LA 0.91 MK 1.02 NL 1.01 ! NY 1.00* OV 1.02 PT 1.05 SC 1.02 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.12 FD 0.97 HN 1.29 KP 1.07 ! LA 0.86 NY 1.00* OV 1.09 PT 1.11 ! SC 0.90 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 FD 0.98 HN 1.29 KP 1.04 ! LA 0.87 NL 1.05 NY 1.00* OV 1.07 ! PT 1.10 SC 0.98 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.03 FD 0.95 HN 1.38 KK 1.00* ! KP 1.02 LA 0.85 MK 1.05 NL 1.12 ! NY 1.00* OV 1.08 PT 1.08 SC 0.95 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.26 FD 0.96 HN 1.19 KP 0.95 ! LA 0.73 NY 1.00* OV 0.87 PT 0.98 ! SC 0.93 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.05 FD 0.94 HN 1.18 KP 0.95 ! LA 0.72 NL 0.97 NY 1.00* OV 0.84 ! PT 0.97 SC 1.03 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.89 FD 0.93 HN 1.25 KK 1.00* ! KP 0.91 LA 0.73 MK 1.16 NL 1.01 ! NY 1.00* OV 0.86 PT 0.94 SC 0.98 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Fit before self-cal, rms=0.160628Jy sigma=1.111791 ! Fit after self-cal, rms=0.100370Jy sigma=0.732820 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.541 mas, bmaj=7.552 mas, bpa=-24.39 degrees ! Estimated noise=0.943979 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.0160271 Jy ! Component: 100 - total flux cleaned = -0.0189663 Jy ! Component: 150 - total flux cleaned = -0.021268 Jy ! Component: 200 - total flux cleaned = -0.0226001 Jy ! Total flux subtracted in 200 components = -0.0226001 Jy ! Clean residual min=-0.009653 max=0.006872 Jy/beam ! Clean residual mean=0.000030 rms=0.001873 Jy/beam ! Combined flux in latest and established models = 0.884922 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 77 components and 0.884922 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.096634Jy sigma=0.706361 ! Fit after self-cal, rms=0.096573Jy sigma=0.705832 ! 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.096573Jy sigma=0.705832 ! Fit after self-cal, rms=0.096375Jy sigma=0.698928 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.551 mas, bmaj=7.567 mas, bpa=-24.23 degrees ! Estimated noise=0.943504 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0111307 Jy ! Component: 100 - total flux cleaned = -0.0128158 Jy ! Component: 150 - total flux cleaned = -0.0130628 Jy ! Component: 200 - total flux cleaned = -0.0132633 Jy ! Total flux subtracted in 200 components = -0.0132633 Jy ! Clean residual min=-0.007593 max=0.005559 Jy/beam ! Clean residual mean=0.000024 rms=0.001506 Jy/beam ! Combined flux in latest and established models = 0.871659 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 90 components and 0.871659 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.094496Jy sigma=0.686429 ! Fit after self-cal, rms=0.094492Jy sigma=0.686373 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 = -1.78275e-06 Jy ! Component: 100 - total flux cleaned = -6.84598e-06 Jy ! Component: 150 - total flux cleaned = -0.00014405 Jy ! Component: 200 - total flux cleaned = -0.000144984 Jy ! Total flux subtracted in 200 components = -0.000144984 Jy ! Clean residual min=-0.006996 max=0.005357 Jy/beam ! Clean residual mean=0.000025 rms=0.001452 Jy/beam ! Combined flux in latest and established models = 0.871514 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 106 components and 0.871514 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.094074Jy sigma=0.683421 ! Fit after self-cal, rms=0.094066Jy sigma=0.683384 ! 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.094066Jy sigma=0.683384 ! Fit after self-cal, rms=0.094111Jy sigma=0.678504 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.555 mas, bmaj=7.58 mas, bpa=-24.03 degrees ! Estimated noise=0.943631 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00771988 Jy ! Component: 100 - total flux cleaned = -0.00891179 Jy ! Component: 150 - total flux cleaned = -0.00971757 Jy ! Component: 200 - total flux cleaned = -0.00984734 Jy ! Total flux subtracted in 200 components = -0.00984734 Jy ! Clean residual min=-0.005707 max=0.004406 Jy/beam ! Clean residual mean=0.000023 rms=0.001201 Jy/beam ! Combined flux in latest and established models = 0.861667 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 113 components and 0.861667 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093184Jy sigma=0.672536 ! Fit after self-cal, rms=0.093189Jy sigma=0.672522 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 = -4.84773e-06 Jy ! Component: 100 - total flux cleaned = -1.00756e-05 Jy ! Component: 150 - total flux cleaned = -1.64308e-05 Jy ! Component: 200 - total flux cleaned = -1.43212e-05 Jy ! Total flux subtracted in 200 components = -1.43212e-05 Jy ! Clean residual min=-0.005382 max=0.004401 Jy/beam ! Clean residual mean=0.000023 rms=0.001170 Jy/beam ! Combined flux in latest and established models = 0.861652 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 123 components and 0.861652 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092981Jy sigma=0.671180 ! Fit after self-cal, rms=0.092975Jy sigma=0.671161 ! 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.092975Jy sigma=0.671161 ! Fit after self-cal, rms=0.093098Jy sigma=0.668497 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.267 mas, bmaj=12.85 mas, bpa=-30.65 degrees ! Estimated noise=1.11848 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00447036 Jy ! Component: 100 - total flux cleaned = -0.00514113 Jy ! Component: 150 - total flux cleaned = -0.00528269 Jy ! Component: 200 - total flux cleaned = -0.0053465 Jy ! Total flux subtracted in 200 components = -0.0053465 Jy ! Clean residual min=-0.004275 max=0.003905 Jy/beam ! Clean residual mean=0.000031 rms=0.001238 Jy/beam ! Combined flux in latest and established models = 0.856306 Jy selfcal ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 133 components and 0.856306 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092834Jy sigma=0.666573 ! Fit after self-cal, rms=0.092840Jy sigma=0.666527 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 = 7.63549e-05 Jy ! Component: 100 - total flux cleaned = 2.66539e-05 Jy ! Component: 150 - total flux cleaned = -2.2261e-05 Jy ! Component: 200 - total flux cleaned = -2.33107e-05 Jy ! Total flux subtracted in 200 components = -2.33107e-05 Jy ! Clean residual min=-0.004295 max=0.003666 Jy/beam ! Clean residual mean=0.000031 rms=0.001226 Jy/beam ! Combined flux in latest and established models = 0.856282 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 142 components and 0.856282 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.092790Jy sigma=0.666144 ! Fit after self-cal, rms=0.092789Jy sigma=0.666136 ! 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: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:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 6:OV,true ! 6: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. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 54 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 240 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.187469Jy sigma=1.521194 ! Fit after self-cal, rms=0.140570Jy sigma=0.752983 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:NY 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:ZC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC 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 40 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.140570Jy sigma=0.752983 ! Fit after self-cal, rms=0.137951Jy sigma=0.733966 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.929 mas, bmaj=2.352 mas, bpa=82.03 degrees ! Estimated noise=1.53938 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.00968195 Jy ! Component: 100 - total flux cleaned = 0.00783034 Jy ! Component: 150 - total flux cleaned = 0.00979559 Jy ! Component: 200 - total flux cleaned = 0.0114991 Jy ! Total flux subtracted in 200 components = 0.0114991 Jy ! Clean residual min=-0.007963 max=0.009432 Jy/beam ! Clean residual mean=0.000004 rms=0.001761 Jy/beam ! Combined flux in latest and established models = 0.867782 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 157 components and 0.867781 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134118Jy sigma=0.719263 ! Fit after self-cal, rms=0.133325Jy sigma=0.715510 ! Inverting map ! Added new window around map position (-6, 13.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00297661 Jy ! Component: 100 - total flux cleaned = 0.00366997 Jy ! Component: 150 - total flux cleaned = 0.00333132 Jy ! Component: 200 - total flux cleaned = 0.00393964 Jy ! Total flux subtracted in 200 components = 0.00393964 Jy ! Clean residual min=-0.008061 max=0.008331 Jy/beam ! Clean residual mean=0.000003 rms=0.001646 Jy/beam ! Combined flux in latest and established models = 0.871721 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 176 components and 0.871721 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.131961Jy sigma=0.711788 ! Fit after self-cal, rms=0.131842Jy sigma=0.711345 ! 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.654 mas, bmaj=3.003 mas, bpa=86.48 degrees ! Estimated noise=0.810876 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.00752358 Jy ! Component: 100 - total flux cleaned = -0.0110249 Jy ! Component: 150 - total flux cleaned = -0.0115344 Jy ! Component: 200 - total flux cleaned = -0.0118453 Jy ! Component: 250 - total flux cleaned = -0.0121359 Jy ! Component: 300 - total flux cleaned = -0.0124133 Jy ! Component: 350 - total flux cleaned = -0.0125459 Jy ! Component: 400 - total flux cleaned = -0.0126742 Jy ! Total flux subtracted in 400 components = -0.0126742 Jy ! Clean residual min=-0.006229 max=0.006346 Jy/beam ! Clean residual mean=-0.000005 rms=0.001323 Jy/beam ! Combined flux in latest and established models = 0.859047 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 209 components and 0.859047 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.130578Jy sigma=0.706701 ! Fit after self-cal, rms=0.130518Jy sigma=0.706410 ! Inverting map ! Added new window around map position (2, 62.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00215363 Jy ! Component: 100 - total flux cleaned = 0.00139677 Jy ! Component: 150 - total flux cleaned = 0.00126338 Jy ! Component: 200 - total flux cleaned = 0.000875936 Jy ! Component: 250 - total flux cleaned = 0.000632217 Jy ! Component: 300 - total flux cleaned = 0.000513381 Jy ! Component: 350 - total flux cleaned = 0.000283766 Jy ! Component: 400 - total flux cleaned = 0.000172298 Jy ! Total flux subtracted in 400 components = 0.000172298 Jy ! Clean residual min=-0.006076 max=0.005949 Jy/beam ! Clean residual mean=-0.000005 rms=0.001275 Jy/beam ! Combined flux in latest and established models = 0.859219 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 234 components and 0.859219 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129565Jy sigma=0.703675 ! Fit after self-cal, rms=0.129540Jy sigma=0.703524 ! 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.000505429 Jy ! Component: 100 - total flux cleaned = -0.000613696 Jy ! Component: 150 - total flux cleaned = -0.000721065 Jy ! Component: 200 - total flux cleaned = -0.00092839 Jy ! Component: 250 - total flux cleaned = -0.00103152 Jy ! Component: 300 - total flux cleaned = -0.001229 Jy ! Component: 350 - total flux cleaned = -0.00132422 Jy ! Component: 400 - total flux cleaned = -0.00151406 Jy ! Total flux subtracted in 400 components = -0.00151406 Jy ! Clean residual min=-0.006090 max=0.005906 Jy/beam ! Clean residual mean=-0.000004 rms=0.001257 Jy/beam ! Combined flux in latest and established models = 0.857705 Jy ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 253 components and 0.857705 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129051Jy sigma=0.702166 ! Fit after self-cal, rms=0.129032Jy sigma=0.702112 ! 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. ! 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 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 77 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 240 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 56 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 156 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.151362Jy sigma=0.830351 ! Fit after self-cal, rms=0.120547Jy sigma=0.692370 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.608 mas, bmaj=2.992 mas, bpa=87.4 degrees ! Estimated noise=0.808805 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00456698 Jy ! Component: 100 - total flux cleaned = 0.00430218 Jy ! Component: 150 - total flux cleaned = 0.00409453 Jy ! Component: 200 - total flux cleaned = 0.00341164 Jy ! Component: 250 - total flux cleaned = 0.00296162 Jy ! Component: 300 - total flux cleaned = 0.0025597 Jy ! Component: 350 - total flux cleaned = 0.00218782 Jy ! Component: 400 - total flux cleaned = 0.00196048 Jy ! Total flux subtracted in 400 components = 0.00196048 Jy ! Clean residual min=-0.004815 max=0.004918 Jy/beam ! Clean residual mean=0.000003 rms=0.001011 Jy/beam ! Combined flux in latest and established models = 0.859666 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 263 components and 0.859666 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118767Jy sigma=0.684539 ! Fit after self-cal, rms=0.118724Jy sigma=0.684457 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 (14, -46). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00287223 Jy ! Component: 100 - total flux cleaned = 0.00250074 Jy ! Component: 150 - total flux cleaned = 0.00204986 Jy ! Component: 200 - total flux cleaned = 0.00194559 Jy ! Component: 250 - total flux cleaned = 0.00184612 Jy ! Component: 300 - total flux cleaned = 0.00165059 Jy ! Component: 350 - total flux cleaned = 0.00146281 Jy ! Component: 400 - total flux cleaned = 0.00137254 Jy ! Total flux subtracted in 400 components = 0.00137254 Jy ! Clean residual min=-0.004315 max=0.004348 Jy/beam ! Clean residual mean=0.000001 rms=0.000976 Jy/beam ! Combined flux in latest and established models = 0.861038 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 291 components and 0.861038 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118134Jy sigma=0.682359 ! Fit after self-cal, rms=0.118123Jy sigma=0.682273 ! Inverting map ! Added new window around map position (37.2, 22.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00156016 Jy ! Component: 100 - total flux cleaned = 0.00143272 Jy ! Component: 150 - total flux cleaned = 0.0014331 Jy ! Component: 200 - total flux cleaned = 0.0014328 Jy ! Component: 250 - total flux cleaned = 0.00124632 Jy ! Component: 300 - total flux cleaned = 0.00106612 Jy ! Component: 350 - total flux cleaned = 0.000978281 Jy ! Component: 400 - total flux cleaned = 0.000808461 Jy ! Total flux subtracted in 400 components = 0.000808461 Jy ! Clean residual min=-0.004125 max=0.004143 Jy/beam ! Clean residual mean=0.000004 rms=0.000943 Jy/beam ! Combined flux in latest and established models = 0.861847 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 323 components and 0.861847 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117582Jy sigma=0.680604 ! Fit after self-cal, rms=0.117574Jy sigma=0.680550 ! Inverting map ! Added new window around map position (13.2, 17.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00149373 Jy ! Component: 100 - total flux cleaned = 0.00115755 Jy ! Component: 150 - total flux cleaned = 0.000859351 Jy ! Component: 200 - total flux cleaned = 0.000485076 Jy ! Component: 250 - total flux cleaned = 0.00030742 Jy ! Component: 300 - total flux cleaned = -0.000120865 Jy ! Component: 350 - total flux cleaned = -0.000371827 Jy ! Component: 400 - total flux cleaned = -0.000698604 Jy ! Total flux subtracted in 400 components = -0.000698604 Jy ! Clean residual min=-0.003951 max=0.003850 Jy/beam ! Clean residual mean=0.000004 rms=0.000913 Jy/beam ! Combined flux in latest and established models = 0.861148 Jy ! Performing phase self-cal ! Adding 71 model components to the UV plane model. ! The established model now contains 365 components and 0.861148 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117110Jy sigma=0.678986 ! Fit after self-cal, rms=0.117103Jy sigma=0.678936 ! 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.117103Jy sigma=0.678936 ! Fit after self-cal, rms=0.113572Jy sigma=0.674386 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.577 mas, bmaj=2.99 mas, bpa=87.88 degrees ! Estimated noise=0.807743 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00326161 Jy ! Component: 100 - total flux cleaned = 0.00287644 Jy ! Component: 150 - total flux cleaned = 0.00273327 Jy ! Component: 200 - total flux cleaned = 0.00235856 Jy ! Component: 250 - total flux cleaned = 0.00193797 Jy ! Component: 300 - total flux cleaned = 0.00147308 Jy ! Component: 350 - total flux cleaned = 0.00138859 Jy ! Component: 400 - total flux cleaned = 0.00122734 Jy ! Total flux subtracted in 400 components = 0.00122734 Jy ! Clean residual min=-0.003424 max=0.003215 Jy/beam ! Clean residual mean=0.000005 rms=0.000773 Jy/beam ! Combined flux in latest and established models = 0.862375 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 375 components and 0.862375 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112658Jy sigma=0.670154 ! Fit after self-cal, rms=0.112649Jy sigma=0.670124 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.000399027 Jy ! Component: 100 - total flux cleaned = -0.000547048 Jy ! Component: 150 - total flux cleaned = -0.000690443 Jy ! Component: 200 - total flux cleaned = -0.000830347 Jy ! Component: 250 - total flux cleaned = -0.000898894 Jy ! Component: 300 - total flux cleaned = -0.00103212 Jy ! Component: 350 - total flux cleaned = -0.00122757 Jy ! Component: 400 - total flux cleaned = -0.00135426 Jy ! Total flux subtracted in 400 components = -0.00135426 Jy ! Clean residual min=-0.003365 max=0.003108 Jy/beam ! Clean residual mean=0.000005 rms=0.000758 Jy/beam ! Combined flux in latest and established models = 0.861021 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 399 components and 0.861021 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112447Jy sigma=0.669254 ! Fit after self-cal, rms=0.112445Jy sigma=0.669237 ! 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.112445Jy sigma=0.669237 ! Fit after self-cal, rms=0.109456Jy sigma=0.656923 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.081 mas, bmaj=11.62 mas, bpa=-27.73 degrees ! Estimated noise=1.11963 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00105088 Jy ! Component: 100 - total flux cleaned = 0.000886916 Jy ! Component: 150 - total flux cleaned = 0.000764549 Jy ! Component: 200 - total flux cleaned = 0.000592035 Jy ! Component: 250 - total flux cleaned = 0.000538258 Jy ! Component: 300 - total flux cleaned = 0.00038789 Jy ! Component: 350 - total flux cleaned = 0.000294394 Jy ! Component: 400 - total flux cleaned = 0.000249849 Jy ! Total flux subtracted in 400 components = 0.000249849 Jy ! Clean residual min=-0.003561 max=0.003695 Jy/beam ! Clean residual mean=0.000015 rms=0.000879 Jy/beam ! Combined flux in latest and established models = 0.861271 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 432 components and 0.861271 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109477Jy sigma=0.657035 ! Fit after self-cal, rms=0.109474Jy sigma=0.656976 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 (-50.8, -102). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00211987 Jy ! Component: 100 - total flux cleaned = 0.00212187 Jy ! Component: 150 - total flux cleaned = 0.00211862 Jy ! Component: 200 - total flux cleaned = 0.00217781 Jy ! Component: 250 - total flux cleaned = 0.00233691 Jy ! Component: 300 - total flux cleaned = 0.00243349 Jy ! Component: 350 - total flux cleaned = 0.00247828 Jy ! Component: 400 - total flux cleaned = 0.00256584 Jy ! Total flux subtracted in 400 components = 0.00256584 Jy ! Clean residual min=-0.002879 max=0.002582 Jy/beam ! Clean residual mean=0.000010 rms=0.000742 Jy/beam ! Combined flux in latest and established models = 0.863837 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 482 components and 0.863837 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109467Jy sigma=0.656782 ! Fit after self-cal, rms=0.109466Jy sigma=0.656727 ! 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.555 mas, bmaj=2.989 mas, bpa=88.25 degrees ! Estimated noise=0.807024 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.555 x 2.989 at 88.25 degrees (North through East) ! Clean map min=-0.0060654 max=0.69775 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.109466Jy sigma=0.656727 ! Fit after self-cal, rms=0.106411Jy sigma=0.620416 ! Inverting map and beam ! Estimated beam: bmin=2.556 mas, bmaj=2.991 mas, bpa=88.24 degrees ! Estimated noise=0.809716 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00260784 Jy ! Component: 100 - total flux cleaned = 0.00233238 Jy ! Component: 150 - total flux cleaned = 0.00243774 Jy ! Component: 200 - total flux cleaned = 0.00253531 Jy ! Component: 250 - total flux cleaned = 0.00296365 Jy ! Component: 300 - total flux cleaned = 0.00312135 Jy ! Component: 350 - total flux cleaned = 0.00312149 Jy ! Component: 400 - total flux cleaned = 0.00311965 Jy ! Total flux subtracted in 400 components = 0.00311965 Jy ! Clean residual min=-0.003104 max=0.002476 Jy/beam ! Clean residual mean=0.000003 rms=0.000625 Jy/beam ! Combined flux in latest and established models = 0.866956 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 512 components and 0.866956 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105866Jy sigma=0.617770 ! Fit after self-cal, rms=0.105859Jy sigma=0.617754 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000216435 Jy ! Component: 100 - total flux cleaned = 0.000146978 Jy ! Component: 150 - total flux cleaned = 0.000281261 Jy ! Component: 200 - total flux cleaned = 0.000215671 Jy ! Component: 250 - total flux cleaned = 0.000216773 Jy ! Component: 300 - total flux cleaned = 0.00015369 Jy ! Component: 350 - total flux cleaned = 0.000154299 Jy ! Component: 400 - total flux cleaned = 3.29204e-05 Jy ! Total flux subtracted in 400 components = 3.29204e-05 Jy ! Clean residual min=-0.003025 max=0.002477 Jy/beam ! Clean residual mean=0.000003 rms=0.000611 Jy/beam ! Combined flux in latest and established models = 0.866989 Jy ! Performing phase self-cal ! Adding 96 model components to the UV plane model. ! The established model now contains 558 components and 0.866989 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.105686Jy sigma=0.616897 ! Fit after self-cal, rms=0.105684Jy sigma=0.616886 ! 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.105684Jy sigma=0.616886 ! Fit after self-cal, rms=0.104846Jy sigma=0.615854 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.86 mas, bmaj=2.351 mas, bpa=85.09 degrees ! Estimated noise=1.49828 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.506901 Jy ! Component: 100 - total flux cleaned = 0.648962 Jy ! Component: 150 - total flux cleaned = 0.71657 Jy ! Component: 200 - total flux cleaned = 0.752038 Jy ! Component: 250 - total flux cleaned = 0.775552 Jy ! Component: 300 - total flux cleaned = 0.79292 Jy ! Component: 350 - total flux cleaned = 0.806319 Jy ! Component: 400 - total flux cleaned = 0.816949 Jy ! Total flux subtracted in 400 components = 0.816949 Jy ! Clean residual min=-0.005088 max=0.006333 Jy/beam ! Clean residual mean=0.000032 rms=0.001295 Jy/beam ! Combined flux in latest and established models = 0.816949 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 38 model components to the UV plane model. ! The established model now contains 38 components and 0.816949 Jy ! Inverting map and beam ! Estimated beam: bmin=2.541 mas, bmaj=2.991 mas, bpa=88.49 degrees ! Estimated noise=0.809127 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.0152773 Jy ! Component: 100 - total flux cleaned = 0.0243786 Jy ! Component: 150 - total flux cleaned = 0.0308359 Jy ! Component: 200 - total flux cleaned = 0.0357369 Jy ! Component: 250 - total flux cleaned = 0.0390321 Jy ! Component: 300 - total flux cleaned = 0.0410338 Jy ! Component: 350 - total flux cleaned = 0.0425811 Jy ! Component: 400 - total flux cleaned = 0.0436334 Jy ! Total flux subtracted in 400 components = 0.0436334 Jy ! Clean residual min=-0.002874 max=0.003145 Jy/beam ! Clean residual mean=0.000007 rms=0.000776 Jy/beam ! Combined flux in latest and established models = 0.860582 Jy ! Adding 81 model components to the UV plane model. ! The established model now contains 106 components and 0.860583 Jy ! Inverting map ! Added new window around map position (6, -13.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00200533 Jy ! Component: 100 - total flux cleaned = 0.00248615 Jy ! Component: 150 - total flux cleaned = 0.0029438 Jy ! Component: 200 - total flux cleaned = 0.00327571 Jy ! Component: 250 - total flux cleaned = 0.00390967 Jy ! Component: 300 - total flux cleaned = 0.00421304 Jy ! Component: 350 - total flux cleaned = 0.00441016 Jy ! Component: 400 - total flux cleaned = 0.00488934 Jy ! Total flux subtracted in 400 components = 0.00488934 Jy ! Clean residual min=-0.002655 max=0.002773 Jy/beam ! Clean residual mean=0.000003 rms=0.000703 Jy/beam ! Combined flux in latest and established models = 0.865472 Jy ! Adding 99 model components to the UV plane model. ! The established model now contains 169 components and 0.865472 Jy ! Inverting map ! Added new window around map position (-21.2, 36.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00174817 Jy ! Component: 100 - total flux cleaned = 0.00211387 Jy ! Component: 150 - total flux cleaned = 0.00255733 Jy ! Component: 200 - total flux cleaned = 0.00264424 Jy ! Component: 250 - total flux cleaned = 0.0028147 Jy ! Component: 300 - total flux cleaned = 0.00306459 Jy ! Component: 350 - total flux cleaned = 0.00322601 Jy ! Component: 400 - total flux cleaned = 0.0033852 Jy ! Total flux subtracted in 400 components = 0.0033852 Jy ! Clean residual min=-0.002630 max=0.002482 Jy/beam ! Clean residual mean=0.000000 rms=0.000663 Jy/beam ! Combined flux in latest and established models = 0.868857 Jy ! Adding 121 model components to the UV plane model. ! The established model now contains 223 components and 0.868857 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.084 mas, bmaj=11.63 mas, bpa=-27.65 degrees ! Estimated noise=1.12487 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00134832 Jy ! Component: 100 - total flux cleaned = 0.00228048 Jy ! Component: 150 - total flux cleaned = 0.00274245 Jy ! Component: 200 - total flux cleaned = 0.00299019 Jy ! Component: 250 - total flux cleaned = 0.00291049 Jy ! Component: 300 - total flux cleaned = 0.00283912 Jy ! Component: 350 - total flux cleaned = 0.00236299 Jy ! Component: 400 - total flux cleaned = 0.00230043 Jy ! Total flux subtracted in 400 components = 0.00230043 Jy ! Clean residual min=-0.002802 max=0.002787 Jy/beam ! Clean residual mean=-0.000001 rms=0.000760 Jy/beam ! Combined flux in latest and established models = 0.871158 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.000122455 Jy ! Component: 100 - total flux cleaned = -0.000183417 Jy ! Component: 150 - total flux cleaned = -0.000414047 Jy ! Component: 200 - total flux cleaned = -0.000467836 Jy ! Component: 250 - total flux cleaned = -0.00057649 Jy ! Component: 300 - total flux cleaned = -0.000734451 Jy ! Component: 350 - total flux cleaned = -0.000684198 Jy ! Component: 400 - total flux cleaned = -0.000734568 Jy ! Total flux subtracted in 400 components = -0.000734568 Jy ! Clean residual min=-0.002815 max=0.002559 Jy/beam ! Clean residual mean=0.000000 rms=0.000718 Jy/beam ! Combined flux in latest and established models = 0.870423 Jy ! Adding 97 model components to the UV plane model. ! The established model now contains 313 components and 0.870423 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.541 mas, bmaj=2.991 mas, bpa=88.49 degrees ! Estimated noise=0.809127 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000661507 Jy ! Component: 100 - total flux cleaned = -0.000821364 Jy ! Component: 150 - total flux cleaned = -0.00120033 Jy ! Component: 200 - total flux cleaned = -0.00141751 Jy ! Component: 250 - total flux cleaned = -0.00155871 Jy ! Component: 300 - total flux cleaned = -0.00176723 Jy ! Component: 350 - total flux cleaned = -0.0018349 Jy ! Component: 400 - total flux cleaned = -0.00196965 Jy ! Total flux subtracted in 400 components = -0.00196965 Jy ! Clean residual min=-0.002354 max=0.002290 Jy/beam ! Clean residual mean=0.000002 rms=0.000593 Jy/beam ! Combined flux in latest and established models = 0.868453 Jy selfcal ! Performing phase self-cal ! Adding 89 model components to the UV plane model. ! The established model now contains 367 components and 0.868453 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104604Jy sigma=0.614917 ! Fit after self-cal, rms=0.104529Jy sigma=0.614542 wmodel /scr/pima/rdv73_uvs/J0607+6720_S_map.mod ! Writing 367 model components to file: /scr/pima/rdv73_uvs/J0607+6720_S_map.mod wobs /scr/pima/rdv73_uvs/J0607+6720_S_uvs.fits ! Applying 25 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J0607+6720_S_uvs.fits wwins /scr/pima/rdv73_uvs/J0607+6720_S_map.win ! wwins: Wrote 10 windows to /scr/pima/rdv73_uvs/J0607+6720_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.000309679 Jy ! Component: 100 - total flux cleaned = -0.000394572 Jy ! Component: 150 - total flux cleaned = -0.00039537 Jy ! Component: 200 - total flux cleaned = -0.000355101 Jy ! Component: 250 - total flux cleaned = -0.000237219 Jy ! Component: 300 - total flux cleaned = -0.000198428 Jy ! Component: 350 - total flux cleaned = -4.78788e-05 Jy ! Component: 400 - total flux cleaned = -4.83555e-05 Jy ! Component: 450 - total flux cleaned = 2.43228e-05 Jy ! Component: 500 - total flux cleaned = 2.43841e-05 Jy ! Component: 550 - total flux cleaned = 0.00013073 Jy ! Component: 600 - total flux cleaned = 0.000306397 Jy ! Component: 650 - total flux cleaned = 0.000445209 Jy ! Component: 700 - total flux cleaned = 0.00058325 Jy ! Component: 750 - total flux cleaned = 0.000821395 Jy ! Component: 800 - total flux cleaned = 0.000888815 Jy ! Component: 850 - total flux cleaned = 0.0010564 Jy ! Component: 900 - total flux cleaned = 0.00128873 Jy ! Component: 950 - total flux cleaned = 0.00155172 Jy ! Component: 1000 - total flux cleaned = 0.0017477 Jy ! Component: 1050 - total flux cleaned = 0.0019744 Jy ! Component: 1100 - total flux cleaned = 0.00219937 Jy ! Component: 1150 - total flux cleaned = 0.00229534 Jy ! Component: 1200 - total flux cleaned = 0.00261284 Jy ! Component: 1250 - total flux cleaned = 0.00283331 Jy ! Component: 1300 - total flux cleaned = 0.00311453 Jy ! Component: 1350 - total flux cleaned = 0.00326959 Jy ! Component: 1400 - total flux cleaned = 0.00339337 Jy ! Component: 1450 - total flux cleaned = 0.0036083 Jy ! Component: 1500 - total flux cleaned = 0.00391344 Jy ! Component: 1550 - total flux cleaned = 0.00391327 Jy ! Component: 1600 - total flux cleaned = 0.00412466 Jy ! Component: 1650 - total flux cleaned = 0.00421474 Jy ! Component: 1700 - total flux cleaned = 0.00436375 Jy ! Component: 1750 - total flux cleaned = 0.00454196 Jy ! Component: 1800 - total flux cleaned = 0.0046895 Jy ! Component: 1850 - total flux cleaned = 0.00468937 Jy ! Component: 1900 - total flux cleaned = 0.00480651 Jy ! Component: 1950 - total flux cleaned = 0.00495185 Jy ! Component: 2000 - total flux cleaned = 0.00506754 Jy ! Component: 2050 - total flux cleaned = 0.0052403 Jy ! Component: 2100 - total flux cleaned = 0.00526889 Jy ! Component: 2150 - total flux cleaned = 0.00538296 Jy ! Component: 2200 - total flux cleaned = 0.00549643 Jy ! Component: 2250 - total flux cleaned = 0.00566588 Jy ! Component: 2300 - total flux cleaned = 0.0058344 Jy ! Component: 2350 - total flux cleaned = 0.00594624 Jy ! Component: 2400 - total flux cleaned = 0.00602987 Jy ! Component: 2450 - total flux cleaned = 0.00611315 Jy ! Component: 2500 - total flux cleaned = 0.00625105 Jy ! Component: 2550 - total flux cleaned = 0.00644342 Jy ! Component: 2600 - total flux cleaned = 0.00655256 Jy ! Component: 2650 - total flux cleaned = 0.00668889 Jy ! Component: 2700 - total flux cleaned = 0.00679723 Jy ! Component: 2750 - total flux cleaned = 0.00687826 Jy ! Component: 2800 - total flux cleaned = 0.00709326 Jy ! Component: 2850 - total flux cleaned = 0.00711981 Jy ! Component: 2900 - total flux cleaned = 0.0072264 Jy ! Component: 2950 - total flux cleaned = 0.00735923 Jy ! Component: 3000 - total flux cleaned = 0.00743841 Jy ! Component: 3050 - total flux cleaned = 0.00762276 Jy ! Component: 3100 - total flux cleaned = 0.00778019 Jy ! Component: 3150 - total flux cleaned = 0.00772796 Jy ! Component: 3200 - total flux cleaned = 0.00791016 Jy ! Component: 3250 - total flux cleaned = 0.00809147 Jy ! Component: 3300 - total flux cleaned = 0.00814316 Jy ! Component: 3350 - total flux cleaned = 0.00819467 Jy ! Component: 3400 - total flux cleaned = 0.00822024 Jy ! Component: 3450 - total flux cleaned = 0.00852652 Jy ! Component: 3500 - total flux cleaned = 0.00857732 Jy ! Component: 3550 - total flux cleaned = 0.00870398 Jy ! Component: 3600 - total flux cleaned = 0.00865372 Jy ! Component: 3650 - total flux cleaned = 0.00888016 Jy ! Component: 3700 - total flux cleaned = 0.00898041 Jy ! Component: 3750 - total flux cleaned = 0.0091552 Jy ! Component: 3800 - total flux cleaned = 0.00920494 Jy ! Component: 3850 - total flux cleaned = 0.00920505 Jy ! Component: 3900 - total flux cleaned = 0.00935322 Jy ! Component: 3950 - total flux cleaned = 0.00942705 Jy ! Component: 4000 - total flux cleaned = 0.00954972 Jy ! Component: 4050 - total flux cleaned = 0.00962303 Jy ! Total flux subtracted in 4096 components = 0.00967176 Jy ! Clean residual min=-0.001478 max=0.001252 Jy/beam ! Clean residual mean=0.000002 rms=0.000444 Jy/beam ! Combined flux in latest and established models = 0.878125 Jy keep ! Adding 310 model components to the UV plane model. ! The established model now contains 675 components and 0.878125 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J0607+6720_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.541 x 2.991 at 88.49 degrees (North through East) ! Clean map min=-0.0036709 max=0.69733 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J0607+6720_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.6973 rms= 0.00044 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J0607+6720_S_uvs.log closed on Sat Jan 9 17:46:06 2016