! Started logfile: /scr/pima/rdv73_uvs/J0921+6215_S_uvs.log on Sat Jan 9 17:47:56 2016 obs /scr/pima/rdv73_uvs/J0921+6215_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J0921+6215_S_uva.fits ! AN table 1: 5 integrations on 10 of 10 possible baselines. ! AN table 2: 3 integrations on 21 of 21 possible baselines. ! AN table 3: 10 integrations on 3 of 3 possible baselines. ! AN table 4: 5 integrations on 6 of 6 possible baselines. ! AN table 5: 6 integrations on 66 of 66 possible baselines. ! AN table 6: 4 integrations on 10 of 10 possible baselines. ! AN table 7: 10 integrations on 91 of 91 possible baselines. ! AN table 8: 2 integrations on 1 of 1 possible baselines. ! AN table 9: 10 integrations on 66 of 66 possible baselines. ! Apparent sampling: 0.644658 visibilities/baseline/integration-bin. ! Found source: J0921+6215 ! ! 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 5624 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/J0921+6215_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:HN self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR 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:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:BR 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 ! 6:ZC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:OV self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC self-cal status=correctable weight=1 ! 8:MK self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KK self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 9: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 4 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 63 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 12 telescope corrections were flagged in sub-array 6. ! A total of 73 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 64 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 63 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 63 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=1.440466Jy sigma=4.625568 ! Fit after self-cal, rms=0.340650Jy sigma=1.075509 ! 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:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 3:ZC ! Flagging baselines of antenna 3:ZC 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:NY ! Flagging baselines of antenna 5:NY 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 6:ZC ! Flagging baselines of antenna 6:ZC 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. flag 8:ZC ! Flagging baselines of antenna 8:ZC in the currently selected channels. flag 9:KK ! Flagging baselines of antenna 9:KK 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 9:ZC ! Flagging baselines of antenna 9: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=3.295 mas, bmaj=5.55 mas, bpa=58.58 degrees ! Estimated noise=5.916 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.754988 Jy ! Component: 100 - total flux cleaned = 0.939329 Jy ! Total flux subtracted in 100 components = 0.939329 Jy ! Clean residual min=-0.035595 max=0.101720 Jy/beam ! Clean residual mean=-0.000050 rms=0.011783 Jy/beam ! Combined flux in latest and established models = 0.939329 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.939329 Jy ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 25 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! A total of 15 telescope corrections were flagged in sub-array 7. ! A total of 12 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 25 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 25 telescope corrections were flagged in sub-array 7. ! A total of 18 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.328787Jy sigma=0.927009 ! Fit after self-cal, rms=0.328063Jy sigma=0.923807 ! Inverting map and beam ! Estimated beam: bmin=3.295 mas, bmaj=5.55 mas, bpa=58.58 degrees ! Estimated noise=5.916 mJy/beam. ! Added new window around map position (-1.6, 5.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.108608 Jy ! Component: 100 - total flux cleaned = 0.184255 Jy ! Total flux subtracted in 100 components = 0.184255 Jy ! Clean residual min=-0.031010 max=0.060390 Jy/beam ! Clean residual mean=-0.000004 rms=0.008737 Jy/beam ! Combined flux in latest and established models = 1.12358 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 39 components and 1.12358 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.259250Jy sigma=0.709938 ! Fit after self-cal, rms=0.258214Jy sigma=0.705329 ! 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.813 mas, bmaj=6.331 mas, bpa=60.45 degrees ! Estimated noise=4.59878 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.0792859 Jy ! Component: 100 - total flux cleaned = 0.127689 Jy ! Component: 150 - total flux cleaned = 0.159644 Jy ! Component: 200 - total flux cleaned = 0.18113 Jy ! Total flux subtracted in 200 components = 0.18113 Jy ! Clean residual min=-0.029507 max=0.033317 Jy/beam ! Clean residual mean=-0.000414 rms=0.007640 Jy/beam ! Combined flux in latest and established models = 1.30471 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 61 components and 1.30471 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.218872Jy sigma=0.622449 ! Fit after self-cal, rms=0.217876Jy sigma=0.619151 ! 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 15 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 HN 1.10 KP 1.00 LA 0.99 ! OV 1.04 PT 1.07 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HN 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.96 FD 0.87 HN 1.30 KP 0.96 ! LA 0.98 NL 1.04 NY 1.00* OV 1.00 ! PT 1.07 SC 0.94 WZ 1.00* ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 0.99 FD 0.90 HN 1.11 KK 1.00* ! KP 0.97 LA 0.94 MK 0.84 NL 0.96 ! NY 1.00* OV 1.00 PT 1.03 SC 0.93 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! BR 1.00 FD 0.90 HN 1.17 KK 1.00* ! KP 1.02 LA 0.97 MK 0.87 NY 1.00* ! OV 0.99 PT 1.09 WZ 1.00* ZC 1.00* ! ! ! Correcting IF 2. ! A total of 15 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 HN 1.16 KP 1.03 LA 0.94 ! OV 1.05 PT 1.11 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HN 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 0.95 HN 1.33 KP 1.05 ! LA 0.90 NL 1.09 NY 1.00* OV 1.04 ! PT 1.03 SC 1.02 WZ 1.00* ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 1.01 FD 0.98 HN 1.18 KK 1.00* ! KP 0.99 LA 0.89 MK 0.90 NL 0.97 ! NY 1.00* OV 1.02 PT 1.02 SC 1.02 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! BR 1.03 FD 0.97 HN 1.27 KK 1.00* ! KP 1.05 LA 0.93 MK 0.97 NY 1.00* ! OV 1.04 PT 1.07 WZ 1.00* ZC 1.00* ! ! ! Correcting IF 3. ! A total of 15 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 HN 1.17 KP 1.05 LA 0.88 ! OV 1.14 PT 1.13 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HN 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 FD 0.96 HN 1.37 KP 1.05 ! LA 0.84 NL 1.18 NY 1.00* OV 1.06 ! PT 1.07 SC 0.93 WZ 1.00* ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 1.04 FD 0.98 HN 1.23 KK 1.00* ! KP 0.96 LA 0.82 MK 0.95 NL 1.07 ! NY 1.00* OV 1.09 PT 1.06 SC 0.92 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! BR 1.08 FD 0.96 HN 1.36 KK 1.00* ! KP 1.04 LA 0.83 MK 0.99 NY 1.00* ! OV 1.09 PT 1.11 WZ 1.00* ZC 1.00* ! ! ! Correcting IF 4. ! A total of 15 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* HN 0.00* NL 0.00* NY 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 HN 1.05 KP 0.94 LA 0.74 ! OV 0.90 PT 0.96 ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HN 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.96 FD 0.94 HN 1.21 KP 0.92 ! LA 0.71 NL 1.08 NY 1.00* OV 0.84 ! PT 0.97 SC 0.96 WZ 1.00* ZC 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 6. ! Telescope amplitude corrections in sub-array 6: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ZC 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! BR 0.93 FD 0.94 HN 1.14 KK 1.00* ! KP 0.90 LA 0.72 MK 1.04 NL 0.95 ! NY 1.00* OV 0.86 PT 0.91 SC 0.96 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! MK 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! BR 1.03 FD 0.96 HN 1.23 KK 1.00* ! KP 0.94 LA 0.73 MK 1.08 NY 1.00* ! OV 0.87 PT 0.97 WZ 1.00* ZC 1.00* ! ! ! Fit before self-cal, rms=0.217876Jy sigma=0.619151 ! Fit after self-cal, rms=0.114182Jy sigma=0.343307 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.78 mas, bmaj=6.381 mas, bpa=62.37 degrees ! Estimated noise=4.73 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 (-6, 20.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.014111 Jy ! Component: 100 - total flux cleaned = 0.0174139 Jy ! Component: 150 - total flux cleaned = 0.0194539 Jy ! Component: 200 - total flux cleaned = 0.0216194 Jy ! Total flux subtracted in 200 components = 0.0216194 Jy ! Clean residual min=-0.012835 max=0.012867 Jy/beam ! Clean residual mean=-0.000251 rms=0.002881 Jy/beam ! Combined flux in latest and established models = 1.32633 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 90 components and 1.32633 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107699Jy sigma=0.319261 ! Fit after self-cal, rms=0.107059Jy sigma=0.317647 ! 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.107059Jy sigma=0.317647 ! Fit after self-cal, rms=0.106195Jy sigma=0.315668 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.782 mas, bmaj=6.391 mas, bpa=62.41 degrees ! Estimated noise=4.71584 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00244884 Jy ! Component: 100 - total flux cleaned = -0.0020648 Jy ! Component: 150 - total flux cleaned = -0.000575729 Jy ! Component: 200 - total flux cleaned = 0.000740593 Jy ! Total flux subtracted in 200 components = 0.000740593 Jy ! Clean residual min=-0.010729 max=0.011046 Jy/beam ! Clean residual mean=-0.000232 rms=0.002536 Jy/beam ! Combined flux in latest and established models = 1.32707 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 100 components and 1.32707 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103066Jy sigma=0.304985 ! Fit after self-cal, rms=0.102826Jy sigma=0.304463 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 (-11.6, 29.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00861362 Jy ! Component: 100 - total flux cleaned = 0.0123127 Jy ! Component: 150 - total flux cleaned = 0.0140209 Jy ! Component: 200 - total flux cleaned = 0.0148012 Jy ! Total flux subtracted in 200 components = 0.0148012 Jy ! Clean residual min=-0.009944 max=0.008249 Jy/beam ! Clean residual mean=-0.000190 rms=0.002283 Jy/beam ! Combined flux in latest and established models = 1.34188 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 124 components and 1.34188 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100728Jy sigma=0.298000 ! Fit after self-cal, rms=0.100582Jy sigma=0.297552 ! 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.100582Jy sigma=0.297552 ! Fit after self-cal, rms=0.099888Jy sigma=0.295654 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.784 mas, bmaj=6.403 mas, bpa=62.51 degrees ! Estimated noise=4.70967 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0060021 Jy ! Component: 100 - total flux cleaned = -0.00511032 Jy ! Component: 150 - total flux cleaned = -0.00484012 Jy ! Component: 200 - total flux cleaned = -0.0041052 Jy ! Total flux subtracted in 200 components = -0.0041052 Jy ! Clean residual min=-0.009044 max=0.007374 Jy/beam ! Clean residual mean=-0.000174 rms=0.002072 Jy/beam ! Combined flux in latest and established models = 1.33777 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 131 components and 1.33777 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098150Jy sigma=0.290333 ! Fit after self-cal, rms=0.098075Jy sigma=0.290140 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.000684429 Jy ! Component: 100 - total flux cleaned = 0.0017918 Jy ! Component: 150 - total flux cleaned = 0.00260241 Jy ! Component: 200 - total flux cleaned = 0.00316793 Jy ! Total flux subtracted in 200 components = 0.00316793 Jy ! Clean residual min=-0.008519 max=0.006730 Jy/beam ! Clean residual mean=-0.000166 rms=0.001974 Jy/beam ! Combined flux in latest and established models = 1.34094 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 144 components and 1.34094 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097193Jy sigma=0.287427 ! Fit after self-cal, rms=0.097145Jy sigma=0.287318 ! 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.097145Jy sigma=0.287318 ! Fit after self-cal, rms=0.096686Jy sigma=0.285978 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.148 mas, bmaj=10.71 mas, bpa=70.92 degrees ! Estimated noise=5.73531 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00465266 Jy ! Component: 100 - total flux cleaned = -0.00380195 Jy ! Component: 150 - total flux cleaned = -0.00299571 Jy ! Component: 200 - total flux cleaned = -0.00224323 Jy ! Total flux subtracted in 200 components = -0.00224323 Jy ! Clean residual min=-0.008450 max=0.006491 Jy/beam ! Clean residual mean=-0.000273 rms=0.002159 Jy/beam ! Combined flux in latest and established models = 1.3387 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 173 components and 1.3387 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095946Jy sigma=0.284160 ! Fit after self-cal, rms=0.095912Jy sigma=0.284026 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.000866351 Jy ! Component: 100 - total flux cleaned = 0.00143386 Jy ! Component: 150 - total flux cleaned = 0.00153493 Jy ! Component: 200 - total flux cleaned = 0.00162834 Jy ! Total flux subtracted in 200 components = 0.00162834 Jy ! Clean residual min=-0.008120 max=0.006677 Jy/beam ! Clean residual mean=-0.000264 rms=0.002129 Jy/beam ! Combined flux in latest and established models = 1.34032 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 182 components and 1.34032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.095706Jy sigma=0.283530 ! Fit after self-cal, rms=0.095694Jy sigma=0.283483 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR 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:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN 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:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA 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 5:SC,true ! 5:SC self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR 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 selfant 7:BR,true ! 7:BR self-cal status=fixed weight=1 selfant 7:FD,true ! 7:FD self-cal status=fixed weight=1 selfant 7:HN,true ! 7:HN self-cal status=fixed weight=1 selfant 7:KP,true ! 7:KP self-cal status=fixed weight=1 selfant 7:LA,true ! 7:LA self-cal status=fixed weight=1 selfant 7:MK,true ! 7:MK self-cal status=fixed weight=1 selfant 7:NL,true ! 7:NL self-cal status=fixed weight=1 selfant 7:OV,true ! 7:OV self-cal status=fixed weight=1 selfant 7:PT,true ! 7:PT self-cal status=fixed weight=1 selfant 7:SC,true ! 7:SC self-cal status=fixed weight=1 selfant 8:MK,true ! 8:MK self-cal status=fixed weight=1 selfant 9:BR,true ! 9:BR 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:MK,true ! 9:MK self-cal status=fixed weight=1 selfant 9:OV,true ! 9:OV self-cal status=fixed weight=1 selfant 9:PT,true ! 9:PT 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 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 22 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 25 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 28 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 22 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 4 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! A total of 19 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 22 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.306044Jy sigma=1.099853 ! Fit after self-cal, rms=0.190750Jy sigma=0.675020 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:ZC self-cal status=correctable weight=1 ! 2:BR 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:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:SC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:BR 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 ! 6:ZC self-cal status=correctable weight=1 ! 7:BR self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:OV self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC self-cal status=correctable weight=1 ! 8:MK self-cal status=correctable weight=1 ! 8:ZC self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:KK self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 9: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 3 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 36 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 48 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 36 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 36 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 2 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 2 telescope corrections were flagged in sub-array 8. ! A total of 36 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.190750Jy sigma=0.675020 ! Fit after self-cal, rms=0.126949Jy sigma=0.339625 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.812 mas, bmaj=2.651 mas, bpa=2.076 degrees ! Estimated noise=4.8165 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.00975868 Jy ! Component: 100 - total flux cleaned = -0.0154992 Jy ! Component: 150 - total flux cleaned = -0.0195036 Jy ! Component: 200 - total flux cleaned = -0.0226414 Jy ! Total flux subtracted in 200 components = -0.0226414 Jy ! Clean residual min=-0.010747 max=0.010313 Jy/beam ! Clean residual mean=-0.000013 rms=0.002437 Jy/beam ! Combined flux in latest and established models = 1.31768 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 195 components and 1.31768 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122161Jy sigma=0.326156 ! Fit after self-cal, rms=0.121074Jy sigma=0.323611 ! 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.286 mas, bmaj=3.19 mas, bpa=7.638 degrees ! Estimated noise=3.58751 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-9.6, 7.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000142182 Jy ! Component: 100 - total flux cleaned = 0.000130782 Jy ! Component: 150 - total flux cleaned = 0.00113328 Jy ! Component: 200 - total flux cleaned = 0.0014413 Jy ! Component: 250 - total flux cleaned = 0.0023225 Jy ! Component: 300 - total flux cleaned = 0.00287253 Jy ! Component: 350 - total flux cleaned = 0.00287324 Jy ! Component: 400 - total flux cleaned = 0.00312163 Jy ! Total flux subtracted in 400 components = 0.00312163 Jy ! Clean residual min=-0.006525 max=0.006802 Jy/beam ! Clean residual mean=-0.000084 rms=0.001622 Jy/beam ! Combined flux in latest and established models = 1.3208 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 222 components and 1.3208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117461Jy sigma=0.314733 ! Fit after self-cal, rms=0.117142Jy sigma=0.314355 ! Inverting map ! Added new window around map position (-0.8, 10). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00259735 Jy ! Component: 100 - total flux cleaned = 0.00210114 Jy ! Component: 150 - total flux cleaned = 0.00233088 Jy ! Component: 200 - total flux cleaned = 0.00233298 Jy ! Component: 250 - total flux cleaned = 0.00254288 Jy ! Component: 300 - total flux cleaned = 0.0027413 Jy ! Component: 350 - total flux cleaned = 0.00293527 Jy ! Component: 400 - total flux cleaned = 0.00293387 Jy ! Total flux subtracted in 400 components = 0.00293387 Jy ! Clean residual min=-0.006140 max=0.006416 Jy/beam ! Clean residual mean=-0.000081 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 1.32374 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 251 components and 1.32374 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115377Jy sigma=0.309951 ! Fit after self-cal, rms=0.115207Jy sigma=0.309738 ! Inverting map ! Added new window around map position (21.6, -27.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00583356 Jy ! Component: 100 - total flux cleaned = 0.0070367 Jy ! Component: 150 - total flux cleaned = 0.00760308 Jy ! Component: 200 - total flux cleaned = 0.00832273 Jy ! Component: 250 - total flux cleaned = 0.0083209 Jy ! Component: 300 - total flux cleaned = 0.00865799 Jy ! Component: 350 - total flux cleaned = 0.00865733 Jy ! Component: 400 - total flux cleaned = 0.00865948 Jy ! Total flux subtracted in 400 components = 0.00865948 Jy ! Clean residual min=-0.006079 max=0.005799 Jy/beam ! Clean residual mean=-0.000064 rms=0.001464 Jy/beam ! Combined flux in latest and established models = 1.3324 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 284 components and 1.3324 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.113916Jy sigma=0.306406 ! Fit after self-cal, rms=0.113813Jy sigma=0.306246 ! 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.00161327 Jy ! Component: 100 - total flux cleaned = 0.00209747 Jy ! Component: 150 - total flux cleaned = 0.00209393 Jy ! Component: 200 - total flux cleaned = 0.00194411 Jy ! Component: 250 - total flux cleaned = 0.00209413 Jy ! Component: 300 - total flux cleaned = 0.00224151 Jy ! Component: 350 - total flux cleaned = 0.00238191 Jy ! Component: 400 - total flux cleaned = 0.00238125 Jy ! Total flux subtracted in 400 components = 0.00238125 Jy ! Clean residual min=-0.005967 max=0.005705 Jy/beam ! Clean residual mean=-0.000060 rms=0.001422 Jy/beam ! Combined flux in latest and established models = 1.33478 Jy ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 303 components and 1.33478 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112948Jy sigma=0.304139 ! Fit after self-cal, rms=0.112899Jy sigma=0.304075 ! 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. 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 7 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 66 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 12 telescope corrections were flagged in sub-array 6. ! A total of 73 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 64 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 66 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 30 telescope corrections were flagged in sub-array 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 14 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 66 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 58 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.236370Jy sigma=0.849864 ! Fit after self-cal, rms=0.112216Jy sigma=0.300316 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.29 mas, bmaj=3.189 mas, bpa=7.894 degrees ! Estimated noise=3.58701 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00808186 Jy ! Component: 100 - total flux cleaned = -0.0106126 Jy ! Component: 150 - total flux cleaned = -0.011994 Jy ! Component: 200 - total flux cleaned = -0.0128837 Jy ! Component: 250 - total flux cleaned = -0.0132104 Jy ! Component: 300 - total flux cleaned = -0.0132148 Jy ! Component: 350 - total flux cleaned = -0.0133607 Jy ! Component: 400 - total flux cleaned = -0.0132207 Jy ! Total flux subtracted in 400 components = -0.0132207 Jy ! Clean residual min=-0.005269 max=0.004904 Jy/beam ! Clean residual mean=-0.000053 rms=0.001277 Jy/beam ! Combined flux in latest and established models = 1.32156 Jy selfcal ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 312 components and 1.32156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110664Jy sigma=0.296835 ! Fit after self-cal, rms=0.110587Jy sigma=0.296737 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 = -3.59093e-05 Jy ! Component: 100 - total flux cleaned = -0.000176808 Jy ! Component: 150 - total flux cleaned = -0.000174929 Jy ! Component: 200 - total flux cleaned = 8.14667e-05 Jy ! Component: 250 - total flux cleaned = 7.90741e-05 Jy ! Component: 300 - total flux cleaned = -4.13978e-05 Jy ! Component: 350 - total flux cleaned = 7.52406e-05 Jy ! Component: 400 - total flux cleaned = 7.64674e-05 Jy ! Total flux subtracted in 400 components = 7.64674e-05 Jy ! Clean residual min=-0.005230 max=0.004762 Jy/beam ! Clean residual mean=-0.000053 rms=0.001249 Jy/beam ! Combined flux in latest and established models = 1.32163 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 323 components and 1.32163 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109900Jy sigma=0.295356 ! Fit after self-cal, rms=0.109862Jy sigma=0.295307 ! 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.109862Jy sigma=0.295307 ! Fit after self-cal, rms=0.109816Jy sigma=0.294213 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.294 mas, bmaj=3.191 mas, bpa=8.072 degrees ! Estimated noise=3.58506 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00593002 Jy ! Component: 100 - total flux cleaned = -0.00798641 Jy ! Component: 150 - total flux cleaned = -0.00901564 Jy ! Component: 200 - total flux cleaned = -0.0100692 Jy ! Component: 250 - total flux cleaned = -0.0109233 Jy ! Component: 300 - total flux cleaned = -0.0111519 Jy ! Component: 350 - total flux cleaned = -0.0112659 Jy ! Component: 400 - total flux cleaned = -0.0114857 Jy ! Total flux subtracted in 400 components = -0.0114857 Jy ! Clean residual min=-0.004865 max=0.004458 Jy/beam ! Clean residual mean=-0.000047 rms=0.001176 Jy/beam ! Combined flux in latest and established models = 1.31015 Jy selfcal ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 341 components and 1.31015 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108920Jy sigma=0.292277 ! Fit after self-cal, rms=0.108889Jy sigma=0.292223 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.00024753 Jy ! Component: 100 - total flux cleaned = -0.000243978 Jy ! Component: 150 - total flux cleaned = -0.000344507 Jy ! Component: 200 - total flux cleaned = -0.000142376 Jy ! Component: 250 - total flux cleaned = -0.000141326 Jy ! Component: 300 - total flux cleaned = -4.57658e-05 Jy ! Component: 350 - total flux cleaned = -4.60849e-05 Jy ! Component: 400 - total flux cleaned = -4.53814e-05 Jy ! Total flux subtracted in 400 components = -4.53814e-05 Jy ! Clean residual min=-0.004824 max=0.004445 Jy/beam ! Clean residual mean=-0.000046 rms=0.001158 Jy/beam ! Combined flux in latest and established models = 1.3101 Jy ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 355 components and 1.3101 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108448Jy sigma=0.291357 ! Fit after self-cal, rms=0.108431Jy sigma=0.291330 ! 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.108431Jy sigma=0.291330 ! Fit after self-cal, rms=0.107835Jy sigma=0.288825 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.4 mas, bmaj=10.08 mas, bpa=59.36 degrees ! Estimated noise=5.2438 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00482105 Jy ! Component: 100 - total flux cleaned = -0.00530051 Jy ! Component: 150 - total flux cleaned = -0.00530193 Jy ! Component: 200 - total flux cleaned = -0.00508254 Jy ! Component: 250 - total flux cleaned = -0.00503512 Jy ! Component: 300 - total flux cleaned = -0.00507722 Jy ! Component: 350 - total flux cleaned = -0.00507826 Jy ! Component: 400 - total flux cleaned = -0.00511772 Jy ! Total flux subtracted in 400 components = -0.00511772 Jy ! Clean residual min=-0.005843 max=0.006050 Jy/beam ! Clean residual mean=-0.000122 rms=0.001726 Jy/beam ! Combined flux in latest and established models = 1.30498 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 400 components and 1.30498 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107760Jy sigma=0.288527 ! Fit after self-cal, rms=0.107743Jy sigma=0.288470 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 (-45.6, 77.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00618319 Jy ! Component: 100 - total flux cleaned = 0.00797751 Jy ! Component: 150 - total flux cleaned = 0.0086488 Jy ! Component: 200 - total flux cleaned = 0.00945284 Jy ! Component: 250 - total flux cleaned = 0.00997161 Jy ! Component: 300 - total flux cleaned = 0.0103514 Jy ! Component: 350 - total flux cleaned = 0.0108006 Jy ! Component: 400 - total flux cleaned = 0.0113096 Jy ! Total flux subtracted in 400 components = 0.0113096 Jy ! Clean residual min=-0.004983 max=0.005136 Jy/beam ! Clean residual mean=-0.000106 rms=0.001507 Jy/beam ! Combined flux in latest and established models = 1.31629 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 468 components and 1.31629 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107565Jy sigma=0.287724 ! Fit after self-cal, rms=0.107558Jy sigma=0.287627 ! 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.297 mas, bmaj=3.192 mas, bpa=8.193 degrees ! Estimated noise=3.58347 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.297 x 3.192 at 8.193 degrees (North through East) ! Clean map min=-0.0056349 max=0.90162 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.107558Jy sigma=0.287627 ! Fit after self-cal, rms=0.103872Jy sigma=0.272647 ! Inverting map and beam ! Estimated beam: bmin=2.299 mas, bmaj=3.195 mas, bpa=8.189 degrees ! Estimated noise=3.57896 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00157672 Jy ! Component: 100 - total flux cleaned = -0.00158981 Jy ! Component: 150 - total flux cleaned = -0.00159978 Jy ! Component: 200 - total flux cleaned = -0.00160032 Jy ! Component: 250 - total flux cleaned = -0.00160014 Jy ! Component: 300 - total flux cleaned = -0.00160478 Jy ! Component: 350 - total flux cleaned = -0.0017348 Jy ! Component: 400 - total flux cleaned = -0.00221964 Jy ! Total flux subtracted in 400 components = -0.00221964 Jy ! Clean residual min=-0.004976 max=0.004407 Jy/beam ! Clean residual mean=-0.000028 rms=0.000995 Jy/beam ! Combined flux in latest and established models = 1.31407 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 534 components and 1.31407 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103009Jy sigma=0.270231 ! Fit after self-cal, rms=0.102976Jy sigma=0.270177 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000117855 Jy ! Component: 100 - total flux cleaned = -4.39559e-07 Jy ! Component: 150 - total flux cleaned = -1.91472e-06 Jy ! Component: 200 - total flux cleaned = -3.63025e-06 Jy ! Component: 250 - total flux cleaned = -0.000115443 Jy ! Component: 300 - total flux cleaned = 0.000100588 Jy ! Component: 350 - total flux cleaned = -0.000111117 Jy ! Component: 400 - total flux cleaned = 0.000201621 Jy ! Total flux subtracted in 400 components = 0.000201621 Jy ! Clean residual min=-0.004890 max=0.004346 Jy/beam ! Clean residual mean=-0.000028 rms=0.000967 Jy/beam ! Combined flux in latest and established models = 1.31428 Jy ! Performing phase self-cal ! Adding 112 model components to the UV plane model. ! The established model now contains 577 components and 1.31428 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102467Jy sigma=0.268895 ! Fit after self-cal, rms=0.102462Jy sigma=0.268882 ! 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.102462Jy sigma=0.268882 ! Fit after self-cal, rms=0.102381Jy sigma=0.268631 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.816 mas, bmaj=2.648 mas, bpa=2.748 degrees ! Estimated noise=4.82742 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.707434 Jy ! Component: 100 - total flux cleaned = 0.867629 Jy ! Component: 150 - total flux cleaned = 0.930147 Jy ! Component: 200 - total flux cleaned = 0.973681 Jy ! Component: 250 - total flux cleaned = 1.00772 Jy ! Component: 300 - total flux cleaned = 1.03531 Jy ! Component: 350 - total flux cleaned = 1.05815 Jy ! Component: 400 - total flux cleaned = 1.07729 Jy ! Total flux subtracted in 400 components = 1.07729 Jy ! Clean residual min=-0.008171 max=0.011676 Jy/beam ! Clean residual mean=-0.000030 rms=0.002096 Jy/beam ! Combined flux in latest and established models = 1.07729 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 73 model components to the UV plane model. ! The established model now contains 73 components and 1.07729 Jy ! Inverting map and beam ! Estimated beam: bmin=2.301 mas, bmaj=3.196 mas, bpa=8.228 degrees ! Estimated noise=3.57506 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.0340161 Jy ! Component: 100 - total flux cleaned = 0.0571473 Jy ! Component: 150 - total flux cleaned = 0.0749521 Jy ! Component: 200 - total flux cleaned = 0.089822 Jy ! Component: 250 - total flux cleaned = 0.102591 Jy ! Component: 300 - total flux cleaned = 0.113731 Jy ! Component: 350 - total flux cleaned = 0.123531 Jy ! Component: 400 - total flux cleaned = 0.132239 Jy ! Total flux subtracted in 400 components = 0.132239 Jy ! Clean residual min=-0.005808 max=0.006966 Jy/beam ! Clean residual mean=-0.000165 rms=0.001529 Jy/beam ! Combined flux in latest and established models = 1.20953 Jy ! Adding 143 model components to the UV plane model. ! The established model now contains 192 components and 1.20953 Jy ! Inverting map ! Added new window around map position (6.4, -0.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00810806 Jy ! Component: 100 - total flux cleaned = 0.0152093 Jy ! Component: 150 - total flux cleaned = 0.0217079 Jy ! Component: 200 - total flux cleaned = 0.027701 Jy ! Component: 250 - total flux cleaned = 0.0330612 Jy ! Component: 300 - total flux cleaned = 0.0368312 Jy ! Component: 350 - total flux cleaned = 0.0401935 Jy ! Component: 400 - total flux cleaned = 0.0431865 Jy ! Total flux subtracted in 400 components = 0.0431865 Jy ! Clean residual min=-0.005079 max=0.005576 Jy/beam ! Clean residual mean=-0.000101 rms=0.001177 Jy/beam ! Combined flux in latest and established models = 1.25272 Jy ! Adding 116 model components to the UV plane model. ! The established model now contains 266 components and 1.25272 Jy ! Inverting map ! Added new window around map position (-30, 26). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00559762 Jy ! Component: 100 - total flux cleaned = 0.00859523 Jy ! Component: 150 - total flux cleaned = 0.0109309 Jy ! Component: 200 - total flux cleaned = 0.0133323 Jy ! Component: 250 - total flux cleaned = 0.0151799 Jy ! Component: 300 - total flux cleaned = 0.0171137 Jy ! Component: 350 - total flux cleaned = 0.0188432 Jy ! Component: 400 - total flux cleaned = 0.0203802 Jy ! Total flux subtracted in 400 components = 0.0203802 Jy ! Clean residual min=-0.004519 max=0.004275 Jy/beam ! Clean residual mean=-0.000064 rms=0.001036 Jy/beam ! Combined flux in latest and established models = 1.2731 Jy ! Adding 105 model components to the UV plane model. ! The established model now contains 329 components and 1.2731 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.416 mas, bmaj=10.1 mas, bpa=59.58 degrees ! Estimated noise=5.20876 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00692321 Jy ! Component: 100 - total flux cleaned = 0.0122883 Jy ! Component: 150 - total flux cleaned = 0.0165402 Jy ! Component: 200 - total flux cleaned = 0.0195845 Jy ! Component: 250 - total flux cleaned = 0.0217609 Jy ! Component: 300 - total flux cleaned = 0.023624 Jy ! Component: 350 - total flux cleaned = 0.0254481 Jy ! Component: 400 - total flux cleaned = 0.0270393 Jy ! Total flux subtracted in 400 components = 0.0270393 Jy ! Clean residual min=-0.004796 max=0.005387 Jy/beam ! Clean residual mean=-0.000087 rms=0.001370 Jy/beam ! Combined flux in latest and established models = 1.30014 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) ! Added new window around map position (50, 46). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00326011 Jy ! Component: 100 - total flux cleaned = 0.00364709 Jy ! Component: 150 - total flux cleaned = 0.00502921 Jy ! Component: 200 - total flux cleaned = 0.00640339 Jy ! Component: 250 - total flux cleaned = 0.00748903 Jy ! Component: 300 - total flux cleaned = 0.00842578 Jy ! Component: 350 - total flux cleaned = 0.00913987 Jy ! Component: 400 - total flux cleaned = 0.010081 Jy ! Total flux subtracted in 400 components = 0.010081 Jy ! Clean residual min=-0.004129 max=0.003517 Jy/beam ! Clean residual mean=-0.000093 rms=0.001179 Jy/beam ! Combined flux in latest and established models = 1.31022 Jy ! Adding 181 model components to the UV plane model. ! The established model now contains 500 components and 1.31022 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.301 mas, bmaj=3.196 mas, bpa=8.228 degrees ! Estimated noise=3.57506 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000570443 Jy ! Component: 100 - total flux cleaned = 0.000424265 Jy ! Component: 150 - total flux cleaned = 0.000562368 Jy ! Component: 200 - total flux cleaned = 0.000819158 Jy ! Component: 250 - total flux cleaned = 0.000817895 Jy ! Component: 300 - total flux cleaned = 0.000576162 Jy ! Component: 350 - total flux cleaned = 0.000695299 Jy ! Component: 400 - total flux cleaned = 0.000695225 Jy ! Total flux subtracted in 400 components = 0.000695225 Jy ! Clean residual min=-0.003564 max=0.003768 Jy/beam ! Clean residual mean=-0.000028 rms=0.000864 Jy/beam ! Combined flux in latest and established models = 1.31091 Jy selfcal ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 571 components and 1.31091 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100830Jy sigma=0.265257 ! Fit after self-cal, rms=0.100731Jy sigma=0.264757 wmodel /scr/pima/rdv73_uvs/J0921+6215_S_map.mod ! Writing 571 model components to file: /scr/pima/rdv73_uvs/J0921+6215_S_map.mod wobs /scr/pima/rdv73_uvs/J0921+6215_S_uvs.fits ! Applying 1968 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J0921+6215_S_uvs.fits wwins /scr/pima/rdv73_uvs/J0921+6215_S_map.win ! wwins: Wrote 11 windows to /scr/pima/rdv73_uvs/J0921+6215_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.000999169 Jy ! Component: 100 - total flux cleaned = -0.00220936 Jy ! Component: 150 - total flux cleaned = -0.00306875 Jy ! Component: 200 - total flux cleaned = -0.0034873 Jy ! Component: 250 - total flux cleaned = -0.00383901 Jy ! Component: 300 - total flux cleaned = -0.00412725 Jy ! Component: 350 - total flux cleaned = -0.00441098 Jy ! Component: 400 - total flux cleaned = -0.00480239 Jy ! Component: 450 - total flux cleaned = -0.00513344 Jy ! Component: 500 - total flux cleaned = -0.00529754 Jy ! Component: 550 - total flux cleaned = -0.00572879 Jy ! Component: 600 - total flux cleaned = -0.00578233 Jy ! Component: 650 - total flux cleaned = -0.00615211 Jy ! Component: 700 - total flux cleaned = -0.00625587 Jy ! Component: 750 - total flux cleaned = -0.00635917 Jy ! Component: 800 - total flux cleaned = -0.00641057 Jy ! Component: 850 - total flux cleaned = -0.00666597 Jy ! Component: 900 - total flux cleaned = -0.00676678 Jy ! Component: 950 - total flux cleaned = -0.00686779 Jy ! Component: 1000 - total flux cleaned = -0.00711691 Jy ! Component: 1050 - total flux cleaned = -0.00716679 Jy ! Component: 1100 - total flux cleaned = -0.00746208 Jy ! Component: 1150 - total flux cleaned = -0.00760843 Jy ! Component: 1200 - total flux cleaned = -0.00770529 Jy ! Component: 1250 - total flux cleaned = -0.0077536 Jy ! Component: 1300 - total flux cleaned = -0.00799279 Jy ! Component: 1350 - total flux cleaned = -0.00794544 Jy ! Component: 1400 - total flux cleaned = -0.00822958 Jy ! Component: 1450 - total flux cleaned = -0.00827676 Jy ! Component: 1500 - total flux cleaned = -0.00846366 Jy ! Component: 1550 - total flux cleaned = -0.00860263 Jy ! Component: 1600 - total flux cleaned = -0.00864902 Jy ! Component: 1650 - total flux cleaned = -0.00901601 Jy ! Component: 1700 - total flux cleaned = -0.00910729 Jy ! Component: 1750 - total flux cleaned = -0.00915225 Jy ! Component: 1800 - total flux cleaned = -0.00933274 Jy ! Component: 1850 - total flux cleaned = -0.00951241 Jy ! Component: 1900 - total flux cleaned = -0.00969104 Jy ! Component: 1950 - total flux cleaned = -0.00986881 Jy ! Component: 2000 - total flux cleaned = -0.0100896 Jy ! Component: 2050 - total flux cleaned = -0.0102215 Jy ! Component: 2100 - total flux cleaned = -0.0105276 Jy ! Component: 2150 - total flux cleaned = -0.0106148 Jy ! Component: 2200 - total flux cleaned = -0.0107884 Jy ! Component: 2250 - total flux cleaned = -0.0110468 Jy ! Component: 2300 - total flux cleaned = -0.0111754 Jy ! Component: 2350 - total flux cleaned = -0.0112606 Jy ! Component: 2400 - total flux cleaned = -0.0115158 Jy ! Component: 2450 - total flux cleaned = -0.011389 Jy ! Component: 2500 - total flux cleaned = -0.0117255 Jy ! Component: 2550 - total flux cleaned = -0.0118089 Jy ! Component: 2600 - total flux cleaned = -0.0118925 Jy ! Component: 2650 - total flux cleaned = -0.0121 Jy ! Component: 2700 - total flux cleaned = -0.0120997 Jy ! Component: 2750 - total flux cleaned = -0.0122646 Jy ! Component: 2800 - total flux cleaned = -0.0124283 Jy ! Component: 2850 - total flux cleaned = -0.0126732 Jy ! Component: 2900 - total flux cleaned = -0.0125917 Jy ! Component: 2950 - total flux cleaned = -0.0127941 Jy ! Component: 3000 - total flux cleaned = -0.0129148 Jy ! Component: 3050 - total flux cleaned = -0.0131152 Jy ! Component: 3100 - total flux cleaned = -0.0131156 Jy ! Component: 3150 - total flux cleaned = -0.0134741 Jy ! Component: 3200 - total flux cleaned = -0.0136723 Jy ! Component: 3250 - total flux cleaned = -0.0136722 Jy ! Component: 3300 - total flux cleaned = -0.0138295 Jy ! Component: 3350 - total flux cleaned = -0.014143 Jy ! Component: 3400 - total flux cleaned = -0.0141817 Jy ! Component: 3450 - total flux cleaned = -0.0144151 Jy ! Component: 3500 - total flux cleaned = -0.014415 Jy ! Component: 3550 - total flux cleaned = -0.0144151 Jy ! Component: 3600 - total flux cleaned = -0.0145304 Jy ! Component: 3650 - total flux cleaned = -0.0147219 Jy ! Component: 3700 - total flux cleaned = -0.014989 Jy ! Component: 3750 - total flux cleaned = -0.0149506 Jy ! Component: 3800 - total flux cleaned = -0.0151779 Jy ! Component: 3850 - total flux cleaned = -0.0152911 Jy ! Component: 3900 - total flux cleaned = -0.0154791 Jy ! Component: 3950 - total flux cleaned = -0.0156287 Jy ! Component: 4000 - total flux cleaned = -0.0156286 Jy ! Component: 4050 - total flux cleaned = -0.0157772 Jy ! Total flux subtracted in 4096 components = -0.0160367 Jy ! Clean residual min=-0.002409 max=0.002001 Jy/beam ! Clean residual mean=-0.000003 rms=0.000647 Jy/beam ! Combined flux in latest and established models = 1.29488 Jy keep ! Adding 292 model components to the UV plane model. ! The established model now contains 859 components and 1.29488 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J0921+6215_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.301 x 3.196 at 8.228 degrees (North through East) ! Clean map min=-0.003167 max=0.90722 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J0921+6215_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.9072 rms= 0.00064 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J0921+6215_S_uvs.log closed on Sat Jan 9 17:48:04 2016