! Started logfile: /scr/pima/rdv73_uvs/J1037-2934_S_uvs.log on Sat Jan 9 17:48:20 2016 obs /scr/pima/rdv73_uvs/J1037-2934_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J1037-2934_S_uva.fits ! AN table 1: 3 integrations on 15 of 15 possible baselines. ! AN table 2: 7 integrations on 21 of 21 possible baselines. ! AN table 3: 3 integrations on 36 of 36 possible baselines. ! AN table 4: 3 integrations on 3 of 3 possible baselines. ! AN table 5: 3 integrations on 28 of 28 possible baselines. ! AN table 6: 2 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.929114 visibilities/baseline/integration-bin. ! Found source: J1037-2934 ! ! 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 1468 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/J1037-2934_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:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC 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: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:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC 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 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6: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 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.902917Jy sigma=5.177825 ! Fit after self-cal, rms=0.757794Jy sigma=1.854134 ! 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 4:KK ! Flagging baselines of antenna 4:KK in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6: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=2.682 mas, bmaj=8.523 mas, bpa=-9.088 degrees ! Estimated noise=9.5565 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.98961 Jy ! Component: 100 - total flux cleaned = 1.31146 Jy ! Total flux subtracted in 100 components = 1.31146 Jy ! Clean residual min=-0.070882 max=0.177479 Jy/beam ! Clean residual mean=0.001416 rms=0.023366 Jy/beam ! Combined flux in latest and established models = 1.31146 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 12 components and 1.31146 Jy ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.580674Jy sigma=1.368724 ! Fit after self-cal, rms=0.581526Jy sigma=1.367482 ! Inverting map and beam ! Estimated beam: bmin=2.681 mas, bmaj=8.505 mas, bpa=-8.761 degrees ! Estimated noise=9.66872 mJy/beam. ! Added new window around map position (-3.2, -6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.20431 Jy ! Component: 100 - total flux cleaned = 0.33983 Jy ! Total flux subtracted in 100 components = 0.33983 Jy ! Clean residual min=-0.055437 max=0.092512 Jy/beam ! Clean residual mean=0.000856 rms=0.015755 Jy/beam ! Combined flux in latest and established models = 1.65129 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 30 components and 1.65129 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.415398Jy sigma=0.937077 ! Fit after self-cal, rms=0.415028Jy sigma=0.936374 ! 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=3.377 mas, bmaj=10.6 mas, bpa=-6.245 degrees ! Estimated noise=7.91593 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, 9.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.171244 Jy ! Component: 100 - total flux cleaned = 0.272591 Jy ! Component: 150 - total flux cleaned = 0.33453 Jy ! Component: 200 - total flux cleaned = 0.365086 Jy ! Total flux subtracted in 200 components = 0.365086 Jy ! Clean residual min=-0.052360 max=0.033274 Jy/beam ! Clean residual mean=-0.000001 rms=0.011770 Jy/beam ! Combined flux in latest and established models = 2.01638 Jy ! Performing phase self-cal ! Adding 65 model components to the UV plane model. ! The established model now contains 91 components and 2.01638 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.310923Jy sigma=0.705597 ! Fit after self-cal, rms=0.310687Jy sigma=0.704996 ! 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: ! FD 0.85 HN 1.17 LA 1.01 NL 1.09 ! PT 1.08 SC 0.89 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.91 HN 1.04 KP 1.01 LA 0.96 ! OV 1.03 PT 1.08 SC 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.06 FD 0.90 HN 1.21 KP 1.02 ! LA 0.95 MK 0.93 OV 0.97 PT 1.08 ! SC 1.09 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.94 FD 0.92 KK 1.00* KP 1.02 ! LA 1.01 MK 0.77 OV 1.01 PT 1.08 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 0.90 HN 1.27 LA 0.94 NL 1.17 ! PT 1.05 SC 0.95 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.98 HN 1.07 KP 1.05 LA 0.93 ! OV 1.09 PT 1.09 SC 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.08 FD 0.99 HN 1.20 KP 1.03 ! LA 0.90 MK 1.03 OV 1.04 PT 1.10 ! SC 1.12 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.95 FD 1.01 KK 1.00* KP 1.04 ! LA 0.95 MK 0.82 OV 1.04 PT 1.08 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 0.94 HN 1.30 LA 0.87 NL 1.23 ! PT 1.11 SC 0.82 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.98 HN 1.12 KP 1.06 LA 0.84 ! OV 1.15 PT 1.09 SC 0.95 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.19 FD 0.98 HN 1.33 KP 1.02 ! LA 0.85 MK 1.01 OV 1.09 PT 1.11 ! SC 1.01 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.11 FD 0.99 KK 1.00* KP 1.07 ! LA 0.88 MK 0.83 OV 1.06 PT 1.12 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 0.92 HN 1.17 LA 0.75 NL 1.13 ! PT 0.94 SC 0.89 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.95 HN 1.00 KP 0.94 LA 0.70 ! OV 0.92 PT 0.97 SC 0.95 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 0.95 HN 1.14 KP 0.89 ! LA 0.73 MK 1.10 OV 0.86 PT 0.95 ! SC 0.96 ! ! Telescope amplitude corrections in sub-array 4: ! BR 0.00* KK 0.00* MK 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 1.01 KK 1.00* KP 0.94 ! LA 0.74 MK 0.92 OV 0.87 PT 0.97 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Fit before self-cal, rms=0.310687Jy sigma=0.704996 ! Fit after self-cal, rms=0.151550Jy sigma=0.352676 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.394 mas, bmaj=10.77 mas, bpa=-5.915 degrees ! Estimated noise=7.89742 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.00316351 Jy ! Component: 100 - total flux cleaned = -0.00465932 Jy ! Component: 150 - total flux cleaned = -0.00597141 Jy ! Component: 200 - total flux cleaned = -0.00676459 Jy ! Total flux subtracted in 200 components = -0.00676459 Jy ! Clean residual min=-0.018871 max=0.017843 Jy/beam ! Clean residual mean=-0.000289 rms=0.004229 Jy/beam ! Combined flux in latest and established models = 2.00961 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 111 components and 2.00961 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147851Jy sigma=0.345145 ! Fit after self-cal, rms=0.147885Jy sigma=0.344703 ! 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.147885Jy sigma=0.344703 ! Fit after self-cal, rms=0.147327Jy sigma=0.343338 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.39 mas, bmaj=10.76 mas, bpa=-5.749 degrees ! Estimated noise=7.89108 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00452211 Jy ! Component: 100 - total flux cleaned = -0.00538342 Jy ! Component: 150 - total flux cleaned = -0.00611449 Jy ! Component: 200 - total flux cleaned = -0.00709994 Jy ! Total flux subtracted in 200 components = -0.00709994 Jy ! Clean residual min=-0.017546 max=0.016583 Jy/beam ! Clean residual mean=-0.000277 rms=0.004001 Jy/beam ! Combined flux in latest and established models = 2.00251 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 129 components and 2.00251 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.144558Jy sigma=0.338126 ! Fit after self-cal, rms=0.144511Jy sigma=0.337905 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.000691232 Jy ! Component: 100 - total flux cleaned = -0.000981078 Jy ! Component: 150 - total flux cleaned = -0.00177784 Jy ! Component: 200 - total flux cleaned = -0.00251732 Jy ! Total flux subtracted in 200 components = -0.00251732 Jy ! Clean residual min=-0.016704 max=0.016156 Jy/beam ! Clean residual mean=-0.000270 rms=0.003905 Jy/beam ! Combined flux in latest and established models = 2 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 148 components and 2 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143256Jy sigma=0.335486 ! Fit after self-cal, rms=0.143231Jy sigma=0.335345 ! 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.143231Jy sigma=0.335345 ! Fit after self-cal, rms=0.141457Jy sigma=0.330143 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.384 mas, bmaj=10.76 mas, bpa=-5.584 degrees ! Estimated noise=7.88184 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0054593 Jy ! Component: 100 - total flux cleaned = -0.006157 Jy ! Component: 150 - total flux cleaned = -0.00674076 Jy ! Component: 200 - total flux cleaned = -0.00750232 Jy ! Total flux subtracted in 200 components = -0.00750232 Jy ! Clean residual min=-0.015565 max=0.016065 Jy/beam ! Clean residual mean=-0.000271 rms=0.003592 Jy/beam ! Combined flux in latest and established models = 1.99249 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 158 components and 1.99249 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.139387Jy sigma=0.326813 ! Fit after self-cal, rms=0.139364Jy sigma=0.326709 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 (19.2, -20.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0100676 Jy ! Component: 100 - total flux cleaned = 0.0113224 Jy ! Component: 150 - total flux cleaned = 0.0120187 Jy ! Component: 200 - total flux cleaned = 0.0117384 Jy ! Total flux subtracted in 200 components = 0.0117384 Jy ! Clean residual min=-0.014071 max=0.010931 Jy/beam ! Clean residual mean=-0.000303 rms=0.003273 Jy/beam ! Combined flux in latest and established models = 2.00423 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 180 components and 2.00423 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.135874Jy sigma=0.321447 ! Fit after self-cal, rms=0.135734Jy sigma=0.321090 ! 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.135734Jy sigma=0.321090 ! Fit after self-cal, rms=0.135536Jy sigma=0.319816 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.061 mas, bmaj=18.07 mas, bpa=-2.041 degrees ! Estimated noise=10.5376 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000475967 Jy ! Component: 100 - total flux cleaned = -0.00119876 Jy ! Component: 150 - total flux cleaned = -0.00264221 Jy ! Component: 200 - total flux cleaned = -0.00368104 Jy ! Total flux subtracted in 200 components = -0.00368104 Jy ! Clean residual min=-0.015336 max=0.012093 Jy/beam ! Clean residual mean=-0.000497 rms=0.003973 Jy/beam ! Combined flux in latest and established models = 2.00055 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 200 components and 2.00055 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133639Jy sigma=0.319003 ! Fit after self-cal, rms=0.133602Jy sigma=0.318709 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.000804309 Jy ! Component: 100 - total flux cleaned = 0.00034984 Jy ! Component: 150 - total flux cleaned = -7.01018e-05 Jy ! Component: 200 - total flux cleaned = -0.000474122 Jy ! Total flux subtracted in 200 components = -0.000474122 Jy ! Clean residual min=-0.014663 max=0.012054 Jy/beam ! Clean residual mean=-0.000493 rms=0.003904 Jy/beam ! Combined flux in latest and established models = 2.00008 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 211 components and 2.00008 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132590Jy sigma=0.318987 ! Fit after self-cal, rms=0.132586Jy sigma=0.318866 ! 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:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC 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: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:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC 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 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:MK,true ! 5:MK self-cal status=fixed weight=1 selfant 5:OV,true ! 5:OV self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:MK,true ! 6: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 3 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 1 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.397940Jy sigma=0.895195 ! Fit after self-cal, rms=0.381523Jy sigma=0.839560 selfant "",false ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC 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: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:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC 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 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:KK self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:MK self-cal status=correctable weight=1 ! 5:OV self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6: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 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 6 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! A total of 2 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.381523Jy sigma=0.839560 ! Fit after self-cal, rms=0.136715Jy sigma=0.337394 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.649 mas, bmaj=8.743 mas, bpa=-7.615 degrees ! Estimated noise=9.38017 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.00975062 Jy ! Component: 100 - total flux cleaned = -0.00933089 Jy ! Component: 150 - total flux cleaned = -0.00846203 Jy ! Component: 200 - total flux cleaned = -0.00804743 Jy ! Total flux subtracted in 200 components = -0.00804743 Jy ! Clean residual min=-0.014527 max=0.012427 Jy/beam ! Clean residual mean=-0.000069 rms=0.003416 Jy/beam ! Combined flux in latest and established models = 1.99203 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 240 components and 1.99203 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133640Jy sigma=0.329661 ! Fit after self-cal, rms=0.133497Jy sigma=0.329228 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.32 mas, bmaj=10.9 mas, bpa=-5.662 degrees ! Estimated noise=7.683 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.000826691 Jy ! Component: 100 - total flux cleaned = -0.00014084 Jy ! Component: 150 - total flux cleaned = -0.00013529 Jy ! Component: 200 - total flux cleaned = -0.000503864 Jy ! Component: 250 - total flux cleaned = -0.00102199 Jy ! Component: 300 - total flux cleaned = -0.00134271 Jy ! Component: 350 - total flux cleaned = -0.00164392 Jy ! Component: 400 - total flux cleaned = -0.00193343 Jy ! Total flux subtracted in 400 components = -0.00193343 Jy ! Clean residual min=-0.013642 max=0.010857 Jy/beam ! Clean residual mean=-0.000217 rms=0.003067 Jy/beam ! Combined flux in latest and established models = 1.9901 Jy ! Performing phase self-cal ! Adding 40 model components to the UV plane model. ! The established model now contains 272 components and 1.9901 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132355Jy sigma=0.326369 ! Fit after self-cal, rms=0.132332Jy sigma=0.326267 ! 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.00074992 Jy ! Component: 100 - total flux cleaned = 0.000751199 Jy ! Component: 150 - total flux cleaned = 0.000487403 Jy ! Component: 200 - total flux cleaned = 0.000240955 Jy ! Component: 250 - total flux cleaned = 0.000119582 Jy ! Component: 300 - total flux cleaned = -0.000226067 Jy ! Component: 350 - total flux cleaned = -0.000558368 Jy ! Component: 400 - total flux cleaned = -0.000666349 Jy ! Total flux subtracted in 400 components = -0.000666349 Jy ! Clean residual min=-0.012706 max=0.010565 Jy/beam ! Clean residual mean=-0.000216 rms=0.003022 Jy/beam ! Combined flux in latest and established models = 1.98943 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 297 components and 1.98943 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.131773Jy sigma=0.325264 ! Fit after self-cal, rms=0.131784Jy sigma=0.325217 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 1 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 8 telescope corrections were flagged in sub-array 5. ! A total of 4 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.393925Jy sigma=0.887323 ! Fit after self-cal, rms=0.130884Jy sigma=0.322000 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.323 mas, bmaj=10.89 mas, bpa=-5.534 degrees ! Estimated noise=7.69098 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0104571 Jy ! Component: 100 - total flux cleaned = -0.0104951 Jy ! Component: 150 - total flux cleaned = -0.010694 Jy ! Component: 200 - total flux cleaned = -0.0100534 Jy ! Component: 250 - total flux cleaned = -0.00977024 Jy ! Component: 300 - total flux cleaned = -0.00989938 Jy ! Component: 350 - total flux cleaned = -0.0101388 Jy ! Component: 400 - total flux cleaned = -0.0101341 Jy ! Total flux subtracted in 400 components = -0.0101341 Jy ! Clean residual min=-0.011108 max=0.009513 Jy/beam ! Clean residual mean=-0.000134 rms=0.002801 Jy/beam ! Combined flux in latest and established models = 1.9793 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 326 components and 1.9793 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129446Jy sigma=0.318316 ! Fit after self-cal, rms=0.129452Jy sigma=0.318282 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.000388675 Jy ! Component: 100 - total flux cleaned = 0.000605587 Jy ! Component: 150 - total flux cleaned = 0.000505388 Jy ! Component: 200 - total flux cleaned = 0.000600479 Jy ! Component: 250 - total flux cleaned = 0.000601103 Jy ! Component: 300 - total flux cleaned = 0.000600161 Jy ! Component: 350 - total flux cleaned = 0.000600119 Jy ! Component: 400 - total flux cleaned = 0.000514446 Jy ! Total flux subtracted in 400 components = 0.000514446 Jy ! Clean residual min=-0.010556 max=0.009434 Jy/beam ! Clean residual mean=-0.000136 rms=0.002774 Jy/beam ! Combined flux in latest and established models = 1.97981 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 340 components and 1.97981 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129128Jy sigma=0.317691 ! Fit after self-cal, rms=0.129146Jy sigma=0.317668 ! 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.129146Jy sigma=0.317668 ! Fit after self-cal, rms=0.128506Jy sigma=0.315479 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.324 mas, bmaj=10.88 mas, bpa=-5.426 degrees ! Estimated noise=7.69553 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00839161 Jy ! Component: 100 - total flux cleaned = -0.00840273 Jy ! Component: 150 - total flux cleaned = -0.00825764 Jy ! Component: 200 - total flux cleaned = -0.00812757 Jy ! Component: 250 - total flux cleaned = -0.00776835 Jy ! Component: 300 - total flux cleaned = -0.00755573 Jy ! Component: 350 - total flux cleaned = -0.00745625 Jy ! Component: 400 - total flux cleaned = -0.00745655 Jy ! Total flux subtracted in 400 components = -0.00745655 Jy ! Clean residual min=-0.009911 max=0.008957 Jy/beam ! Clean residual mean=-0.000104 rms=0.002620 Jy/beam ! Combined flux in latest and established models = 1.97235 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 351 components and 1.97235 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127383Jy sigma=0.313093 ! Fit after self-cal, rms=0.127386Jy sigma=0.313077 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.000406124 Jy ! Component: 100 - total flux cleaned = 0.000410012 Jy ! Component: 150 - total flux cleaned = 0.000327557 Jy ! Component: 200 - total flux cleaned = 0.000328122 Jy ! Component: 250 - total flux cleaned = 0.000249459 Jy ! Component: 300 - total flux cleaned = 0.000249309 Jy ! Component: 350 - total flux cleaned = 0.000248486 Jy ! Component: 400 - total flux cleaned = 0.00024973 Jy ! Total flux subtracted in 400 components = 0.00024973 Jy ! Clean residual min=-0.009794 max=0.009009 Jy/beam ! Clean residual mean=-0.000105 rms=0.002597 Jy/beam ! Combined flux in latest and established models = 1.9726 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 367 components and 1.9726 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127186Jy sigma=0.312633 ! Fit after self-cal, rms=0.127192Jy sigma=0.312624 ! 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.127192Jy sigma=0.312624 ! Fit after self-cal, rms=0.126009Jy sigma=0.307829 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.093 mas, bmaj=18.37 mas, bpa=-1.776 degrees ! Estimated noise=10.3428 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00505554 Jy ! Component: 100 - total flux cleaned = -0.00505975 Jy ! Component: 150 - total flux cleaned = -0.00506597 Jy ! Component: 200 - total flux cleaned = -0.00528911 Jy ! Component: 250 - total flux cleaned = -0.0054989 Jy ! Component: 300 - total flux cleaned = -0.00559757 Jy ! Component: 350 - total flux cleaned = -0.00569164 Jy ! Component: 400 - total flux cleaned = -0.00586595 Jy ! Total flux subtracted in 400 components = -0.00586595 Jy ! Clean residual min=-0.010920 max=0.011413 Jy/beam ! Clean residual mean=-0.000175 rms=0.003347 Jy/beam ! Combined flux in latest and established models = 1.96674 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 380 components and 1.96674 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125399Jy sigma=0.306863 ! Fit after self-cal, rms=0.125401Jy sigma=0.306795 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 (47.6, -50). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.011868 Jy ! Component: 100 - total flux cleaned = 0.0164804 Jy ! Component: 150 - total flux cleaned = 0.0169696 Jy ! Component: 200 - total flux cleaned = 0.0169715 Jy ! Component: 250 - total flux cleaned = 0.0171767 Jy ! Component: 300 - total flux cleaned = 0.0173686 Jy ! Component: 350 - total flux cleaned = 0.0177332 Jy ! Component: 400 - total flux cleaned = 0.0177348 Jy ! Total flux subtracted in 400 components = 0.0177348 Jy ! Clean residual min=-0.008961 max=0.010823 Jy/beam ! Clean residual mean=-0.000294 rms=0.002663 Jy/beam ! Combined flux in latest and established models = 1.98447 Jy ! Performing phase self-cal ! Adding 98 model components to the UV plane model. ! The established model now contains 478 components and 1.98447 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123511Jy sigma=0.303649 ! Fit after self-cal, rms=0.123388Jy sigma=0.303457 ! Inverting map ! Added new window around map position (-6.8, 102). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00515134 Jy ! Component: 100 - total flux cleaned = 0.00585714 Jy ! Component: 150 - total flux cleaned = 0.00606398 Jy ! Component: 200 - total flux cleaned = 0.0060637 Jy ! Component: 250 - total flux cleaned = 0.00642885 Jy ! Component: 300 - total flux cleaned = 0.00659515 Jy ! Component: 350 - total flux cleaned = 0.00691952 Jy ! Component: 400 - total flux cleaned = 0.00722806 Jy ! Total flux subtracted in 400 components = 0.00722806 Jy ! Clean residual min=-0.007896 max=0.007394 Jy/beam ! Clean residual mean=-0.000296 rms=0.002295 Jy/beam ! Combined flux in latest and established models = 1.9917 Jy ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 539 components and 1.9917 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121875Jy sigma=0.301898 ! Fit after self-cal, rms=0.121787Jy sigma=0.301736 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=3.322 mas, bmaj=10.87 mas, bpa=-5.358 degrees ! Estimated noise=7.7011 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.322 x 10.87 at -5.358 degrees (North through East) ! Clean map min=-0.02068 max=1.3472 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121787Jy sigma=0.301736 ! Fit after self-cal, rms=0.115866Jy sigma=0.282254 ! Inverting map and beam ! Estimated beam: bmin=3.322 mas, bmaj=10.86 mas, bpa=-5.336 degrees ! Estimated noise=7.69265 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000946187 Jy ! Component: 100 - total flux cleaned = 0.000515245 Jy ! Component: 150 - total flux cleaned = 0.00074555 Jy ! Component: 200 - total flux cleaned = 0.000746486 Jy ! Component: 250 - total flux cleaned = 0.000348906 Jy ! Component: 300 - total flux cleaned = -2.98733e-05 Jy ! Component: 350 - total flux cleaned = -0.000208089 Jy ! Component: 400 - total flux cleaned = -0.000206375 Jy ! Total flux subtracted in 400 components = -0.000206375 Jy ! Clean residual min=-0.007185 max=0.008329 Jy/beam ! Clean residual mean=-0.000074 rms=0.001912 Jy/beam ! Combined flux in latest and established models = 1.99149 Jy ! Performing phase self-cal ! Adding 79 model components to the UV plane model. ! The established model now contains 613 components and 1.99149 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114239Jy sigma=0.278476 ! Fit after self-cal, rms=0.114183Jy sigma=0.278336 ! Inverting map ! Added new window around map position (-74.4, 76.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00660101 Jy ! Component: 100 - total flux cleaned = 0.00805605 Jy ! Component: 150 - total flux cleaned = 0.00975947 Jy ! Component: 200 - total flux cleaned = 0.0102975 Jy ! Component: 250 - total flux cleaned = 0.0116683 Jy ! Component: 300 - total flux cleaned = 0.0121566 Jy ! Component: 350 - total flux cleaned = 0.0124702 Jy ! Component: 400 - total flux cleaned = 0.0124702 Jy ! Total flux subtracted in 400 components = 0.0124702 Jy ! Clean residual min=-0.007447 max=0.006265 Jy/beam ! Clean residual mean=-0.000108 rms=0.001715 Jy/beam ! Combined flux in latest and established models = 2.00396 Jy ! Performing phase self-cal ! Adding 81 model components to the UV plane model. ! The established model now contains 671 components and 2.00396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112885Jy sigma=0.274676 ! Fit after self-cal, rms=0.112851Jy sigma=0.274514 ! Inverting map ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112851Jy sigma=0.274514 ! Fit after self-cal, rms=0.112658Jy sigma=0.274159 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.646 mas, bmaj=8.677 mas, bpa=-7.311 degrees ! Estimated noise=9.40421 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.976013 Jy ! Component: 100 - total flux cleaned = 1.29709 Jy ! Component: 150 - total flux cleaned = 1.47908 Jy ! Component: 200 - total flux cleaned = 1.59485 Jy ! Component: 250 - total flux cleaned = 1.67419 Jy ! Component: 300 - total flux cleaned = 1.73358 Jy ! Component: 350 - total flux cleaned = 1.7792 Jy ! Component: 400 - total flux cleaned = 1.81462 Jy ! Total flux subtracted in 400 components = 1.81462 Jy ! Clean residual min=-0.016876 max=0.026661 Jy/beam ! Clean residual mean=0.000346 rms=0.004528 Jy/beam ! Combined flux in latest and established models = 1.81462 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 69 model components to the UV plane model. ! The established model now contains 69 components and 1.81462 Jy ! Inverting map and beam ! Estimated beam: bmin=3.323 mas, bmaj=10.86 mas, bpa=-5.303 degrees ! Estimated noise=7.68215 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (-5.6, 6.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0514207 Jy ! Component: 100 - total flux cleaned = 0.0824705 Jy ! Component: 150 - total flux cleaned = 0.104291 Jy ! Component: 200 - total flux cleaned = 0.120991 Jy ! Component: 250 - total flux cleaned = 0.134546 Jy ! Component: 300 - total flux cleaned = 0.144801 Jy ! Component: 350 - total flux cleaned = 0.150196 Jy ! Component: 400 - total flux cleaned = 0.155533 Jy ! Total flux subtracted in 400 components = 0.155533 Jy ! Clean residual min=-0.012063 max=0.009121 Jy/beam ! Clean residual mean=0.000012 rms=0.002689 Jy/beam ! Combined flux in latest and established models = 1.97016 Jy ! Adding 108 model components to the UV plane model. ! The established model now contains 171 components and 1.97016 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.109 mas, bmaj=18.41 mas, bpa=-1.674 degrees ! Estimated noise=10.3023 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0119633 Jy ! Component: 100 - total flux cleaned = 0.01778 Jy ! Component: 150 - total flux cleaned = 0.0218646 Jy ! Component: 200 - total flux cleaned = 0.0242313 Jy ! Component: 250 - total flux cleaned = 0.0257957 Jy ! Component: 300 - total flux cleaned = 0.0275453 Jy ! Component: 350 - total flux cleaned = 0.0289102 Jy ! Component: 400 - total flux cleaned = 0.0299452 Jy ! Total flux subtracted in 400 components = 0.0299452 Jy ! Clean residual min=-0.010020 max=0.007977 Jy/beam ! Clean residual mean=-0.000157 rms=0.002579 Jy/beam ! Combined flux in latest and established models = 2.0001 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.00123279 Jy ! Component: 100 - total flux cleaned = 0.00262988 Jy ! Component: 150 - total flux cleaned = 0.00352241 Jy ! Component: 200 - total flux cleaned = 0.004379 Jy ! Component: 250 - total flux cleaned = 0.00479579 Jy ! Component: 300 - total flux cleaned = 0.00559517 Jy ! Component: 350 - total flux cleaned = 0.00617569 Jy ! Component: 400 - total flux cleaned = 0.00673464 Jy ! Total flux subtracted in 400 components = 0.00673464 Jy ! Clean residual min=-0.008213 max=0.006833 Jy/beam ! Clean residual mean=-0.000201 rms=0.002271 Jy/beam ! Combined flux in latest and established models = 2.00684 Jy ! Adding 183 model components to the UV plane model. ! The established model now contains 346 components and 2.00684 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.323 mas, bmaj=10.86 mas, bpa=-5.303 degrees ! Estimated noise=7.68215 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00129311 Jy ! Component: 100 - total flux cleaned = 0.00155389 Jy ! Component: 150 - total flux cleaned = 0.0015391 Jy ! Component: 200 - total flux cleaned = 0.00129374 Jy ! Component: 250 - total flux cleaned = 0.00129569 Jy ! Component: 300 - total flux cleaned = 0.00152273 Jy ! Component: 350 - total flux cleaned = 0.00174305 Jy ! Component: 400 - total flux cleaned = 0.00236307 Jy ! Total flux subtracted in 400 components = 0.00236307 Jy ! Clean residual min=-0.008743 max=0.007178 Jy/beam ! Clean residual mean=-0.000088 rms=0.001870 Jy/beam ! Combined flux in latest and established models = 2.0092 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 412 components and 2.0092 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114476Jy sigma=0.277587 ! Fit after self-cal, rms=0.114299Jy sigma=0.277107 wmodel /scr/pima/rdv73_uvs/J1037-2934_S_map.mod ! Writing 412 model components to file: /scr/pima/rdv73_uvs/J1037-2934_S_map.mod wobs /scr/pima/rdv73_uvs/J1037-2934_S_uvs.fits ! Applying 798 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J1037-2934_S_uvs.fits wwins /scr/pima/rdv73_uvs/J1037-2934_S_map.win ! wwins: Wrote 8 windows to /scr/pima/rdv73_uvs/J1037-2934_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.00156047 Jy ! Component: 100 - total flux cleaned = -0.00241834 Jy ! Component: 150 - total flux cleaned = -0.00322809 Jy ! Component: 200 - total flux cleaned = -0.00400371 Jy ! Component: 250 - total flux cleaned = -0.00400381 Jy ! Component: 300 - total flux cleaned = -0.00400368 Jy ! Component: 350 - total flux cleaned = -0.00450586 Jy ! Component: 400 - total flux cleaned = -0.00479972 Jy ! Component: 450 - total flux cleaned = -0.00537626 Jy ! Component: 500 - total flux cleaned = -0.0058455 Jy ! Component: 550 - total flux cleaned = -0.00603117 Jy ! Component: 600 - total flux cleaned = -0.00639631 Jy ! Component: 650 - total flux cleaned = -0.0066657 Jy ! Component: 700 - total flux cleaned = -0.00701914 Jy ! Component: 750 - total flux cleaned = -0.00728023 Jy ! Component: 800 - total flux cleaned = -0.0073664 Jy ! Component: 850 - total flux cleaned = -0.00778825 Jy ! Component: 900 - total flux cleaned = -0.00787189 Jy ! Component: 950 - total flux cleaned = -0.00819999 Jy ! Component: 1000 - total flux cleaned = -0.00852401 Jy ! Component: 1050 - total flux cleaned = -0.0086044 Jy ! Component: 1100 - total flux cleaned = -0.00892044 Jy ! Component: 1150 - total flux cleaned = -0.00923185 Jy ! Component: 1200 - total flux cleaned = -0.00953852 Jy ! Component: 1250 - total flux cleaned = -0.00969015 Jy ! Component: 1300 - total flux cleaned = -0.00999128 Jy ! Component: 1350 - total flux cleaned = -0.00999189 Jy ! Component: 1400 - total flux cleaned = -0.0104352 Jy ! Component: 1450 - total flux cleaned = -0.0105069 Jy ! Component: 1500 - total flux cleaned = -0.0105794 Jy ! Component: 1550 - total flux cleaned = -0.010723 Jy ! Component: 1600 - total flux cleaned = -0.0107233 Jy ! Component: 1650 - total flux cleaned = -0.0109343 Jy ! Component: 1700 - total flux cleaned = -0.0111432 Jy ! Component: 1750 - total flux cleaned = -0.0113493 Jy ! Component: 1800 - total flux cleaned = -0.0114861 Jy ! Component: 1850 - total flux cleaned = -0.011622 Jy ! Component: 1900 - total flux cleaned = -0.0119582 Jy ! Component: 1950 - total flux cleaned = -0.0120917 Jy ! Component: 2000 - total flux cleaned = -0.0122905 Jy ! Component: 2050 - total flux cleaned = -0.0124222 Jy ! Component: 2100 - total flux cleaned = -0.0126833 Jy ! Component: 2150 - total flux cleaned = -0.0130071 Jy ! Component: 2200 - total flux cleaned = -0.0130073 Jy ! Component: 2250 - total flux cleaned = -0.0132622 Jy ! Component: 2300 - total flux cleaned = -0.0135783 Jy ! Component: 2350 - total flux cleaned = -0.0136413 Jy ! Component: 2400 - total flux cleaned = -0.0138281 Jy ! Component: 2450 - total flux cleaned = -0.0138898 Jy ! Component: 2500 - total flux cleaned = -0.0140126 Jy ! Component: 2550 - total flux cleaned = -0.0140132 Jy ! Component: 2600 - total flux cleaned = -0.0144378 Jy ! Component: 2650 - total flux cleaned = -0.0143177 Jy ! Component: 2700 - total flux cleaned = -0.014617 Jy ! Component: 2750 - total flux cleaned = -0.014795 Jy ! Component: 2800 - total flux cleaned = -0.0149135 Jy ! Component: 2850 - total flux cleaned = -0.0150307 Jy ! Component: 2900 - total flux cleaned = -0.0153218 Jy ! Component: 2950 - total flux cleaned = -0.0155528 Jy ! Component: 3000 - total flux cleaned = -0.0156103 Jy ! Component: 3050 - total flux cleaned = -0.0157817 Jy ! Component: 3100 - total flux cleaned = -0.0161786 Jy ! Component: 3150 - total flux cleaned = -0.0163478 Jy ! Component: 3200 - total flux cleaned = -0.0167401 Jy ! Component: 3250 - total flux cleaned = -0.0166846 Jy ! Component: 3300 - total flux cleaned = -0.0171828 Jy ! Component: 3350 - total flux cleaned = -0.0172929 Jy ! Component: 3400 - total flux cleaned = -0.0176213 Jy ! Component: 3450 - total flux cleaned = -0.0178932 Jy ! Component: 3500 - total flux cleaned = -0.0180557 Jy ! Component: 3550 - total flux cleaned = -0.0182708 Jy ! Component: 3600 - total flux cleaned = -0.0185381 Jy ! Component: 3650 - total flux cleaned = -0.0186448 Jy ! Component: 3700 - total flux cleaned = -0.0190679 Jy ! Component: 3750 - total flux cleaned = -0.0191723 Jy ! Component: 3800 - total flux cleaned = -0.0192244 Jy ! Component: 3850 - total flux cleaned = -0.0197961 Jy ! Component: 3900 - total flux cleaned = -0.0198474 Jy ! Component: 3950 - total flux cleaned = -0.0201561 Jy ! Component: 4000 - total flux cleaned = -0.0202588 Jy ! Component: 4050 - total flux cleaned = -0.0202588 Jy ! Total flux subtracted in 4096 components = -0.0205125 Jy ! Clean residual min=-0.002750 max=0.003538 Jy/beam ! Clean residual mean=-0.000032 rms=0.001061 Jy/beam ! Combined flux in latest and established models = 1.98869 Jy keep ! Adding 312 model components to the UV plane model. ! The established model now contains 723 components and 1.98869 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J1037-2934_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.323 x 10.86 at -5.303 degrees (North through East) ! Clean map min=-0.0086952 max=1.3474 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J1037-2934_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= 1.3473 rms= 0.00105 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J1037-2934_S_uvs.log closed on Sat Jan 9 17:48:26 2016