! Started logfile: /scr/pima/rdv62_uvs/J2331-1556_S_uvs.log on Fri Nov 22 22:42:52 2013 obs /scr/pima/rdv62_uvs/J2331-1556_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J2331-1556_S_uva.fits ! AN table 1: 8 integrations on 6 of 6 possible baselines. ! AN table 2: 11 integrations on 21 of 21 possible baselines. ! AN table 3: 13 integrations on 6 of 6 possible baselines. ! AN table 4: 11 integrations on 1 of 1 possible baselines. ! AN table 5: 28 integrations on 3 of 3 possible baselines. ! AN table 6: 13 integrations on 3 of 3 possible baselines. ! AN table 7: 12 integrations on 15 of 15 possible baselines. ! AN table 8: 16 integrations on 21 of 21 possible baselines. ! AN table 9: 6 integrations on 21 of 21 possible baselines. ! AN table 10: 11 integrations on 36 of 36 possible baselines. ! AN table 11: 35 integrations on 3 of 3 possible baselines. ! AN table 12: 17 integrations on 36 of 36 possible baselines. ! Apparent sampling: 0.567676 visibilities/baseline/integration-bin. ! Found source: J2331-1556 ! ! 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 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 5100 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/rdv62_uvs/J2331-1556_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:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:SC self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:NL self-cal status=correctable weight=1 ! 11:SC self-cal status=correctable weight=1 ! 12:FD self-cal status=correctable weight=1 ! 12:HN self-cal status=correctable weight=1 ! 12:KP self-cal status=correctable weight=1 ! 12:LA self-cal status=correctable weight=1 ! 12:MK self-cal status=correctable weight=1 ! 12:NL self-cal status=correctable weight=1 ! 12:OV self-cal status=correctable weight=1 ! 12:PT self-cal status=correctable weight=1 ! 12:SC 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 9 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 22 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 23 telescope corrections were flagged in sub-array 10. ! A total of 45 telescope corrections were flagged in sub-array 11. ! A total of 76 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 22 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 23 telescope corrections were flagged in sub-array 10. ! A total of 45 telescope corrections were flagged in sub-array 11. ! A total of 76 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 22 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 23 telescope corrections were flagged in sub-array 10. ! A total of 45 telescope corrections were flagged in sub-array 11. ! A total of 76 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 30 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 22 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 23 telescope corrections were flagged in sub-array 10. ! A total of 45 telescope corrections were flagged in sub-array 11. ! A total of 76 telescope corrections were flagged in sub-array 12. ! ! Fit before self-cal, rms=1.388616Jy sigma=11.854734 ! Fit after self-cal, rms=0.331709Jy sigma=2.098630 ! 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 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 5:HH ! Flagging baselines of antenna 5:HH 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:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.539 mas, bmaj=8.727 mas, bpa=-0.4043 degrees ! Estimated noise=2.39007 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.666618 Jy ! Component: 100 - total flux cleaned = 0.879306 Jy ! Total flux subtracted in 100 components = 0.879306 Jy ! Clean residual min=-0.030385 max=0.097357 Jy/beam ! Clean residual mean=-0.000230 rms=0.009708 Jy/beam ! Combined flux in latest and established models = 0.879306 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.879306 Jy ! ! Correcting IF 1. ! A total of 11 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 11 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 11 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=0.254523Jy sigma=2.311932 ! Fit after self-cal, rms=0.255021Jy sigma=2.311495 ! Inverting map and beam ! Estimated beam: bmin=3.539 mas, bmaj=8.727 mas, bpa=-0.4043 degrees ! Estimated noise=2.39007 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.106738 Jy ! Component: 100 - total flux cleaned = 0.161633 Jy ! Total flux subtracted in 100 components = 0.161633 Jy ! Clean residual min=-0.021346 max=0.026587 Jy/beam ! Clean residual mean=-0.000059 rms=0.005779 Jy/beam ! Combined flux in latest and established models = 1.04094 Jy ! Performing phase self-cal ! Adding 11 model components to the UV plane model. ! The established model now contains 21 components and 1.04094 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.191651Jy sigma=1.667528 ! Fit after self-cal, rms=0.191758Jy sigma=1.667279 ! 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=4.696 mas, bmaj=9.922 mas, bpa=-1.723 degrees ! Estimated noise=1.34193 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.0535407 Jy ! Component: 100 - total flux cleaned = 0.0796645 Jy ! Component: 150 - total flux cleaned = 0.0937909 Jy ! Component: 200 - total flux cleaned = 0.1017 Jy ! Total flux subtracted in 200 components = 0.1017 Jy ! Clean residual min=-0.024351 max=0.036502 Jy/beam ! Clean residual mean=0.000236 rms=0.007758 Jy/beam ! Combined flux in latest and established models = 1.14264 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 77 components and 1.14264 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.179701Jy sigma=1.527335 ! Fit after self-cal, rms=0.179633Jy sigma=1.527230 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 HN 0.95 NL 1.00 SC 0.98 ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 1.00 KP 1.01 LA 1.11 ! MK 0.94 OV 1.02 PT 1.04 ! ! A total of 33 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* WZ 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.94 HN 0.96 KP 1.16 LA 1.25 ! PT 0.97 SC 0.98 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.91 HN 0.93 KP 1.07 LA 1.17 ! NL 1.02 PT 1.00 SC 0.93 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.93 KP 0.98 LA 1.13 NL 1.06 ! OV 0.98 PT 1.04 SC 0.93 ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 0.97 FD 0.96 HN 0.95 KP 1.04 ! LA 1.11 NL 1.05 OV 1.03 PT 1.02 ! SC 0.88 ! ! A total of 60 telescope corrections were flagged in sub-array 11. ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NL 0.00* SC 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 12. ! Telescope amplitude corrections in sub-array 12: ! FD 1.05 HN 0.92 KP 1.00 LA 1.16 ! MK 1.01 NL 1.01 OV 1.00 PT 1.11 ! SC 0.90 ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 HN 1.08 NL 1.05 SC 1.03 ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 1.08 KP 1.05 LA 1.03 ! MK 1.02 OV 1.04 PT 1.05 ! ! A total of 33 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* WZ 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 1.00 HN 1.12 KP 1.14 LA 1.12 ! PT 1.01 SC 1.07 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.96 HN 1.05 KP 1.11 LA 1.06 ! NL 1.05 PT 1.01 SC 0.98 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.99 KP 1.03 LA 1.02 NL 1.04 ! OV 1.05 PT 1.04 SC 1.01 ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.01 FD 1.04 HN 1.07 KP 1.06 ! LA 1.02 NL 1.05 OV 1.10 PT 1.01 ! SC 0.97 ! ! A total of 60 telescope corrections were flagged in sub-array 11. ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NL 0.00* SC 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 12. ! Telescope amplitude corrections in sub-array 12: ! FD 1.10 HN 1.03 KP 1.04 LA 1.06 ! MK 1.05 NL 1.03 OV 1.06 PT 1.12 ! SC 1.01 ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 HN 1.01 NL 1.09 SC 1.04 ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 1.11 KP 1.03 LA 0.83 ! MK 1.04 OV 1.10 PT 1.04 ! ! A total of 33 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* WZ 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 1.04 HN 1.06 KP 1.07 LA 0.92 ! PT 1.02 SC 1.04 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.97 HN 1.02 KP 1.04 LA 0.87 ! NL 1.11 PT 1.03 SC 1.03 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.99 KP 1.01 LA 0.83 NL 1.16 ! OV 1.11 PT 1.09 SC 1.01 ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.10 FD 1.03 HN 1.03 KP 1.06 ! LA 0.84 NL 1.13 OV 1.13 PT 1.05 ! SC 1.00 ! ! A total of 60 telescope corrections were flagged in sub-array 11. ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NL 0.00* SC 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 12. ! Telescope amplitude corrections in sub-array 12: ! FD 1.15 HN 0.99 KP 1.03 LA 0.83 ! MK 1.12 NL 1.11 OV 1.11 PT 1.13 ! SC 1.00 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 HN 0.98 NL 0.99 SC 1.16 ! ! A total of 3 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 1.07 KP 0.97 LA 0.71 ! MK 1.08 OV 0.88 PT 0.99 ! ! A total of 33 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! KK 0.00* KP 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HN 0.00* SC 0.00* WZ 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.98 HN 0.97 KP 1.05 LA 0.78 ! PT 1.02 SC 1.20 ! ! Telescope amplitude corrections in sub-array 8: ! FD 0.96 HN 0.94 KP 0.99 LA 0.72 ! NL 1.01 PT 0.99 SC 1.11 ! ! Telescope amplitude corrections in sub-array 9: ! FD 0.96 KP 0.95 LA 0.74 NL 0.99 ! OV 0.96 PT 1.01 SC 1.19 ! ! A total of 6 telescope corrections were flagged in sub-array 10. ! Telescope amplitude corrections in sub-array 10: ! BR 1.00 FD 0.98 HN 0.95 KP 0.95 ! LA 0.73 NL 0.99 OV 0.96 PT 1.01 ! SC 1.10 ! ! A total of 60 telescope corrections were flagged in sub-array 11. ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NL 0.00* SC 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 12. ! Telescope amplitude corrections in sub-array 12: ! FD 1.13 HN 0.94 KP 0.94 LA 0.74 ! MK 1.09 NL 1.02 OV 0.88 PT 1.10 ! SC 1.11 ! ! ! Fit before self-cal, rms=0.179633Jy sigma=1.527230 ! Fit after self-cal, rms=0.132268Jy sigma=0.973205 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.788 mas, bmaj=10.88 mas, bpa=-4.025 degrees ! Estimated noise=1.42205 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 (52.8, -1.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0354016 Jy ! Component: 100 - total flux cleaned = 0.0522371 Jy ! Component: 150 - total flux cleaned = 0.0610826 Jy ! Component: 200 - total flux cleaned = 0.0642673 Jy ! Total flux subtracted in 200 components = 0.0642673 Jy ! Clean residual min=-0.011524 max=0.018532 Jy/beam ! Clean residual mean=0.000098 rms=0.003678 Jy/beam ! Combined flux in latest and established models = 1.20691 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 156 components and 1.20691 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123437Jy sigma=0.872724 ! Fit after self-cal, rms=0.121576Jy sigma=0.854353 ! Inverting map ! Added new window around map position (12.4, 3.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0194783 Jy ! Component: 100 - total flux cleaned = 0.0313968 Jy ! Component: 150 - total flux cleaned = 0.0371129 Jy ! Component: 200 - total flux cleaned = 0.0406161 Jy ! Total flux subtracted in 200 components = 0.0406161 Jy ! Clean residual min=-0.008862 max=0.011293 Jy/beam ! Clean residual mean=0.000142 rms=0.002702 Jy/beam ! Combined flux in latest and established models = 1.24752 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 209 components and 1.24752 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117946Jy sigma=0.811280 ! Fit after self-cal, rms=0.117032Jy sigma=0.805226 ! 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.117032Jy sigma=0.805226 ! Fit after self-cal, rms=0.115113Jy sigma=0.788921 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.788 mas, bmaj=10.88 mas, bpa=-3.977 degrees ! Estimated noise=1.42148 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0128371 Jy ! Component: 100 - total flux cleaned = 0.0208957 Jy ! Component: 150 - total flux cleaned = 0.024548 Jy ! Component: 200 - total flux cleaned = 0.026938 Jy ! Total flux subtracted in 200 components = 0.026938 Jy ! Clean residual min=-0.006746 max=0.007923 Jy/beam ! Clean residual mean=0.000088 rms=0.002054 Jy/beam ! Combined flux in latest and established models = 1.27446 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 252 components and 1.27446 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112747Jy sigma=0.767363 ! Fit after self-cal, rms=0.112167Jy sigma=0.763693 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.00728999 Jy ! Component: 100 - total flux cleaned = 0.00962686 Jy ! Component: 150 - total flux cleaned = 0.0116021 Jy ! Component: 200 - total flux cleaned = 0.0126417 Jy ! Total flux subtracted in 200 components = 0.0126417 Jy ! Clean residual min=-0.006334 max=0.007207 Jy/beam ! Clean residual mean=0.000095 rms=0.001849 Jy/beam ! Combined flux in latest and established models = 1.2871 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 281 components and 1.2871 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.110898Jy sigma=0.755348 ! Fit after self-cal, rms=0.110702Jy sigma=0.754048 ! 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.110702Jy sigma=0.754048 ! Fit after self-cal, rms=0.110376Jy sigma=0.749969 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.789 mas, bmaj=10.88 mas, bpa=-3.942 degrees ! Estimated noise=1.4217 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00757921 Jy ! Component: 100 - total flux cleaned = 0.0101846 Jy ! Component: 150 - total flux cleaned = 0.0118254 Jy ! Component: 200 - total flux cleaned = 0.0126293 Jy ! Total flux subtracted in 200 components = 0.0126293 Jy ! Clean residual min=-0.005536 max=0.006354 Jy/beam ! Clean residual mean=0.000061 rms=0.001575 Jy/beam ! Combined flux in latest and established models = 1.29973 Jy selfcal ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 304 components and 1.29973 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109233Jy sigma=0.742349 ! Fit after self-cal, rms=0.109066Jy sigma=0.741132 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.00418937 Jy ! Component: 100 - total flux cleaned = 0.00579198 Jy ! Component: 150 - total flux cleaned = 0.00706747 Jy ! Component: 200 - total flux cleaned = 0.00783974 Jy ! Total flux subtracted in 200 components = 0.00783974 Jy ! Clean residual min=-0.005476 max=0.006001 Jy/beam ! Clean residual mean=0.000065 rms=0.001472 Jy/beam ! Combined flux in latest and established models = 1.30757 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 314 components and 1.30757 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.108299Jy sigma=0.737705 ! Fit after self-cal, rms=0.108217Jy sigma=0.737167 ! 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.108217Jy sigma=0.737167 ! Fit after self-cal, rms=0.108088Jy sigma=0.735603 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.804 mas, bmaj=15.46 mas, bpa=-6.743 degrees ! Estimated noise=1.43409 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00621671 Jy ! Component: 100 - total flux cleaned = 0.00883983 Jy ! Component: 150 - total flux cleaned = 0.0108275 Jy ! Component: 200 - total flux cleaned = 0.0121823 Jy ! Total flux subtracted in 200 components = 0.0121823 Jy ! Clean residual min=-0.004639 max=0.005251 Jy/beam ! Clean residual mean=0.000077 rms=0.001490 Jy/beam ! Combined flux in latest and established models = 1.31975 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 345 components and 1.31975 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107890Jy sigma=0.732076 ! Fit after self-cal, rms=0.107735Jy sigma=0.731427 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.00310018 Jy ! Component: 100 - total flux cleaned = 0.00422734 Jy ! Component: 150 - total flux cleaned = 0.00504872 Jy ! Component: 200 - total flux cleaned = 0.00560683 Jy ! Total flux subtracted in 200 components = 0.00560683 Jy ! Clean residual min=-0.004551 max=0.005021 Jy/beam ! Clean residual mean=0.000080 rms=0.001429 Jy/beam ! Combined flux in latest and established models = 1.32536 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 371 components and 1.32536 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107697Jy sigma=0.730311 ! Fit after self-cal, rms=0.107647Jy sigma=0.730098 ! 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 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:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:SC,true ! 6:SC 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:PT,true ! 7:PT self-cal status=fixed weight=1 selfant 7:SC,true ! 7:SC self-cal status=fixed weight=1 selfant 8:FD,true ! 8:FD self-cal status=fixed weight=1 selfant 8:HN,true ! 8:HN self-cal status=fixed weight=1 selfant 8:KP,true ! 8:KP self-cal status=fixed weight=1 selfant 8:LA,true ! 8:LA self-cal status=fixed weight=1 selfant 8:NL,true ! 8:NL self-cal status=fixed weight=1 selfant 8:PT,true ! 8:PT self-cal status=fixed weight=1 selfant 8:SC,true ! 8:SC self-cal status=fixed weight=1 selfant 9:FD,true ! 9:FD self-cal status=fixed weight=1 selfant 9:KP,true ! 9:KP self-cal status=fixed weight=1 selfant 9:LA,true ! 9:LA self-cal status=fixed weight=1 selfant 9:NL,true ! 9:NL self-cal status=fixed weight=1 selfant 9:OV,true ! 9:OV self-cal status=fixed weight=1 selfant 9:PT,true ! 9:PT self-cal status=fixed weight=1 selfant 9:SC,true ! 9:SC self-cal status=fixed weight=1 selfant 10:BR,true ! 10:BR self-cal status=fixed weight=1 selfant 10:FD,true ! 10:FD self-cal status=fixed weight=1 selfant 10:HN,true ! 10:HN self-cal status=fixed weight=1 selfant 10:KP,true ! 10:KP self-cal status=fixed weight=1 selfant 10:LA,true ! 10:LA self-cal status=fixed weight=1 selfant 10:NL,true ! 10:NL self-cal status=fixed weight=1 selfant 10:OV,true ! 10:OV self-cal status=fixed weight=1 selfant 10:PT,true ! 10:PT self-cal status=fixed weight=1 selfant 10:SC,true ! 10:SC self-cal status=fixed weight=1 selfant 11:HN,true ! 11:HN self-cal status=fixed weight=1 selfant 11:NL,true ! 11:NL self-cal status=fixed weight=1 selfant 11:SC,true ! 11:SC self-cal status=fixed weight=1 selfant 12:FD,true ! 12:FD self-cal status=fixed weight=1 selfant 12:HN,true ! 12:HN self-cal status=fixed weight=1 selfant 12:KP,true ! 12:KP self-cal status=fixed weight=1 selfant 12:LA,true ! 12:LA self-cal status=fixed weight=1 selfant 12:MK,true ! 12:MK self-cal status=fixed weight=1 selfant 12:NL,true ! 12:NL self-cal status=fixed weight=1 selfant 12:OV,true ! 12:OV self-cal status=fixed weight=1 selfant 12:PT,true ! 12:PT self-cal status=fixed weight=1 selfant 12:SC,true ! 12:SC self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.303147Jy sigma=1.846349 ! Fit after self-cal, rms=0.253209Jy sigma=1.651821 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:SC self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:MK self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:KP self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:PT self-cal status=correctable weight=1 ! 7:SC self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:PT self-cal status=correctable weight=1 ! 8:SC self-cal status=correctable weight=1 ! 9:FD self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:LA self-cal status=correctable weight=1 ! 9:NL self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 ! 9:PT self-cal status=correctable weight=1 ! 9:SC self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:HN self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:NL self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 ! 10:SC self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:NL self-cal status=correctable weight=1 ! 11:SC self-cal status=correctable weight=1 ! 12:FD self-cal status=correctable weight=1 ! 12:HN self-cal status=correctable weight=1 ! 12:KP self-cal status=correctable weight=1 ! 12:LA self-cal status=correctable weight=1 ! 12:MK self-cal status=correctable weight=1 ! 12:NL self-cal status=correctable weight=1 ! 12:OV self-cal status=correctable weight=1 ! 12:PT self-cal status=correctable weight=1 ! 12:SC 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 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Fit before self-cal, rms=0.253209Jy sigma=1.651821 ! Fit after self-cal, rms=0.122714Jy sigma=0.821645 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.528 mas, bmaj=8.831 mas, bpa=-0.3912 degrees ! Estimated noise=2.47331 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.00088299 Jy ! Component: 100 - total flux cleaned = 0.00186916 Jy ! Component: 150 - total flux cleaned = 0.00300437 Jy ! Component: 200 - total flux cleaned = 0.00364762 Jy ! Total flux subtracted in 200 components = 0.00364762 Jy ! Clean residual min=-0.007812 max=0.007932 Jy/beam ! Clean residual mean=0.000009 rms=0.001996 Jy/beam ! Combined flux in latest and established models = 1.32901 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 401 components and 1.32901 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121934Jy sigma=0.821109 ! Fit after self-cal, rms=0.122005Jy sigma=0.819978 ! 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=4.677 mas, bmaj=11.1 mas, bpa=-4.737 degrees ! Estimated noise=1.3863 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 (28.4, -1.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00399887 Jy ! Component: 100 - total flux cleaned = 0.0056362 Jy ! Component: 150 - total flux cleaned = 0.00679301 Jy ! Component: 200 - total flux cleaned = 0.007615 Jy ! Component: 250 - total flux cleaned = 0.00827506 Jy ! Component: 300 - total flux cleaned = 0.00879255 Jy ! Component: 350 - total flux cleaned = 0.0091862 Jy ! Component: 400 - total flux cleaned = 0.00946831 Jy ! Total flux subtracted in 400 components = 0.00946831 Jy ! Clean residual min=-0.003802 max=0.003306 Jy/beam ! Clean residual mean=0.000048 rms=0.000993 Jy/beam ! Combined flux in latest and established models = 1.33848 Jy ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 466 components and 1.33848 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121006Jy sigma=0.816397 ! Fit after self-cal, rms=0.120995Jy sigma=0.816107 ! 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.0024172 Jy ! Component: 100 - total flux cleaned = 0.00347867 Jy ! Component: 150 - total flux cleaned = 0.00409847 Jy ! Component: 200 - total flux cleaned = 0.00443709 Jy ! Component: 250 - total flux cleaned = 0.00468194 Jy ! Component: 300 - total flux cleaned = 0.00492201 Jy ! Component: 350 - total flux cleaned = 0.00500187 Jy ! Component: 400 - total flux cleaned = 0.00522846 Jy ! Total flux subtracted in 400 components = 0.00522846 Jy ! Clean residual min=-0.003751 max=0.003244 Jy/beam ! Clean residual mean=0.000053 rms=0.000939 Jy/beam ! Combined flux in latest and established models = 1.3437 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 505 components and 1.3437 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120368Jy sigma=0.814815 ! Fit after self-cal, rms=0.120361Jy sigma=0.814716 ! 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 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 33 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 84 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 7. ! A total of 17 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! A total of 29 telescope corrections were flagged in sub-array 10. ! A total of 105 telescope corrections were flagged in sub-array 11. ! A total of 91 telescope corrections were flagged in sub-array 12. ! ! Fit before self-cal, rms=0.301872Jy sigma=1.817217 ! Fit after self-cal, rms=0.119882Jy sigma=0.813953 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.675 mas, bmaj=11.1 mas, bpa=-4.731 degrees ! Estimated noise=1.38513 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00271913 Jy ! Component: 100 - total flux cleaned = 0.00402499 Jy ! Component: 150 - total flux cleaned = 0.00462306 Jy ! Component: 200 - total flux cleaned = 0.00494546 Jy ! Component: 250 - total flux cleaned = 0.00525335 Jy ! Component: 300 - total flux cleaned = 0.00532626 Jy ! Component: 350 - total flux cleaned = 0.00539891 Jy ! Component: 400 - total flux cleaned = 0.0055386 Jy ! Total flux subtracted in 400 components = 0.0055386 Jy ! Clean residual min=-0.003364 max=0.002998 Jy/beam ! Clean residual mean=0.000038 rms=0.000855 Jy/beam ! Combined flux in latest and established models = 1.34924 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 529 components and 1.34924 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119329Jy sigma=0.812709 ! Fit after self-cal, rms=0.119314Jy sigma=0.812606 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.00154473 Jy ! Component: 100 - total flux cleaned = 0.00218816 Jy ! Component: 150 - total flux cleaned = 0.00239301 Jy ! Component: 200 - total flux cleaned = 0.00259101 Jy ! Component: 250 - total flux cleaned = 0.00271944 Jy ! Component: 300 - total flux cleaned = 0.00278196 Jy ! Component: 350 - total flux cleaned = 0.00278258 Jy ! Component: 400 - total flux cleaned = 0.00284585 Jy ! Total flux subtracted in 400 components = 0.00284585 Jy ! Clean residual min=-0.003262 max=0.002929 Jy/beam ! Clean residual mean=0.000040 rms=0.000831 Jy/beam ! Combined flux in latest and established models = 1.35209 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 552 components and 1.35209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118810Jy sigma=0.812032 ! Fit after self-cal, rms=0.118808Jy sigma=0.811990 ! 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.118808Jy sigma=0.811990 ! Fit after self-cal, rms=0.118666Jy sigma=0.811776 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.675 mas, bmaj=11.1 mas, bpa=-4.722 degrees ! Estimated noise=1.38488 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00168537 Jy ! Component: 100 - total flux cleaned = 0.00233298 Jy ! Component: 150 - total flux cleaned = 0.00266854 Jy ! Component: 200 - total flux cleaned = 0.00286143 Jy ! Component: 250 - total flux cleaned = 0.00298572 Jy ! Component: 300 - total flux cleaned = 0.00304564 Jy ! Component: 350 - total flux cleaned = 0.00304684 Jy ! Component: 400 - total flux cleaned = 0.00316655 Jy ! Total flux subtracted in 400 components = 0.00316655 Jy ! Clean residual min=-0.003114 max=0.002913 Jy/beam ! Clean residual mean=0.000034 rms=0.000796 Jy/beam ! Combined flux in latest and established models = 1.35526 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 565 components and 1.35526 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.118201Jy sigma=0.811175 ! Fit after self-cal, rms=0.118204Jy sigma=0.811127 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.00107249 Jy ! Component: 100 - total flux cleaned = 0.00137306 Jy ! Component: 150 - total flux cleaned = 0.00148704 Jy ! Component: 200 - total flux cleaned = 0.00148591 Jy ! Component: 250 - total flux cleaned = 0.00148442 Jy ! Component: 300 - total flux cleaned = 0.00148443 Jy ! Component: 350 - total flux cleaned = 0.00159546 Jy ! Component: 400 - total flux cleaned = 0.00159571 Jy ! Total flux subtracted in 400 components = 0.00159571 Jy ! Clean residual min=-0.003026 max=0.002839 Jy/beam ! Clean residual mean=0.000035 rms=0.000781 Jy/beam ! Combined flux in latest and established models = 1.35685 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 583 components and 1.35685 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117776Jy sigma=0.810773 ! Fit after self-cal, rms=0.117784Jy sigma=0.810751 ! 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.117784Jy sigma=0.810751 ! Fit after self-cal, rms=0.119781Jy sigma=0.739179 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.699 mas, bmaj=15.59 mas, bpa=-6.571 degrees ! Estimated noise=1.40751 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00164399 Jy ! Component: 100 - total flux cleaned = 0.00200487 Jy ! Component: 150 - total flux cleaned = 0.00224536 Jy ! Component: 200 - total flux cleaned = 0.00248713 Jy ! Component: 250 - total flux cleaned = 0.00259356 Jy ! Component: 300 - total flux cleaned = 0.00272278 Jy ! Component: 350 - total flux cleaned = 0.00281152 Jy ! Component: 400 - total flux cleaned = 0.0029738 Jy ! Total flux subtracted in 400 components = 0.0029738 Jy ! Clean residual min=-0.002884 max=0.002826 Jy/beam ! Clean residual mean=0.000056 rms=0.000932 Jy/beam ! Combined flux in latest and established models = 1.35982 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 642 components and 1.35982 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119798Jy sigma=0.738549 ! Fit after self-cal, rms=0.119842Jy sigma=0.738470 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (88.4, -61.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00140989 Jy ! Component: 100 - total flux cleaned = 0.00213345 Jy ! Component: 150 - total flux cleaned = 0.00320676 Jy ! Component: 200 - total flux cleaned = 0.00381141 Jy ! Component: 250 - total flux cleaned = 0.00423234 Jy ! Component: 300 - total flux cleaned = 0.0045581 Jy ! Component: 350 - total flux cleaned = 0.00481021 Jy ! Component: 400 - total flux cleaned = 0.00500075 Jy ! Total flux subtracted in 400 components = 0.00500075 Jy ! Clean residual min=-0.003031 max=0.002815 Jy/beam ! Clean residual mean=0.000050 rms=0.000853 Jy/beam ! Combined flux in latest and established models = 1.36483 Jy ! Performing phase self-cal ! Adding 74 model components to the UV plane model. ! The established model now contains 715 components and 1.36483 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119682Jy sigma=0.736768 ! Fit after self-cal, rms=0.119676Jy sigma=0.736611 ! Inverting map ! Added new window around map position (30.8, 90.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00149407 Jy ! Component: 100 - total flux cleaned = 0.00246467 Jy ! Component: 150 - total flux cleaned = 0.00328174 Jy ! Component: 200 - total flux cleaned = 0.00400675 Jy ! Component: 250 - total flux cleaned = 0.00462535 Jy ! Component: 300 - total flux cleaned = 0.00510771 Jy ! Component: 350 - total flux cleaned = 0.00567221 Jy ! Component: 400 - total flux cleaned = 0.00611503 Jy ! Total flux subtracted in 400 components = 0.00611503 Jy ! Clean residual min=-0.002436 max=0.002487 Jy/beam ! Clean residual mean=0.000027 rms=0.000732 Jy/beam ! Combined flux in latest and established models = 1.37094 Jy ! Performing phase self-cal ! Adding 84 model components to the UV plane model. ! The established model now contains 788 components and 1.37094 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119576Jy sigma=0.735084 ! Fit after self-cal, rms=0.119575Jy sigma=0.735009 ! 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=4.681 mas, bmaj=11.1 mas, bpa=-4.684 degrees ! Estimated noise=1.38719 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.681 x 11.1 at -4.684 degrees (North through East) ! Clean map min=-0.0040249 max=0.93494 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.119575Jy sigma=0.735009 ! Fit after self-cal, rms=0.134469Jy sigma=0.615776 ! Inverting map and beam ! Estimated beam: bmin=4.691 mas, bmaj=11.11 mas, bpa=-4.66 degrees ! Estimated noise=1.38646 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00195879 Jy ! Component: 100 - total flux cleaned = 0.00232791 Jy ! Component: 150 - total flux cleaned = 0.00258027 Jy ! Component: 200 - total flux cleaned = 0.00305472 Jy ! Component: 250 - total flux cleaned = 0.0034299 Jy ! Component: 300 - total flux cleaned = 0.00371721 Jy ! Component: 350 - total flux cleaned = 0.00406453 Jy ! Component: 400 - total flux cleaned = 0.00420078 Jy ! Total flux subtracted in 400 components = 0.00420078 Jy ! Clean residual min=-0.002774 max=0.002939 Jy/beam ! Clean residual mean=0.000024 rms=0.000641 Jy/beam ! Combined flux in latest and established models = 1.37514 Jy ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 841 components and 1.37514 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133949Jy sigma=0.614292 ! Fit after self-cal, rms=0.133962Jy sigma=0.614225 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00117481 Jy ! Component: 100 - total flux cleaned = 0.00158029 Jy ! Component: 150 - total flux cleaned = 0.00177802 Jy ! Component: 200 - total flux cleaned = 0.00197004 Jy ! Component: 250 - total flux cleaned = 0.00215786 Jy ! Component: 300 - total flux cleaned = 0.00215885 Jy ! Component: 350 - total flux cleaned = 0.00234069 Jy ! Component: 400 - total flux cleaned = 0.00251942 Jy ! Total flux subtracted in 400 components = 0.00251942 Jy ! Clean residual min=-0.002728 max=0.002743 Jy/beam ! Clean residual mean=0.000021 rms=0.000613 Jy/beam ! Combined flux in latest and established models = 1.37766 Jy ! Performing phase self-cal ! Adding 91 model components to the UV plane model. ! The established model now contains 884 components and 1.37766 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133610Jy sigma=0.613367 ! Fit after self-cal, rms=0.133621Jy sigma=0.613345 ! 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.133621Jy sigma=0.613345 ! Fit after self-cal, rms=0.133220Jy sigma=0.613203 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=3.542 mas, bmaj=8.833 mas, bpa=-0.3114 degrees ! Estimated noise=2.47989 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.671768 Jy ! Component: 100 - total flux cleaned = 0.883411 Jy ! Component: 150 - total flux cleaned = 0.989208 Jy ! Component: 200 - total flux cleaned = 1.04444 Jy ! Component: 250 - total flux cleaned = 1.08214 Jy ! Component: 300 - total flux cleaned = 1.11204 Jy ! Component: 350 - total flux cleaned = 1.13655 Jy ! Component: 400 - total flux cleaned = 1.15724 Jy ! Total flux subtracted in 400 components = 1.15724 Jy ! Clean residual min=-0.009480 max=0.014364 Jy/beam ! Clean residual mean=-0.000077 rms=0.002935 Jy/beam ! Combined flux in latest and established models = 1.15724 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 81 model components to the UV plane model. ! The established model now contains 81 components and 1.15724 Jy ! Inverting map and beam ! Estimated beam: bmin=4.687 mas, bmaj=11.11 mas, bpa=-4.678 degrees ! Estimated noise=1.38621 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.0366659 Jy ! Component: 100 - total flux cleaned = 0.063726 Jy ! Component: 150 - total flux cleaned = 0.0850431 Jy ! Component: 200 - total flux cleaned = 0.102332 Jy ! Component: 250 - total flux cleaned = 0.116613 Jy ! Component: 300 - total flux cleaned = 0.128523 Jy ! Component: 350 - total flux cleaned = 0.138619 Jy ! Component: 400 - total flux cleaned = 0.147306 Jy ! Total flux subtracted in 400 components = 0.147306 Jy ! Clean residual min=-0.005444 max=0.006516 Jy/beam ! Clean residual mean=0.000012 rms=0.001967 Jy/beam ! Combined flux in latest and established models = 1.30455 Jy ! Adding 165 model components to the UV plane model. ! The established model now contains 230 components and 1.30454 Jy ! Inverting map ! Added new window around map position (42.4, 4.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00776404 Jy ! Component: 100 - total flux cleaned = 0.0139143 Jy ! Component: 150 - total flux cleaned = 0.0193735 Jy ! Component: 200 - total flux cleaned = 0.0236433 Jy ! Component: 250 - total flux cleaned = 0.0272984 Jy ! Component: 300 - total flux cleaned = 0.0306237 Jy ! Component: 350 - total flux cleaned = 0.0337073 Jy ! Component: 400 - total flux cleaned = 0.0364243 Jy ! Total flux subtracted in 400 components = 0.0364243 Jy ! Clean residual min=-0.003689 max=0.003226 Jy/beam ! Clean residual mean=0.000015 rms=0.001048 Jy/beam ! Combined flux in latest and established models = 1.34097 Jy ! Adding 112 model components to the UV plane model. ! The established model now contains 317 components and 1.34097 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.697 mas, bmaj=15.59 mas, bpa=-6.579 degrees ! Estimated noise=1.40646 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00332387 Jy ! Component: 100 - total flux cleaned = 0.0059009 Jy ! Component: 150 - total flux cleaned = 0.00770215 Jy ! Component: 200 - total flux cleaned = 0.00934058 Jy ! Component: 250 - total flux cleaned = 0.0111478 Jy ! Component: 300 - total flux cleaned = 0.0127454 Jy ! Component: 350 - total flux cleaned = 0.0140765 Jy ! Component: 400 - total flux cleaned = 0.0153297 Jy ! Total flux subtracted in 400 components = 0.0153297 Jy ! Clean residual min=-0.002637 max=0.002333 Jy/beam ! Clean residual mean=0.000018 rms=0.000769 Jy/beam ! Combined flux in latest and established models = 1.3563 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.000818291 Jy ! Component: 100 - total flux cleaned = 0.00152796 Jy ! Component: 150 - total flux cleaned = 0.00227164 Jy ! Component: 200 - total flux cleaned = 0.00292295 Jy ! Component: 250 - total flux cleaned = 0.00342427 Jy ! Component: 300 - total flux cleaned = 0.00378962 Jy ! Component: 350 - total flux cleaned = 0.00420189 Jy ! Component: 400 - total flux cleaned = 0.0045444 Jy ! Total flux subtracted in 400 components = 0.0045444 Jy ! Clean residual min=-0.002448 max=0.002127 Jy/beam ! Clean residual mean=0.000018 rms=0.000661 Jy/beam ! Combined flux in latest and established models = 1.36084 Jy ! Adding 204 model components to the UV plane model. ! The established model now contains 515 components and 1.36084 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.687 mas, bmaj=11.11 mas, bpa=-4.678 degrees ! Estimated noise=1.38621 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00182509 Jy ! Component: 100 - total flux cleaned = -0.00194649 Jy ! Component: 150 - total flux cleaned = -0.00217052 Jy ! Component: 200 - total flux cleaned = -0.00173959 Jy ! Component: 250 - total flux cleaned = -0.00163758 Jy ! Component: 300 - total flux cleaned = -0.00143741 Jy ! Component: 350 - total flux cleaned = -0.0012448 Jy ! Component: 400 - total flux cleaned = -0.000963894 Jy ! Total flux subtracted in 400 components = -0.000963894 Jy ! Clean residual min=-0.002850 max=0.002730 Jy/beam ! Clean residual mean=0.000015 rms=0.000709 Jy/beam ! Combined flux in latest and established models = 1.35988 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 543 components and 1.35988 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136853Jy sigma=0.616803 ! Fit after self-cal, rms=0.136758Jy sigma=0.616258 wmodel /scr/pima/rdv62_uvs/J2331-1556_S_map.mod ! Writing 543 model components to file: /scr/pima/rdv62_uvs/J2331-1556_S_map.mod wobs /scr/pima/rdv62_uvs/J2331-1556_S_uvs.fits ! Applying 2053 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J2331-1556_S_uvs.fits wwins /scr/pima/rdv62_uvs/J2331-1556_S_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv62_uvs/J2331-1556_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.00116205 Jy ! Component: 100 - total flux cleaned = -0.00177788 Jy ! Component: 150 - total flux cleaned = -0.00205071 Jy ! Component: 200 - total flux cleaned = -0.0023139 Jy ! Component: 250 - total flux cleaned = -0.00256866 Jy ! Component: 300 - total flux cleaned = -0.0029003 Jy ! Component: 350 - total flux cleaned = -0.00318306 Jy ! Component: 400 - total flux cleaned = -0.00330226 Jy ! Component: 450 - total flux cleaned = -0.00345895 Jy ! Component: 500 - total flux cleaned = -0.0036904 Jy ! Component: 550 - total flux cleaned = -0.00388098 Jy ! Component: 600 - total flux cleaned = -0.00410719 Jy ! Component: 650 - total flux cleaned = -0.0043305 Jy ! Component: 700 - total flux cleaned = -0.00455146 Jy ! Component: 750 - total flux cleaned = -0.00477007 Jy ! Component: 800 - total flux cleaned = -0.00491454 Jy ! Component: 850 - total flux cleaned = -0.0050936 Jy ! Component: 900 - total flux cleaned = -0.00520036 Jy ! Component: 950 - total flux cleaned = -0.00534182 Jy ! Component: 1000 - total flux cleaned = -0.0054825 Jy ! Component: 1050 - total flux cleaned = -0.00572613 Jy ! Component: 1100 - total flux cleaned = -0.00579515 Jy ! Component: 1150 - total flux cleaned = -0.00600067 Jy ! Component: 1200 - total flux cleaned = -0.00610276 Jy ! Component: 1250 - total flux cleaned = -0.00623789 Jy ! Component: 1300 - total flux cleaned = -0.00630515 Jy ! Component: 1350 - total flux cleaned = -0.00653904 Jy ! Component: 1400 - total flux cleaned = -0.00663826 Jy ! Component: 1450 - total flux cleaned = -0.00673702 Jy ! Component: 1500 - total flux cleaned = -0.00680268 Jy ! Component: 1550 - total flux cleaned = -0.0069977 Jy ! Component: 1600 - total flux cleaned = -0.00709444 Jy ! Component: 1650 - total flux cleaned = -0.0070947 Jy ! Component: 1700 - total flux cleaned = -0.00725438 Jy ! Component: 1750 - total flux cleaned = -0.00731801 Jy ! Component: 1800 - total flux cleaned = -0.00744427 Jy ! Component: 1850 - total flux cleaned = -0.00756976 Jy ! Component: 1900 - total flux cleaned = -0.0076009 Jy ! Component: 1950 - total flux cleaned = -0.00775615 Jy ! Component: 2000 - total flux cleaned = -0.0078487 Jy ! Component: 2050 - total flux cleaned = -0.00794055 Jy ! Component: 2100 - total flux cleaned = -0.00809304 Jy ! Component: 2150 - total flux cleaned = -0.00806273 Jy ! Component: 2200 - total flux cleaned = -0.00821366 Jy ! Component: 2250 - total flux cleaned = -0.00818366 Jy ! Component: 2300 - total flux cleaned = -0.00827344 Jy ! Component: 2350 - total flux cleaned = -0.0084217 Jy ! Component: 2400 - total flux cleaned = -0.00839216 Jy ! Component: 2450 - total flux cleaned = -0.00845114 Jy ! Component: 2500 - total flux cleaned = -0.00848052 Jy ! Component: 2550 - total flux cleaned = -0.00853896 Jy ! Component: 2600 - total flux cleaned = -0.00859703 Jy ! Component: 2650 - total flux cleaned = -0.00862615 Jy ! Component: 2700 - total flux cleaned = -0.00868336 Jy ! Component: 2750 - total flux cleaned = -0.00874069 Jy ! Component: 2800 - total flux cleaned = -0.00879756 Jy ! Component: 2850 - total flux cleaned = -0.00876914 Jy ! Component: 2900 - total flux cleaned = -0.00882546 Jy ! Component: 2950 - total flux cleaned = -0.00890974 Jy ! Component: 3000 - total flux cleaned = -0.00882596 Jy ! Component: 3050 - total flux cleaned = -0.00890942 Jy ! Component: 3100 - total flux cleaned = -0.0088817 Jy ! Component: 3150 - total flux cleaned = -0.00896441 Jy ! Component: 3200 - total flux cleaned = -0.00890924 Jy ! Component: 3250 - total flux cleaned = -0.00896417 Jy ! Component: 3300 - total flux cleaned = -0.00899118 Jy ! Component: 3350 - total flux cleaned = -0.00896395 Jy ! Component: 3400 - total flux cleaned = -0.00907217 Jy ! Component: 3450 - total flux cleaned = -0.00899126 Jy ! Component: 3500 - total flux cleaned = -0.00907185 Jy ! Component: 3550 - total flux cleaned = -0.00909859 Jy ! Component: 3600 - total flux cleaned = -0.00912517 Jy ! Component: 3650 - total flux cleaned = -0.00904537 Jy ! Component: 3700 - total flux cleaned = -0.00912471 Jy ! Component: 3750 - total flux cleaned = -0.00920391 Jy ! Component: 3800 - total flux cleaned = -0.0092564 Jy ! Component: 3850 - total flux cleaned = -0.00930877 Jy ! Component: 3900 - total flux cleaned = -0.00936092 Jy ! Component: 3950 - total flux cleaned = -0.00938684 Jy ! Component: 4000 - total flux cleaned = -0.00941283 Jy ! Component: 4050 - total flux cleaned = -0.00930972 Jy ! Total flux subtracted in 4096 components = -0.00943857 Jy ! Clean residual min=-0.001317 max=0.001308 Jy/beam ! Clean residual mean=0.000009 rms=0.000495 Jy/beam ! Combined flux in latest and established models = 1.35044 Jy keep ! Adding 231 model components to the UV plane model. ! The established model now contains 768 components and 1.35044 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J2331-1556_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.687 x 11.11 at -4.678 degrees (North through East) ! Clean map min=-0.003469 max=0.93803 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J2331-1556_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.938 rms= 0.00049 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J2331-1556_S_uvs.log closed on Fri Nov 22 22:43:00 2013