! Started logfile: /scr/pima/rdv94_uvs/J0309+1029_X_uvs.log on Tue Oct 8 12:10:36 2013 obs /scr/pima/rdv94_uvs/J0309+1029_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J0309+1029_X_uva.fits ! AN table 1: 24 integrations on 3 of 3 possible baselines. ! AN table 2: 156 integrations on 6 of 6 possible baselines. ! AN table 3: 10 integrations on 1 of 1 possible baselines. ! AN table 4: 22 integrations on 15 of 15 possible baselines. ! AN table 5: 22 integrations on 10 of 10 possible baselines. ! AN table 6: 22 integrations on 21 of 21 possible baselines. ! AN table 7: 22 integrations on 36 of 36 possible baselines. ! AN table 8: 5 integrations on 36 of 36 possible baselines. ! AN table 9: 5 integrations on 6 of 6 possible baselines. ! AN table 10: 3 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.729887 visibilities/baseline/integration-bin. ! Found source: J0309+1029 ! ! 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 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 9036 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.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv94_uvs/J0309+1029_X] 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:NY self-cal status=correctable weight=1 ! 1:ON self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:ON self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:ON self-cal status=correctable weight=1 ! 4:WF self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:ON self-cal status=correctable weight=1 ! 5:WF self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:ON self-cal status=correctable weight=1 ! 6:WF self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HH self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WF self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:BR 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:MK self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:WF self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:WF self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MK self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 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 46 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 25 telescope corrections were flagged in sub-array 4. ! A total of 15 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! A total of 58 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 46 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 25 telescope corrections were flagged in sub-array 4. ! A total of 15 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! A total of 58 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 294 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 52 telescope corrections were flagged in sub-array 4. ! A total of 44 telescope corrections were flagged in sub-array 5. ! A total of 66 telescope corrections were flagged in sub-array 6. ! A total of 90 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 294 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 52 telescope corrections were flagged in sub-array 4. ! A total of 44 telescope corrections were flagged in sub-array 5. ! A total of 66 telescope corrections were flagged in sub-array 6. ! A total of 90 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=1.486861Jy sigma=15.218708 ! Fit after self-cal, rms=0.297687Jy sigma=2.356399 ! 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:WZ ! Flagging baselines of antenna 1:WZ in the currently selected channels. flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:ON ! Flagging baselines of antenna 2:ON in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 4:HH ! Flagging baselines of antenna 4:HH in the currently selected channels. flag 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 4:ON ! Flagging baselines of antenna 4:ON in the currently selected channels. flag 4:WF ! Flagging baselines of antenna 4:WF in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 5:HH ! Flagging baselines of antenna 5:HH in the currently selected channels. flag 5:NY ! Flagging baselines of antenna 5:NY in the currently selected channels. flag 5:ON ! Flagging baselines of antenna 5:ON in the currently selected channels. flag 5:WF ! Flagging baselines of antenna 5:WF in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 6:HH ! Flagging baselines of antenna 6:HH in the currently selected channels. flag 6:NY ! Flagging baselines of antenna 6:NY in the currently selected channels. flag 6:ON ! Flagging baselines of antenna 6:ON in the currently selected channels. flag 6:WF ! Flagging baselines of antenna 6:WF in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. flag 7:HH ! Flagging baselines of antenna 7:HH 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:WF ! Flagging baselines of antenna 7:WF in the currently selected channels. flag 7:WZ ! Flagging baselines of antenna 7:WZ in the currently selected channels. flag 8:WF ! Flagging baselines of antenna 8:WF in the currently selected channels. flag 9:WF ! Flagging baselines of antenna 9:WF 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=0.9569 mas, bmaj=2.598 mas, bpa=-21.4 degrees ! Estimated noise=1.61418 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.818233 Jy ! Component: 100 - total flux cleaned = 0.9999 Jy ! Total flux subtracted in 100 components = 0.9999 Jy ! Clean residual min=-0.034320 max=0.060716 Jy/beam ! Clean residual mean=0.000136 rms=0.009649 Jy/beam ! Combined flux in latest and established models = 0.9999 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 9 components and 0.9999 Jy ! ! Correcting IF 1. ! A total of 63 telescope corrections were flagged in sub-array 1. ! A total of 578 telescope corrections were flagged in sub-array 2. ! A total of 107 telescope corrections were flagged in sub-array 4. ! A total of 95 telescope corrections were flagged in sub-array 5. ! A total of 120 telescope corrections were flagged in sub-array 6. ! A total of 101 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 7 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 63 telescope corrections were flagged in sub-array 1. ! A total of 578 telescope corrections were flagged in sub-array 2. ! A total of 107 telescope corrections were flagged in sub-array 4. ! A total of 95 telescope corrections were flagged in sub-array 5. ! A total of 120 telescope corrections were flagged in sub-array 6. ! A total of 101 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 7 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 330 telescope corrections were flagged in sub-array 2. ! A total of 80 telescope corrections were flagged in sub-array 4. ! A total of 66 telescope corrections were flagged in sub-array 5. ! A total of 88 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 7 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 330 telescope corrections were flagged in sub-array 2. ! A total of 80 telescope corrections were flagged in sub-array 4. ! A total of 66 telescope corrections were flagged in sub-array 5. ! A total of 88 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 4 telescope corrections were flagged in sub-array 8. ! A total of 7 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.172062Jy sigma=3.121243 ! Fit after self-cal, rms=0.169709Jy sigma=3.168729 ! Inverting map and beam ! Estimated beam: bmin=0.9546 mas, bmaj=2.592 mas, bpa=-21.4 degrees ! Estimated noise=1.60938 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0630994 Jy ! Component: 100 - total flux cleaned = 0.0974357 Jy ! Total flux subtracted in 100 components = 0.0974357 Jy ! Clean residual min=-0.032052 max=0.035215 Jy/beam ! Clean residual mean=0.000026 rms=0.007594 Jy/beam ! Combined flux in latest and established models = 1.09734 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 25 components and 1.09734 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137913Jy sigma=2.573748 ! Fit after self-cal, rms=0.137914Jy sigma=2.573593 ! 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=1.162 mas, bmaj=2.923 mas, bpa=-20.13 degrees ! Estimated noise=1.24016 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.0283857 Jy ! Component: 100 - total flux cleaned = 0.0430994 Jy ! Component: 150 - total flux cleaned = 0.0487019 Jy ! Component: 200 - total flux cleaned = 0.0508822 Jy ! Total flux subtracted in 200 components = 0.0508822 Jy ! Clean residual min=-0.031779 max=0.028560 Jy/beam ! Clean residual mean=-0.000148 rms=0.007556 Jy/beam ! Combined flux in latest and established models = 1.14822 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 47 components and 1.14822 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132113Jy sigma=2.466759 ! Fit after self-cal, rms=0.132082Jy sigma=2.466296 ! 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: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.95 HH 1.00* HN 0.94 LA 1.09 ! NL 0.99 NY 1.00* ON 1.00* WF 1.00* ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 0.97 FD 0.92 HN 0.97 KP 0.99 ! LA 1.08 MK 1.05 NL 1.04 OV 0.97 ! WF 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.00* HN 0.00* MK 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! BR 0.94 FD 0.93 KP 0.97 LA 1.08 ! MK 0.99 OV 0.99 PT 1.05 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.95 HH 1.00* HN 0.96 LA 1.03 ! NL 0.99 NY 1.00* ON 1.00* WF 1.00* ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 1.35 FD 0.91 HN 0.96 KP 0.92 ! LA 1.04 MK 1.09 NL 1.04 OV 0.96 ! WF 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.00* HN 0.00* MK 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! BR 1.36 FD 0.94 KP 0.94 LA 1.05 ! MK 1.00 OV 0.97 PT 1.06 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 0.96 HH 1.00* HN 0.95 LA 1.11 ! NL 1.02 NY 1.00* ON 1.00* WF 1.00* ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 1.02 FD 0.90 HN 0.98 KP 0.84 ! LA 1.09 MK 1.05 NL 1.06 OV 1.01 ! WF 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.00* HN 0.00* MK 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! BR 0.99 FD 0.92 KP 0.83 LA 1.10 ! MK 0.97 OV 1.03 PT 1.04 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* NY 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* NY 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* NY 0.00* ! ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 7. ! Telescope amplitude corrections in sub-array 7: ! FD 1.00 HH 1.00* HN 0.93 LA 1.10 ! NL 1.02 NY 1.00* ON 1.00* WF 1.00* ! WZ 1.00* ! ! Telescope amplitude corrections in sub-array 8: ! BR 1.00 FD 0.96 HN 0.98 KP 0.88 ! LA 1.07 MK 1.07 NL 1.08 OV 0.99 ! WF 1.00* ! ! A total of 12 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! BR 0.00* HN 0.00* MK 0.00* WF 0.00* ! ! Telescope amplitude corrections in sub-array 10: ! BR 0.96 FD 0.96 KP 0.88 LA 1.09 ! MK 1.01 OV 0.99 PT 1.07 ! ! ! Fit before self-cal, rms=0.132082Jy sigma=2.466296 ! Fit after self-cal, rms=0.070586Jy sigma=1.276472 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.204 mas, bmaj=3.092 mas, bpa=-17.31 degrees ! Estimated noise=1.25576 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.00317757 Jy ! Component: 100 - total flux cleaned = -0.00412397 Jy ! Component: 150 - total flux cleaned = -0.00440485 Jy ! Component: 200 - total flux cleaned = -0.00447248 Jy ! Total flux subtracted in 200 components = -0.00447248 Jy ! Clean residual min=-0.008657 max=0.008509 Jy/beam ! Clean residual mean=0.000009 rms=0.001997 Jy/beam ! Combined flux in latest and established models = 1.14375 Jy ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 53 components and 1.14375 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070167Jy sigma=1.265631 ! Fit after self-cal, rms=0.070230Jy sigma=1.264772 ! 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.070230Jy sigma=1.264772 ! Fit after self-cal, rms=0.069969Jy sigma=1.261692 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.204 mas, bmaj=3.1 mas, bpa=-17.32 degrees ! Estimated noise=1.25491 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00347761 Jy ! Component: 100 - total flux cleaned = -0.00387093 Jy ! Component: 150 - total flux cleaned = -0.00409084 Jy ! Component: 200 - total flux cleaned = -0.00428161 Jy ! Total flux subtracted in 200 components = -0.00428161 Jy ! Clean residual min=-0.008167 max=0.008359 Jy/beam ! Clean residual mean=0.000012 rms=0.001890 Jy/beam ! Combined flux in latest and established models = 1.13946 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 58 components and 1.13946 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.069724Jy sigma=1.255594 ! Fit after self-cal, rms=0.069720Jy sigma=1.255445 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 (35.6, 49.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00506554 Jy ! Component: 100 - total flux cleaned = 0.00621709 Jy ! Component: 150 - total flux cleaned = 0.00680007 Jy ! Component: 200 - total flux cleaned = 0.00698415 Jy ! Total flux subtracted in 200 components = 0.00698415 Jy ! Clean residual min=-0.008031 max=0.006073 Jy/beam ! Clean residual mean=0.000005 rms=0.001754 Jy/beam ! Combined flux in latest and established models = 1.14645 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 73 components and 1.14645 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068910Jy sigma=1.238534 ! Fit after self-cal, rms=0.068856Jy sigma=1.237611 ! 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.068856Jy sigma=1.237611 ! Fit after self-cal, rms=0.068638Jy sigma=1.234699 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.204 mas, bmaj=3.107 mas, bpa=-17.33 degrees ! Estimated noise=1.25361 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00119279 Jy ! Component: 100 - total flux cleaned = -0.00100525 Jy ! Component: 150 - total flux cleaned = -0.0010002 Jy ! Component: 200 - total flux cleaned = -0.000998505 Jy ! Total flux subtracted in 200 components = -0.000998505 Jy ! Clean residual min=-0.007619 max=0.006107 Jy/beam ! Clean residual mean=-0.000002 rms=0.001654 Jy/beam ! Combined flux in latest and established models = 1.14545 Jy selfcal ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 80 components and 1.14545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068382Jy sigma=1.229145 ! Fit after self-cal, rms=0.068379Jy sigma=1.228975 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 = 8.09638e-05 Jy ! Component: 100 - total flux cleaned = -6.04451e-05 Jy ! Component: 150 - total flux cleaned = -0.000196453 Jy ! Component: 200 - total flux cleaned = -0.000262433 Jy ! Total flux subtracted in 200 components = -0.000262433 Jy ! Clean residual min=-0.007441 max=0.006004 Jy/beam ! Clean residual mean=-0.000002 rms=0.001634 Jy/beam ! Combined flux in latest and established models = 1.14519 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 83 components and 1.14519 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068260Jy sigma=1.226553 ! Fit after self-cal, rms=0.068261Jy sigma=1.226483 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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.068261Jy sigma=1.226483 ! Fit after self-cal, rms=0.068162Jy sigma=1.225304 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.279 mas, bmaj=3.625 mas, bpa=-8.58 degrees ! Estimated noise=1.42434 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000183447 Jy ! Component: 100 - total flux cleaned = -0.000435172 Jy ! Component: 150 - total flux cleaned = -0.000437138 Jy ! Component: 200 - total flux cleaned = -0.000507649 Jy ! Total flux subtracted in 200 components = -0.000507649 Jy ! Clean residual min=-0.007522 max=0.006568 Jy/beam ! Clean residual mean=-0.000009 rms=0.001816 Jy/beam ! Combined flux in latest and established models = 1.14468 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 94 components and 1.14468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068080Jy sigma=1.222910 ! Fit after self-cal, rms=0.068081Jy sigma=1.222704 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 4.48819e-05 Jy ! Component: 100 - total flux cleaned = 4.6515e-05 Jy ! Component: 150 - total flux cleaned = 9.94958e-05 Jy ! Component: 200 - total flux cleaned = 4.97514e-06 Jy ! Total flux subtracted in 200 components = 4.97514e-06 Jy ! Clean residual min=-0.007440 max=0.006577 Jy/beam ! Clean residual mean=-0.000009 rms=0.001788 Jy/beam ! Combined flux in latest and established models = 1.14468 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 97 components and 1.14468 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068069Jy sigma=1.222112 ! Fit after self-cal, rms=0.068073Jy sigma=1.221989 ! 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 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:NL,true ! 6:NL 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:LA,true ! 7:LA self-cal status=fixed weight=1 selfant 7:NL,true ! 7:NL self-cal status=fixed weight=1 selfant 8:BR,true ! 8:BR 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:MK,true ! 8:MK self-cal status=fixed weight=1 selfant 8:NL,true ! 8:NL self-cal status=fixed weight=1 selfant 8:OV,true ! 8:OV self-cal status=fixed weight=1 selfant 9:BR,true ! 9:BR self-cal status=fixed weight=1 selfant 9:HN,true ! 9:HN self-cal status=fixed weight=1 selfant 9:MK,true ! 9:MK 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:KP,true ! 10:KP self-cal status=fixed weight=1 selfant 10:LA,true ! 10:LA self-cal status=fixed weight=1 selfant 10:MK,true ! 10:MK 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 ![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 72 telescope corrections were flagged in sub-array 1. ! A total of 184 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 27 telescope corrections were flagged in sub-array 6. ! A total of 27 telescope corrections were flagged in sub-array 7. ! A total of 1 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 184 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 27 telescope corrections were flagged in sub-array 6. ! A total of 27 telescope corrections were flagged in sub-array 7. ! A total of 1 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 624 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 50 telescope corrections were flagged in sub-array 4. ! A total of 50 telescope corrections were flagged in sub-array 5. ! A total of 50 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 1 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 624 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 50 telescope corrections were flagged in sub-array 4. ! A total of 50 telescope corrections were flagged in sub-array 5. ! A total of 50 telescope corrections were flagged in sub-array 6. ! A total of 50 telescope corrections were flagged in sub-array 7. ! A total of 1 telescope corrections were flagged in sub-array 8. ! A total of 1 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.366913Jy sigma=2.482051 ! Fit after self-cal, rms=0.156422Jy sigma=1.277816 selfant "",false ! 1:NY self-cal status=correctable weight=1 ! 1:ON self-cal status=correctable weight=1 ! 1:WZ self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:ON self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:ON self-cal status=correctable weight=1 ! 4:WF self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:NY self-cal status=correctable weight=1 ! 5:ON self-cal status=correctable weight=1 ! 5:WF self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:NY self-cal status=correctable weight=1 ! 6:ON self-cal status=correctable weight=1 ! 6:WF self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:FD self-cal status=correctable weight=1 ! 7:HH self-cal status=correctable weight=1 ! 7:HN self-cal status=correctable weight=1 ! 7:LA self-cal status=correctable weight=1 ! 7:NL self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WF self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:BR 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:MK self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:WF self-cal status=correctable weight=1 ! 9:BR self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:WF self-cal status=correctable weight=1 ! 10:BR self-cal status=correctable weight=1 ! 10:FD self-cal status=correctable weight=1 ! 10:KP self-cal status=correctable weight=1 ! 10:LA self-cal status=correctable weight=1 ! 10:MK self-cal status=correctable weight=1 ! 10:OV self-cal status=correctable weight=1 ! 10:PT 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 11 telescope corrections were flagged in sub-array 4. ! A total of 22 telescope corrections were flagged in sub-array 6. ! A total of 46 telescope corrections were flagged in sub-array 7. ! A total of 8 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 22 telescope corrections were flagged in sub-array 6. ! A total of 46 telescope corrections were flagged in sub-array 7. ! A total of 8 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 22 telescope corrections were flagged in sub-array 6. ! A total of 46 telescope corrections were flagged in sub-array 7. ! A total of 8 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 11 telescope corrections were flagged in sub-array 4. ! A total of 22 telescope corrections were flagged in sub-array 6. ! A total of 46 telescope corrections were flagged in sub-array 7. ! A total of 8 telescope corrections were flagged in sub-array 8. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.156422Jy sigma=1.277816 ! Fit after self-cal, rms=0.140718Jy sigma=1.175580 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 ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5727 mas, bmaj=0.7854 mas, bpa=45.04 degrees ! Estimated noise=1.24015 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 = 9.51806e-05 Jy ! Component: 100 - total flux cleaned = 0.00039293 Jy ! Component: 150 - total flux cleaned = 0.000180417 Jy ! Component: 200 - total flux cleaned = -3.13462e-05 Jy ! Total flux subtracted in 200 components = -3.13462e-05 Jy ! Clean residual min=-0.007397 max=0.007928 Jy/beam ! Clean residual mean=0.000029 rms=0.001552 Jy/beam ! Combined flux in latest and established models = 1.14465 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 105 components and 1.14465 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139476Jy sigma=1.162864 ! Fit after self-cal, rms=0.139119Jy sigma=1.157048 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Added new window around map position (0.8, -3.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00335597 Jy ! Component: 100 - total flux cleaned = 0.00301852 Jy ! Component: 150 - total flux cleaned = 0.00219007 Jy ! Component: 200 - total flux cleaned = 0.00144341 Jy ! Total flux subtracted in 200 components = 0.00144341 Jy ! Clean residual min=-0.007236 max=0.007002 Jy/beam ! Clean residual mean=0.000024 rms=0.001468 Jy/beam ! Combined flux in latest and established models = 1.1461 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 115 components and 1.1461 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138706Jy sigma=1.153076 ! Fit after self-cal, rms=0.138670Jy sigma=1.152611 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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 ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5922 mas, bmaj=0.8151 mas, bpa=56.19 degrees ! Estimated noise=0.939328 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 (-36.4, 6.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00562403 Jy ! Component: 100 - total flux cleaned = 0.00854106 Jy ! Component: 150 - total flux cleaned = 0.0100385 Jy ! Component: 200 - total flux cleaned = 0.0109386 Jy ! Component: 250 - total flux cleaned = 0.0113433 Jy ! Component: 300 - total flux cleaned = 0.0114668 Jy ! Component: 350 - total flux cleaned = 0.0115865 Jy ! Component: 400 - total flux cleaned = 0.0115871 Jy ! Total flux subtracted in 400 components = 0.0115871 Jy ! Clean residual min=-0.005630 max=0.006009 Jy/beam ! Clean residual mean=0.000225 rms=0.001430 Jy/beam ! Combined flux in latest and established models = 1.15768 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 152 components and 1.15768 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137868Jy sigma=1.145502 ! Fit after self-cal, rms=0.137851Jy sigma=1.144918 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Added new window around map position (2, 3.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00582564 Jy ! Component: 100 - total flux cleaned = 0.0081378 Jy ! Component: 150 - total flux cleaned = 0.00839203 Jy ! Component: 200 - total flux cleaned = 0.00850773 Jy ! Component: 250 - total flux cleaned = 0.00850923 Jy ! Component: 300 - total flux cleaned = 0.00850857 Jy ! Component: 350 - total flux cleaned = 0.00860558 Jy ! Component: 400 - total flux cleaned = 0.00860466 Jy ! Total flux subtracted in 400 components = 0.00860466 Jy ! Clean residual min=-0.005071 max=0.005731 Jy/beam ! Clean residual mean=0.000099 rms=0.001357 Jy/beam ! Combined flux in latest and established models = 1.16629 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 177 components and 1.16629 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137372Jy sigma=1.140872 ! Fit after self-cal, rms=0.137336Jy sigma=1.140350 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Added new window around map position (8.6, 7.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00544568 Jy ! Component: 100 - total flux cleaned = 0.00794063 Jy ! Component: 150 - total flux cleaned = 0.00826894 Jy ! Component: 200 - total flux cleaned = 0.00827449 Jy ! Component: 250 - total flux cleaned = 0.00846285 Jy ! Component: 300 - total flux cleaned = 0.00855014 Jy ! Component: 350 - total flux cleaned = 0.00872023 Jy ! Component: 400 - total flux cleaned = 0.00880142 Jy ! Total flux subtracted in 400 components = 0.00880142 Jy ! Clean residual min=-0.005045 max=0.005216 Jy/beam ! Clean residual mean=-0.000035 rms=0.001293 Jy/beam ! Combined flux in latest and established models = 1.17509 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 206 components and 1.17509 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.137004Jy sigma=1.136324 ! Fit after self-cal, rms=0.136965Jy sigma=1.135836 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Added new window around map position (25.8, -13.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00494979 Jy ! Component: 100 - total flux cleaned = 0.00737975 Jy ! Component: 150 - total flux cleaned = 0.00787552 Jy ! Component: 200 - total flux cleaned = 0.00823295 Jy ! Component: 250 - total flux cleaned = 0.00831737 Jy ! Component: 300 - total flux cleaned = 0.00848037 Jy ! Component: 350 - total flux cleaned = 0.00855886 Jy ! Component: 400 - total flux cleaned = 0.00855888 Jy ! Total flux subtracted in 400 components = 0.00855888 Jy ! Clean residual min=-0.005115 max=0.004653 Jy/beam ! Clean residual mean=-0.000132 rms=0.001226 Jy/beam ! Combined flux in latest and established models = 1.18365 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 227 components and 1.18365 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136605Jy sigma=1.132564 ! Fit after self-cal, rms=0.136583Jy sigma=1.132290 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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.00112596 Jy ! Component: 100 - total flux cleaned = 0.00120228 Jy ! Component: 150 - total flux cleaned = 0.00112683 Jy ! Component: 200 - total flux cleaned = 0.0010573 Jy ! Component: 250 - total flux cleaned = 0.00112481 Jy ! Component: 300 - total flux cleaned = 0.00105797 Jy ! Component: 350 - total flux cleaned = 0.00112289 Jy ! Component: 400 - total flux cleaned = 0.00105886 Jy ! Total flux subtracted in 400 components = 0.00105886 Jy ! Clean residual min=-0.005019 max=0.004588 Jy/beam ! Clean residual mean=-0.000148 rms=0.001201 Jy/beam ! Combined flux in latest and established models = 1.18471 Jy ! Performing phase self-cal ! Adding 69 model components to the UV plane model. ! The established model now contains 240 components and 1.18471 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.136404Jy sigma=1.131109 ! Fit after self-cal, rms=0.136400Jy sigma=1.131051 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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 72 telescope corrections were flagged in sub-array 1. ! A total of 184 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 34 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 49 telescope corrections were flagged in sub-array 6. ! A total of 73 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 4 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 184 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 34 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 49 telescope corrections were flagged in sub-array 6. ! A total of 73 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 4 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 624 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 61 telescope corrections were flagged in sub-array 4. ! A total of 50 telescope corrections were flagged in sub-array 5. ! A total of 72 telescope corrections were flagged in sub-array 6. ! A total of 96 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 4 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 72 telescope corrections were flagged in sub-array 1. ! A total of 624 telescope corrections were flagged in sub-array 2. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 61 telescope corrections were flagged in sub-array 4. ! A total of 50 telescope corrections were flagged in sub-array 5. ! A total of 72 telescope corrections were flagged in sub-array 6. ! A total of 96 telescope corrections were flagged in sub-array 7. ! A total of 9 telescope corrections were flagged in sub-array 8. ! A total of 4 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.253688Jy sigma=1.748141 ! Fit after self-cal, rms=0.136174Jy sigma=1.127045 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5936 mas, bmaj=0.8142 mas, bpa=55.92 degrees ! Estimated noise=0.940085 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00193439 Jy ! Component: 100 - total flux cleaned = -0.00224325 Jy ! Component: 150 - total flux cleaned = -0.00215456 Jy ! Component: 200 - total flux cleaned = -0.00198797 Jy ! Component: 250 - total flux cleaned = -0.0018332 Jy ! Component: 300 - total flux cleaned = -0.00183411 Jy ! Component: 350 - total flux cleaned = -0.00169287 Jy ! Component: 400 - total flux cleaned = -0.00162571 Jy ! Total flux subtracted in 400 components = -0.00162571 Jy ! Clean residual min=-0.004937 max=0.004374 Jy/beam ! Clean residual mean=-0.000165 rms=0.001147 Jy/beam ! Combined flux in latest and established models = 1.18308 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 263 components and 1.18308 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135881Jy sigma=1.124687 ! Fit after self-cal, rms=0.135880Jy sigma=1.124634 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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.000122297 Jy ! Component: 100 - total flux cleaned = 0.000188303 Jy ! Component: 150 - total flux cleaned = 0.00025107 Jy ! Component: 200 - total flux cleaned = 0.000188378 Jy ! Component: 250 - total flux cleaned = 0.000248987 Jy ! Component: 300 - total flux cleaned = 0.000248103 Jy ! Component: 350 - total flux cleaned = 0.000247037 Jy ! Component: 400 - total flux cleaned = 0.000190381 Jy ! Total flux subtracted in 400 components = 0.000190381 Jy ! Clean residual min=-0.004868 max=0.004389 Jy/beam ! Clean residual mean=-0.000168 rms=0.001133 Jy/beam ! Combined flux in latest and established models = 1.18327 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 278 components and 1.18327 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135761Jy sigma=1.123872 ! Fit after self-cal, rms=0.135761Jy sigma=1.123838 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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.135761Jy sigma=1.123838 ! Fit after self-cal, rms=0.134339Jy sigma=1.107147 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5943 mas, bmaj=0.8144 mas, bpa=55.73 degrees ! Estimated noise=0.94125 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00239693 Jy ! Component: 100 - total flux cleaned = -0.00310269 Jy ! Component: 150 - total flux cleaned = -0.00354552 Jy ! Component: 200 - total flux cleaned = -0.00379029 Jy ! Component: 250 - total flux cleaned = -0.00394799 Jy ! Component: 300 - total flux cleaned = -0.00424271 Jy ! Component: 350 - total flux cleaned = -0.00452054 Jy ! Component: 400 - total flux cleaned = -0.00471778 Jy ! Total flux subtracted in 400 components = -0.00471778 Jy ! Clean residual min=-0.003928 max=0.003344 Jy/beam ! Clean residual mean=-0.000163 rms=0.000864 Jy/beam ! Combined flux in latest and established models = 1.17855 Jy selfcal ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 304 components and 1.17855 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134109Jy sigma=1.105046 ! Fit after self-cal, rms=0.134108Jy sigma=1.104942 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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.000230221 Jy ! Component: 100 - total flux cleaned = -0.000292886 Jy ! Component: 150 - total flux cleaned = -0.000414024 Jy ! Component: 200 - total flux cleaned = -0.000530709 Jy ! Component: 250 - total flux cleaned = -0.000642956 Jy ! Component: 300 - total flux cleaned = -0.000750036 Jy ! Component: 350 - total flux cleaned = -0.000855051 Jy ! Component: 400 - total flux cleaned = -0.000956926 Jy ! Total flux subtracted in 400 components = -0.000956926 Jy ! Clean residual min=-0.003950 max=0.003266 Jy/beam ! Clean residual mean=-0.000154 rms=0.000850 Jy/beam ! Combined flux in latest and established models = 1.1776 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 319 components and 1.1776 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.134008Jy sigma=1.104072 ! Fit after self-cal, rms=0.134007Jy sigma=1.104041 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 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.134007Jy sigma=1.104041 ! Fit after self-cal, rms=0.129790Jy sigma=1.062653 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=2.427 mas, bmaj=3.193 mas, bpa=-22.02 degrees ! Estimated noise=1.30611 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000770011 Jy ! Component: 100 - total flux cleaned = 0.000690802 Jy ! Component: 150 - total flux cleaned = 0.000626084 Jy ! Component: 200 - total flux cleaned = 0.00038865 Jy ! Component: 250 - total flux cleaned = 0.000278107 Jy ! Component: 300 - total flux cleaned = 0.000123029 Jy ! Component: 350 - total flux cleaned = 2.67738e-05 Jy ! Component: 400 - total flux cleaned = -0.000107993 Jy ! Total flux subtracted in 400 components = -0.000107993 Jy ! Clean residual min=-0.004363 max=0.003644 Jy/beam ! Clean residual mean=-0.000266 rms=0.001148 Jy/beam ! Combined flux in latest and established models = 1.17749 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 346 components and 1.17749 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129783Jy sigma=1.062249 ! Fit after self-cal, rms=0.129781Jy sigma=1.062105 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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 = -8.87255e-05 Jy ! Component: 100 - total flux cleaned = -0.00027352 Jy ! Component: 150 - total flux cleaned = -0.000402103 Jy ! Component: 200 - total flux cleaned = -0.000441766 Jy ! Component: 250 - total flux cleaned = -0.000442123 Jy ! Component: 300 - total flux cleaned = -0.000478823 Jy ! Component: 350 - total flux cleaned = -0.000549591 Jy ! Component: 400 - total flux cleaned = -0.000582857 Jy ! Total flux subtracted in 400 components = -0.000582857 Jy ! Clean residual min=-0.004367 max=0.003693 Jy/beam ! Clean residual mean=-0.000253 rms=0.001136 Jy/beam ! Combined flux in latest and established models = 1.17691 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 371 components and 1.17691 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129792Jy sigma=1.062177 ! Fit after self-cal, rms=0.129792Jy sigma=1.062145 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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 ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5946 mas, bmaj=0.8143 mas, bpa=55.61 degrees ! Estimated noise=0.941668 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5946 x 0.8143 at 55.61 degrees (North through East) ! Clean map min=-0.0076934 max=0.97918 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.129792Jy sigma=1.062145 ! Fit after self-cal, rms=0.115440Jy sigma=0.906320 ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5945 mas, bmaj=0.814 mas, bpa=55.59 degrees ! Estimated noise=0.941091 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000883876 Jy ! Component: 100 - total flux cleaned = -0.00111045 Jy ! Component: 150 - total flux cleaned = -0.00118369 Jy ! Component: 200 - total flux cleaned = -0.00124716 Jy ! Component: 250 - total flux cleaned = -0.00124601 Jy ! Component: 300 - total flux cleaned = -0.00124522 Jy ! Component: 350 - total flux cleaned = -0.00113961 Jy ! Component: 400 - total flux cleaned = -0.00113916 Jy ! Total flux subtracted in 400 components = -0.00113916 Jy ! Clean residual min=-0.003816 max=0.003174 Jy/beam ! Clean residual mean=-0.000082 rms=0.000773 Jy/beam ! Combined flux in latest and established models = 1.17577 Jy ! Performing phase self-cal ! Adding 51 model components to the UV plane model. ! The established model now contains 394 components and 1.17577 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115306Jy sigma=0.905135 ! Fit after self-cal, rms=0.115311Jy sigma=0.905078 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -3.40369e-06 Jy ! Component: 100 - total flux cleaned = -5.51937e-05 Jy ! Component: 150 - total flux cleaned = 4.28605e-05 Jy ! Component: 200 - total flux cleaned = 8.95223e-05 Jy ! Component: 250 - total flux cleaned = 0.000134783 Jy ! Component: 300 - total flux cleaned = 0.0001794 Jy ! Component: 350 - total flux cleaned = 0.000179666 Jy ! Component: 400 - total flux cleaned = 0.000179341 Jy ! Total flux subtracted in 400 components = 0.000179341 Jy ! Clean residual min=-0.003746 max=0.003090 Jy/beam ! Clean residual mean=-0.000090 rms=0.000763 Jy/beam ! Combined flux in latest and established models = 1.17595 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 409 components and 1.17595 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115242Jy sigma=0.904585 ! Fit after self-cal, rms=0.115244Jy sigma=0.904553 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115244Jy sigma=0.904553 ! Fit after self-cal, rms=0.115206Jy sigma=0.904278 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 ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5738 mas, bmaj=0.7853 mas, bpa=44.28 degrees ! Estimated noise=1.24502 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.800812 Jy ! Component: 100 - total flux cleaned = 0.975932 Jy ! Component: 150 - total flux cleaned = 1.03472 Jy ! Component: 200 - total flux cleaned = 1.06795 Jy ! Component: 250 - total flux cleaned = 1.08793 Jy ! Component: 300 - total flux cleaned = 1.10131 Jy ! Component: 350 - total flux cleaned = 1.11163 Jy ! Component: 400 - total flux cleaned = 1.11945 Jy ! Total flux subtracted in 400 components = 1.11945 Jy ! Clean residual min=-0.006466 max=0.006345 Jy/beam ! Clean residual mean=0.000087 rms=0.001444 Jy/beam ! Combined flux in latest and established models = 1.11945 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 33 model components to the UV plane model. ! The established model now contains 33 components and 1.11945 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5947 mas, bmaj=0.8138 mas, bpa=55.52 degrees ! Estimated noise=0.941214 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.00933708 Jy ! Component: 100 - total flux cleaned = 0.0165592 Jy ! Component: 150 - total flux cleaned = 0.0225639 Jy ! Component: 200 - total flux cleaned = 0.0269612 Jy ! Component: 250 - total flux cleaned = 0.0308901 Jy ! Component: 300 - total flux cleaned = 0.0341143 Jy ! Component: 350 - total flux cleaned = 0.0370684 Jy ! Component: 400 - total flux cleaned = 0.0399298 Jy ! Total flux subtracted in 400 components = 0.0399298 Jy ! Clean residual min=-0.004301 max=0.003987 Jy/beam ! Clean residual mean=0.000103 rms=0.000947 Jy/beam ! Combined flux in latest and established models = 1.15938 Jy ! Adding 60 model components to the UV plane model. ! The established model now contains 84 components and 1.15938 Jy ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=2.428 mas, bmaj=3.195 mas, bpa=-21.93 degrees ! Estimated noise=1.30579 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00394272 Jy ! Component: 100 - total flux cleaned = 0.00652469 Jy ! Component: 150 - total flux cleaned = 0.0080788 Jy ! Component: 200 - total flux cleaned = 0.00939577 Jy ! Component: 250 - total flux cleaned = 0.0106331 Jy ! Component: 300 - total flux cleaned = 0.0117282 Jy ! Component: 350 - total flux cleaned = 0.0125377 Jy ! Component: 400 - total flux cleaned = 0.0132629 Jy ! Total flux subtracted in 400 components = 0.0132629 Jy ! Clean residual min=-0.004494 max=0.003827 Jy/beam ! Clean residual mean=-0.000173 rms=0.001181 Jy/beam ! Combined flux in latest and established models = 1.17264 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.000336998 Jy ! Component: 100 - total flux cleaned = 0.000713315 Jy ! Component: 150 - total flux cleaned = 0.00100992 Jy ! Component: 200 - total flux cleaned = 0.0012354 Jy ! Component: 250 - total flux cleaned = 0.00134322 Jy ! Component: 300 - total flux cleaned = 0.00149899 Jy ! Component: 350 - total flux cleaned = 0.00164795 Jy ! Component: 400 - total flux cleaned = 0.00174387 Jy ! Total flux subtracted in 400 components = 0.00174387 Jy ! Clean residual min=-0.004214 max=0.003758 Jy/beam ! Clean residual mean=-0.000227 rms=0.001141 Jy/beam ! Combined flux in latest and established models = 1.17439 Jy ! Adding 74 model components to the UV plane model. ! The established model now contains 152 components and 1.17439 Jy ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! Estimated beam: bmin=0.5947 mas, bmaj=0.8138 mas, bpa=55.52 degrees ! Estimated noise=0.941214 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000371967 Jy ! Component: 100 - total flux cleaned = -3.36057e-05 Jy ! Component: 150 - total flux cleaned = 0.00027379 Jy ! Component: 200 - total flux cleaned = 0.000270888 Jy ! Component: 250 - total flux cleaned = 0.000180445 Jy ! Component: 300 - total flux cleaned = 0.000180802 Jy ! Component: 350 - total flux cleaned = 9.85919e-05 Jy ! Component: 400 - total flux cleaned = 0.00018402 Jy ! Total flux subtracted in 400 components = 0.00018402 Jy ! Clean residual min=-0.003919 max=0.003406 Jy/beam ! Clean residual mean=-0.000068 rms=0.000847 Jy/beam ! Combined flux in latest and established models = 1.17457 Jy selfcal ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 194 components and 1.17457 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.116269Jy sigma=0.916995 ! Fit after self-cal, rms=0.116123Jy sigma=0.915410 wmodel /scr/pima/rdv94_uvs/J0309+1029_X_map.mod ! Writing 194 model components to file: /scr/pima/rdv94_uvs/J0309+1029_X_map.mod wobs /scr/pima/rdv94_uvs/J0309+1029_X_uvs.fits ! Applying 2236 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J0309+1029_X_uvs.fits wwins /scr/pima/rdv94_uvs/J0309+1029_X_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv94_uvs/J0309+1029_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00177396 Jy ! Component: 100 - total flux cleaned = -0.00307784 Jy ! Component: 150 - total flux cleaned = -0.00419172 Jy ! Component: 200 - total flux cleaned = -0.00520844 Jy ! Component: 250 - total flux cleaned = -0.00567387 Jy ! Component: 300 - total flux cleaned = -0.00612902 Jy ! Component: 350 - total flux cleaned = -0.00640782 Jy ! Component: 400 - total flux cleaned = -0.00679145 Jy ! Component: 450 - total flux cleaned = -0.00706087 Jy ! Component: 500 - total flux cleaned = -0.00732547 Jy ! Component: 550 - total flux cleaned = -0.00763759 Jy ! Component: 600 - total flux cleaned = -0.007843 Jy ! Component: 650 - total flux cleaned = -0.00809696 Jy ! Component: 700 - total flux cleaned = -0.00844817 Jy ! Component: 750 - total flux cleaned = -0.00864664 Jy ! Component: 800 - total flux cleaned = -0.00894054 Jy ! Component: 850 - total flux cleaned = -0.00903675 Jy ! Component: 900 - total flux cleaned = -0.00913245 Jy ! Component: 950 - total flux cleaned = -0.00932132 Jy ! Component: 1000 - total flux cleaned = -0.00932107 Jy ! Component: 1050 - total flux cleaned = -0.00950728 Jy ! Component: 1100 - total flux cleaned = -0.00959984 Jy ! Component: 1150 - total flux cleaned = -0.00973723 Jy ! Component: 1200 - total flux cleaned = -0.00991767 Jy ! Component: 1250 - total flux cleaned = -0.0100074 Jy ! Component: 1300 - total flux cleaned = -0.0101853 Jy ! Component: 1350 - total flux cleaned = -0.0103617 Jy ! Component: 1400 - total flux cleaned = -0.0104925 Jy ! Component: 1450 - total flux cleaned = -0.0104495 Jy ! Component: 1500 - total flux cleaned = -0.0106212 Jy ! Component: 1550 - total flux cleaned = -0.0108334 Jy ! Component: 1600 - total flux cleaned = -0.0109603 Jy ! Component: 1650 - total flux cleaned = -0.0108765 Jy ! Component: 1700 - total flux cleaned = -0.0111676 Jy ! Component: 1750 - total flux cleaned = -0.0110852 Jy ! Component: 1800 - total flux cleaned = -0.011167 Jy ! Component: 1850 - total flux cleaned = -0.0113297 Jy ! Component: 1900 - total flux cleaned = -0.0114103 Jy ! Component: 1950 - total flux cleaned = -0.0114108 Jy ! Component: 2000 - total flux cleaned = -0.0115696 Jy ! Component: 2050 - total flux cleaned = -0.0115697 Jy ! Component: 2100 - total flux cleaned = -0.0117265 Jy ! Component: 2150 - total flux cleaned = -0.0117265 Jy ! Component: 2200 - total flux cleaned = -0.0116106 Jy ! Component: 2250 - total flux cleaned = -0.0117644 Jy ! Component: 2300 - total flux cleaned = -0.0117261 Jy ! Component: 2350 - total flux cleaned = -0.0118399 Jy ! Component: 2400 - total flux cleaned = -0.011915 Jy ! Component: 2450 - total flux cleaned = -0.0118775 Jy ! Component: 2500 - total flux cleaned = -0.0119146 Jy ! Component: 2550 - total flux cleaned = -0.0119887 Jy ! Component: 2600 - total flux cleaned = -0.0119888 Jy ! Component: 2650 - total flux cleaned = -0.0120984 Jy ! Component: 2700 - total flux cleaned = -0.012062 Jy ! Component: 2750 - total flux cleaned = -0.0122064 Jy ! Component: 2800 - total flux cleaned = -0.0122064 Jy ! Component: 2850 - total flux cleaned = -0.0122417 Jy ! Component: 2900 - total flux cleaned = -0.0122062 Jy ! Component: 2950 - total flux cleaned = -0.0123119 Jy ! Component: 3000 - total flux cleaned = -0.0123116 Jy ! Component: 3050 - total flux cleaned = -0.012486 Jy ! Component: 3100 - total flux cleaned = -0.0124508 Jy ! Component: 3150 - total flux cleaned = -0.012416 Jy ! Component: 3200 - total flux cleaned = -0.0123817 Jy ! Component: 3250 - total flux cleaned = -0.0125522 Jy ! Component: 3300 - total flux cleaned = -0.0125182 Jy ! Component: 3350 - total flux cleaned = -0.0125182 Jy ! Component: 3400 - total flux cleaned = -0.0125852 Jy ! Component: 3450 - total flux cleaned = -0.0126521 Jy ! Component: 3500 - total flux cleaned = -0.0127518 Jy ! Component: 3550 - total flux cleaned = -0.0128508 Jy ! Component: 3600 - total flux cleaned = -0.0128179 Jy ! Component: 3650 - total flux cleaned = -0.012916 Jy ! Component: 3700 - total flux cleaned = -0.0128508 Jy ! Component: 3750 - total flux cleaned = -0.0129159 Jy ! Component: 3800 - total flux cleaned = -0.01311 Jy ! Component: 3850 - total flux cleaned = -0.0131099 Jy ! Component: 3900 - total flux cleaned = -0.0131739 Jy ! Component: 3950 - total flux cleaned = -0.013174 Jy ! Component: 4000 - total flux cleaned = -0.0132056 Jy ! Component: 4050 - total flux cleaned = -0.0133006 Jy ! Total flux subtracted in 4096 components = -0.0132375 Jy ! Clean residual min=-0.002173 max=0.001935 Jy/beam ! Clean residual mean=0.000002 rms=0.000521 Jy/beam ! Combined flux in latest and established models = 1.16133 Jy keep ! Adding 349 model components to the UV plane model. ! The established model now contains 543 components and 1.16133 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J0309+1029_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.5947 x 0.8138 at 55.52 degrees (North through East) ! Clean map min=-0.004176 max=1.0115 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J0309+1029_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 7.97% of the data. ! The x-axis pixel size should ideally be below 0.1764 milli-arcsec ! The y-axis pixel size should ideally be below 0.1868 milli-arcsec 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= 1.0114 rms= 0.00051 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J0309+1029_X_uvs.log closed on Tue Oct 8 12:10:43 2013