! Started logfile: /scr/pima/rdv73_uvs/J0407-1211_S_uvs.log on Sat Jan 9 17:44:49 2016 obs /scr/pima/rdv73_uvs/J0407-1211_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J0407-1211_S_uva.fits ! AN table 1: 31 integrations on 45 of 45 possible baselines. ! AN table 2: 11 integrations on 36 of 36 possible baselines. ! AN table 3: 3 integrations on 1 of 1 possible baselines. ! AN table 4: 14 integrations on 6 of 6 possible baselines. ! AN table 5: 8 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.920997 visibilities/baseline/integration-bin. ! Found source: J0407-1211 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 1.6e+07 1 1.6e+07 ! 02 2 2.26299e+09 1.6e+07 1 1.6e+07 ! 03 3 2.35299e+09 1.6e+07 1 1.6e+07 ! 04 4 2.37299e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 6948 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv73_uvs/J0407-1211_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:WZ 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:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 38 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.157014Jy sigma=8.844617 ! Fit after self-cal, rms=0.496198Jy sigma=4.026632 ! 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:WZ ! Flagging baselines of antenna 1:WZ 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. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.144 mas, bmaj=9.266 mas, bpa=-1.203 degrees ! Estimated noise=2.13854 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.399734 Jy ! Component: 100 - total flux cleaned = 0.529601 Jy ! Total flux subtracted in 100 components = 0.529601 Jy ! Clean residual min=-0.038504 max=0.073995 Jy/beam ! Clean residual mean=0.000339 rms=0.010714 Jy/beam ! Combined flux in latest and established models = 0.529601 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 16 components and 0.529601 Jy ! ! Correcting IF 1. ! A total of 31 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 31 telescope corrections were flagged in sub-array 1. ! ! Correcting IF 3. ! A total of 31 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 31 telescope corrections were flagged in sub-array 1. ! A total of 24 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.211294Jy sigma=1.535461 ! Fit after self-cal, rms=0.210321Jy sigma=1.528992 ! Inverting map and beam ! Estimated beam: bmin=4.145 mas, bmaj=9.267 mas, bpa=-1.201 degrees ! Estimated noise=2.14603 mJy/beam. ! Added new window around map position (6, -5.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0721342 Jy ! Component: 100 - total flux cleaned = 0.114405 Jy ! Total flux subtracted in 100 components = 0.114405 Jy ! Clean residual min=-0.030793 max=0.037083 Jy/beam ! Clean residual mean=0.000200 rms=0.007622 Jy/beam ! Combined flux in latest and established models = 0.644007 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 34 components and 0.644007 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167880Jy sigma=1.219561 ! Fit after self-cal, rms=0.167537Jy sigma=1.214359 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.924 mas, bmaj=10.57 mas, bpa=-2.9 degrees ! Estimated noise=1.248 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.0504955 Jy ! Component: 100 - total flux cleaned = 0.0837331 Jy ! Component: 150 - total flux cleaned = 0.106335 Jy ! Component: 200 - total flux cleaned = 0.121789 Jy ! Total flux subtracted in 200 components = 0.121789 Jy ! Clean residual min=-0.025109 max=0.036884 Jy/beam ! Clean residual mean=0.000100 rms=0.007030 Jy/beam ! Combined flux in latest and established models = 0.765795 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 66 components and 0.765795 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146189Jy sigma=1.063764 ! Fit after self-cal, rms=0.145850Jy sigma=1.061371 ! 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 1.11 FD 0.90 HN 1.04 KP 0.98 ! LA 1.00 NL 0.99 OV 0.96 PT 1.05 ! SC 1.02 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 0.88 HN 1.05 KP 1.00 ! LA 1.02 NL 1.01 OV 0.97 PT 1.07 ! SC 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* OV 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.12 FD 0.95 HN 1.07 KP 1.02 ! LA 0.95 NL 1.02 OV 1.01 PT 1.05 ! SC 1.09 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.00 FD 0.96 HN 1.15 KP 1.04 ! LA 0.94 NL 1.01 OV 1.02 PT 1.06 ! SC 1.12 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* OV 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.15 FD 0.97 HN 1.16 KP 1.02 ! LA 0.87 NL 1.10 OV 1.05 PT 1.07 ! SC 0.99 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.08 FD 0.97 HN 1.17 KP 1.01 ! LA 0.89 NL 1.07 OV 1.05 PT 1.11 ! SC 1.01 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* OV 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.10 FD 0.94 HN 1.03 KP 0.92 ! LA 0.74 NL 1.01 OV 0.84 PT 0.92 ! SC 1.02 WZ 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.78 FD 0.94 HN 1.05 KP 0.92 ! LA 0.77 NL 0.99 OV 0.84 PT 0.95 ! SC 1.02 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* OV 0.00* ! ! A total of 18 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* OV 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! KK 0.00* MK 0.00* ! ! ! Fit before self-cal, rms=0.145850Jy sigma=1.061371 ! Fit after self-cal, rms=0.114455Jy sigma=0.842819 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.982 mas, bmaj=10.85 mas, bpa=-2.106 degrees ! Estimated noise=1.27527 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 (8, 39.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0326341 Jy ! Component: 100 - total flux cleaned = 0.0450561 Jy ! Component: 150 - total flux cleaned = 0.047077 Jy ! Component: 200 - total flux cleaned = 0.0486787 Jy ! Total flux subtracted in 200 components = 0.0486787 Jy ! Clean residual min=-0.010101 max=0.011386 Jy/beam ! Clean residual mean=-0.000028 rms=0.002548 Jy/beam ! Combined flux in latest and established models = 0.814474 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 99 components and 0.814474 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103530Jy sigma=0.770430 ! Fit after self-cal, rms=0.103184Jy sigma=0.767518 ! 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.103184Jy sigma=0.767518 ! Fit after self-cal, rms=0.101593Jy sigma=0.764273 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.967 mas, bmaj=10.82 mas, bpa=-1.967 degrees ! Estimated noise=1.26901 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00174804 Jy ! Component: 100 - total flux cleaned = 0.00150752 Jy ! Component: 150 - total flux cleaned = 0.00110385 Jy ! Component: 200 - total flux cleaned = 0.000171255 Jy ! Total flux subtracted in 200 components = 0.000171255 Jy ! Clean residual min=-0.009335 max=0.009392 Jy/beam ! Clean residual mean=-0.000012 rms=0.002187 Jy/beam ! Combined flux in latest and established models = 0.814645 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 125 components and 0.814645 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100351Jy sigma=0.757417 ! Fit after self-cal, rms=0.100336Jy sigma=0.757193 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.000347996 Jy ! Component: 100 - total flux cleaned = -0.000680993 Jy ! Component: 150 - total flux cleaned = -0.000982804 Jy ! Component: 200 - total flux cleaned = -0.00126838 Jy ! Total flux subtracted in 200 components = -0.00126838 Jy ! Clean residual min=-0.008992 max=0.009081 Jy/beam ! Clean residual mean=-0.000008 rms=0.002081 Jy/beam ! Combined flux in latest and established models = 0.813377 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 142 components and 0.813377 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099806Jy sigma=0.754337 ! Fit after self-cal, rms=0.099794Jy sigma=0.754266 ! Inverting map ! Added new window around map position (17.2, -60.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00417264 Jy ! Component: 100 - total flux cleaned = 0.00538697 Jy ! Component: 150 - total flux cleaned = 0.00598449 Jy ! Component: 200 - total flux cleaned = 0.00653057 Jy ! Total flux subtracted in 200 components = 0.00653057 Jy ! Clean residual min=-0.007781 max=0.006526 Jy/beam ! Clean residual mean=-0.000010 rms=0.001852 Jy/beam ! Combined flux in latest and established models = 0.819908 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 163 components and 0.819908 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098822Jy sigma=0.749095 ! Fit after self-cal, rms=0.098773Jy sigma=0.748917 ! 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.098773Jy sigma=0.748917 ! Fit after self-cal, rms=0.098392Jy sigma=0.748214 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.967 mas, bmaj=10.82 mas, bpa=-1.827 degrees ! Estimated noise=1.26783 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000992239 Jy ! Component: 100 - total flux cleaned = -0.000860235 Jy ! Component: 150 - total flux cleaned = -0.000859451 Jy ! Component: 200 - total flux cleaned = -0.0009735 Jy ! Total flux subtracted in 200 components = -0.0009735 Jy ! Clean residual min=-0.007359 max=0.006588 Jy/beam ! Clean residual mean=-0.000014 rms=0.001728 Jy/beam ! Combined flux in latest and established models = 0.818934 Jy selfcal ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 189 components and 0.818934 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098055Jy sigma=0.746117 ! Fit after self-cal, rms=0.098053Jy sigma=0.746093 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.000477236 Jy ! Component: 100 - total flux cleaned = 0.000477653 Jy ! Component: 150 - total flux cleaned = 0.000686036 Jy ! Component: 200 - total flux cleaned = 0.000787877 Jy ! Total flux subtracted in 200 components = 0.000787877 Jy ! Clean residual min=-0.007070 max=0.006536 Jy/beam ! Clean residual mean=-0.000014 rms=0.001673 Jy/beam ! Combined flux in latest and established models = 0.819722 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 202 components and 0.819722 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097808Jy sigma=0.744805 ! Fit after self-cal, rms=0.097808Jy sigma=0.744791 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097808Jy sigma=0.744791 ! Fit after self-cal, rms=0.097859Jy sigma=0.744616 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.088 mas, bmaj=14.48 mas, bpa=-14.63 degrees ! Estimated noise=1.45791 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00183662 Jy ! Component: 100 - total flux cleaned = -0.00161033 Jy ! Component: 150 - total flux cleaned = -0.00180994 Jy ! Component: 200 - total flux cleaned = -0.00180623 Jy ! Total flux subtracted in 200 components = -0.00180623 Jy ! Clean residual min=-0.006586 max=0.007107 Jy/beam ! Clean residual mean=-0.000025 rms=0.001828 Jy/beam ! Combined flux in latest and established models = 0.817916 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 211 components and 0.817916 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097682Jy sigma=0.743598 ! Fit after self-cal, rms=0.097658Jy sigma=0.743488 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.000376751 Jy ! Component: 100 - total flux cleaned = 0.000287783 Jy ! Component: 150 - total flux cleaned = 0.000286152 Jy ! Component: 200 - total flux cleaned = 0.000360251 Jy ! Total flux subtracted in 200 components = 0.000360251 Jy ! Clean residual min=-0.006304 max=0.007058 Jy/beam ! Clean residual mean=-0.000023 rms=0.001774 Jy/beam ! Combined flux in latest and established models = 0.818276 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 222 components and 0.818276 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097553Jy sigma=0.742922 ! Fit after self-cal, rms=0.097534Jy sigma=0.742851 ! 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: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:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2: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:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK 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 8 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 14 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 8 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.166060Jy sigma=1.110602 ! Fit after self-cal, rms=0.123433Jy sigma=0.927468 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:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 1:WZ 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:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:MK 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 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 24 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 42 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 24 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 24 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.123433Jy sigma=0.927468 ! Fit after self-cal, rms=0.109482Jy sigma=0.803005 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.538 mas, bmaj=9.683 mas, bpa=-4.332 degrees ! Estimated noise=1.9802 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.015093 Jy ! Component: 100 - total flux cleaned = -0.0156856 Jy ! Component: 150 - total flux cleaned = -0.0152817 Jy ! Component: 200 - total flux cleaned = -0.0144696 Jy ! Total flux subtracted in 200 components = -0.0144696 Jy ! Clean residual min=-0.013929 max=0.011086 Jy/beam ! Clean residual mean=0.000031 rms=0.002717 Jy/beam ! Combined flux in latest and established models = 0.803806 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 247 components and 0.803806 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.106433Jy sigma=0.786178 ! Fit after self-cal, rms=0.105669Jy sigma=0.781993 ! 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=2.844 mas, bmaj=11.13 mas, bpa=-4.876 degrees ! Estimated noise=1.17455 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 (3.2, 12.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00909618 Jy ! Component: 100 - total flux cleaned = 0.0101406 Jy ! Component: 150 - total flux cleaned = 0.00991585 Jy ! Component: 200 - total flux cleaned = 0.010128 Jy ! Component: 250 - total flux cleaned = 0.0103213 Jy ! Component: 300 - total flux cleaned = 0.0101359 Jy ! Component: 350 - total flux cleaned = 0.00979485 Jy ! Component: 400 - total flux cleaned = 0.0096319 Jy ! Total flux subtracted in 400 components = 0.0096319 Jy ! Clean residual min=-0.007083 max=0.005483 Jy/beam ! Clean residual mean=0.000003 rms=0.001417 Jy/beam ! Combined flux in latest and established models = 0.813438 Jy ! Performing phase self-cal ! Adding 82 model components to the UV plane model. ! The established model now contains 316 components and 0.813438 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.103462Jy sigma=0.768080 ! Fit after self-cal, rms=0.103365Jy sigma=0.767654 ! 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.000235887 Jy ! Component: 100 - total flux cleaned = 0.000237425 Jy ! Component: 150 - total flux cleaned = -4.37568e-05 Jy ! Component: 200 - total flux cleaned = -0.000308966 Jy ! Component: 250 - total flux cleaned = -0.000309033 Jy ! Component: 300 - total flux cleaned = -0.000556247 Jy ! Component: 350 - total flux cleaned = -0.00067516 Jy ! Component: 400 - total flux cleaned = -0.000787817 Jy ! Total flux subtracted in 400 components = -0.000787817 Jy ! Clean residual min=-0.006567 max=0.005050 Jy/beam ! Clean residual mean=-0.000000 rms=0.001277 Jy/beam ! Combined flux in latest and established models = 0.812651 Jy ! Performing phase self-cal ! Adding 86 model components to the UV plane model. ! The established model now contains 359 components and 0.81265 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102675Jy sigma=0.763395 ! Fit after self-cal, rms=0.102647Jy sigma=0.763289 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 56 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 4. ! A total of 16 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.126639Jy sigma=0.997647 ! Fit after self-cal, rms=0.102600Jy sigma=0.762244 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.839 mas, bmaj=11.15 mas, bpa=-4.87 degrees ! Estimated noise=1.17404 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00210955 Jy ! Component: 100 - total flux cleaned = -0.00246102 Jy ! Component: 150 - total flux cleaned = -0.00289579 Jy ! Component: 200 - total flux cleaned = -0.00310339 Jy ! Component: 250 - total flux cleaned = -0.00349631 Jy ! Component: 300 - total flux cleaned = -0.0034011 Jy ! Component: 350 - total flux cleaned = -0.00358399 Jy ! Component: 400 - total flux cleaned = -0.00375664 Jy ! Total flux subtracted in 400 components = -0.00375664 Jy ! Clean residual min=-0.006133 max=0.004691 Jy/beam ! Clean residual mean=-0.000007 rms=0.001131 Jy/beam ! Combined flux in latest and established models = 0.808894 Jy selfcal ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 401 components and 0.808894 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.102173Jy sigma=0.759499 ! Fit after self-cal, rms=0.102159Jy sigma=0.759434 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.000209302 Jy ! Component: 100 - total flux cleaned = -4.08323e-05 Jy ! Component: 150 - total flux cleaned = -0.000198885 Jy ! Component: 200 - total flux cleaned = -0.000274918 Jy ! Component: 250 - total flux cleaned = -0.000276161 Jy ! Component: 300 - total flux cleaned = -0.000421871 Jy ! Component: 350 - total flux cleaned = -0.000492195 Jy ! Component: 400 - total flux cleaned = -0.000560672 Jy ! Total flux subtracted in 400 components = -0.000560672 Jy ! Clean residual min=-0.005849 max=0.004579 Jy/beam ! Clean residual mean=-0.000008 rms=0.001084 Jy/beam ! Combined flux in latest and established models = 0.808333 Jy ! Performing phase self-cal ! Adding 78 model components to the UV plane model. ! The established model now contains 438 components and 0.808333 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101931Jy sigma=0.758086 ! Fit after self-cal, rms=0.101928Jy sigma=0.758058 ! 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.101928Jy sigma=0.758058 ! Fit after self-cal, rms=0.101960Jy sigma=0.757837 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.839 mas, bmaj=11.16 mas, bpa=-4.867 degrees ! Estimated noise=1.17432 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00132647 Jy ! Component: 100 - total flux cleaned = -0.00161417 Jy ! Component: 150 - total flux cleaned = -0.00168072 Jy ! Component: 200 - total flux cleaned = -0.00180925 Jy ! Component: 250 - total flux cleaned = -0.00180867 Jy ! Component: 300 - total flux cleaned = -0.00187014 Jy ! Component: 350 - total flux cleaned = -0.00187062 Jy ! Component: 400 - total flux cleaned = -0.00198567 Jy ! Total flux subtracted in 400 components = -0.00198567 Jy ! Clean residual min=-0.005452 max=0.004409 Jy/beam ! Clean residual mean=-0.000009 rms=0.001017 Jy/beam ! Combined flux in latest and established models = 0.806347 Jy selfcal ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 459 components and 0.806347 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101786Jy sigma=0.756759 ! Fit after self-cal, rms=0.101784Jy sigma=0.756735 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.00019284 Jy ! Component: 100 - total flux cleaned = 0.000247966 Jy ! Component: 150 - total flux cleaned = 0.000193604 Jy ! Component: 200 - total flux cleaned = 0.000192592 Jy ! Component: 250 - total flux cleaned = 0.000139908 Jy ! Component: 300 - total flux cleaned = 3.95816e-05 Jy ! Component: 350 - total flux cleaned = 4.01258e-05 Jy ! Component: 400 - total flux cleaned = 4.00463e-05 Jy ! Total flux subtracted in 400 components = 4.00463e-05 Jy ! Clean residual min=-0.005255 max=0.004358 Jy/beam ! Clean residual mean=-0.000009 rms=0.000992 Jy/beam ! Combined flux in latest and established models = 0.806388 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 479 components and 0.806388 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101669Jy sigma=0.756132 ! Fit after self-cal, rms=0.101669Jy sigma=0.756117 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101669Jy sigma=0.756117 ! Fit after self-cal, rms=0.100030Jy sigma=0.739849 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.061 mas, bmaj=14.48 mas, bpa=-15.28 degrees ! Estimated noise=1.42576 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000321518 Jy ! Component: 100 - total flux cleaned = -0.000320617 Jy ! Component: 150 - total flux cleaned = -0.000316696 Jy ! Component: 200 - total flux cleaned = -0.000358739 Jy ! Component: 250 - total flux cleaned = -0.00040058 Jy ! Component: 300 - total flux cleaned = -0.000479191 Jy ! Component: 350 - total flux cleaned = -0.000554625 Jy ! Component: 400 - total flux cleaned = -0.000590475 Jy ! Total flux subtracted in 400 components = -0.000590475 Jy ! Clean residual min=-0.003640 max=0.003879 Jy/beam ! Clean residual mean=-0.000029 rms=0.001078 Jy/beam ! Combined flux in latest and established models = 0.805797 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 503 components and 0.805797 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099995Jy sigma=0.739416 ! Fit after self-cal, rms=0.099993Jy sigma=0.739381 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-22, 98). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00141205 Jy ! Component: 100 - total flux cleaned = 0.000550632 Jy ! Component: 150 - total flux cleaned = 0.000397931 Jy ! Component: 200 - total flux cleaned = 0.000697173 Jy ! Component: 250 - total flux cleaned = 0.000749423 Jy ! Component: 300 - total flux cleaned = 0.000797757 Jy ! Component: 350 - total flux cleaned = 0.000703122 Jy ! Component: 400 - total flux cleaned = 0.000613351 Jy ! Total flux subtracted in 400 components = 0.000613351 Jy ! Clean residual min=-0.002975 max=0.003001 Jy/beam ! Clean residual mean=-0.000002 rms=0.000886 Jy/beam ! Combined flux in latest and established models = 0.80641 Jy ! Performing phase self-cal ! Adding 58 model components to the UV plane model. ! The established model now contains 558 components and 0.80641 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099780Jy sigma=0.738162 ! Fit after self-cal, rms=0.099774Jy sigma=0.738112 ! Inverting map ! Added new window around map position (79.6, 57.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 6.20622e-05 Jy ! Component: 100 - total flux cleaned = -3.95656e-05 Jy ! Component: 150 - total flux cleaned = -0.000121291 Jy ! Component: 200 - total flux cleaned = 1.196e-05 Jy ! Component: 250 - total flux cleaned = 0.000193077 Jy ! Component: 300 - total flux cleaned = 0.000360042 Jy ! Component: 350 - total flux cleaned = 0.000411102 Jy ! Component: 400 - total flux cleaned = 0.000411888 Jy ! Total flux subtracted in 400 components = 0.000411888 Jy ! Clean residual min=-0.002559 max=0.002230 Jy/beam ! Clean residual mean=-0.000012 rms=0.000643 Jy/beam ! Combined flux in latest and established models = 0.806822 Jy ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 633 components and 0.806822 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099583Jy sigma=0.736623 ! Fit after self-cal, rms=0.099569Jy sigma=0.736537 ! 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=2.841 mas, bmaj=11.17 mas, bpa=-4.867 degrees ! Estimated noise=1.17613 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.841 x 11.17 at -4.867 degrees (North through East) ! Clean map min=-0.0093457 max=0.47676 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=2.532 mas, bmaj=9.708 mas, bpa=-4.301 degrees ! Estimated noise=1.99309 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.385127 Jy ! Component: 100 - total flux cleaned = 0.502355 Jy ! Component: 150 - total flux cleaned = 0.568398 Jy ! Component: 200 - total flux cleaned = 0.612288 Jy ! Component: 250 - total flux cleaned = 0.644132 Jy ! Component: 300 - total flux cleaned = 0.66885 Jy ! Component: 350 - total flux cleaned = 0.689269 Jy ! Component: 400 - total flux cleaned = 0.706619 Jy ! Total flux subtracted in 400 components = 0.706619 Jy ! Clean residual min=-0.013493 max=0.014744 Jy/beam ! Clean residual mean=0.000074 rms=0.003117 Jy/beam ! Combined flux in latest and established models = 0.706619 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 111 model components to the UV plane model. ! The established model now contains 111 components and 0.706619 Jy ! Inverting map and beam ! Estimated beam: bmin=2.841 mas, bmaj=11.17 mas, bpa=-4.867 degrees ! Estimated noise=1.17613 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.0234955 Jy ! Component: 100 - total flux cleaned = 0.0394254 Jy ! Component: 150 - total flux cleaned = 0.0516012 Jy ! Component: 200 - total flux cleaned = 0.0612643 Jy ! Component: 250 - total flux cleaned = 0.0679607 Jy ! Component: 300 - total flux cleaned = 0.072083 Jy ! Component: 350 - total flux cleaned = 0.075021 Jy ! Component: 400 - total flux cleaned = 0.0774826 Jy ! Total flux subtracted in 400 components = 0.0774826 Jy ! Clean residual min=-0.006229 max=0.006116 Jy/beam ! Clean residual mean=0.000032 rms=0.001569 Jy/beam ! Combined flux in latest and established models = 0.784101 Jy ! Adding 113 model components to the UV plane model. ! The established model now contains 215 components and 0.784101 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.061 mas, bmaj=14.48 mas, bpa=-15.28 degrees ! Estimated noise=1.42576 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0074282 Jy ! Component: 100 - total flux cleaned = 0.0101363 Jy ! Component: 150 - total flux cleaned = 0.0127875 Jy ! Component: 200 - total flux cleaned = 0.0143866 Jy ! Component: 250 - total flux cleaned = 0.0161469 Jy ! Component: 300 - total flux cleaned = 0.017768 Jy ! Component: 350 - total flux cleaned = 0.018776 Jy ! Component: 400 - total flux cleaned = 0.019607 Jy ! Total flux subtracted in 400 components = 0.019607 Jy ! Clean residual min=-0.004017 max=0.004013 Jy/beam ! Clean residual mean=-0.000003 rms=0.001182 Jy/beam ! Combined flux in latest and established models = 0.803708 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.000450137 Jy ! Component: 100 - total flux cleaned = 0.00056079 Jy ! Component: 150 - total flux cleaned = 0.000665733 Jy ! Component: 200 - total flux cleaned = 0.000961088 Jy ! Component: 250 - total flux cleaned = 0.000959379 Jy ! Component: 300 - total flux cleaned = 0.00113627 Jy ! Component: 350 - total flux cleaned = 0.00113667 Jy ! Component: 400 - total flux cleaned = 0.00122284 Jy ! Total flux subtracted in 400 components = 0.00122284 Jy ! Clean residual min=-0.003378 max=0.003334 Jy/beam ! Clean residual mean=-0.000004 rms=0.000930 Jy/beam ! Combined flux in latest and established models = 0.804931 Jy ! Adding 156 model components to the UV plane model. ! The established model now contains 363 components and 0.804931 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.841 mas, bmaj=11.17 mas, bpa=-4.867 degrees ! Estimated noise=1.17613 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00108359 Jy ! Component: 100 - total flux cleaned = 0.00215732 Jy ! Component: 150 - total flux cleaned = 0.0031498 Jy ! Component: 200 - total flux cleaned = 0.00396296 Jy ! Component: 250 - total flux cleaned = 0.00484146 Jy ! Component: 300 - total flux cleaned = 0.00567992 Jy ! Component: 350 - total flux cleaned = 0.00638325 Jy ! Component: 400 - total flux cleaned = 0.00677236 Jy ! Total flux subtracted in 400 components = 0.00677236 Jy ! Clean residual min=-0.004936 max=0.003391 Jy/beam ! Clean residual mean=0.000000 rms=0.000886 Jy/beam ! Combined flux in latest and established models = 0.811703 Jy selfcal ! Performing phase self-cal ! Adding 82 model components to the UV plane model. ! The established model now contains 442 components and 0.811704 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100044Jy sigma=0.738584 ! Fit after self-cal, rms=0.099941Jy sigma=0.737975 wmodel /scr/pima/rdv73_uvs/J0407-1211_S_map.mod ! Writing 442 model components to file: /scr/pima/rdv73_uvs/J0407-1211_S_map.mod wobs /scr/pima/rdv73_uvs/J0407-1211_S_uvs.fits ! Applying 1265 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J0407-1211_S_uvs.fits wwins /scr/pima/rdv73_uvs/J0407-1211_S_map.win ! wwins: Wrote 7 windows to /scr/pima/rdv73_uvs/J0407-1211_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00153998 Jy ! Component: 100 - total flux cleaned = -0.00244247 Jy ! Component: 150 - total flux cleaned = -0.00244568 Jy ! Component: 200 - total flux cleaned = -0.00233589 Jy ! Component: 250 - total flux cleaned = -0.00223014 Jy ! Component: 300 - total flux cleaned = -0.00207431 Jy ! Component: 350 - total flux cleaned = -0.00187575 Jy ! Component: 400 - total flux cleaned = -0.00168427 Jy ! Component: 450 - total flux cleaned = -0.00145146 Jy ! Component: 500 - total flux cleaned = -0.00126981 Jy ! Component: 550 - total flux cleaned = -0.00104878 Jy ! Component: 600 - total flux cleaned = -0.00083353 Jy ! Component: 650 - total flux cleaned = -0.000707071 Jy ! Component: 700 - total flux cleaned = -0.000499305 Jy ! Component: 750 - total flux cleaned = -0.000254949 Jy ! Component: 800 - total flux cleaned = -0.000255544 Jy ! Component: 850 - total flux cleaned = -0.000175991 Jy ! Component: 900 - total flux cleaned = -9.83966e-05 Jy ! Component: 950 - total flux cleaned = -2.29995e-05 Jy ! Component: 1000 - total flux cleaned = -6.05753e-05 Jy ! Component: 1050 - total flux cleaned = 1.30207e-05 Jy ! Component: 1100 - total flux cleaned = -2.33495e-05 Jy ! Component: 1150 - total flux cleaned = 0.000119383 Jy ! Component: 1200 - total flux cleaned = 0.000118894 Jy ! Component: 1250 - total flux cleaned = 0.000188337 Jy ! Component: 1300 - total flux cleaned = 0.000290705 Jy ! Component: 1350 - total flux cleaned = 0.000358151 Jy ! Component: 1400 - total flux cleaned = 0.000324814 Jy ! Component: 1450 - total flux cleaned = 0.000424045 Jy ! Component: 1500 - total flux cleaned = 0.00039143 Jy ! Component: 1550 - total flux cleaned = 0.000488622 Jy ! Component: 1600 - total flux cleaned = 0.000393304 Jy ! Component: 1650 - total flux cleaned = 0.000520164 Jy ! Component: 1700 - total flux cleaned = 0.000551337 Jy ! Component: 1750 - total flux cleaned = 0.000520939 Jy ! Component: 1800 - total flux cleaned = 0.000521135 Jy ! Component: 1850 - total flux cleaned = 0.000460861 Jy ! Component: 1900 - total flux cleaned = 0.00052056 Jy ! Component: 1950 - total flux cleaned = 0.000490846 Jy ! Component: 2000 - total flux cleaned = 0.000461924 Jy ! Component: 2050 - total flux cleaned = 0.000462175 Jy ! Component: 2100 - total flux cleaned = 0.000405337 Jy ! Component: 2150 - total flux cleaned = 0.000490733 Jy ! Component: 2200 - total flux cleaned = 0.000462487 Jy ! Component: 2250 - total flux cleaned = 0.00051794 Jy ! Component: 2300 - total flux cleaned = 0.000490179 Jy ! Component: 2350 - total flux cleaned = 0.000408694 Jy ! Component: 2400 - total flux cleaned = 0.00043554 Jy ! Component: 2450 - total flux cleaned = 0.000382078 Jy ! Component: 2500 - total flux cleaned = 0.000329243 Jy ! Component: 2550 - total flux cleaned = 0.000328655 Jy ! Component: 2600 - total flux cleaned = 0.000276195 Jy ! Component: 2650 - total flux cleaned = 9.49254e-05 Jy ! Component: 2700 - total flux cleaned = 9.51175e-05 Jy ! Component: 2750 - total flux cleaned = 9.4963e-05 Jy ! Component: 2800 - total flux cleaned = -3.14967e-05 Jy ! Component: 2850 - total flux cleaned = -8.1495e-05 Jy ! Component: 2900 - total flux cleaned = -0.000131205 Jy ! Component: 2950 - total flux cleaned = -0.000303924 Jy ! Component: 3000 - total flux cleaned = -0.000328623 Jy ! Component: 3050 - total flux cleaned = -0.000450116 Jy ! Component: 3100 - total flux cleaned = -0.000425692 Jy ! Component: 3150 - total flux cleaned = -0.000545342 Jy ! Component: 3200 - total flux cleaned = -0.000592647 Jy ! Component: 3250 - total flux cleaned = -0.000616159 Jy ! Component: 3300 - total flux cleaned = -0.000639468 Jy ! Component: 3350 - total flux cleaned = -0.000662637 Jy ! Component: 3400 - total flux cleaned = -0.000731592 Jy ! Component: 3450 - total flux cleaned = -0.000823112 Jy ! Component: 3500 - total flux cleaned = -0.000754861 Jy ! Component: 3550 - total flux cleaned = -0.000709608 Jy ! Component: 3600 - total flux cleaned = -0.000754478 Jy ! Component: 3650 - total flux cleaned = -0.000754178 Jy ! Component: 3700 - total flux cleaned = -0.000842856 Jy ! Component: 3750 - total flux cleaned = -0.000798523 Jy ! Component: 3800 - total flux cleaned = -0.000732936 Jy ! Component: 3850 - total flux cleaned = -0.000754899 Jy ! Component: 3900 - total flux cleaned = -0.000754756 Jy ! Component: 3950 - total flux cleaned = -0.000711575 Jy ! Component: 4000 - total flux cleaned = -0.000711503 Jy ! Component: 4050 - total flux cleaned = -0.000647504 Jy ! Total flux subtracted in 4096 components = -0.000605247 Jy ! Clean residual min=-0.001407 max=0.001463 Jy/beam ! Clean residual mean=0.000002 rms=0.000404 Jy/beam ! Combined flux in latest and established models = 0.811098 Jy keep ! Adding 251 model components to the UV plane model. ! The established model now contains 689 components and 0.811099 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J0407-1211_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.841 x 11.17 at -4.867 degrees (North through East) ! Clean map min=-0.006037 max=0.49307 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J0407-1211_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.493 rms= 0.00039 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J0407-1211_S_uvs.log closed on Sat Jan 9 17:44:56 2016