! Started logfile: /scr/pima/rdv62_uvs/J0559+2353_X_uvs.log on Fri Nov 22 22:34:20 2013 obs /scr/pima/rdv62_uvs/J0559+2353_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J0559+2353_X_uva.fits ! AN table 1: 6 integrations on 10 of 10 possible baselines. ! AN table 2: 9 integrations on 21 of 21 possible baselines. ! AN table 3: 7 integrations on 36 of 36 possible baselines. ! AN table 4: 14 integrations on 1 of 1 possible baselines. ! AN table 5: 15 integrations on 15 of 15 possible baselines. ! AN table 6: 2 integrations on 21 of 21 possible baselines. ! AN table 7: 47 integrations on 3 of 3 possible baselines. ! AN table 8: 34 integrations on 1 of 1 possible baselines. ! AN table 9: 3 integrations on 3 of 3 possible baselines. ! Apparent sampling: 0.819876 visibilities/baseline/integration-bin. ! Found source: J0559+2353 ! ! 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 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3168 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/rdv62_uvs/J0559+2353_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:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:KK 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:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9:WZ 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 2 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 28 telescope corrections were flagged in sub-array 4. ! A total of 21 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 2 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 28 telescope corrections were flagged in sub-array 4. ! A total of 21 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 2 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 28 telescope corrections were flagged in sub-array 4. ! A total of 21 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 2 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 28 telescope corrections were flagged in sub-array 4. ! A total of 21 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=1.057384Jy sigma=4.024227 ! Fit after self-cal, rms=0.655922Jy sigma=2.538035 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 3:KK ! Flagging baselines of antenna 3:KK in the currently selected channels. flag 4:KK ! Flagging baselines of antenna 4:KK in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 7:KK ! Flagging baselines of antenna 7:KK in the currently selected channels. flag 7:NY ! Flagging baselines of antenna 7:NY in the currently selected channels. flag 8:NY ! Flagging baselines of antenna 8:NY in the currently selected channels. flag 8:WZ ! Flagging baselines of antenna 8:WZ in the currently selected channels. flag 9:NY ! Flagging baselines of antenna 9:NY in the currently selected channels. flag 9:WZ ! Flagging baselines of antenna 9:WZ in the currently selected channels. ![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=1.466 mas, bmaj=2.376 mas, bpa=-21.46 degrees ! Estimated noise=5.87307 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.268973 Jy ! Component: 100 - total flux cleaned = 0.32766 Jy ! Total flux subtracted in 100 components = 0.32766 Jy ! Clean residual min=-0.010276 max=0.016814 Jy/beam ! Clean residual mean=-0.000015 rms=0.002470 Jy/beam ! Combined flux in latest and established models = 0.32766 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.32766 Jy ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 5. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 5. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 5. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 11 telescope corrections were flagged in sub-array 5. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.084931Jy sigma=0.342533 ! Fit after self-cal, rms=0.084815Jy sigma=0.342143 ! Inverting map and beam ! Estimated beam: bmin=1.466 mas, bmaj=2.376 mas, bpa=-21.46 degrees ! Estimated noise=5.87378 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0155076 Jy ! Component: 100 - total flux cleaned = 0.0220251 Jy ! Total flux subtracted in 100 components = 0.0220251 Jy ! Clean residual min=-0.010105 max=0.008709 Jy/beam ! Clean residual mean=-0.000020 rms=0.002116 Jy/beam ! Combined flux in latest and established models = 0.349685 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 19 components and 0.349685 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079901Jy sigma=0.318013 ! Fit after self-cal, rms=0.079912Jy sigma=0.317974 ! 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.737 mas, bmaj=2.788 mas, bpa=-17.9 degrees ! Estimated noise=3.86192 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.013061 Jy ! Component: 100 - total flux cleaned = 0.0185371 Jy ! Component: 150 - total flux cleaned = 0.0211688 Jy ! Component: 200 - total flux cleaned = 0.0221773 Jy ! Total flux subtracted in 200 components = 0.0221773 Jy ! Clean residual min=-0.008668 max=0.007022 Jy/beam ! Clean residual mean=-0.000014 rms=0.001953 Jy/beam ! Combined flux in latest and established models = 0.371862 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 41 components and 0.371862 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078236Jy sigma=0.309043 ! Fit after self-cal, rms=0.078252Jy sigma=0.309007 ! 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.94 FD 1.03 LA 0.98 NL 1.05 ! PT 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 0.96 HN 0.98 KP 0.99 ! LA 0.92 NL 1.01 PT 1.00 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 0.98 KK 1.00* KP 1.00 ! LA 0.94 MK 1.01 NL 0.98 OV 0.99 ! PT 0.94 ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.93 KK 1.00* KP 0.96 MK 0.95 ! OV 1.00 PT 0.98 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.06 FD 1.10 KP 1.10 LA 1.00 ! NL 0.99 OV 1.01 PT 1.01 ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HN 0.00* NY 0.00* WZ 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.96 FD 1.03 LA 0.98 NL 0.93 ! PT 1.11 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 0.99 HN 0.99 KP 0.94 ! LA 0.99 NL 0.97 PT 1.07 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.98 FD 0.98 KK 1.00* KP 0.93 ! LA 0.93 MK 0.96 NL 0.96 OV 0.98 ! PT 1.02 ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.96 KK 1.00* KP 0.93 MK 0.99 ! OV 0.96 PT 1.01 ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.93 FD 1.14 KP 1.04 LA 0.94 ! NL 1.04 OV 1.04 PT 1.16 ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HN 0.00* NY 0.00* WZ 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.14 FD 1.08 LA 1.03 NL 0.98 ! PT 1.09 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.09 FD 1.03 HN 0.99 KP 0.86 ! LA 0.99 NL 1.02 PT 1.09 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.02 FD 1.06 KK 1.00* KP 0.92 ! LA 0.98 MK 1.04 NL 0.98 OV 1.05 ! PT 1.03 ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.06 KK 1.00* KP 0.88 MK 1.10 ! OV 1.01 PT 1.02 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.06 FD 1.19 KP 0.89 LA 1.08 ! NL 1.14 OV 1.06 PT 1.14 ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HN 0.00* NY 0.00* WZ 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.08 FD 1.09 LA 1.06 NL 1.03 ! PT 1.18 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 FD 1.09 HN 1.06 KP 0.93 ! LA 1.01 NL 1.07 PT 1.07 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 1.05 KK 1.00* KP 0.98 ! LA 1.00 MK 1.02 NL 1.04 OV 1.02 ! PT 1.05 ! ! Telescope amplitude corrections in sub-array 4: ! KK 0.00* MK 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 KK 1.00* KP 0.92 MK 1.00 ! OV 1.07 PT 1.03 ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.19 FD 1.10 KP 1.04 LA 1.02 ! NL 1.00 OV 1.07 PT 1.12 ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* MK 0.00* NY 0.00* ! ! Telescope amplitude corrections in sub-array 8: ! NY 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 9: ! HN 0.00* NY 0.00* WZ 0.00* ! ! ! Fit before self-cal, rms=0.078252Jy sigma=0.309007 ! Fit after self-cal, rms=0.072697Jy sigma=0.281400 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.735 mas, bmaj=2.778 mas, bpa=-17.56 degrees ! Estimated noise=3.87901 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 (16.6, -44.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00422982 Jy ! Component: 100 - total flux cleaned = 0.00512974 Jy ! Component: 150 - total flux cleaned = 0.00527826 Jy ! Component: 200 - total flux cleaned = 0.00547106 Jy ! Total flux subtracted in 200 components = 0.00547106 Jy ! Clean residual min=-0.004226 max=0.003857 Jy/beam ! Clean residual mean=-0.000020 rms=0.000944 Jy/beam ! Combined flux in latest and established models = 0.377333 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 58 components and 0.377333 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072402Jy sigma=0.280229 ! Fit after self-cal, rms=0.072417Jy sigma=0.280107 ! 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.072417Jy sigma=0.280107 ! Fit after self-cal, rms=0.072277Jy sigma=0.280043 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.734 mas, bmaj=2.778 mas, bpa=-17.59 degrees ! Estimated noise=3.87229 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0011752 Jy ! Component: 100 - total flux cleaned = 0.00128931 Jy ! Component: 150 - total flux cleaned = 0.00134033 Jy ! Component: 200 - total flux cleaned = 0.00138734 Jy ! Total flux subtracted in 200 components = 0.00138734 Jy ! Clean residual min=-0.004057 max=0.003818 Jy/beam ! Clean residual mean=-0.000020 rms=0.000920 Jy/beam ! Combined flux in latest and established models = 0.378721 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 78 components and 0.37872 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072194Jy sigma=0.279728 ! Fit after self-cal, rms=0.072188Jy sigma=0.279716 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.000337673 Jy ! Component: 100 - total flux cleaned = 0.000420168 Jy ! Component: 150 - total flux cleaned = 0.000458405 Jy ! Component: 200 - total flux cleaned = 0.000494836 Jy ! Total flux subtracted in 200 components = 0.000494836 Jy ! Clean residual min=-0.004012 max=0.003777 Jy/beam ! Clean residual mean=-0.000020 rms=0.000911 Jy/beam ! Combined flux in latest and established models = 0.379215 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 89 components and 0.379215 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072148Jy sigma=0.279561 ! Fit after self-cal, rms=0.072145Jy sigma=0.279558 ! 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.072145Jy sigma=0.279558 ! Fit after self-cal, rms=0.072071Jy sigma=0.279524 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.735 mas, bmaj=2.778 mas, bpa=-17.57 degrees ! Estimated noise=3.86768 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000178408 Jy ! Component: 100 - total flux cleaned = -0.000141974 Jy ! Component: 150 - total flux cleaned = -0.000141404 Jy ! Component: 200 - total flux cleaned = -0.000107027 Jy ! Total flux subtracted in 200 components = -0.000107027 Jy ! Clean residual min=-0.003911 max=0.003745 Jy/beam ! Clean residual mean=-0.000020 rms=0.000902 Jy/beam ! Combined flux in latest and established models = 0.379108 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 95 components and 0.379108 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072036Jy sigma=0.279393 ! Fit after self-cal, rms=0.072035Jy sigma=0.279390 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 = 3.16688e-05 Jy ! Component: 100 - total flux cleaned = 6.30548e-05 Jy ! Component: 150 - total flux cleaned = 6.32448e-05 Jy ! Component: 200 - total flux cleaned = 6.29351e-05 Jy ! Total flux subtracted in 200 components = 6.29351e-05 Jy ! Clean residual min=-0.003888 max=0.003707 Jy/beam ! Clean residual mean=-0.000020 rms=0.000897 Jy/beam ! Combined flux in latest and established models = 0.379171 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 99 components and 0.379171 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072011Jy sigma=0.279306 ! Fit after self-cal, rms=0.072010Jy sigma=0.279304 ! 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.072010Jy sigma=0.279304 ! Fit after self-cal, rms=0.071976Jy sigma=0.279289 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.724 mas, bmaj=3.475 mas, bpa=6.098 degrees ! Estimated noise=4.09297 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -8.37994e-05 Jy ! Component: 100 - total flux cleaned = -0.000106436 Jy ! Component: 150 - total flux cleaned = -0.000147701 Jy ! Component: 200 - total flux cleaned = -0.000147746 Jy ! Total flux subtracted in 200 components = -0.000147746 Jy ! Clean residual min=-0.003707 max=0.003747 Jy/beam ! Clean residual mean=-0.000029 rms=0.000951 Jy/beam ! Combined flux in latest and established models = 0.379023 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 113 components and 0.379023 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.071965Jy sigma=0.279238 ! Fit after self-cal, rms=0.071966Jy sigma=0.279232 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 = -2.64877e-05 Jy ! Component: 100 - total flux cleaned = -8.84592e-06 Jy ! Component: 150 - total flux cleaned = 7.81926e-06 Jy ! Component: 200 - total flux cleaned = 8.66907e-06 Jy ! Total flux subtracted in 200 components = 8.66907e-06 Jy ! Clean residual min=-0.003699 max=0.003713 Jy/beam ! Clean residual mean=-0.000029 rms=0.000950 Jy/beam ! Combined flux in latest and established models = 0.379032 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 121 components and 0.379032 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.071962Jy sigma=0.279211 ! Fit after self-cal, rms=0.071963Jy sigma=0.279209 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL 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:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT 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: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 3:PT,true ! 3:PT 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:KP,true ! 5:KP 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 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:FD,true ! 6:FD self-cal status=fixed weight=1 selfant 6:KP,true ! 6:KP self-cal status=fixed weight=1 selfant 6:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:NL,true ! 6:NL self-cal status=fixed weight=1 selfant 6:OV,true ! 6:OV self-cal status=fixed weight=1 selfant 6:PT,true ! 6:PT self-cal status=fixed weight=1 selfant 7:MK,true ! 7:MK self-cal status=fixed weight=1 selfant 9:HN,true ! 9:HN 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 1 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 94 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 6 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 94 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 6 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 94 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 6 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 94 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 6 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.148041Jy sigma=0.407005 ! Fit after self-cal, rms=0.082761Jy sigma=0.285680 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:KK 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:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:KP self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:MK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 8:NY self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:HN self-cal status=correctable weight=1 ! 9:NY self-cal status=correctable weight=1 ! 9: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 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 8 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 47 telescope corrections were flagged in sub-array 7. ! A total of 3 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.082761Jy sigma=0.285680 ! Fit after self-cal, rms=0.081967Jy sigma=0.281942 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.217 mas, bmaj=2.336 mas, bpa=-20.51 degrees ! Estimated noise=5.84239 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 (44.8, 12.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00149542 Jy ! Component: 100 - total flux cleaned = 0.00208543 Jy ! Component: 150 - total flux cleaned = 0.00259232 Jy ! Component: 200 - total flux cleaned = 0.00319745 Jy ! Total flux subtracted in 200 components = 0.00319745 Jy ! Clean residual min=-0.005220 max=0.005723 Jy/beam ! Clean residual mean=-0.000014 rms=0.001377 Jy/beam ! Combined flux in latest and established models = 0.38223 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 147 components and 0.38223 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.081752Jy sigma=0.281949 ! Fit after self-cal, rms=0.081352Jy sigma=0.281275 ! 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=1.455 mas, bmaj=2.716 mas, bpa=-18.63 degrees ! Estimated noise=3.76839 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 (-34.8, 48.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00319026 Jy ! Component: 100 - total flux cleaned = 0.00265725 Jy ! Component: 150 - total flux cleaned = 0.00227993 Jy ! Component: 200 - total flux cleaned = 0.00210949 Jy ! Component: 250 - total flux cleaned = 0.00187899 Jy ! Component: 300 - total flux cleaned = 0.00188096 Jy ! Component: 350 - total flux cleaned = 0.00194859 Jy ! Component: 400 - total flux cleaned = 0.00201008 Jy ! Total flux subtracted in 400 components = 0.00201008 Jy ! Clean residual min=-0.003348 max=0.003598 Jy/beam ! Clean residual mean=-0.000022 rms=0.000874 Jy/beam ! Combined flux in latest and established models = 0.38424 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 191 components and 0.38424 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.080940Jy sigma=0.279596 ! Fit after self-cal, rms=0.080937Jy sigma=0.279521 ! Inverting map ! Added new window around map position (-11.4, 7.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0033765 Jy ! Component: 100 - total flux cleaned = 0.00320374 Jy ! Component: 150 - total flux cleaned = 0.00311715 Jy ! Component: 200 - total flux cleaned = 0.0029649 Jy ! Component: 250 - total flux cleaned = 0.00275598 Jy ! Component: 300 - total flux cleaned = 0.00262693 Jy ! Component: 350 - total flux cleaned = 0.00238127 Jy ! Component: 400 - total flux cleaned = 0.00226482 Jy ! Total flux subtracted in 400 components = 0.00226482 Jy ! Clean residual min=-0.003049 max=0.002968 Jy/beam ! Clean residual mean=-0.000018 rms=0.000798 Jy/beam ! Combined flux in latest and established models = 0.386505 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 230 components and 0.386504 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.080558Jy sigma=0.278414 ! Fit after self-cal, rms=0.080539Jy sigma=0.278345 ! 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.000397543 Jy ! Component: 100 - total flux cleaned = 0.000164612 Jy ! Component: 150 - total flux cleaned = -1.38055e-06 Jy ! Component: 200 - total flux cleaned = -0.000108173 Jy ! Component: 250 - total flux cleaned = -0.000263641 Jy ! Component: 300 - total flux cleaned = -0.000313122 Jy ! Component: 350 - total flux cleaned = -0.000411475 Jy ! Component: 400 - total flux cleaned = -0.000505998 Jy ! Total flux subtracted in 400 components = -0.000505998 Jy ! Clean residual min=-0.002922 max=0.002868 Jy/beam ! Clean residual mean=-0.000017 rms=0.000762 Jy/beam ! Combined flux in latest and established models = 0.385998 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 259 components and 0.385998 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.080387Jy sigma=0.277864 ! Fit after self-cal, rms=0.080384Jy sigma=0.277846 ! Inverting map ! Added new window around map position (-40.4, -35.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00237769 Jy ! Component: 100 - total flux cleaned = 0.00291372 Jy ! Component: 150 - total flux cleaned = 0.00304652 Jy ! Component: 200 - total flux cleaned = 0.00292888 Jy ! Component: 250 - total flux cleaned = 0.0027106 Jy ! Component: 300 - total flux cleaned = 0.00255403 Jy ! Component: 350 - total flux cleaned = 0.00250205 Jy ! Component: 400 - total flux cleaned = 0.00230919 Jy ! Total flux subtracted in 400 components = 0.00230919 Jy ! Clean residual min=-0.002828 max=0.002804 Jy/beam ! Clean residual mean=-0.000017 rms=0.000728 Jy/beam ! Combined flux in latest and established models = 0.388308 Jy ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 300 components and 0.388308 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.080145Jy sigma=0.277138 ! Fit after self-cal, rms=0.080142Jy sigma=0.277103 ! Inverting map ! Added new window around map position (31.2, -45). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00117289 Jy ! Component: 100 - total flux cleaned = 0.00182353 Jy ! Component: 150 - total flux cleaned = 0.0021026 Jy ! Component: 200 - total flux cleaned = 0.00216586 Jy ! Component: 250 - total flux cleaned = 0.0022809 Jy ! Component: 300 - total flux cleaned = 0.00244218 Jy ! Component: 350 - total flux cleaned = 0.00264792 Jy ! Component: 400 - total flux cleaned = 0.00254788 Jy ! Total flux subtracted in 400 components = 0.00254788 Jy ! Clean residual min=-0.002936 max=0.002635 Jy/beam ! Clean residual mean=-0.000014 rms=0.000688 Jy/beam ! Combined flux in latest and established models = 0.390856 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 333 components and 0.390856 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079880Jy sigma=0.276301 ! Fit after self-cal, rms=0.079872Jy sigma=0.276261 ! 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. 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 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 6 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 28 telescope corrections were flagged in sub-array 4. ! A total of 27 telescope corrections were flagged in sub-array 5. ! A total of 3 telescope corrections were flagged in sub-array 6. ! A total of 141 telescope corrections were flagged in sub-array 7. ! A total of 68 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.147923Jy sigma=0.403633 ! Fit after self-cal, rms=0.079560Jy sigma=0.275973 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.453 mas, bmaj=2.717 mas, bpa=-18.67 degrees ! Estimated noise=3.76459 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00103708 Jy ! Component: 100 - total flux cleaned = 0.00125861 Jy ! Component: 150 - total flux cleaned = 0.0013125 Jy ! Component: 200 - total flux cleaned = 0.0013629 Jy ! Component: 250 - total flux cleaned = 0.00136249 Jy ! Component: 300 - total flux cleaned = 0.00141081 Jy ! Component: 350 - total flux cleaned = 0.00145684 Jy ! Component: 400 - total flux cleaned = 0.00141127 Jy ! Total flux subtracted in 400 components = 0.00141127 Jy ! Clean residual min=-0.002693 max=0.002474 Jy/beam ! Clean residual mean=-0.000013 rms=0.000641 Jy/beam ! Combined flux in latest and established models = 0.392267 Jy selfcal ! Performing phase self-cal ! Adding 57 model components to the UV plane model. ! The established model now contains 353 components and 0.392267 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079421Jy sigma=0.275542 ! Fit after self-cal, rms=0.079418Jy sigma=0.275531 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.000146494 Jy ! Component: 100 - total flux cleaned = 0.000280154 Jy ! Component: 150 - total flux cleaned = 0.000235479 Jy ! Component: 200 - total flux cleaned = 0.000321893 Jy ! Component: 250 - total flux cleaned = 0.000279097 Jy ! Component: 300 - total flux cleaned = 0.000403418 Jy ! Component: 350 - total flux cleaned = 0.000402257 Jy ! Component: 400 - total flux cleaned = 0.000322216 Jy ! Total flux subtracted in 400 components = 0.000322216 Jy ! Clean residual min=-0.002698 max=0.002433 Jy/beam ! Clean residual mean=-0.000013 rms=0.000626 Jy/beam ! Combined flux in latest and established models = 0.392589 Jy ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 380 components and 0.392589 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079320Jy sigma=0.275243 ! Fit after self-cal, rms=0.079319Jy sigma=0.275238 ! 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.079319Jy sigma=0.275238 ! Fit after self-cal, rms=0.079317Jy sigma=0.275206 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.452 mas, bmaj=2.717 mas, bpa=-18.68 degrees ! Estimated noise=3.76494 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000630805 Jy ! Component: 100 - total flux cleaned = 0.000672148 Jy ! Component: 150 - total flux cleaned = 0.000713435 Jy ! Component: 200 - total flux cleaned = 0.000713633 Jy ! Component: 250 - total flux cleaned = 0.000752989 Jy ! Component: 300 - total flux cleaned = 0.000752849 Jy ! Component: 350 - total flux cleaned = 0.000753223 Jy ! Component: 400 - total flux cleaned = 0.000791992 Jy ! Total flux subtracted in 400 components = 0.000791992 Jy ! Clean residual min=-0.002646 max=0.002397 Jy/beam ! Clean residual mean=-0.000012 rms=0.000609 Jy/beam ! Combined flux in latest and established models = 0.393381 Jy selfcal ! Performing phase self-cal ! Adding 91 model components to the UV plane model. ! The established model now contains 402 components and 0.393381 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079229Jy sigma=0.274942 ! Fit after self-cal, rms=0.079228Jy sigma=0.274937 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.000237023 Jy ! Component: 100 - total flux cleaned = 0.000273881 Jy ! Component: 150 - total flux cleaned = 0.000198975 Jy ! Component: 200 - total flux cleaned = 0.000272582 Jy ! Component: 250 - total flux cleaned = 0.000272625 Jy ! Component: 300 - total flux cleaned = 0.000308293 Jy ! Component: 350 - total flux cleaned = 0.000308146 Jy ! Component: 400 - total flux cleaned = 0.00030751 Jy ! Total flux subtracted in 400 components = 0.00030751 Jy ! Clean residual min=-0.002654 max=0.002360 Jy/beam ! Clean residual mean=-0.000012 rms=0.000599 Jy/beam ! Combined flux in latest and established models = 0.393689 Jy ! Performing phase self-cal ! Adding 112 model components to the UV plane model. ! The established model now contains 435 components and 0.393688 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.079158Jy sigma=0.274727 ! Fit after self-cal, rms=0.079157Jy sigma=0.274724 ! 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, 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.079157Jy sigma=0.274724 ! Fit after self-cal, rms=0.077644Jy sigma=0.269658 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.602 mas, bmaj=3.492 mas, bpa=5.389 degrees ! Estimated noise=4.05451 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000998335 Jy ! Component: 100 - total flux cleaned = 0.000962811 Jy ! Component: 150 - total flux cleaned = 0.000926943 Jy ! Component: 200 - total flux cleaned = 0.000750064 Jy ! Component: 250 - total flux cleaned = 0.000648652 Jy ! Component: 300 - total flux cleaned = 0.000551121 Jy ! Component: 350 - total flux cleaned = 0.000488169 Jy ! Component: 400 - total flux cleaned = 0.000395909 Jy ! Total flux subtracted in 400 components = 0.000395909 Jy ! Clean residual min=-0.002382 max=0.002191 Jy/beam ! Clean residual mean=-0.000018 rms=0.000596 Jy/beam ! Combined flux in latest and established models = 0.394084 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 473 components and 0.394084 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.077593Jy sigma=0.269460 ! Fit after self-cal, rms=0.077595Jy sigma=0.269442 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.00011609 Jy ! Component: 100 - total flux cleaned = 5.65861e-05 Jy ! Component: 150 - total flux cleaned = -1.43209e-06 Jy ! Component: 200 - total flux cleaned = -0.000143302 Jy ! Component: 250 - total flux cleaned = -0.000280823 Jy ! Component: 300 - total flux cleaned = -0.000334576 Jy ! Component: 350 - total flux cleaned = -0.000413803 Jy ! Component: 400 - total flux cleaned = -0.000517173 Jy ! Total flux subtracted in 400 components = -0.000517173 Jy ! Clean residual min=-0.002323 max=0.002192 Jy/beam ! Clean residual mean=-0.000017 rms=0.000581 Jy/beam ! Combined flux in latest and established models = 0.393567 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 507 components and 0.393567 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.077566Jy sigma=0.269328 ! Fit after self-cal, rms=0.077564Jy sigma=0.269324 ! 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=1.451 mas, bmaj=2.717 mas, bpa=-18.68 degrees ! Estimated noise=3.77324 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.451 x 2.717 at -18.68 degrees (North through East) ! Clean map min=-0.0038305 max=0.35032 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 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.213 mas, bmaj=2.338 mas, bpa=-20.58 degrees ! Estimated noise=5.85365 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.272407 Jy ! Component: 100 - total flux cleaned = 0.332881 Jy ! Component: 150 - total flux cleaned = 0.351108 Jy ! Component: 200 - total flux cleaned = 0.356504 Jy ! Component: 250 - total flux cleaned = 0.359879 Jy ! Component: 300 - total flux cleaned = 0.362473 Jy ! Component: 350 - total flux cleaned = 0.365293 Jy ! Component: 400 - total flux cleaned = 0.367303 Jy ! Total flux subtracted in 400 components = 0.367303 Jy ! Clean residual min=-0.005422 max=0.005650 Jy/beam ! Clean residual mean=-0.000013 rms=0.001410 Jy/beam ! Combined flux in latest and established models = 0.367303 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 64 model components to the UV plane model. ! The established model now contains 64 components and 0.367303 Jy ! Inverting map and beam ! Estimated beam: bmin=1.451 mas, bmaj=2.717 mas, bpa=-18.68 degrees ! Estimated noise=3.77324 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.00660191 Jy ! Component: 100 - total flux cleaned = 0.0114607 Jy ! Component: 150 - total flux cleaned = 0.0142354 Jy ! Component: 200 - total flux cleaned = 0.0162379 Jy ! Component: 250 - total flux cleaned = 0.0172778 Jy ! Component: 300 - total flux cleaned = 0.0176394 Jy ! Component: 350 - total flux cleaned = 0.0179828 Jy ! Component: 400 - total flux cleaned = 0.0184178 Jy ! Total flux subtracted in 400 components = 0.0184178 Jy ! Clean residual min=-0.003465 max=0.003128 Jy/beam ! Clean residual mean=-0.000018 rms=0.000847 Jy/beam ! Combined flux in latest and established models = 0.385721 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 127 components and 0.385721 Jy ! Inverting map 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 ! Estimated beam: bmin=2.602 mas, bmaj=3.492 mas, bpa=5.389 degrees ! Estimated noise=4.05451 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000863754 Jy ! Component: 100 - total flux cleaned = 0.00134648 Jy ! Component: 150 - total flux cleaned = 0.00195435 Jy ! Component: 200 - total flux cleaned = 0.00252246 Jy ! Component: 250 - total flux cleaned = 0.00298192 Jy ! Component: 300 - total flux cleaned = 0.00319878 Jy ! Component: 350 - total flux cleaned = 0.00354229 Jy ! Component: 400 - total flux cleaned = 0.00367375 Jy ! Total flux subtracted in 400 components = 0.00367375 Jy ! Clean residual min=-0.002949 max=0.002827 Jy/beam ! Clean residual mean=-0.000022 rms=0.000745 Jy/beam ! Combined flux in latest and established models = 0.389395 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.000256502 Jy ! Component: 100 - total flux cleaned = 0.000625112 Jy ! Component: 150 - total flux cleaned = 0.000743755 Jy ! Component: 200 - total flux cleaned = 0.000918721 Jy ! Component: 250 - total flux cleaned = 0.00103352 Jy ! Component: 300 - total flux cleaned = 0.00125404 Jy ! Component: 350 - total flux cleaned = 0.00141547 Jy ! Component: 400 - total flux cleaned = 0.00157285 Jy ! Total flux subtracted in 400 components = 0.00157285 Jy ! Clean residual min=-0.002739 max=0.002601 Jy/beam ! Clean residual mean=-0.000021 rms=0.000691 Jy/beam ! Combined flux in latest and established models = 0.390968 Jy ! Adding 116 model components to the UV plane model. ! The established model now contains 230 components and 0.390967 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.451 mas, bmaj=2.717 mas, bpa=-18.68 degrees ! Estimated noise=3.77324 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000491475 Jy ! Component: 100 - total flux cleaned = -0.000714182 Jy ! Component: 150 - total flux cleaned = -0.000643818 Jy ! Component: 200 - total flux cleaned = -0.0006433 Jy ! Component: 250 - total flux cleaned = -0.000642754 Jy ! Component: 300 - total flux cleaned = -0.000580297 Jy ! Component: 350 - total flux cleaned = -0.000518609 Jy ! Component: 400 - total flux cleaned = -0.000338126 Jy ! Total flux subtracted in 400 components = -0.000338126 Jy ! Clean residual min=-0.002785 max=0.002592 Jy/beam ! Clean residual mean=-0.000012 rms=0.000670 Jy/beam ! Combined flux in latest and established models = 0.390629 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 280 components and 0.390629 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078168Jy sigma=0.271331 ! Fit after self-cal, rms=0.078122Jy sigma=0.271179 wmodel /scr/pima/rdv62_uvs/J0559+2353_X_map.mod ! Writing 280 model components to file: /scr/pima/rdv62_uvs/J0559+2353_X_map.mod wobs /scr/pima/rdv62_uvs/J0559+2353_X_uvs.fits ! Applying 2136 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J0559+2353_X_uvs.fits wwins /scr/pima/rdv62_uvs/J0559+2353_X_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv62_uvs/J0559+2353_X_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.00074775 Jy ! Component: 100 - total flux cleaned = -0.000839238 Jy ! Component: 150 - total flux cleaned = -0.00105967 Jy ! Component: 200 - total flux cleaned = -0.00135818 Jy ! Component: 250 - total flux cleaned = -0.00148239 Jy ! Component: 300 - total flux cleaned = -0.0014027 Jy ! Component: 350 - total flux cleaned = -0.00140369 Jy ! Component: 400 - total flux cleaned = -0.00144203 Jy ! Component: 450 - total flux cleaned = -0.00140318 Jy ! Component: 500 - total flux cleaned = -0.00132892 Jy ! Component: 550 - total flux cleaned = -0.00140185 Jy ! Component: 600 - total flux cleaned = -0.00140184 Jy ! Component: 650 - total flux cleaned = -0.00133102 Jy ! Component: 700 - total flux cleaned = -0.00136612 Jy ! Component: 750 - total flux cleaned = -0.00136616 Jy ! Component: 800 - total flux cleaned = -0.00126492 Jy ! Component: 850 - total flux cleaned = -0.00133224 Jy ! Component: 900 - total flux cleaned = -0.00133216 Jy ! Component: 950 - total flux cleaned = -0.00133248 Jy ! Component: 1000 - total flux cleaned = -0.00139732 Jy ! Component: 1050 - total flux cleaned = -0.00142946 Jy ! Component: 1100 - total flux cleaned = -0.00139788 Jy ! Component: 1150 - total flux cleaned = -0.00149172 Jy ! Component: 1200 - total flux cleaned = -0.00155388 Jy ! Component: 1250 - total flux cleaned = -0.00155399 Jy ! Component: 1300 - total flux cleaned = -0.00161486 Jy ! Component: 1350 - total flux cleaned = -0.0017056 Jy ! Component: 1400 - total flux cleaned = -0.001795 Jy ! Component: 1450 - total flux cleaned = -0.00188384 Jy ! Component: 1500 - total flux cleaned = -0.00191313 Jy ! Component: 1550 - total flux cleaned = -0.00194264 Jy ! Component: 1600 - total flux cleaned = -0.00214356 Jy ! Component: 1650 - total flux cleaned = -0.00211483 Jy ! Component: 1700 - total flux cleaned = -0.00214314 Jy ! Component: 1750 - total flux cleaned = -0.00217111 Jy ! Component: 1800 - total flux cleaned = -0.0022268 Jy ! Component: 1850 - total flux cleaned = -0.0022542 Jy ! Component: 1900 - total flux cleaned = -0.0022268 Jy ! Component: 1950 - total flux cleaned = -0.0023084 Jy ! Component: 2000 - total flux cleaned = -0.00230837 Jy ! Component: 2050 - total flux cleaned = -0.00228148 Jy ! Component: 2100 - total flux cleaned = -0.0024141 Jy ! Component: 2150 - total flux cleaned = -0.00230834 Jy ! Component: 2200 - total flux cleaned = -0.00230834 Jy ! Component: 2250 - total flux cleaned = -0.00228233 Jy ! Component: 2300 - total flux cleaned = -0.00233428 Jy ! Component: 2350 - total flux cleaned = -0.00230837 Jy ! Component: 2400 - total flux cleaned = -0.00246198 Jy ! Component: 2450 - total flux cleaned = -0.00253811 Jy ! Component: 2500 - total flux cleaned = -0.00251268 Jy ! Component: 2550 - total flux cleaned = -0.00251267 Jy ! Component: 2600 - total flux cleaned = -0.00253793 Jy ! Component: 2650 - total flux cleaned = -0.00258751 Jy ! Component: 2700 - total flux cleaned = -0.0026125 Jy ! Component: 2750 - total flux cleaned = -0.00256323 Jy ! Component: 2800 - total flux cleaned = -0.00261218 Jy ! Component: 2850 - total flux cleaned = -0.00268498 Jy ! Component: 2900 - total flux cleaned = -0.00268503 Jy ! Component: 2950 - total flux cleaned = -0.00266103 Jy ! Component: 3000 - total flux cleaned = -0.00263715 Jy ! Component: 3050 - total flux cleaned = -0.00263698 Jy ! Component: 3100 - total flux cleaned = -0.00249579 Jy ! Component: 3150 - total flux cleaned = -0.00261314 Jy ! Component: 3200 - total flux cleaned = -0.00258977 Jy ! Component: 3250 - total flux cleaned = -0.00261277 Jy ! Component: 3300 - total flux cleaned = -0.00254362 Jy ! Component: 3350 - total flux cleaned = -0.00261242 Jy ! Component: 3400 - total flux cleaned = -0.00245264 Jy ! Component: 3450 - total flux cleaned = -0.00252095 Jy ! Component: 3500 - total flux cleaned = -0.00252066 Jy ! Component: 3550 - total flux cleaned = -0.00245337 Jy ! Component: 3600 - total flux cleaned = -0.00245339 Jy ! Component: 3650 - total flux cleaned = -0.00245361 Jy ! Component: 3700 - total flux cleaned = -0.00256455 Jy ! Component: 3750 - total flux cleaned = -0.00245408 Jy ! Component: 3800 - total flux cleaned = -0.00241011 Jy ! Component: 3850 - total flux cleaned = -0.00247582 Jy ! Component: 3900 - total flux cleaned = -0.0024975 Jy ! Component: 3950 - total flux cleaned = -0.00243241 Jy ! Component: 4000 - total flux cleaned = -0.00234603 Jy ! Component: 4050 - total flux cleaned = -0.00230336 Jy ! Total flux subtracted in 4096 components = -0.00232472 Jy ! Clean residual min=-0.001527 max=0.001468 Jy/beam ! Clean residual mean=-0.000007 rms=0.000390 Jy/beam ! Combined flux in latest and established models = 0.388305 Jy keep ! Adding 285 model components to the UV plane model. ! The established model now contains 564 components and 0.388305 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J0559+2353_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.451 x 2.717 at -18.68 degrees (North through East) ! Clean map min=-0.003246 max=0.3501 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J0559+2353_X_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.35 rms= 0.00038 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J0559+2353_X_uvs.log closed on Fri Nov 22 22:34:27 2013