! Started logfile: /scr/pima/rdv73_uvs/J0607+6720_X_uvs.log on Sat Jan 9 17:31:31 2016 obs /scr/pima/rdv73_uvs/J0607+6720_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J0607+6720_X_uva.fits ! AN table 1: 28 integrations on 36 of 36 possible baselines. ! AN table 2: 25 integrations on 66 of 66 possible baselines. ! AN table 3: 7 integrations on 3 of 3 possible baselines. ! AN table 4: 31 integrations on 91 of 91 possible baselines. ! AN table 5: 2 integrations on 1 of 1 possible baselines. ! AN table 6: 13 integrations on 3 of 3 possible baselines. ! AN table 7: 60 integrations on 6 of 6 possible baselines. ! Apparent sampling: 0.733604 visibilities/baseline/integration-bin. ! Found source: J0607+6720 ! ! 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 17316 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/rdv73_uvs/J0607+6720_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:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:ZC 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:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC 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 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 108 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 108 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 108 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 15 telescope corrections were flagged in sub-array 6. ! A total of 108 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=1.078236Jy sigma=9.167384 ! Fit after self-cal, rms=0.739094Jy sigma=6.156323 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:NY ! Flagging baselines of antenna 1:NY in the currently selected channels. flag 1:ZC ! Flagging baselines of antenna 1:ZC in the currently selected channels. flag 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. 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 4:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 4:ZC ! Flagging baselines of antenna 4:ZC in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 5:ZC ! Flagging baselines of antenna 5:ZC in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. flag 7:KK ! Flagging baselines of antenna 7:KK in the currently selected channels. flag 7:NY ! Flagging baselines of antenna 7:NY in the currently selected channels. flag 7:WZ ! Flagging baselines of antenna 7:WZ in the currently selected channels. flag 7:ZC ! Flagging baselines of antenna 7:ZC in the currently selected channels. ![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.7282 mas, bmaj=1.567 mas, bpa=-19.08 degrees ! Estimated noise=1.6244 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.284202 Jy ! Component: 100 - total flux cleaned = 0.402903 Jy ! Total flux subtracted in 100 components = 0.402903 Jy ! Clean residual min=-0.021348 max=0.058031 Jy/beam ! Clean residual mean=0.000004 rms=0.006063 Jy/beam ! Combined flux in latest and established models = 0.402903 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 14 components and 0.402903 Jy ! ! Correcting IF 1. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 132 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 132 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 132 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 1. ! A total of 58 telescope corrections were flagged in sub-array 2. ! A total of 12 telescope corrections were flagged in sub-array 3. ! A total of 103 telescope corrections were flagged in sub-array 4. ! A total of 24 telescope corrections were flagged in sub-array 6. ! A total of 132 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.195700Jy sigma=1.718000 ! Fit after self-cal, rms=0.195196Jy sigma=1.715033 ! Inverting map and beam ! Estimated beam: bmin=0.7245 mas, bmaj=1.598 mas, bpa=-19.18 degrees ! Estimated noise=1.6137 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0689002 Jy ! Component: 100 - total flux cleaned = 0.11167 Jy ! Total flux subtracted in 100 components = 0.11167 Jy ! Clean residual min=-0.017222 max=0.029089 Jy/beam ! Clean residual mean=0.000006 rms=0.004282 Jy/beam ! Combined flux in latest and established models = 0.514573 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 35 components and 0.514573 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.141931Jy sigma=1.213472 ! Fit after self-cal, rms=0.141558Jy sigma=1.211106 ! Inverting map ! Added new window around map position (-0.4, 2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.033182 Jy ! Component: 100 - total flux cleaned = 0.0561789 Jy ! Total flux subtracted in 100 components = 0.0561789 Jy ! Clean residual min=-0.015834 max=0.021174 Jy/beam ! Clean residual mean=0.000007 rms=0.003693 Jy/beam ! Combined flux in latest and established models = 0.570752 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 41 components and 0.570752 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123324Jy sigma=1.036602 ! Fit after self-cal, rms=0.123051Jy sigma=1.035015 ! 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.021 mas, bmaj=2.085 mas, bpa=-20.86 degrees ! Estimated noise=0.885615 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.4, -1.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.040449 Jy ! Component: 100 - total flux cleaned = 0.0653252 Jy ! Component: 150 - total flux cleaned = 0.0818322 Jy ! Component: 200 - total flux cleaned = 0.0930194 Jy ! Total flux subtracted in 200 components = 0.0930194 Jy ! Clean residual min=-0.016497 max=0.013901 Jy/beam ! Clean residual mean=0.000018 rms=0.003497 Jy/beam ! Combined flux in latest and established models = 0.663772 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 74 components and 0.663771 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.107996Jy sigma=0.892477 ! Fit after self-cal, rms=0.107908Jy sigma=0.891673 ! 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.96 FD 0.94 HN 1.12 KP 0.99 ! LA 1.05 NY 1.00* OV 0.97 PT 0.98 ! ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 0.93 HN 1.83 KP 0.99 ! LA 1.03 NL 1.12 NY 1.00* OV 0.94 ! PT 0.98 SC 0.94 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.87 FD 0.95 HN 2.12 KK 1.00* ! KP 0.94 LA 1.07 MK 0.97 NL 1.19 ! NY 1.00* OV 0.94 PT 0.94 SC 0.96 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.95 FD 0.96 HN 1.13 KP 0.94 ! LA 1.02 NY 1.00* OV 0.94 PT 1.04 ! ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.92 FD 0.93 HN 1.85 KP 0.92 ! LA 1.01 NL 1.09 NY 1.00* OV 0.90 ! PT 1.06 SC 1.09 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.86 FD 0.95 HN 2.21 KK 1.00* ! KP 0.91 LA 1.06 MK 0.97 NL 1.17 ! NY 1.00* OV 0.92 PT 1.00 SC 1.00 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.97 HN 0.99 KP 0.87 ! LA 1.07 NY 1.00* OV 0.99 PT 1.06 ! ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 0.92 HN 1.85 KP 0.85 ! LA 1.07 NL 1.09 NY 1.00* OV 0.95 ! PT 1.05 SC 1.09 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.92 FD 0.95 HN 2.06 KK 1.00* ! KP 0.85 LA 1.13 MK 0.97 NL 1.19 ! NY 1.00* OV 0.96 PT 1.01 SC 1.01 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.99 HN 1.02 KP 0.94 ! LA 1.07 NY 1.00* OV 0.99 PT 1.05 ! ZC 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 0.97 HN 1.86 KP 0.92 ! LA 1.05 NL 1.09 NY 1.00* OV 0.93 ! PT 1.06 SC 1.07 WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.91 FD 0.98 HN 2.06 KK 1.00* ! KP 0.89 LA 1.14 MK 1.00 NL 1.18 ! NY 1.00* OV 0.96 PT 0.99 SC 1.02 ! WZ 1.00* ZC 1.00* ! ! Telescope amplitude corrections in sub-array 5: ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! KK 0.00* NY 0.00* WZ 0.00* ZC 0.00* ! ! ! Fit before self-cal, rms=0.107908Jy sigma=0.891673 ! Fit after self-cal, rms=0.102171Jy sigma=0.709666 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.031 mas, bmaj=2.141 mas, bpa=-18.9 degrees ! Estimated noise=0.920947 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.00524336 Jy ! Component: 100 - total flux cleaned = -0.00707405 Jy ! Component: 150 - total flux cleaned = -0.00777832 Jy ! Component: 200 - total flux cleaned = -0.00857212 Jy ! Total flux subtracted in 200 components = -0.00857212 Jy ! Clean residual min=-0.008440 max=0.005479 Jy/beam ! Clean residual mean=0.000005 rms=0.001277 Jy/beam ! Combined flux in latest and established models = 0.655199 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 91 components and 0.655199 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100151Jy sigma=0.696173 ! Fit after self-cal, rms=0.100175Jy sigma=0.695292 ! 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.100175Jy sigma=0.695292 ! Fit after self-cal, rms=0.102486Jy sigma=0.692679 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.03 mas, bmaj=2.15 mas, bpa=-18.39 degrees ! Estimated noise=0.927876 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0038286 Jy ! Component: 100 - total flux cleaned = -0.00515612 Jy ! Component: 150 - total flux cleaned = -0.00606677 Jy ! Component: 200 - total flux cleaned = -0.00667569 Jy ! Total flux subtracted in 200 components = -0.00667569 Jy ! Clean residual min=-0.007019 max=0.004741 Jy/beam ! Clean residual mean=0.000004 rms=0.001112 Jy/beam ! Combined flux in latest and established models = 0.648524 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 99 components and 0.648524 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101368Jy sigma=0.685618 ! Fit after self-cal, rms=0.101244Jy sigma=0.685187 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.00113914 Jy ! Component: 100 - total flux cleaned = -0.00113635 Jy ! Component: 150 - total flux cleaned = -0.00103446 Jy ! Component: 200 - total flux cleaned = -0.00112827 Jy ! Total flux subtracted in 200 components = -0.00112827 Jy ! Clean residual min=-0.006403 max=0.004924 Jy/beam ! Clean residual mean=0.000003 rms=0.001072 Jy/beam ! Combined flux in latest and established models = 0.647395 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 112 components and 0.647395 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100883Jy sigma=0.683169 ! Fit after self-cal, rms=0.100794Jy sigma=0.682875 ! 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.100794Jy sigma=0.682875 ! Fit after self-cal, rms=0.103381Jy sigma=0.680324 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.029 mas, bmaj=2.157 mas, bpa=-17.95 degrees ! Estimated noise=0.934015 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00314698 Jy ! Component: 100 - total flux cleaned = -0.0040779 Jy ! Component: 150 - total flux cleaned = -0.00459326 Jy ! Component: 200 - total flux cleaned = -0.00492393 Jy ! Total flux subtracted in 200 components = -0.00492393 Jy ! Clean residual min=-0.005681 max=0.004083 Jy/beam ! Clean residual mean=0.000004 rms=0.000951 Jy/beam ! Combined flux in latest and established models = 0.642471 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 115 components and 0.642472 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102811Jy sigma=0.676937 ! Fit after self-cal, rms=0.102735Jy sigma=0.676672 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.000728413 Jy ! Component: 100 - total flux cleaned = -0.000725648 Jy ! Component: 150 - total flux cleaned = -0.0007942 Jy ! Component: 200 - total flux cleaned = -0.000792583 Jy ! Total flux subtracted in 200 components = -0.000792583 Jy ! Clean residual min=-0.005256 max=0.004179 Jy/beam ! Clean residual mean=0.000004 rms=0.000927 Jy/beam ! Combined flux in latest and established models = 0.641679 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 119 components and 0.641679 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102525Jy sigma=0.675612 ! Fit after self-cal, rms=0.102459Jy sigma=0.675411 ! 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.102459Jy sigma=0.675411 ! Fit after self-cal, rms=0.104548Jy sigma=0.674053 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.496 mas, bmaj=3.817 mas, bpa=-26.4 degrees ! Estimated noise=1.06182 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0033975 Jy ! Component: 100 - total flux cleaned = -0.00402869 Jy ! Component: 150 - total flux cleaned = -0.00440428 Jy ! Component: 200 - total flux cleaned = -0.00467263 Jy ! Total flux subtracted in 200 components = -0.00467263 Jy ! Clean residual min=-0.003949 max=0.003313 Jy/beam ! Clean residual mean=0.000007 rms=0.000932 Jy/beam ! Combined flux in latest and established models = 0.637006 Jy selfcal ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 126 components and 0.637006 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104384Jy sigma=0.672861 ! Fit after self-cal, rms=0.104395Jy sigma=0.672738 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.000535063 Jy ! Component: 100 - total flux cleaned = -0.000486666 Jy ! Component: 150 - total flux cleaned = -0.000445546 Jy ! Component: 200 - total flux cleaned = -0.000446227 Jy ! Total flux subtracted in 200 components = -0.000446227 Jy ! Clean residual min=-0.003920 max=0.003282 Jy/beam ! Clean residual mean=0.000006 rms=0.000924 Jy/beam ! Combined flux in latest and established models = 0.63656 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 131 components and 0.63656 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.104304Jy sigma=0.672340 ! Fit after self-cal, rms=0.104284Jy sigma=0.672245 ! 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. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2: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:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 6:OV,true ! 6:OV self-cal status=fixed weight=1 ![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 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 16 telescope corrections were flagged in sub-array 1. ! A total of 17 telescope corrections were flagged in sub-array 2. ! A total of 7 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 13 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.160365Jy sigma=1.203088 ! Fit after self-cal, rms=0.410782Jy sigma=0.841092 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NY self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:ZC 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:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:KK self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 5:ZC self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 7:KK self-cal status=correctable weight=1 ! 7:NY self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 7:ZC self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 35 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 35 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 35 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 35 telescope corrections were flagged in sub-array 1. ! A total of 38 telescope corrections were flagged in sub-array 2. ! A total of 14 telescope corrections were flagged in sub-array 3. ! A total of 22 telescope corrections were flagged in sub-array 4. ! A total of 26 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.410782Jy sigma=0.841092 ! Fit after self-cal, rms=0.411822Jy sigma=0.817675 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 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.5565 mas, bmaj=0.6681 mas, bpa=73.87 degrees ! Estimated noise=1.72268 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.0085322 Jy ! Component: 100 - total flux cleaned = 0.00427215 Jy ! Component: 150 - total flux cleaned = 0.00263352 Jy ! Component: 200 - total flux cleaned = 0.000535439 Jy ! Total flux subtracted in 200 components = 0.000535439 Jy ! Clean residual min=-0.012926 max=0.014763 Jy/beam ! Clean residual mean=0.000006 rms=0.002336 Jy/beam ! Combined flux in latest and established models = 0.637095 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 138 components and 0.637096 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.406616Jy sigma=0.793235 ! Fit after self-cal, rms=0.405397Jy sigma=0.784709 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00311336 Jy ! Component: 100 - total flux cleaned = 0.00436676 Jy ! Component: 150 - total flux cleaned = 0.00210429 Jy ! Component: 200 - total flux cleaned = 1.53587e-05 Jy ! Total flux subtracted in 200 components = 1.53587e-05 Jy ! Clean residual min=-0.011543 max=0.013890 Jy/beam ! Clean residual mean=0.000006 rms=0.002232 Jy/beam ! Combined flux in latest and established models = 0.637111 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 146 components and 0.637111 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.402560Jy sigma=0.778326 ! Fit after self-cal, rms=0.401276Jy sigma=0.777711 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (2.2, 1.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00808676 Jy ! Component: 100 - total flux cleaned = 0.00761827 Jy ! Component: 150 - total flux cleaned = 0.00572194 Jy ! Component: 200 - total flux cleaned = 0.00437143 Jy ! Total flux subtracted in 200 components = 0.00437143 Jy ! Clean residual min=-0.011644 max=0.012568 Jy/beam ! Clean residual mean=0.000006 rms=0.002184 Jy/beam ! Combined flux in latest and established models = 0.641482 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 153 components and 0.641483 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.398854Jy sigma=0.773089 ! Fit after self-cal, rms=0.397928Jy sigma=0.772387 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-3.2, 0.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00898004 Jy ! Component: 100 - total flux cleaned = 0.0103417 Jy ! Component: 150 - total flux cleaned = 0.00950435 Jy ! Component: 200 - total flux cleaned = 0.00992978 Jy ! Total flux subtracted in 200 components = 0.00992978 Jy ! Clean residual min=-0.011752 max=0.011812 Jy/beam ! Clean residual mean=0.000006 rms=0.002151 Jy/beam ! Combined flux in latest and established models = 0.651412 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 162 components and 0.651412 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.396099Jy sigma=0.770827 ! Fit after self-cal, rms=0.395511Jy sigma=0.770073 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.8247 mas, bmaj=0.8859 mas, bpa=64.89 degrees ! Estimated noise=0.882199 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=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00934385 Jy ! Component: 100 - total flux cleaned = -0.0120131 Jy ! Component: 150 - total flux cleaned = -0.0129683 Jy ! Component: 200 - total flux cleaned = -0.0138577 Jy ! Component: 250 - total flux cleaned = -0.0144879 Jy ! Component: 300 - total flux cleaned = -0.0148871 Jy ! Component: 350 - total flux cleaned = -0.01508 Jy ! Component: 400 - total flux cleaned = -0.015269 Jy ! Total flux subtracted in 400 components = -0.015269 Jy ! Clean residual min=-0.007374 max=0.008351 Jy/beam ! Clean residual mean=0.000008 rms=0.001502 Jy/beam ! Combined flux in latest and established models = 0.636144 Jy ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 177 components and 0.636144 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.393952Jy sigma=0.763215 ! Fit after self-cal, rms=0.393448Jy sigma=0.762531 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-6.2, -1.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00113485 Jy ! Component: 100 - total flux cleaned = -0.000738908 Jy ! Component: 150 - total flux cleaned = -0.000974943 Jy ! Component: 200 - total flux cleaned = -0.000983821 Jy ! Component: 250 - total flux cleaned = -0.000990515 Jy ! Component: 300 - total flux cleaned = -0.000990917 Jy ! Component: 350 - total flux cleaned = -0.000802714 Jy ! Component: 400 - total flux cleaned = -0.000623045 Jy ! Total flux subtracted in 400 components = -0.000623045 Jy ! Clean residual min=-0.007144 max=0.008322 Jy/beam ! Clean residual mean=0.000007 rms=0.001492 Jy/beam ! Combined flux in latest and established models = 0.635521 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 198 components and 0.63552 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.390977Jy sigma=0.759986 ! Fit after self-cal, rms=0.390584Jy sigma=0.759486 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-5.6, 15.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000788151 Jy ! Component: 100 - total flux cleaned = -2.74928e-05 Jy ! Component: 150 - total flux cleaned = -0.000745092 Jy ! Component: 200 - total flux cleaned = -0.00118229 Jy ! Component: 250 - total flux cleaned = -0.00138153 Jy ! Component: 300 - total flux cleaned = -0.00157795 Jy ! Component: 350 - total flux cleaned = -0.00195122 Jy ! Component: 400 - total flux cleaned = -0.00231093 Jy ! Total flux subtracted in 400 components = -0.00231093 Jy ! Clean residual min=-0.007333 max=0.007812 Jy/beam ! Clean residual mean=0.000008 rms=0.001467 Jy/beam ! Combined flux in latest and established models = 0.633209 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 222 components and 0.633209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.387910Jy sigma=0.756136 ! Fit after self-cal, rms=0.387642Jy sigma=0.755603 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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.000225759 Jy ! Component: 100 - total flux cleaned = -0.000782036 Jy ! Component: 150 - total flux cleaned = -0.000781858 Jy ! Component: 200 - total flux cleaned = -0.00129018 Jy ! Component: 250 - total flux cleaned = -0.00145935 Jy ! Component: 300 - total flux cleaned = -0.00178789 Jy ! Component: 350 - total flux cleaned = -0.00195031 Jy ! Component: 400 - total flux cleaned = -0.00226662 Jy ! Total flux subtracted in 400 components = -0.00226662 Jy ! Clean residual min=-0.007413 max=0.007778 Jy/beam ! Clean residual mean=0.000008 rms=0.001466 Jy/beam ! Combined flux in latest and established models = 0.630943 Jy ! Performing phase self-cal ! Adding 72 model components to the UV plane model. ! The established model now contains 231 components and 0.630943 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.385970Jy sigma=0.754572 ! Fit after self-cal, rms=0.385738Jy sigma=0.754239 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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. ! 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 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 2. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 3. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Correcting IF 4. ! A total of 51 telescope corrections were flagged in sub-array 1. ! A total of 55 telescope corrections were flagged in sub-array 2. ! A total of 21 telescope corrections were flagged in sub-array 3. ! A total of 45 telescope corrections were flagged in sub-array 4. ! A total of 4 telescope corrections were flagged in sub-array 5. ! A total of 39 telescope corrections were flagged in sub-array 6. ! A total of 168 telescope corrections were flagged in sub-array 7. ! ! Fit before self-cal, rms=0.382266Jy sigma=0.784829 ! Fit after self-cal, rms=0.200211Jy sigma=0.734467 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.7761 mas, bmaj=0.8624 mas, bpa=76.78 degrees ! Estimated noise=0.858048 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00276148 Jy ! Component: 100 - total flux cleaned = 0.00246155 Jy ! Component: 150 - total flux cleaned = 0.0016203 Jy ! Component: 200 - total flux cleaned = 0.000850909 Jy ! Component: 250 - total flux cleaned = -9.42188e-05 Jy ! Component: 300 - total flux cleaned = -0.00052791 Jy ! Component: 350 - total flux cleaned = -0.000954496 Jy ! Component: 400 - total flux cleaned = -0.00154641 Jy ! Total flux subtracted in 400 components = -0.00154641 Jy ! Clean residual min=-0.005344 max=0.006038 Jy/beam ! Clean residual mean=0.000006 rms=0.001145 Jy/beam ! Combined flux in latest and established models = 0.629396 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 240 components and 0.629396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.197629Jy sigma=0.725212 ! Fit after self-cal, rms=0.196743Jy sigma=0.724332 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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) ! Added new window around map position (9.4, 0.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00215288 Jy ! Component: 100 - total flux cleaned = 0.00175992 Jy ! Component: 150 - total flux cleaned = 0.000999611 Jy ! Component: 200 - total flux cleaned = 0.000457931 Jy ! Component: 250 - total flux cleaned = -0.000220894 Jy ! Component: 300 - total flux cleaned = -0.000876737 Jy ! Component: 350 - total flux cleaned = -0.00135314 Jy ! Component: 400 - total flux cleaned = -0.00165609 Jy ! Total flux subtracted in 400 components = -0.00165609 Jy ! Clean residual min=-0.005107 max=0.005986 Jy/beam ! Clean residual mean=0.000006 rms=0.001127 Jy/beam ! Combined flux in latest and established models = 0.62774 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 260 components and 0.62774 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.195291Jy sigma=0.720724 ! Fit after self-cal, rms=0.195024Jy sigma=0.720271 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-3.8, 11.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000924966 Jy ! Component: 100 - total flux cleaned = 0.000349399 Jy ! Component: 150 - total flux cleaned = 0.00104209 Jy ! Component: 200 - total flux cleaned = 0.00121003 Jy ! Component: 250 - total flux cleaned = 0.00120666 Jy ! Component: 300 - total flux cleaned = 0.000752793 Jy ! Component: 350 - total flux cleaned = 0.000753944 Jy ! Component: 400 - total flux cleaned = 0.000320954 Jy ! Total flux subtracted in 400 components = 0.000320954 Jy ! Clean residual min=-0.005132 max=0.005378 Jy/beam ! Clean residual mean=0.000006 rms=0.001099 Jy/beam ! Combined flux in latest and established models = 0.628061 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 287 components and 0.628061 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.193596Jy sigma=0.717385 ! Fit after self-cal, rms=0.193472Jy sigma=0.717075 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (0, -12). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00120972 Jy ! Component: 100 - total flux cleaned = 0.000299185 Jy ! Component: 150 - total flux cleaned = -3.47675e-05 Jy ! Component: 200 - total flux cleaned = -0.000193429 Jy ! Component: 250 - total flux cleaned = -0.000337169 Jy ! Component: 300 - total flux cleaned = -0.000337656 Jy ! Component: 350 - total flux cleaned = -0.000615285 Jy ! Component: 400 - total flux cleaned = -0.000750386 Jy ! Total flux subtracted in 400 components = -0.000750386 Jy ! Clean residual min=-0.004438 max=0.005245 Jy/beam ! Clean residual mean=0.000006 rms=0.001069 Jy/beam ! Combined flux in latest and established models = 0.62731 Jy ! Performing phase self-cal ! Adding 70 model components to the UV plane model. ! The established model now contains 316 components and 0.62731 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.192165Jy sigma=0.714572 ! Fit after self-cal, rms=0.192056Jy sigma=0.714397 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (3.6, -12.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00143155 Jy ! Component: 100 - total flux cleaned = 0.000969521 Jy ! Component: 150 - total flux cleaned = 0.000827081 Jy ! Component: 200 - total flux cleaned = 0.000825714 Jy ! Component: 250 - total flux cleaned = 0.000690041 Jy ! Component: 300 - total flux cleaned = 0.000421891 Jy ! Component: 350 - total flux cleaned = 3.0907e-05 Jy ! Component: 400 - total flux cleaned = -0.000225682 Jy ! Total flux subtracted in 400 components = -0.000225682 Jy ! Clean residual min=-0.004459 max=0.004933 Jy/beam ! Clean residual mean=0.000006 rms=0.001048 Jy/beam ! Combined flux in latest and established models = 0.627085 Jy ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 337 components and 0.627085 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.190940Jy sigma=0.712406 ! Fit after self-cal, rms=0.190869Jy sigma=0.712292 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (0.2, 5.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000455411 Jy ! Component: 100 - total flux cleaned = -0.000676447 Jy ! Component: 150 - total flux cleaned = -0.00140059 Jy ! Component: 200 - total flux cleaned = -0.00153549 Jy ! Component: 250 - total flux cleaned = -0.00205826 Jy ! Component: 300 - total flux cleaned = -0.00205636 Jy ! Component: 350 - total flux cleaned = -0.00205837 Jy ! Component: 400 - total flux cleaned = -0.00218358 Jy ! Total flux subtracted in 400 components = -0.00218358 Jy ! Clean residual min=-0.004304 max=0.004706 Jy/beam ! Clean residual mean=0.000006 rms=0.001030 Jy/beam ! Combined flux in latest and established models = 0.624901 Jy ! Performing phase self-cal ! Adding 81 model components to the UV plane model. ! The established model now contains 370 components and 0.624901 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.189752Jy sigma=0.710807 ! Fit after self-cal, rms=0.189689Jy sigma=0.710666 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (11.2, -20.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0018308 Jy ! Component: 100 - total flux cleaned = 0.00129535 Jy ! Component: 150 - total flux cleaned = 0.00142256 Jy ! Component: 200 - total flux cleaned = 0.00142267 Jy ! Component: 250 - total flux cleaned = 0.00130004 Jy ! Component: 300 - total flux cleaned = 0.00118036 Jy ! Component: 350 - total flux cleaned = 0.00142008 Jy ! Component: 400 - total flux cleaned = 0.00118188 Jy ! Total flux subtracted in 400 components = 0.00118188 Jy ! Clean residual min=-0.004138 max=0.004372 Jy/beam ! Clean residual mean=0.000005 rms=0.001008 Jy/beam ! Combined flux in latest and established models = 0.626083 Jy ! Performing phase self-cal ! Adding 98 model components to the UV plane model. ! The established model now contains 395 components and 0.626083 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.188787Jy sigma=0.709246 ! Fit after self-cal, rms=0.188757Jy sigma=0.709163 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-2.2, -2.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00144757 Jy ! Component: 100 - total flux cleaned = 0.00144166 Jy ! Component: 150 - total flux cleaned = 0.00131524 Jy ! Component: 200 - total flux cleaned = 0.0016842 Jy ! Component: 250 - total flux cleaned = 0.00156503 Jy ! Component: 300 - total flux cleaned = 0.00156561 Jy ! Component: 350 - total flux cleaned = 0.00144973 Jy ! Component: 400 - total flux cleaned = 0.00121745 Jy ! Total flux subtracted in 400 components = 0.00121745 Jy ! Clean residual min=-0.004109 max=0.004394 Jy/beam ! Clean residual mean=0.000005 rms=0.000993 Jy/beam ! Combined flux in latest and established models = 0.627301 Jy ! Performing phase self-cal ! Adding 89 model components to the UV plane model. ! The established model now contains 414 components and 0.6273 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.187776Jy sigma=0.707911 ! Fit after self-cal, rms=0.187715Jy sigma=0.707813 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (3.4, 4.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000648779 Jy ! Component: 100 - total flux cleaned = 0.000105498 Jy ! Component: 150 - total flux cleaned = -0.000144068 Jy ! Component: 200 - total flux cleaned = -0.000387283 Jy ! Component: 250 - total flux cleaned = -0.000504343 Jy ! Component: 300 - total flux cleaned = -0.000620381 Jy ! Component: 350 - total flux cleaned = -0.000506527 Jy ! Component: 400 - total flux cleaned = -0.000506967 Jy ! Total flux subtracted in 400 components = -0.000506967 Jy ! Clean residual min=-0.004097 max=0.004369 Jy/beam ! Clean residual mean=0.000005 rms=0.000980 Jy/beam ! Combined flux in latest and established models = 0.626793 Jy ! Performing phase self-cal ! Adding 86 model components to the UV plane model. ! The established model now contains 442 components and 0.626793 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.186807Jy sigma=0.706116 ! Fit after self-cal, rms=0.186790Jy sigma=0.706039 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-3.4, 5.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000426129 Jy ! Component: 100 - total flux cleaned = 0.000109844 Jy ! Component: 150 - total flux cleaned = -0.000154611 Jy ! Component: 200 - total flux cleaned = -0.000280227 Jy ! Component: 250 - total flux cleaned = -0.000280417 Jy ! Component: 300 - total flux cleaned = -0.00039805 Jy ! Component: 350 - total flux cleaned = -0.000395757 Jy ! Component: 400 - total flux cleaned = -0.000626002 Jy ! Total flux subtracted in 400 components = -0.000626002 Jy ! Clean residual min=-0.003773 max=0.003985 Jy/beam ! Clean residual mean=0.000005 rms=0.000957 Jy/beam ! Combined flux in latest and established models = 0.626167 Jy ! Performing phase self-cal ! Adding 86 model components to the UV plane model. ! The established model now contains 478 components and 0.626167 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.185799Jy sigma=0.704508 ! Fit after self-cal, rms=0.185774Jy sigma=0.704386 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-19.4, 33). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00080996 Jy ! Component: 100 - total flux cleaned = 0.000684981 Jy ! Component: 150 - total flux cleaned = 0.000310643 Jy ! Component: 200 - total flux cleaned = 0.000307467 Jy ! Component: 250 - total flux cleaned = 0.000188822 Jy ! Component: 300 - total flux cleaned = 0.000189522 Jy ! Component: 350 - total flux cleaned = 0.000301424 Jy ! Component: 400 - total flux cleaned = 0.00018972 Jy ! Total flux subtracted in 400 components = 0.00018972 Jy ! Clean residual min=-0.003787 max=0.003828 Jy/beam ! Clean residual mean=0.000005 rms=0.000936 Jy/beam ! Combined flux in latest and established models = 0.626357 Jy ! Performing phase self-cal ! Adding 91 model components to the UV plane model. ! The established model now contains 507 components and 0.626357 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.184874Jy sigma=0.703005 ! Fit after self-cal, rms=0.184864Jy sigma=0.702908 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-7, 12.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0011845 Jy ! Component: 100 - total flux cleaned = 0.000643008 Jy ! Component: 150 - total flux cleaned = 0.000271619 Jy ! Component: 200 - total flux cleaned = 0.000273432 Jy ! Component: 250 - total flux cleaned = 0.000389024 Jy ! Component: 300 - total flux cleaned = 0.000276203 Jy ! Component: 350 - total flux cleaned = 5.2742e-05 Jy ! Component: 400 - total flux cleaned = 5.43781e-05 Jy ! Total flux subtracted in 400 components = 5.43781e-05 Jy ! Clean residual min=-0.003698 max=0.003839 Jy/beam ! Clean residual mean=0.000005 rms=0.000915 Jy/beam ! Combined flux in latest and established models = 0.626411 Jy ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 537 components and 0.626411 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183937Jy sigma=0.701506 ! Fit after self-cal, rms=0.183908Jy sigma=0.701415 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (17, 46.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00167067 Jy ! Component: 100 - total flux cleaned = 0.00155338 Jy ! Component: 150 - total flux cleaned = 0.00133056 Jy ! Component: 200 - total flux cleaned = 0.00122092 Jy ! Component: 250 - total flux cleaned = 0.00121913 Jy ! Component: 300 - total flux cleaned = 0.000899137 Jy ! Component: 350 - total flux cleaned = 0.000792489 Jy ! Component: 400 - total flux cleaned = 0.000583127 Jy ! Total flux subtracted in 400 components = 0.000583127 Jy ! Clean residual min=-0.003808 max=0.003616 Jy/beam ! Clean residual mean=0.000004 rms=0.000900 Jy/beam ! Combined flux in latest and established models = 0.626994 Jy ! Performing phase self-cal ! Adding 108 model components to the UV plane model. ! The established model now contains 550 components and 0.626994 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.183152Jy sigma=0.700354 ! Fit after self-cal, rms=0.183127Jy sigma=0.700292 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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.183127Jy sigma=0.700292 ! Fit after self-cal, rms=0.152458Jy sigma=0.692127 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.7501 mas, bmaj=0.8549 mas, bpa=84.02 degrees ! Estimated noise=0.842973 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0013391 Jy ! Component: 100 - total flux cleaned = 0.00116013 Jy ! Component: 150 - total flux cleaned = 0.00100639 Jy ! Component: 200 - total flux cleaned = 0.000871722 Jy ! Component: 250 - total flux cleaned = 0.00100141 Jy ! Component: 300 - total flux cleaned = 0.00100672 Jy ! Component: 350 - total flux cleaned = 0.000897185 Jy ! Component: 400 - total flux cleaned = 0.000792975 Jy ! Total flux subtracted in 400 components = 0.000792975 Jy ! Clean residual min=-0.003014 max=0.003139 Jy/beam ! Clean residual mean=0.000002 rms=0.000703 Jy/beam ! Combined flux in latest and established models = 0.627787 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 563 components and 0.627788 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151665Jy sigma=0.689134 ! Fit after self-cal, rms=0.151447Jy sigma=0.688975 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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) ! Added new window around map position (7.8, 7.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00150391 Jy ! Component: 100 - total flux cleaned = 0.00213608 Jy ! Component: 150 - total flux cleaned = 0.00272935 Jy ! Component: 200 - total flux cleaned = 0.00301462 Jy ! Component: 250 - total flux cleaned = 0.00301658 Jy ! Component: 300 - total flux cleaned = 0.00319755 Jy ! Component: 350 - total flux cleaned = 0.00310894 Jy ! Component: 400 - total flux cleaned = 0.00336919 Jy ! Total flux subtracted in 400 components = 0.00336919 Jy ! Clean residual min=-0.002981 max=0.002869 Jy/beam ! Clean residual mean=0.000001 rms=0.000682 Jy/beam ! Combined flux in latest and established models = 0.631157 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 593 components and 0.631157 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151077Jy sigma=0.687734 ! Fit after self-cal, rms=0.151029Jy sigma=0.687667 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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.151029Jy sigma=0.687667 ! Fit after self-cal, rms=0.148942Jy sigma=0.674587 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=2.477 mas, bmaj=3.53 mas, bpa=-23.29 degrees ! Estimated noise=1.0736 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00125527 Jy ! Component: 100 - total flux cleaned = 0.00138742 Jy ! Component: 150 - total flux cleaned = 0.00144331 Jy ! Component: 200 - total flux cleaned = 0.00155225 Jy ! Component: 250 - total flux cleaned = 0.00171005 Jy ! Component: 300 - total flux cleaned = 0.00176195 Jy ! Component: 350 - total flux cleaned = 0.00181065 Jy ! Component: 400 - total flux cleaned = 0.00185566 Jy ! Total flux subtracted in 400 components = 0.00185566 Jy ! Clean residual min=-0.003023 max=0.002245 Jy/beam ! Clean residual mean=0.000006 rms=0.000639 Jy/beam ! Combined flux in latest and established models = 0.633012 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 630 components and 0.633013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148869Jy sigma=0.674266 ! Fit after self-cal, rms=0.148883Jy sigma=0.674177 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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 = -5.28852e-05 Jy ! Component: 100 - total flux cleaned = -0.000147293 Jy ! Component: 150 - total flux cleaned = -0.00010193 Jy ! Component: 200 - total flux cleaned = -5.92922e-05 Jy ! Component: 250 - total flux cleaned = -1.9163e-05 Jy ! Component: 300 - total flux cleaned = 2.09645e-05 Jy ! Component: 350 - total flux cleaned = 0.000100014 Jy ! Component: 400 - total flux cleaned = 0.000138219 Jy ! Total flux subtracted in 400 components = 0.000138219 Jy ! Clean residual min=-0.003012 max=0.002140 Jy/beam ! Clean residual mean=0.000005 rms=0.000616 Jy/beam ! Combined flux in latest and established models = 0.633151 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 644 components and 0.633151 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148833Jy sigma=0.674153 ! Fit after self-cal, rms=0.148818Jy sigma=0.674129 ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.7434 mas, bmaj=0.854 mas, bpa=85.31 degrees ! Estimated noise=0.842041 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7434 x 0.854 at 85.31 degrees (North through East) ! Clean map min=-0.015947 max=0.27365 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 ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.5153 mas, bmaj=0.6514 mas, bpa=85.61 degrees ! Estimated noise=1.51398 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.21432 Jy ! Component: 100 - total flux cleaned = 0.323336 Jy ! Component: 150 - total flux cleaned = 0.389802 Jy ! Component: 200 - total flux cleaned = 0.437161 Jy ! Component: 250 - total flux cleaned = 0.473507 Jy ! Component: 300 - total flux cleaned = 0.50245 Jy ! Component: 350 - total flux cleaned = 0.526009 Jy ! Component: 400 - total flux cleaned = 0.545665 Jy ! Total flux subtracted in 400 components = 0.545665 Jy ! Clean residual min=-0.010190 max=0.013204 Jy/beam ! Clean residual mean=0.000003 rms=0.001838 Jy/beam ! Combined flux in latest and established models = 0.545666 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 41 model components to the UV plane model. ! The established model now contains 41 components and 0.545666 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.7434 mas, bmaj=0.854 mas, bpa=85.31 degrees ! Estimated noise=0.842041 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.0371139 Jy ! Component: 100 - total flux cleaned = 0.0583412 Jy ! Component: 150 - total flux cleaned = 0.0691025 Jy ! Component: 200 - total flux cleaned = 0.0720189 Jy ! Component: 250 - total flux cleaned = 0.0752321 Jy ! Component: 300 - total flux cleaned = 0.07745 Jy ! Component: 350 - total flux cleaned = 0.0804085 Jy ! Component: 400 - total flux cleaned = 0.0826148 Jy ! Total flux subtracted in 400 components = 0.0826148 Jy ! Clean residual min=-0.004961 max=0.006741 Jy/beam ! Clean residual mean=0.000001 rms=0.001068 Jy/beam ! Combined flux in latest and established models = 0.62828 Jy ! Adding 64 model components to the UV plane model. ! The established model now contains 81 components and 0.628281 Jy ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (-1.6, 2.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00335584 Jy ! Component: 100 - total flux cleaned = 0.00485447 Jy ! Component: 150 - total flux cleaned = 0.00627972 Jy ! Component: 200 - total flux cleaned = 0.00786625 Jy ! Component: 250 - total flux cleaned = 0.00852205 Jy ! Component: 300 - total flux cleaned = 0.00894752 Jy ! Component: 350 - total flux cleaned = 0.00957002 Jy ! Component: 400 - total flux cleaned = 0.0105694 Jy ! Total flux subtracted in 400 components = 0.0105694 Jy ! Clean residual min=-0.004428 max=0.004498 Jy/beam ! Clean residual mean=0.000001 rms=0.000957 Jy/beam ! Combined flux in latest and established models = 0.63885 Jy ! Adding 71 model components to the UV plane model. ! The established model now contains 128 components and 0.63885 Jy ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (0.2, -6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000953478 Jy ! Component: 100 - total flux cleaned = 0.000765579 Jy ! Component: 150 - total flux cleaned = 0.00132933 Jy ! Component: 200 - total flux cleaned = 0.0016955 Jy ! Component: 250 - total flux cleaned = 0.00169862 Jy ! Component: 300 - total flux cleaned = 0.00222411 Jy ! Component: 350 - total flux cleaned = 0.00239819 Jy ! Component: 400 - total flux cleaned = 0.00256982 Jy ! Total flux subtracted in 400 components = 0.00256982 Jy ! Clean residual min=-0.004151 max=0.004225 Jy/beam ! Clean residual mean=0.000000 rms=0.000897 Jy/beam ! Combined flux in latest and established models = 0.64142 Jy ! Adding 97 model components to the UV plane model. ! The established model now contains 168 components and 0.641419 Jy ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Added new window around map position (9.6, -11.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00162699 Jy ! Component: 100 - total flux cleaned = 0.00162646 Jy ! Component: 150 - total flux cleaned = 0.00178577 Jy ! Component: 200 - total flux cleaned = 0.0019456 Jy ! Component: 250 - total flux cleaned = 0.00194611 Jy ! Component: 300 - total flux cleaned = 0.00194489 Jy ! Component: 350 - total flux cleaned = 0.00179476 Jy ! Component: 400 - total flux cleaned = 0.00164795 Jy ! Total flux subtracted in 400 components = 0.00164795 Jy ! Clean residual min=-0.003953 max=0.003826 Jy/beam ! Clean residual mean=-0.000000 rms=0.000851 Jy/beam ! Combined flux in latest and established models = 0.643067 Jy ! Adding 118 model components to the UV plane model. ! The established model now contains 201 components and 0.643068 Jy ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=2.477 mas, bmaj=3.53 mas, bpa=-23.29 degrees ! Estimated noise=1.0736 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000855788 Jy ! Component: 100 - total flux cleaned = 0.000858513 Jy ! Component: 150 - total flux cleaned = 0.000855083 Jy ! Component: 200 - total flux cleaned = 0.00095664 Jy ! Component: 250 - total flux cleaned = 0.000954993 Jy ! Component: 300 - total flux cleaned = 0.000952812 Jy ! Component: 350 - total flux cleaned = 0.00085967 Jy ! Component: 400 - total flux cleaned = 0.000952371 Jy ! Total flux subtracted in 400 components = 0.000952371 Jy ! Clean residual min=-0.002989 max=0.002738 Jy/beam ! Clean residual mean=-0.000003 rms=0.000768 Jy/beam ! Combined flux in latest and established models = 0.64402 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 = -7.70066e-07 Jy ! Component: 100 - total flux cleaned = 8.50614e-05 Jy ! Component: 150 - total flux cleaned = 0.000170523 Jy ! Component: 200 - total flux cleaned = 0.000167637 Jy ! Component: 250 - total flux cleaned = 8.4169e-05 Jy ! Component: 300 - total flux cleaned = 8.76015e-05 Jy ! Component: 350 - total flux cleaned = 0.000328536 Jy ! Component: 400 - total flux cleaned = 0.000485289 Jy ! Total flux subtracted in 400 components = 0.000485289 Jy ! Clean residual min=-0.002918 max=0.002682 Jy/beam ! Clean residual mean=-0.000003 rms=0.000728 Jy/beam ! Combined flux in latest and established models = 0.644505 Jy ! Adding 64 model components to the UV plane model. ! The established model now contains 257 components and 0.644505 Jy ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! Estimated beam: bmin=0.7434 mas, bmaj=0.854 mas, bpa=85.31 degrees ! Estimated noise=0.842041 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00111243 Jy ! Component: 100 - total flux cleaned = -0.000958509 Jy ! Component: 150 - total flux cleaned = -0.000960625 Jy ! Component: 200 - total flux cleaned = -0.000956577 Jy ! Component: 250 - total flux cleaned = -0.00082097 Jy ! Component: 300 - total flux cleaned = -0.000820576 Jy ! Component: 350 - total flux cleaned = -0.000687528 Jy ! Component: 400 - total flux cleaned = -0.000816727 Jy ! Total flux subtracted in 400 components = -0.000816727 Jy ! Clean residual min=-0.003371 max=0.003339 Jy/beam ! Clean residual mean=-0.000001 rms=0.000767 Jy/beam ! Combined flux in latest and established models = 0.643688 Jy selfcal ! Performing phase self-cal ! Adding 86 model components to the UV plane model. ! The established model now contains 294 components and 0.643689 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.152858Jy sigma=0.686739 ! Fit after self-cal, rms=0.152565Jy sigma=0.683863 wmodel /scr/pima/rdv73_uvs/J0607+6720_X_map.mod ! Writing 294 model components to file: /scr/pima/rdv73_uvs/J0607+6720_X_map.mod wobs /scr/pima/rdv73_uvs/J0607+6720_X_uvs.fits ! Applying 2081 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J0607+6720_X_uvs.fits wwins /scr/pima/rdv73_uvs/J0607+6720_X_map.win ! wwins: Wrote 23 windows to /scr/pima/rdv73_uvs/J0607+6720_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 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 6.47712e-05 Jy ! Component: 100 - total flux cleaned = 0.000417456 Jy ! Component: 150 - total flux cleaned = 0.000646384 Jy ! Component: 200 - total flux cleaned = 0.000815053 Jy ! Component: 250 - total flux cleaned = 0.00109064 Jy ! Component: 300 - total flux cleaned = 0.00125377 Jy ! Component: 350 - total flux cleaned = 0.00141418 Jy ! Component: 400 - total flux cleaned = 0.00162496 Jy ! Component: 450 - total flux cleaned = 0.00167597 Jy ! Component: 500 - total flux cleaned = 0.00167569 Jy ! Component: 550 - total flux cleaned = 0.00182924 Jy ! Component: 600 - total flux cleaned = 0.00193087 Jy ! Component: 650 - total flux cleaned = 0.00203143 Jy ! Component: 700 - total flux cleaned = 0.00208135 Jy ! Component: 750 - total flux cleaned = 0.00227926 Jy ! Component: 800 - total flux cleaned = 0.00232809 Jy ! Component: 850 - total flux cleaned = 0.00237678 Jy ! Component: 900 - total flux cleaned = 0.00252199 Jy ! Component: 950 - total flux cleaned = 0.00252183 Jy ! Component: 1000 - total flux cleaned = 0.00280715 Jy ! Component: 1050 - total flux cleaned = 0.00266513 Jy ! Component: 1100 - total flux cleaned = 0.00271177 Jy ! Component: 1150 - total flux cleaned = 0.00285219 Jy ! Component: 1200 - total flux cleaned = 0.00280572 Jy ! Component: 1250 - total flux cleaned = 0.00275928 Jy ! Component: 1300 - total flux cleaned = 0.00285089 Jy ! Component: 1350 - total flux cleaned = 0.00280492 Jy ! Component: 1400 - total flux cleaned = 0.00275951 Jy ! Component: 1450 - total flux cleaned = 0.00262384 Jy ! Component: 1500 - total flux cleaned = 0.00275858 Jy ! Component: 1550 - total flux cleaned = 0.00275825 Jy ! Component: 1600 - total flux cleaned = 0.00249123 Jy ! Component: 1650 - total flux cleaned = 0.00280146 Jy ! Component: 1700 - total flux cleaned = 0.00271302 Jy ! Component: 1750 - total flux cleaned = 0.00284459 Jy ! Component: 1800 - total flux cleaned = 0.00275701 Jy ! Component: 1850 - total flux cleaned = 0.00280072 Jy ! Component: 1900 - total flux cleaned = 0.00293081 Jy ! Component: 1950 - total flux cleaned = 0.00284468 Jy ! Component: 2000 - total flux cleaned = 0.00297372 Jy ! Component: 2050 - total flux cleaned = 0.00293079 Jy ! Component: 2100 - total flux cleaned = 0.00318702 Jy ! Component: 2150 - total flux cleaned = 0.00339976 Jy ! Component: 2200 - total flux cleaned = 0.00361143 Jy ! Component: 2250 - total flux cleaned = 0.00356908 Jy ! Component: 2300 - total flux cleaned = 0.00356904 Jy ! Component: 2350 - total flux cleaned = 0.00369463 Jy ! Component: 2400 - total flux cleaned = 0.00377817 Jy ! Component: 2450 - total flux cleaned = 0.0039031 Jy ! Component: 2500 - total flux cleaned = 0.00398586 Jy ! Component: 2550 - total flux cleaned = 0.00415128 Jy ! Component: 2600 - total flux cleaned = 0.00419244 Jy ! Component: 2650 - total flux cleaned = 0.00427445 Jy ! Component: 2700 - total flux cleaned = 0.00435601 Jy ! Component: 2750 - total flux cleaned = 0.00447852 Jy ! Component: 2800 - total flux cleaned = 0.00460054 Jy ! Component: 2850 - total flux cleaned = 0.00468148 Jy ! Component: 2900 - total flux cleaned = 0.00476218 Jy ! Component: 2950 - total flux cleaned = 0.00488308 Jy ! Component: 3000 - total flux cleaned = 0.00496314 Jy ! Component: 3050 - total flux cleaned = 0.00508307 Jy ! Component: 3100 - total flux cleaned = 0.0051628 Jy ! Component: 3150 - total flux cleaned = 0.00528223 Jy ! Component: 3200 - total flux cleaned = 0.00544083 Jy ! Component: 3250 - total flux cleaned = 0.00540122 Jy ! Component: 3300 - total flux cleaned = 0.00559851 Jy ! Component: 3350 - total flux cleaned = 0.0057951 Jy ! Component: 3400 - total flux cleaned = 0.0058343 Jy ! Component: 3450 - total flux cleaned = 0.00618642 Jy ! Component: 3500 - total flux cleaned = 0.00606921 Jy ! Component: 3550 - total flux cleaned = 0.00618598 Jy ! Component: 3600 - total flux cleaned = 0.00618581 Jy ! Component: 3650 - total flux cleaned = 0.0063018 Jy ! Component: 3700 - total flux cleaned = 0.00637883 Jy ! Component: 3750 - total flux cleaned = 0.00641725 Jy ! Component: 3800 - total flux cleaned = 0.0066091 Jy ! Component: 3850 - total flux cleaned = 0.00653258 Jy ! Component: 3900 - total flux cleaned = 0.00679982 Jy ! Component: 3950 - total flux cleaned = 0.00683782 Jy ! Component: 4000 - total flux cleaned = 0.0069517 Jy ! Component: 4050 - total flux cleaned = 0.00706544 Jy ! Total flux subtracted in 4096 components = 0.00725456 Jy ! Clean residual min=-0.002307 max=0.002403 Jy/beam ! Clean residual mean=0.000002 rms=0.000625 Jy/beam ! Combined flux in latest and established models = 0.650943 Jy keep ! Adding 356 model components to the UV plane model. ! The established model now contains 641 components and 0.650943 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J0607+6720_X\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7434 x 0.854 at 85.31 degrees (North through East) ! Clean map min=-0.0094884 max=0.27041 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J0607+6720_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map ! Your choice of large map pixels excluded 2.85% of the data. ! The x-axis pixel size should ideally be below 0.1702 milli-arcsec ! The y-axis pixel size should ideally be below 0.184 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= 0.2704 rms= 0.00062 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J0607+6720_X_uvs.log closed on Sat Jan 9 17:31:41 2016