! Started logfile: /scr/pima/rdv94_uvs/J1230+1223_S_uvs.log on Tue Oct 8 12:18:43 2013 obs /scr/pima/rdv94_uvs/J1230+1223_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J1230+1223_S_uva.fits ! AN table 1: 4 integrations on 66 of 66 possible baselines. ! AN table 2: 2 integrations on 3 of 3 possible baselines. ! AN table 3: 3 integrations on 45 of 45 possible baselines. ! AN table 4: 2 integrations on 3 of 3 possible baselines. ! AN table 5: 3 integrations on 15 of 15 possible baselines. ! AN table 6: 3 integrations on 1 of 1 possible baselines. ! AN table 7: 3 integrations on 3 of 3 possible baselines. ! AN table 8: 3 integrations on 1 of 1 possible baselines. ! AN table 9: 63 integrations on 6 of 6 possible baselines. ! AN table 10: 24 integrations on 10 of 10 possible baselines. ! AN table 11: 4 integrations on 3 of 3 possible baselines. ! AN table 12: 2 integrations on 3 of 3 possible baselines. ! AN table 13: 12 integrations on 21 of 21 possible baselines. ! AN table 14: 12 integrations on 15 of 15 possible baselines. ! AN table 15: 12 integrations on 78 of 78 possible baselines. ! Apparent sampling: 0.582626 visibilities/baseline/integration-bin. ! Found source: J1230+1223 ! ! 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 5768 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv94_uvs/J1230+1223_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:ON self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:WF self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:WF self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:WF self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:ON self-cal status=correctable weight=1 ! 9:HH self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:ON self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 10:HH self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:ON self-cal status=correctable weight=1 ! 10:WF self-cal status=correctable weight=1 ! 10:WZ self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:NY self-cal status=correctable weight=1 ! 11:WF self-cal status=correctable weight=1 ! 12:HN self-cal status=correctable weight=1 ! 12:NL self-cal status=correctable weight=1 ! 12:WF self-cal status=correctable weight=1 ! 13:HH self-cal status=correctable weight=1 ! 13:HN self-cal status=correctable weight=1 ! 13:NL self-cal status=correctable weight=1 ! 13:NY self-cal status=correctable weight=1 ! 13:ON self-cal status=correctable weight=1 ! 13:WF self-cal status=correctable weight=1 ! 13:WZ self-cal status=correctable weight=1 ! 14:HH self-cal status=correctable weight=1 ! 14:HN self-cal status=correctable weight=1 ! 14:NY self-cal status=correctable weight=1 ! 14:ON self-cal status=correctable weight=1 ! 14:WF self-cal status=correctable weight=1 ! 14:WZ self-cal status=correctable weight=1 ! 15:BR self-cal status=correctable weight=1 ! 15:FD self-cal status=correctable weight=1 ! 15:HH self-cal status=correctable weight=1 ! 15:HN self-cal status=correctable weight=1 ! 15:KP self-cal status=correctable weight=1 ! 15:LA self-cal status=correctable weight=1 ! 15:NL self-cal status=correctable weight=1 ! 15:NY self-cal status=correctable weight=1 ! 15:ON self-cal status=correctable weight=1 ! 15:OV self-cal status=correctable weight=1 ! 15:PT self-cal status=correctable weight=1 ! 15:WF self-cal status=correctable weight=1 ! 15:WZ self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 119 telescope corrections were flagged in sub-array 9. ! A total of 48 telescope corrections were flagged in sub-array 10. ! A total of 3 telescope corrections were flagged in sub-array 11. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 73 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 119 telescope corrections were flagged in sub-array 9. ! A total of 26 telescope corrections were flagged in sub-array 10. ! A total of 3 telescope corrections were flagged in sub-array 11. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 73 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 119 telescope corrections were flagged in sub-array 9. ! A total of 26 telescope corrections were flagged in sub-array 10. ! A total of 3 telescope corrections were flagged in sub-array 11. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 73 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 119 telescope corrections were flagged in sub-array 9. ! A total of 26 telescope corrections were flagged in sub-array 10. ! A total of 3 telescope corrections were flagged in sub-array 11. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 73 telescope corrections were flagged in sub-array 15. ! ! Fit before self-cal, rms=1.792255Jy sigma=28.128907 ! Fit after self-cal, rms=1.027747Jy sigma=12.852660 ! 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:ON ! Flagging baselines of antenna 1:ON in the currently selected channels. flag 1:WF ! Flagging baselines of antenna 1:WF in the currently selected channels. flag 2:WF ! Flagging baselines of antenna 2:WF in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 3:WF ! Flagging baselines of antenna 3:WF in the currently selected channels. flag 6:NY ! Flagging baselines of antenna 6:NY in the currently selected channels. flag 7:NY ! Flagging baselines of antenna 7:NY in the currently selected channels. flag 7:ON ! Flagging baselines of antenna 7:ON in the currently selected channels. flag 7:WZ ! Flagging baselines of antenna 7:WZ in the currently selected channels. flag 8:NY ! Flagging baselines of antenna 8:NY in the currently selected channels. flag 8:ON ! Flagging baselines of antenna 8:ON in the currently selected channels. flag 9:HH ! Flagging baselines of antenna 9:HH in the currently selected channels. flag 9:NY ! Flagging baselines of antenna 9:NY in the currently selected channels. flag 9:ON ! Flagging baselines of antenna 9:ON in the currently selected channels. flag 9:WZ ! Flagging baselines of antenna 9:WZ in the currently selected channels. flag 10:HH ! Flagging baselines of antenna 10:HH in the currently selected channels. flag 10:NY ! Flagging baselines of antenna 10:NY in the currently selected channels. flag 10:ON ! Flagging baselines of antenna 10:ON in the currently selected channels. flag 10:WF ! Flagging baselines of antenna 10:WF in the currently selected channels. flag 10:WZ ! Flagging baselines of antenna 10:WZ in the currently selected channels. flag 11:NY ! Flagging baselines of antenna 11:NY in the currently selected channels. flag 11:WF ! Flagging baselines of antenna 11:WF in the currently selected channels. flag 12:WF ! Flagging baselines of antenna 12:WF in the currently selected channels. flag 13:HH ! Flagging baselines of antenna 13:HH in the currently selected channels. flag 13:NY ! Flagging baselines of antenna 13:NY in the currently selected channels. flag 13:ON ! Flagging baselines of antenna 13:ON in the currently selected channels. flag 13:WF ! Flagging baselines of antenna 13:WF in the currently selected channels. flag 13:WZ ! Flagging baselines of antenna 13:WZ in the currently selected channels. flag 14:HH ! Flagging baselines of antenna 14:HH in the currently selected channels. flag 14:NY ! Flagging baselines of antenna 14:NY in the currently selected channels. flag 14:ON ! Flagging baselines of antenna 14:ON in the currently selected channels. flag 14:WF ! Flagging baselines of antenna 14:WF in the currently selected channels. flag 14:WZ ! Flagging baselines of antenna 14:WZ in the currently selected channels. flag 15:HH ! Flagging baselines of antenna 15:HH in the currently selected channels. flag 15:NY ! Flagging baselines of antenna 15:NY in the currently selected channels. flag 15:ON ! Flagging baselines of antenna 15:ON in the currently selected channels. flag 15:WF ! Flagging baselines of antenna 15:WF in the currently selected channels. flag 15:WZ ! Flagging baselines of antenna 15:WZ in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.42 mas, bmaj=8.958 mas, bpa=-19.52 degrees ! Estimated noise=1.07423 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.890255 Jy ! Component: 100 - total flux cleaned = 1.18082 Jy ! Total flux subtracted in 100 components = 1.18082 Jy ! Clean residual min=-0.071421 max=0.132883 Jy/beam ! Clean residual mean=0.000138 rms=0.023896 Jy/beam ! Combined flux in latest and established models = 1.18082 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 12 components and 1.18082 Jy ! ! Correcting IF 1. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 133 telescope corrections were flagged in sub-array 9. ! A total of 72 telescope corrections were flagged in sub-array 10. ! A total of 9 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 41 telescope corrections were flagged in sub-array 13. ! A total of 37 telescope corrections were flagged in sub-array 14. ! A total of 42 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 2. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 133 telescope corrections were flagged in sub-array 9. ! A total of 94 telescope corrections were flagged in sub-array 10. ! A total of 9 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 41 telescope corrections were flagged in sub-array 13. ! A total of 37 telescope corrections were flagged in sub-array 14. ! A total of 42 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 3. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 133 telescope corrections were flagged in sub-array 9. ! A total of 94 telescope corrections were flagged in sub-array 10. ! A total of 9 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 41 telescope corrections were flagged in sub-array 13. ! A total of 37 telescope corrections were flagged in sub-array 14. ! A total of 42 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 133 telescope corrections were flagged in sub-array 9. ! A total of 94 telescope corrections were flagged in sub-array 10. ! A total of 9 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 41 telescope corrections were flagged in sub-array 13. ! A total of 37 telescope corrections were flagged in sub-array 14. ! A total of 42 telescope corrections were flagged in sub-array 15. ! ! Fit before self-cal, rms=0.407045Jy sigma=13.261963 ! Fit after self-cal, rms=0.408973Jy sigma=13.374019 ! Inverting map and beam ! Estimated beam: bmin=4.436 mas, bmaj=9.009 mas, bpa=-20.35 degrees ! Estimated noise=1.06817 mJy/beam. ! Added new window around map position (-6.5, 1.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.15144 Jy ! Component: 100 - total flux cleaned = 0.239365 Jy ! Total flux subtracted in 100 components = 0.239365 Jy ! Clean residual min=-0.058582 max=0.129265 Jy/beam ! Clean residual mean=0.000277 rms=0.018333 Jy/beam ! Combined flux in latest and established models = 1.42019 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 34 components and 1.42019 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.345381Jy sigma=10.789225 ! Fit after self-cal, rms=0.342377Jy sigma=10.759915 ! Inverting map ! Added new window around map position (-19, 6.5). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.11458 Jy ! Component: 100 - total flux cleaned = 0.177739 Jy ! Total flux subtracted in 100 components = 0.177739 Jy ! Clean residual min=-0.044249 max=0.062592 Jy/beam ! Clean residual mean=0.000344 rms=0.014173 Jy/beam ! Combined flux in latest and established models = 1.59793 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 57 components and 1.59793 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.298515Jy sigma=9.118164 ! Fit after self-cal, rms=0.291141Jy sigma=9.008536 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.651 mas, bmaj=11.18 mas, bpa=-18.78 degrees ! Estimated noise=0.755852 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.0935931 Jy ! Component: 100 - total flux cleaned = 0.148184 Jy ! Component: 150 - total flux cleaned = 0.185752 Jy ! Component: 200 - total flux cleaned = 0.212481 Jy ! Total flux subtracted in 200 components = 0.212481 Jy ! Clean residual min=-0.046567 max=0.054124 Jy/beam ! Clean residual mean=0.000980 rms=0.018140 Jy/beam ! Combined flux in latest and established models = 1.81041 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 112 components and 1.81041 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272415Jy sigma=7.955519 ! Fit after self-cal, rms=0.267654Jy sigma=7.875139 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.09 HN 1.07 KP 1.19 ! LA 1.21 MK 0.84 NL 0.92 NY 1.00* ! ON 1.00* OV 0.95 PT 1.02 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 1.10 HN 1.10 KP 1.11 ! LA 1.27 MK 0.87 NL 0.96 NY 1.00* ! OV 0.95 WF 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.18 KP 1.30 LA 1.36 ! MK 0.87 OV 0.97 ! ! Telescope amplitude corrections in sub-array 6: ! MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* ON 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NY 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 12: ! HN 0.00* NL 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 14: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 15: ! BR 0.95 FD 0.99 HH 1.00* HN 1.13 ! KP 1.08 LA 1.10 NL 0.98 NY 1.00* ! ON 1.00* OV 0.95 PT 1.00 WF 1.00* ! WZ 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.94 FD 1.16 HN 1.13 KP 1.20 ! LA 1.10 MK 0.92 NL 0.98 NY 1.00* ! ON 1.00* OV 0.99 PT 1.01 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.95 FD 1.17 HN 1.14 KP 1.13 ! LA 1.18 MK 0.92 NL 1.02 NY 1.00* ! OV 0.97 WF 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.23 KP 1.30 LA 1.22 ! MK 0.95 OV 1.01 ! ! Telescope amplitude corrections in sub-array 6: ! MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* ON 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NY 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 12: ! HN 0.00* NL 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 14: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 15: ! BR 1.00 FD 1.05 HH 1.00* HN 1.18 ! KP 1.08 LA 1.02 NL 1.04 NY 1.00* ! ON 1.00* OV 0.95 PT 1.00 WF 1.00* ! WZ 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.18 HN 1.06 KP 1.03 ! LA 0.89 MK 0.93 NL 1.02 NY 1.00* ! ON 1.00* OV 0.98 PT 1.02 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.99 FD 1.16 HN 1.10 KP 0.97 ! LA 0.95 MK 0.94 NL 1.05 NY 1.00* ! OV 1.00 WF 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.01 FD 1.31 KP 1.09 LA 1.04 ! MK 0.96 OV 1.04 ! ! Telescope amplitude corrections in sub-array 6: ! MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* ON 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NY 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 12: ! HN 0.00* NL 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 14: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 15: ! BR 1.27 FD 1.04 HH 1.00* HN 1.14 ! KP 0.96 LA 0.85 NL 1.07 NY 1.00* ! ON 1.00* OV 0.95 PT 1.01 WF 1.00* ! WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.87 FD 1.12 HN 0.92 KP 0.99 ! LA 0.73 MK 1.00 NL 0.96 NY 1.00* ! ON 1.00* OV 0.85 PT 0.91 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* HN 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.87 FD 1.14 HN 0.98 KP 0.91 ! LA 0.78 MK 1.03 NL 0.96 NY 1.00* ! OV 0.83 WF 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* FD 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 1.23 KP 1.05 LA 0.80 ! MK 1.04 OV 0.86 ! ! Telescope amplitude corrections in sub-array 6: ! MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* ON 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 11: ! HN 0.00* NY 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 12: ! HN 0.00* NL 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 13: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 14: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 15: ! BR 1.47 FD 0.99 HH 1.00* HN 1.03 ! KP 0.89 LA 0.70 NL 0.98 NY 1.00* ! ON 1.00* OV 0.81 PT 0.92 WF 1.00* ! WZ 1.00* ! ! ! Fit before self-cal, rms=0.267654Jy sigma=7.875139 ! Fit after self-cal, rms=0.127009Jy sigma=3.588654 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=5.665 mas, bmaj=11.66 mas, bpa=-17.86 degrees ! Estimated noise=0.754452 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.00623344 Jy ! Component: 100 - total flux cleaned = -0.0086677 Jy ! Component: 150 - total flux cleaned = -0.010007 Jy ! Component: 200 - total flux cleaned = -0.0116329 Jy ! Total flux subtracted in 200 components = -0.0116329 Jy ! Clean residual min=-0.024988 max=0.044460 Jy/beam ! Clean residual mean=0.000048 rms=0.011644 Jy/beam ! Combined flux in latest and established models = 1.79877 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 131 components and 1.79877 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124848Jy sigma=3.529221 ! Fit after self-cal, rms=0.124767Jy sigma=3.523605 ! 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.124767Jy sigma=3.523605 ! Fit after self-cal, rms=0.124514Jy sigma=3.512408 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.673 mas, bmaj=11.67 mas, bpa=-17.85 degrees ! Estimated noise=0.751406 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0082466 Jy ! Component: 100 - total flux cleaned = -0.0101748 Jy ! Component: 150 - total flux cleaned = -0.0117533 Jy ! Component: 200 - total flux cleaned = -0.0126263 Jy ! Total flux subtracted in 200 components = -0.0126263 Jy ! Clean residual min=-0.024512 max=0.043857 Jy/beam ! Clean residual mean=0.000048 rms=0.011425 Jy/beam ! Combined flux in latest and established models = 1.78615 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 147 components and 1.78615 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123123Jy sigma=3.482881 ! Fit after self-cal, rms=0.122940Jy sigma=3.481641 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.000358031 Jy ! Component: 100 - total flux cleaned = -0.00070818 Jy ! Component: 150 - total flux cleaned = -0.00103043 Jy ! Component: 200 - total flux cleaned = -0.00118057 Jy ! Total flux subtracted in 200 components = -0.00118057 Jy ! Clean residual min=-0.024423 max=0.043890 Jy/beam ! Clean residual mean=0.000049 rms=0.011374 Jy/beam ! Combined flux in latest and established models = 1.78497 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 152 components and 1.78497 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122478Jy sigma=3.470369 ! Fit after self-cal, rms=0.122383Jy sigma=3.469965 ! 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.122383Jy sigma=3.469965 ! Fit after self-cal, rms=0.122109Jy sigma=3.457496 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=5.683 mas, bmaj=11.69 mas, bpa=-17.85 degrees ! Estimated noise=0.747787 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00890239 Jy ! Component: 100 - total flux cleaned = -0.0124616 Jy ! Component: 150 - total flux cleaned = -0.0137559 Jy ! Component: 200 - total flux cleaned = -0.0146734 Jy ! Total flux subtracted in 200 components = -0.0146734 Jy ! Clean residual min=-0.024697 max=0.043531 Jy/beam ! Clean residual mean=0.000047 rms=0.011229 Jy/beam ! Combined flux in latest and established models = 1.77029 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 162 components and 1.77029 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121071Jy sigma=3.447370 ! Fit after self-cal, rms=0.121018Jy sigma=3.446867 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.00057482 Jy ! Component: 100 - total flux cleaned = -0.000703128 Jy ! Component: 150 - total flux cleaned = -0.000831918 Jy ! Component: 200 - total flux cleaned = -0.000960529 Jy ! Total flux subtracted in 200 components = -0.000960529 Jy ! Clean residual min=-0.024703 max=0.043655 Jy/beam ! Clean residual mean=0.000047 rms=0.011214 Jy/beam ! Combined flux in latest and established models = 1.76933 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 166 components and 1.76933 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120680Jy sigma=3.441970 ! Fit after self-cal, rms=0.120688Jy sigma=3.441792 ! 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.120688Jy sigma=3.441792 ! Fit after self-cal, rms=0.120567Jy sigma=3.432179 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.65 mas, bmaj=13.28 mas, bpa=-23.47 degrees ! Estimated noise=0.737786 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0044278 Jy ! Component: 100 - total flux cleaned = -0.00746765 Jy ! Component: 150 - total flux cleaned = -0.00911508 Jy ! Component: 200 - total flux cleaned = -0.00992987 Jy ! Total flux subtracted in 200 components = -0.00992987 Jy ! Clean residual min=-0.030521 max=0.056378 Jy/beam ! Clean residual mean=0.000074 rms=0.014787 Jy/beam ! Combined flux in latest and established models = 1.7594 Jy selfcal ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 188 components and 1.7594 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120502Jy sigma=3.413736 ! Fit after self-cal, rms=0.120310Jy sigma=3.412064 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.000750025 Jy ! Component: 100 - total flux cleaned = 0.000220659 Jy ! Component: 150 - total flux cleaned = -0.000273408 Jy ! Component: 200 - total flux cleaned = -0.000431883 Jy ! Total flux subtracted in 200 components = -0.000431883 Jy ! Clean residual min=-0.030842 max=0.055756 Jy/beam ! Clean residual mean=0.000074 rms=0.014713 Jy/beam ! Combined flux in latest and established models = 1.75897 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 194 components and 1.75897 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120702Jy sigma=3.405435 ! Fit after self-cal, rms=0.120565Jy sigma=3.404868 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:MK,true ! 1:MK self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT 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 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 11:HN,true ! 11:HN self-cal status=fixed weight=1 selfant 12:HN,true ! 12:HN self-cal status=fixed weight=1 selfant 12:NL,true ! 12:NL self-cal status=fixed weight=1 selfant 13:HN,true ! 13:HN self-cal status=fixed weight=1 selfant 13:NL,true ! 13:NL self-cal status=fixed weight=1 selfant 14:HN,true ! 14:HN self-cal status=fixed weight=1 selfant 15:BR,true ! 15:BR self-cal status=fixed weight=1 selfant 15:FD,true ! 15:FD self-cal status=fixed weight=1 selfant 15:HN,true ! 15:HN self-cal status=fixed weight=1 selfant 15:KP,true ! 15:KP self-cal status=fixed weight=1 selfant 15:LA,true ! 15:LA self-cal status=fixed weight=1 selfant 15:NL,true ! 15:NL self-cal status=fixed weight=1 selfant 15:OV,true ! 15:OV self-cal status=fixed weight=1 selfant 15:PT,true ! 15: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 2 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 63 telescope corrections were flagged in sub-array 10. ! A total of 8 telescope corrections were flagged in sub-array 11. ! A total of 2 telescope corrections were flagged in sub-array 12. ! A total of 27 telescope corrections were flagged in sub-array 13. ! A total of 27 telescope corrections were flagged in sub-array 14. ! A total of 27 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 8 telescope corrections were flagged in sub-array 11. ! A total of 2 telescope corrections were flagged in sub-array 12. ! A total of 27 telescope corrections were flagged in sub-array 13. ! A total of 27 telescope corrections were flagged in sub-array 14. ! A total of 27 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 8 telescope corrections were flagged in sub-array 11. ! A total of 2 telescope corrections were flagged in sub-array 12. ! A total of 27 telescope corrections were flagged in sub-array 13. ! A total of 27 telescope corrections were flagged in sub-array 14. ! A total of 27 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 8 telescope corrections were flagged in sub-array 11. ! A total of 2 telescope corrections were flagged in sub-array 12. ! A total of 27 telescope corrections were flagged in sub-array 13. ! A total of 27 telescope corrections were flagged in sub-array 14. ! A total of 27 telescope corrections were flagged in sub-array 15. ! ! Fit before self-cal, rms=0.884963Jy sigma=8.074081 ! Fit after self-cal, rms=0.273090Jy sigma=2.937471 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:MK self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:ON self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:WF self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:WF self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:WF self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:ON self-cal status=correctable weight=1 ! 9:HH self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:ON self-cal status=correctable weight=1 ! 9:WZ self-cal status=correctable weight=1 ! 10:HH self-cal status=correctable weight=1 ! 10:NY self-cal status=correctable weight=1 ! 10:ON self-cal status=correctable weight=1 ! 10:WF self-cal status=correctable weight=1 ! 10:WZ self-cal status=correctable weight=1 ! 11:HN self-cal status=correctable weight=1 ! 11:NY self-cal status=correctable weight=1 ! 11:WF self-cal status=correctable weight=1 ! 12:HN self-cal status=correctable weight=1 ! 12:NL self-cal status=correctable weight=1 ! 12:WF self-cal status=correctable weight=1 ! 13:HH self-cal status=correctable weight=1 ! 13:HN self-cal status=correctable weight=1 ! 13:NL self-cal status=correctable weight=1 ! 13:NY self-cal status=correctable weight=1 ! 13:ON self-cal status=correctable weight=1 ! 13:WF self-cal status=correctable weight=1 ! 13:WZ self-cal status=correctable weight=1 ! 14:HH self-cal status=correctable weight=1 ! 14:HN self-cal status=correctable weight=1 ! 14:NY self-cal status=correctable weight=1 ! 14:ON self-cal status=correctable weight=1 ! 14:WF self-cal status=correctable weight=1 ! 14:WZ self-cal status=correctable weight=1 ! 15:BR self-cal status=correctable weight=1 ! 15:FD self-cal status=correctable weight=1 ! 15:HH self-cal status=correctable weight=1 ! 15:HN self-cal status=correctable weight=1 ! 15:KP self-cal status=correctable weight=1 ! 15:LA self-cal status=correctable weight=1 ! 15:NL self-cal status=correctable weight=1 ! 15:NY self-cal status=correctable weight=1 ! 15:ON self-cal status=correctable weight=1 ! 15:OV self-cal status=correctable weight=1 ! 15:PT self-cal status=correctable weight=1 ! 15:WF self-cal status=correctable weight=1 ! 15:WZ self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 4 telescope corrections were flagged in sub-array 11. ! A total of 4 telescope corrections were flagged in sub-array 12. ! A total of 16 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! A total of 55 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 2. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 4 telescope corrections were flagged in sub-array 11. ! A total of 4 telescope corrections were flagged in sub-array 12. ! A total of 16 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! A total of 55 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 3. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 4 telescope corrections were flagged in sub-array 11. ! A total of 4 telescope corrections were flagged in sub-array 12. ! A total of 16 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! A total of 55 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 4. ! A total of 4 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 4 telescope corrections were flagged in sub-array 11. ! A total of 4 telescope corrections were flagged in sub-array 12. ! A total of 16 telescope corrections were flagged in sub-array 13. ! A total of 8 telescope corrections were flagged in sub-array 14. ! A total of 55 telescope corrections were flagged in sub-array 15. ! ! Fit before self-cal, rms=0.273090Jy sigma=2.937471 ! Fit after self-cal, rms=0.272214Jy sigma=2.900446 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.878 mas, bmaj=4.619 mas, bpa=-24.47 degrees ! Estimated noise=1.19008 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.0271773 Jy ! Component: 100 - total flux cleaned = -0.0347808 Jy ! Component: 150 - total flux cleaned = -0.0356427 Jy ! Component: 200 - total flux cleaned = -0.0371859 Jy ! Total flux subtracted in 200 components = -0.0371859 Jy ! Clean residual min=-0.022672 max=0.024561 Jy/beam ! Clean residual mean=0.000071 rms=0.005109 Jy/beam ! Combined flux in latest and established models = 1.72179 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 214 components and 1.72179 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.263242Jy sigma=2.891449 ! Fit after self-cal, rms=0.259827Jy sigma=2.676064 ! Inverting map ! Added new window around map position (-32, 11.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0241664 Jy ! Component: 100 - total flux cleaned = 0.0307337 Jy ! Component: 150 - total flux cleaned = 0.0320794 Jy ! Component: 200 - total flux cleaned = 0.0333483 Jy ! Total flux subtracted in 200 components = 0.0333483 Jy ! Clean residual min=-0.017683 max=0.018790 Jy/beam ! Clean residual mean=0.000031 rms=0.004257 Jy/beam ! Combined flux in latest and established models = 1.75513 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 232 components and 1.75513 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.253349Jy sigma=2.538581 ! Fit after self-cal, rms=0.253551Jy sigma=2.527817 ! Inverting map ! Added new window around map position (80.5, 17.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0153841 Jy ! Component: 100 - total flux cleaned = 0.0214644 Jy ! Component: 150 - total flux cleaned = 0.0237405 Jy ! Component: 200 - total flux cleaned = 0.0259156 Jy ! Total flux subtracted in 200 components = 0.0259156 Jy ! Clean residual min=-0.016228 max=0.016666 Jy/beam ! Clean residual mean=-0.000001 rms=0.003971 Jy/beam ! Combined flux in latest and established models = 1.78105 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 250 components and 1.78105 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.249040Jy sigma=2.436428 ! Fit after self-cal, rms=0.249676Jy sigma=2.432940 ! Inverting map ! Added new window around map position (-50, 21.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0159352 Jy ! Component: 100 - total flux cleaned = 0.0231513 Jy ! Component: 150 - total flux cleaned = 0.0262379 Jy ! Component: 200 - total flux cleaned = 0.0282149 Jy ! Total flux subtracted in 200 components = 0.0282149 Jy ! Clean residual min=-0.015405 max=0.015474 Jy/beam ! Clean residual mean=-0.000028 rms=0.003705 Jy/beam ! Combined flux in latest and established models = 1.80926 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 261 components and 1.80926 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.245560Jy sigma=2.343401 ! Fit after self-cal, rms=0.246335Jy sigma=2.339203 ! Inverting map ! Added new window around map position (-61.5, 20). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0152641 Jy ! Component: 100 - total flux cleaned = 0.0225142 Jy ! Component: 150 - total flux cleaned = 0.0258938 Jy ! Component: 200 - total flux cleaned = 0.0277611 Jy ! Total flux subtracted in 200 components = 0.0277611 Jy ! Clean residual min=-0.015295 max=0.014792 Jy/beam ! Clean residual mean=-0.000053 rms=0.003467 Jy/beam ! Combined flux in latest and established models = 1.83703 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 276 components and 1.83703 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.242940Jy sigma=2.260235 ! Fit after self-cal, rms=0.243815Jy sigma=2.256855 ! Inverting map ! Added new window around map position (-12, 1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00443777 Jy ! Component: 100 - total flux cleaned = 0.00734362 Jy ! Component: 150 - total flux cleaned = 0.0100891 Jy ! Component: 200 - total flux cleaned = 0.0118439 Jy ! Total flux subtracted in 200 components = 0.0118439 Jy ! Clean residual min=-0.014720 max=0.013452 Jy/beam ! Clean residual mean=-0.000065 rms=0.003348 Jy/beam ! Combined flux in latest and established models = 1.84887 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 292 components and 1.84887 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.240525Jy sigma=2.222340 ! Fit after self-cal, rms=0.241197Jy sigma=2.221094 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.545 mas, bmaj=5.817 mas, bpa=-21.54 degrees ! Estimated noise=0.754681 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (50, 72). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0194523 Jy ! Component: 100 - total flux cleaned = 0.0257041 Jy ! Component: 150 - total flux cleaned = 0.0282978 Jy ! Component: 200 - total flux cleaned = 0.0298219 Jy ! Component: 250 - total flux cleaned = 0.0301701 Jy ! Component: 300 - total flux cleaned = 0.03017 Jy ! Component: 350 - total flux cleaned = 0.0301675 Jy ! Component: 400 - total flux cleaned = 0.0301662 Jy ! Total flux subtracted in 400 components = 0.0301662 Jy ! Clean residual min=-0.017679 max=0.012418 Jy/beam ! Clean residual mean=-0.001694 rms=0.003975 Jy/beam ! Combined flux in latest and established models = 1.87904 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 344 components and 1.87904 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.237456Jy sigma=2.115281 ! Fit after self-cal, rms=0.237340Jy sigma=2.100089 ! 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.00782272 Jy ! Component: 100 - total flux cleaned = 0.00874369 Jy ! Component: 150 - total flux cleaned = 0.00902391 Jy ! Component: 200 - total flux cleaned = 0.00902687 Jy ! Component: 250 - total flux cleaned = 0.00875691 Jy ! Component: 300 - total flux cleaned = 0.00850999 Jy ! Component: 350 - total flux cleaned = 0.00875813 Jy ! Component: 400 - total flux cleaned = 0.00875674 Jy ! Total flux subtracted in 400 components = 0.00875674 Jy ! Clean residual min=-0.016753 max=0.010468 Jy/beam ! Clean residual mean=-0.001892 rms=0.003644 Jy/beam ! Combined flux in latest and established models = 1.88779 Jy ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 363 components and 1.88779 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.234485Jy sigma=2.066966 ! Fit after self-cal, rms=0.234766Jy sigma=2.063645 ! 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 6 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 63 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 82 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 82 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 82 telescope corrections were flagged in sub-array 15. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 6 telescope corrections were flagged in sub-array 6. ! A total of 9 telescope corrections were flagged in sub-array 7. ! A total of 6 telescope corrections were flagged in sub-array 8. ! A total of 140 telescope corrections were flagged in sub-array 9. ! A total of 53 telescope corrections were flagged in sub-array 10. ! A total of 12 telescope corrections were flagged in sub-array 11. ! A total of 6 telescope corrections were flagged in sub-array 12. ! A total of 43 telescope corrections were flagged in sub-array 13. ! A total of 35 telescope corrections were flagged in sub-array 14. ! A total of 82 telescope corrections were flagged in sub-array 15. ! ! Fit before self-cal, rms=0.375264Jy sigma=2.321436 ! Fit after self-cal, rms=0.225789Jy sigma=1.998333 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.54 mas, bmaj=5.782 mas, bpa=-21.54 degrees ! Estimated noise=0.750401 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00588255 Jy ! Component: 100 - total flux cleaned = 0.0109944 Jy ! Component: 150 - total flux cleaned = 0.0147857 Jy ! Component: 200 - total flux cleaned = 0.0170543 Jy ! Component: 250 - total flux cleaned = 0.0191549 Jy ! Component: 300 - total flux cleaned = 0.0211126 Jy ! Component: 350 - total flux cleaned = 0.0220389 Jy ! Component: 400 - total flux cleaned = 0.0229324 Jy ! Total flux subtracted in 400 components = 0.0229324 Jy ! Clean residual min=-0.011953 max=0.010939 Jy/beam ! Clean residual mean=-0.000309 rms=0.003046 Jy/beam ! Combined flux in latest and established models = 1.91072 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 382 components and 1.91072 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.222434Jy sigma=1.952422 ! Fit after self-cal, rms=0.222491Jy sigma=1.946687 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.00485026 Jy ! Component: 100 - total flux cleaned = 0.00675429 Jy ! Component: 150 - total flux cleaned = 0.00741787 Jy ! Component: 200 - total flux cleaned = 0.00805551 Jy ! Component: 250 - total flux cleaned = 0.00846234 Jy ! Component: 300 - total flux cleaned = 0.00865905 Jy ! Component: 350 - total flux cleaned = 0.00905164 Jy ! Component: 400 - total flux cleaned = 0.00924388 Jy ! Total flux subtracted in 400 components = 0.00924388 Jy ! Clean residual min=-0.011721 max=0.010053 Jy/beam ! Clean residual mean=-0.000494 rms=0.002877 Jy/beam ! Combined flux in latest and established models = 1.91997 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 388 components and 1.91997 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.219853Jy sigma=1.930060 ! Fit after self-cal, rms=0.220288Jy sigma=1.928595 ! 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.220288Jy sigma=1.928595 ! Fit after self-cal, rms=0.218291Jy sigma=1.862344 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.534 mas, bmaj=5.773 mas, bpa=-21.57 degrees ! Estimated noise=0.748109 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00321832 Jy ! Component: 100 - total flux cleaned = 0.00621913 Jy ! Component: 150 - total flux cleaned = 0.00775795 Jy ! Component: 200 - total flux cleaned = 0.00857488 Jy ! Component: 250 - total flux cleaned = 0.00915293 Jy ! Component: 300 - total flux cleaned = 0.00952892 Jy ! Component: 350 - total flux cleaned = 0.00988521 Jy ! Component: 400 - total flux cleaned = 0.00988203 Jy ! Total flux subtracted in 400 components = 0.00988203 Jy ! Clean residual min=-0.007639 max=0.009590 Jy/beam ! Clean residual mean=-0.000106 rms=0.002343 Jy/beam ! Combined flux in latest and established models = 1.92985 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 404 components and 1.92985 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.215570Jy sigma=1.845635 ! Fit after self-cal, rms=0.216002Jy sigma=1.843885 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 (-41.5, 15). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00871697 Jy ! Component: 100 - total flux cleaned = 0.00876588 Jy ! Component: 150 - total flux cleaned = 0.00838928 Jy ! Component: 200 - total flux cleaned = 0.00821879 Jy ! Component: 250 - total flux cleaned = 0.00804734 Jy ! Component: 300 - total flux cleaned = 0.00788306 Jy ! Component: 350 - total flux cleaned = 0.00772215 Jy ! Component: 400 - total flux cleaned = 0.00788094 Jy ! Total flux subtracted in 400 components = 0.00788094 Jy ! Clean residual min=-0.007632 max=0.009430 Jy/beam ! Clean residual mean=-0.000257 rms=0.002175 Jy/beam ! Combined flux in latest and established models = 1.93773 Jy ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 424 components and 1.93773 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.213857Jy sigma=1.827893 ! Fit after self-cal, rms=0.214249Jy sigma=1.825654 ! Inverting map ! Added new window around map position (-54.5, -69). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00847905 Jy ! Component: 100 - total flux cleaned = 0.011036 Jy ! Component: 150 - total flux cleaned = 0.0117855 Jy ! Component: 200 - total flux cleaned = 0.0117833 Jy ! Component: 250 - total flux cleaned = 0.0117823 Jy ! Component: 300 - total flux cleaned = 0.0119447 Jy ! Component: 350 - total flux cleaned = 0.0116254 Jy ! Component: 400 - total flux cleaned = 0.011777 Jy ! Total flux subtracted in 400 components = 0.011777 Jy ! Clean residual min=-0.007478 max=0.007592 Jy/beam ! Clean residual mean=-0.000474 rms=0.001902 Jy/beam ! Combined flux in latest and established models = 1.94951 Jy ! Performing phase self-cal ! Adding 76 model components to the UV plane model. ! The established model now contains 460 components and 1.94951 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.212599Jy sigma=1.804061 ! Fit after self-cal, rms=0.212825Jy sigma=1.802341 ! Inverting map ! Added new window around map position (1, 53). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0060565 Jy ! Component: 100 - total flux cleaned = 0.00589063 Jy ! Component: 150 - total flux cleaned = 0.0058893 Jy ! Component: 200 - total flux cleaned = 0.00588986 Jy ! Component: 250 - total flux cleaned = 0.00574081 Jy ! Component: 300 - total flux cleaned = 0.00559614 Jy ! Component: 350 - total flux cleaned = 0.00545692 Jy ! Component: 400 - total flux cleaned = 0.00546155 Jy ! Total flux subtracted in 400 components = 0.00546155 Jy ! Clean residual min=-0.007434 max=0.006400 Jy/beam ! Clean residual mean=-0.000581 rms=0.001750 Jy/beam ! Combined flux in latest and established models = 1.95497 Jy ! Performing phase self-cal ! Adding 82 model components to the UV plane model. ! The established model now contains 489 components and 1.95497 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.211269Jy sigma=1.790912 ! Fit after self-cal, rms=0.211540Jy sigma=1.789894 ! 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.211540Jy sigma=1.789894 ! Fit after self-cal, rms=0.208521Jy sigma=1.663422 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.306 mas, bmaj=12.47 mas, bpa=-28.95 degrees ! Estimated noise=0.618622 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00739881 Jy ! Component: 100 - total flux cleaned = 0.0111625 Jy ! Component: 150 - total flux cleaned = 0.0123412 Jy ! Component: 200 - total flux cleaned = 0.013087 Jy ! Component: 250 - total flux cleaned = 0.0136286 Jy ! Component: 300 - total flux cleaned = 0.0137529 Jy ! Component: 350 - total flux cleaned = 0.0136393 Jy ! Component: 400 - total flux cleaned = 0.0136372 Jy ! Total flux subtracted in 400 components = 0.0136372 Jy ! Clean residual min=-0.006314 max=0.006207 Jy/beam ! Clean residual mean=-0.000409 rms=0.001996 Jy/beam ! Combined flux in latest and established models = 1.96861 Jy selfcal ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 558 components and 1.96861 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.208381Jy sigma=1.649088 ! Fit after self-cal, rms=0.208286Jy sigma=1.647741 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.00140834 Jy ! Component: 100 - total flux cleaned = 0.001633 Jy ! Component: 150 - total flux cleaned = 0.00163008 Jy ! Component: 200 - total flux cleaned = 0.00153468 Jy ! Component: 250 - total flux cleaned = 0.00135276 Jy ! Component: 300 - total flux cleaned = 0.00117656 Jy ! Component: 350 - total flux cleaned = 0.00109198 Jy ! Component: 400 - total flux cleaned = 0.00101079 Jy ! Total flux subtracted in 400 components = 0.00101079 Jy ! Clean residual min=-0.006132 max=0.005595 Jy/beam ! Clean residual mean=-0.000455 rms=0.001888 Jy/beam ! Combined flux in latest and established models = 1.96962 Jy ! Performing phase self-cal ! Adding 92 model components to the UV plane model. ! The established model now contains 611 components and 1.96962 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.208330Jy sigma=1.643397 ! Fit after self-cal, rms=0.208301Jy sigma=1.643046 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=3.533 mas, bmaj=5.781 mas, bpa=-21.6 degrees ! Estimated noise=0.747617 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.533 x 5.781 at -21.6 degrees (North through East) ! Clean map min=-0.01803 max=0.92636 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.208301Jy sigma=1.643046 ! Fit after self-cal, rms=0.209004Jy sigma=1.145644 ! Inverting map and beam ! Estimated beam: bmin=3.53 mas, bmaj=5.767 mas, bpa=-21.63 degrees ! Estimated noise=0.747071 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0020706 Jy ! Component: 100 - total flux cleaned = 0.00316462 Jy ! Component: 150 - total flux cleaned = 0.00359252 Jy ! Component: 200 - total flux cleaned = 0.0037285 Jy ! Component: 250 - total flux cleaned = 0.00386553 Jy ! Component: 300 - total flux cleaned = 0.00399679 Jy ! Component: 350 - total flux cleaned = 0.00438927 Jy ! Component: 400 - total flux cleaned = 0.00438878 Jy ! Total flux subtracted in 400 components = 0.00438878 Jy ! Clean residual min=-0.006129 max=0.005907 Jy/beam ! Clean residual mean=-0.000133 rms=0.001388 Jy/beam ! Combined flux in latest and established models = 1.97401 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 627 components and 1.97401 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207388Jy sigma=1.138716 ! Fit after self-cal, rms=0.207592Jy sigma=1.138284 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00054786 Jy ! Component: 100 - total flux cleaned = 0.000678199 Jy ! Component: 150 - total flux cleaned = 0.000677997 Jy ! Component: 200 - total flux cleaned = 0.000803433 Jy ! Component: 250 - total flux cleaned = 0.00104602 Jy ! Component: 300 - total flux cleaned = 0.00116493 Jy ! Component: 350 - total flux cleaned = 0.00104601 Jy ! Component: 400 - total flux cleaned = 0.00116488 Jy ! Total flux subtracted in 400 components = 0.00116488 Jy ! Clean residual min=-0.005954 max=0.005924 Jy/beam ! Clean residual mean=-0.000157 rms=0.001368 Jy/beam ! Combined flux in latest and established models = 1.97517 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 637 components and 1.97517 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206078Jy sigma=1.134255 ! Fit after self-cal, rms=0.206380Jy sigma=1.134056 ! Inverting map ! Added new window around map position (-83, 30.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00541697 Jy ! Component: 100 - total flux cleaned = 0.00654593 Jy ! Component: 150 - total flux cleaned = 0.00679431 Jy ! Component: 200 - total flux cleaned = 0.00691338 Jy ! Component: 250 - total flux cleaned = 0.00702922 Jy ! Component: 300 - total flux cleaned = 0.00680154 Jy ! Component: 350 - total flux cleaned = 0.00668978 Jy ! Component: 400 - total flux cleaned = 0.006578 Jy ! Total flux subtracted in 400 components = 0.006578 Jy ! Clean residual min=-0.005571 max=0.004902 Jy/beam ! Clean residual mean=-0.000256 rms=0.001299 Jy/beam ! Combined flux in latest and established models = 1.98175 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 662 components and 1.98175 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.205185Jy sigma=1.124264 ! Fit after self-cal, rms=0.205462Jy sigma=1.123135 ! 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.205462Jy sigma=1.123135 ! Fit after self-cal, rms=0.205512Jy sigma=1.120066 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.868 mas, bmaj=4.604 mas, bpa=-24.43 degrees ! Estimated noise=1.1711 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.791885 Jy ! Component: 100 - total flux cleaned = 1.07772 Jy ! Component: 150 - total flux cleaned = 1.24197 Jy ! Component: 200 - total flux cleaned = 1.35704 Jy ! Component: 250 - total flux cleaned = 1.44219 Jy ! Component: 300 - total flux cleaned = 1.5069 Jy ! Component: 350 - total flux cleaned = 1.55822 Jy ! Component: 400 - total flux cleaned = 1.59984 Jy ! Total flux subtracted in 400 components = 1.59984 Jy ! Clean residual min=-0.017226 max=0.029042 Jy/beam ! Clean residual mean=0.000387 rms=0.005643 Jy/beam ! Combined flux in latest and established models = 1.59984 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 113 model components to the UV plane model. ! The established model now contains 113 components and 1.59984 Jy ! Inverting map and beam ! Estimated beam: bmin=3.53 mas, bmaj=5.775 mas, bpa=-21.62 degrees ! Estimated noise=0.74665 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) ! Added new window around map position (-25.5, 9.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0692911 Jy ! Component: 100 - total flux cleaned = 0.121487 Jy ! Component: 150 - total flux cleaned = 0.16344 Jy ! Component: 200 - total flux cleaned = 0.198222 Jy ! Component: 250 - total flux cleaned = 0.227368 Jy ! Component: 300 - total flux cleaned = 0.252023 Jy ! Component: 350 - total flux cleaned = 0.273259 Jy ! Component: 400 - total flux cleaned = 0.291753 Jy ! Total flux subtracted in 400 components = 0.291753 Jy ! Clean residual min=-0.012200 max=0.013923 Jy/beam ! Clean residual mean=0.001578 rms=0.003771 Jy/beam ! Combined flux in latest and established models = 1.8916 Jy ! Adding 197 model components to the UV plane model. ! The established model now contains 289 components and 1.8916 Jy ! Inverting map ! Added new window around map position (-3, -50.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0148367 Jy ! Component: 100 - total flux cleaned = 0.0261046 Jy ! Component: 150 - total flux cleaned = 0.0358249 Jy ! Component: 200 - total flux cleaned = 0.0442462 Jy ! Component: 250 - total flux cleaned = 0.0519862 Jy ! Component: 300 - total flux cleaned = 0.0587015 Jy ! Component: 350 - total flux cleaned = 0.064546 Jy ! Component: 400 - total flux cleaned = 0.0682059 Jy ! Total flux subtracted in 400 components = 0.0682059 Jy ! Clean residual min=-0.008620 max=0.009422 Jy/beam ! Clean residual mean=0.000322 rms=0.002470 Jy/beam ! Combined flux in latest and established models = 1.9598 Jy ! Adding 126 model components to the UV plane model. ! The established model now contains 377 components and 1.9598 Jy ! Inverting map ! Added new window around map position (-112.5, 85.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00769952 Jy ! Component: 100 - total flux cleaned = 0.0109647 Jy ! Component: 150 - total flux cleaned = 0.0137225 Jy ! Component: 200 - total flux cleaned = 0.0169203 Jy ! Component: 250 - total flux cleaned = 0.0188583 Jy ! Component: 300 - total flux cleaned = 0.0204469 Jy ! Component: 350 - total flux cleaned = 0.0222104 Jy ! Component: 400 - total flux cleaned = 0.0236558 Jy ! Total flux subtracted in 400 components = 0.0236558 Jy ! Clean residual min=-0.007170 max=0.006497 Jy/beam ! Clean residual mean=-0.000038 rms=0.001898 Jy/beam ! Combined flux in latest and established models = 1.98346 Jy ! Adding 117 model components to the UV plane model. ! The established model now contains 435 components and 1.98346 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=9.31 mas, bmaj=12.47 mas, bpa=-29.02 degrees ! Estimated noise=0.618052 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00492163 Jy ! Component: 100 - total flux cleaned = 0.00788597 Jy ! Component: 150 - total flux cleaned = 0.00907811 Jy ! Component: 200 - total flux cleaned = 0.00973542 Jy ! Component: 250 - total flux cleaned = 0.010348 Jy ! Component: 300 - total flux cleaned = 0.0111043 Jy ! Component: 350 - total flux cleaned = 0.0114569 Jy ! Component: 400 - total flux cleaned = 0.0117941 Jy ! Total flux subtracted in 400 components = 0.0117941 Jy ! Clean residual min=-0.007280 max=0.005911 Jy/beam ! Clean residual mean=-0.000401 rms=0.001913 Jy/beam ! Combined flux in latest and established models = 1.99525 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.000317942 Jy ! Component: 100 - total flux cleaned = 0.00046716 Jy ! Component: 150 - total flux cleaned = 0.000470387 Jy ! Component: 200 - total flux cleaned = 0.000741205 Jy ! Component: 250 - total flux cleaned = 0.000870256 Jy ! Component: 300 - total flux cleaned = 0.0012322 Jy ! Component: 350 - total flux cleaned = 0.00134876 Jy ! Component: 400 - total flux cleaned = 0.00156875 Jy ! Total flux subtracted in 400 components = 0.00156875 Jy ! Clean residual min=-0.007071 max=0.005423 Jy/beam ! Clean residual mean=-0.000424 rms=0.001648 Jy/beam ! Combined flux in latest and established models = 1.99682 Jy ! Adding 164 model components to the UV plane model. ! The established model now contains 576 components and 1.99682 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.53 mas, bmaj=5.775 mas, bpa=-21.62 degrees ! Estimated noise=0.74665 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000792124 Jy ! Component: 100 - total flux cleaned = -0.000792907 Jy ! Component: 150 - total flux cleaned = -0.000643678 Jy ! Component: 200 - total flux cleaned = -0.000644018 Jy ! Component: 250 - total flux cleaned = -0.000369267 Jy ! Component: 300 - total flux cleaned = -0.000236241 Jy ! Component: 350 - total flux cleaned = -0.00010821 Jy ! Component: 400 - total flux cleaned = -0.000106718 Jy ! Total flux subtracted in 400 components = -0.000106718 Jy ! Clean residual min=-0.005886 max=0.004993 Jy/beam ! Clean residual mean=-0.000242 rms=0.001264 Jy/beam ! Combined flux in latest and established models = 1.99671 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 627 components and 1.99671 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.186538Jy sigma=1.117000 ! Fit after self-cal, rms=0.185710Jy sigma=1.106783 wmodel /scr/pima/rdv94_uvs/J1230+1223_S_map.mod ! Writing 627 model components to file: /scr/pima/rdv94_uvs/J1230+1223_S_map.mod wobs /scr/pima/rdv94_uvs/J1230+1223_S_uvs.fits ! Applying 430 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J1230+1223_S_uvs.fits wwins /scr/pima/rdv94_uvs/J1230+1223_S_map.win ! wwins: Wrote 16 windows to /scr/pima/rdv94_uvs/J1230+1223_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.00230032 Jy ! Component: 100 - total flux cleaned = -0.00345952 Jy ! Component: 150 - total flux cleaned = -0.00446776 Jy ! Component: 200 - total flux cleaned = -0.00518936 Jy ! Component: 250 - total flux cleaned = -0.00573049 Jy ! Component: 300 - total flux cleaned = -0.00603043 Jy ! Component: 350 - total flux cleaned = -0.00639738 Jy ! Component: 400 - total flux cleaned = -0.00661231 Jy ! Component: 450 - total flux cleaned = -0.00668236 Jy ! Component: 500 - total flux cleaned = -0.00682156 Jy ! Component: 550 - total flux cleaned = -0.00716249 Jy ! Component: 600 - total flux cleaned = -0.00729679 Jy ! Component: 650 - total flux cleaned = -0.00756195 Jy ! Component: 700 - total flux cleaned = -0.00782255 Jy ! Component: 750 - total flux cleaned = -0.00795052 Jy ! Component: 800 - total flux cleaned = -0.00801342 Jy ! Component: 850 - total flux cleaned = -0.00795193 Jy ! Component: 900 - total flux cleaned = -0.00819955 Jy ! Component: 950 - total flux cleaned = -0.00826052 Jy ! Component: 1000 - total flux cleaned = -0.00838139 Jy ! Component: 1050 - total flux cleaned = -0.0084404 Jy ! Component: 1100 - total flux cleaned = -0.00855741 Jy ! Component: 1150 - total flux cleaned = -0.00861508 Jy ! Component: 1200 - total flux cleaned = -0.00855665 Jy ! Component: 1250 - total flux cleaned = -0.00861348 Jy ! Component: 1300 - total flux cleaned = -0.00861309 Jy ! Component: 1350 - total flux cleaned = -0.00861302 Jy ! Component: 1400 - total flux cleaned = -0.00872415 Jy ! Component: 1450 - total flux cleaned = -0.00872406 Jy ! Component: 1500 - total flux cleaned = -0.00872461 Jy ! Component: 1550 - total flux cleaned = -0.00888721 Jy ! Component: 1600 - total flux cleaned = -0.00888739 Jy ! Component: 1650 - total flux cleaned = -0.00904734 Jy ! Component: 1700 - total flux cleaned = -0.00910017 Jy ! Component: 1750 - total flux cleaned = -0.00920535 Jy ! Component: 1800 - total flux cleaned = -0.00925757 Jy ! Component: 1850 - total flux cleaned = -0.00930934 Jy ! Component: 1900 - total flux cleaned = -0.00941234 Jy ! Component: 1950 - total flux cleaned = -0.0095656 Jy ! Component: 2000 - total flux cleaned = -0.00951488 Jy ! Component: 2050 - total flux cleaned = -0.00956533 Jy ! Component: 2100 - total flux cleaned = -0.00971559 Jy ! Component: 2150 - total flux cleaned = -0.00976545 Jy ! Component: 2200 - total flux cleaned = -0.00981486 Jy ! Component: 2250 - total flux cleaned = -0.00991302 Jy ! Component: 2300 - total flux cleaned = -0.0100104 Jy ! Component: 2350 - total flux cleaned = -0.0100102 Jy ! Component: 2400 - total flux cleaned = -0.0100103 Jy ! Component: 2450 - total flux cleaned = -0.0102493 Jy ! Component: 2500 - total flux cleaned = -0.0102488 Jy ! Component: 2550 - total flux cleaned = -0.0102487 Jy ! Component: 2600 - total flux cleaned = -0.0103428 Jy ! Component: 2650 - total flux cleaned = -0.0104362 Jy ! Component: 2700 - total flux cleaned = -0.0104362 Jy ! Component: 2750 - total flux cleaned = -0.0104366 Jy ! Component: 2800 - total flux cleaned = -0.0105285 Jy ! Component: 2850 - total flux cleaned = -0.0105284 Jy ! Component: 2900 - total flux cleaned = -0.0106192 Jy ! Component: 2950 - total flux cleaned = -0.0105743 Jy ! Component: 3000 - total flux cleaned = -0.0106643 Jy ! Component: 3050 - total flux cleaned = -0.0106197 Jy ! Component: 3100 - total flux cleaned = -0.0106198 Jy ! Component: 3150 - total flux cleaned = -0.0107527 Jy ! Component: 3200 - total flux cleaned = -0.0107525 Jy ! Component: 3250 - total flux cleaned = -0.0107527 Jy ! Component: 3300 - total flux cleaned = -0.0107964 Jy ! Component: 3350 - total flux cleaned = -0.0108397 Jy ! Component: 3400 - total flux cleaned = -0.0108827 Jy ! Component: 3450 - total flux cleaned = -0.0109255 Jy ! Component: 3500 - total flux cleaned = -0.0109256 Jy ! Component: 3550 - total flux cleaned = -0.0109255 Jy ! Component: 3600 - total flux cleaned = -0.0108833 Jy ! Component: 3650 - total flux cleaned = -0.0109256 Jy ! Component: 3700 - total flux cleaned = -0.0110096 Jy ! Component: 3750 - total flux cleaned = -0.0110515 Jy ! Component: 3800 - total flux cleaned = -0.0111345 Jy ! Component: 3850 - total flux cleaned = -0.0111346 Jy ! Component: 3900 - total flux cleaned = -0.0111345 Jy ! Component: 3950 - total flux cleaned = -0.0111755 Jy ! Component: 4000 - total flux cleaned = -0.0111346 Jy ! Component: 4050 - total flux cleaned = -0.0112567 Jy ! Total flux subtracted in 4096 components = -0.0112568 Jy ! Clean residual min=-0.002205 max=0.002480 Jy/beam ! Clean residual mean=-0.000058 rms=0.000746 Jy/beam ! Combined flux in latest and established models = 1.98546 Jy keep ! Adding 286 model components to the UV plane model. ! The established model now contains 911 components and 1.98546 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J1230+1223_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.53 x 5.775 at -21.62 degrees (North through East) ! Clean map min=-0.013428 max=0.92673 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J1230+1223_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.9267 rms= 0.00074 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J1230+1223_S_uvs.log closed on Tue Oct 8 12:18:51 2013