! Started logfile: /scr/pima/rdv73_uvs/J1257-3155_S_uvs.log on Sat Jan 9 17:49:50 2016 obs /scr/pima/rdv73_uvs/J1257-3155_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv73_uvs/J1257-3155_S_uva.fits ! AN table 1: 3 integrations on 15 of 15 possible baselines. ! AN table 2: 3 integrations on 21 of 21 possible baselines. ! AN table 3: 7 integrations on 3 of 3 possible baselines. ! AN table 4: 3 integrations on 6 of 6 possible baselines. ! AN table 5: 10 integrations on 28 of 28 possible baselines. ! AN table 6: 4 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.87471 visibilities/baseline/integration-bin. ! Found source: J1257-3155 ! ! 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 1508 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/J1257-3155_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:KP self-cal status=correctable weight=1 ! 1:LA 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:KK self-cal status=correctable weight=1 ! 3:MK 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: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 6 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 5. ! A total of 8 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 21 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! A total of 8 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 6 telescope corrections were flagged in sub-array 3. ! A total of 23 telescope corrections were flagged in sub-array 5. ! A total of 8 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 15 telescope corrections were flagged in sub-array 3. ! A total of 29 telescope corrections were flagged in sub-array 5. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.587945Jy sigma=4.953412 ! Fit after self-cal, rms=0.493411Jy sigma=1.450044 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 3:KK ! Flagging baselines of antenna 3:KK in the currently selected channels. flag 4:KK ! Flagging baselines of antenna 4:KK in the currently selected channels. flag 5:KK ! Flagging baselines of antenna 5:KK in the currently selected channels. flag 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=3.111 mas, bmaj=11.81 mas, bpa=-8.137 degrees ! Estimated noise=10.9154 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.847547 Jy ! Component: 100 - total flux cleaned = 1.1456 Jy ! Total flux subtracted in 100 components = 1.1456 Jy ! Clean residual min=-0.064959 max=0.146386 Jy/beam ! Clean residual mean=0.000960 rms=0.023048 Jy/beam ! Combined flux in latest and established models = 1.1456 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 21 components and 1.1456 Jy ! ! Correcting IF 1. ! A total of 15 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! ! Correcting IF 3. ! A total of 15 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 13 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.427533Jy sigma=1.206201 ! Fit after self-cal, rms=0.421499Jy sigma=1.187080 ! Inverting map and beam ! Estimated beam: bmin=3.13 mas, bmaj=12.29 mas, bpa=-7.28 degrees ! Estimated noise=11.0812 mJy/beam. ! Added new window around map position (3.2, 7.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.177141 Jy ! Component: 100 - total flux cleaned = 0.284287 Jy ! Total flux subtracted in 100 components = 0.284287 Jy ! Clean residual min=-0.041713 max=0.073553 Jy/beam ! Clean residual mean=0.000218 rms=0.013452 Jy/beam ! Combined flux in latest and established models = 1.42989 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 41 components and 1.42989 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.285695Jy sigma=0.754542 ! Fit after self-cal, rms=0.284863Jy sigma=0.751849 ! 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.799 mas, bmaj=15.35 mas, bpa=-7.575 degrees ! Estimated noise=8.42192 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.0890393 Jy ! Component: 100 - total flux cleaned = 0.130846 Jy ! Component: 150 - total flux cleaned = 0.152344 Jy ! Component: 200 - total flux cleaned = 0.156297 Jy ! Total flux subtracted in 200 components = 0.156297 Jy ! Clean residual min=-0.040222 max=0.036232 Jy/beam ! Clean residual mean=0.000363 rms=0.010618 Jy/beam ! Combined flux in latest and established models = 1.58619 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 86 components and 1.58619 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.250440Jy sigma=0.680321 ! Fit after self-cal, rms=0.250296Jy sigma=0.679824 ! 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.88 HN 1.24 KP 1.12 LA 1.04 ! PT 1.13 SC 0.88 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.94 HN 1.16 KP 1.00 LA 1.03 ! OV 1.07 PT 1.08 SC 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 9 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: ! BR 1.05 FD 0.90 KK 1.00* KP 1.00 ! LA 1.02 MK 0.85 OV 1.01 PT 1.06 ! ! 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.98 HN 1.32 KP 1.14 LA 1.00 ! PT 1.09 SC 0.96 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.99 HN 1.21 KP 1.09 LA 0.93 ! OV 1.11 PT 1.09 SC 1.05 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 9 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: ! BR 1.06 FD 0.98 KK 1.00* KP 1.03 ! LA 0.95 MK 0.92 OV 1.06 PT 1.06 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 1.01 HN 1.39 KP 1.15 LA 0.90 ! PT 1.15 SC 0.86 ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.00 HN 1.31 KP 1.06 LA 0.89 ! OV 1.20 PT 1.15 SC 0.93 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 9 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: ! BR 0.82 FD 0.99 KK 1.00* KP 1.04 ! LA 0.88 MK 0.94 OV 1.12 PT 1.09 ! ! 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.94 HN 1.29 KP 1.04 LA 0.75 ! PT 1.00 SC 0.89 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.99 HN 1.24 KP 0.98 LA 0.76 ! OV 0.99 PT 0.97 SC 0.98 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.00* KK 0.00* MK 0.00* ! ! A total of 9 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: ! BR 1.00* FD 0.95 KK 1.00* KP 0.96 ! LA 0.72 MK 1.06 OV 0.91 PT 0.95 ! ! Telescope amplitude corrections in sub-array 6: ! KK 0.00* MK 0.00* ! ! ! Fit before self-cal, rms=0.250296Jy sigma=0.679824 ! Fit after self-cal, rms=0.131171Jy sigma=0.370601 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.812 mas, bmaj=17.04 mas, bpa=-5.61 degrees ! Estimated noise=8.57137 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.00272245 Jy ! Component: 100 - total flux cleaned = 0.00365847 Jy ! Component: 150 - total flux cleaned = 0.00391344 Jy ! Component: 200 - total flux cleaned = 0.0041419 Jy ! Total flux subtracted in 200 components = 0.0041419 Jy ! Clean residual min=-0.012953 max=0.014851 Jy/beam ! Clean residual mean=-0.000099 rms=0.003832 Jy/beam ! Combined flux in latest and established models = 1.59033 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 101 components and 1.59033 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129906Jy sigma=0.366580 ! Fit after self-cal, rms=0.129742Jy sigma=0.366186 ! 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.129742Jy sigma=0.366186 ! Fit after self-cal, rms=0.130188Jy sigma=0.365804 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.811 mas, bmaj=17.05 mas, bpa=-5.547 degrees ! Estimated noise=8.58642 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00191946 Jy ! Component: 100 - total flux cleaned = 0.00240745 Jy ! Component: 150 - total flux cleaned = 0.00262259 Jy ! Component: 200 - total flux cleaned = 0.00280599 Jy ! Total flux subtracted in 200 components = 0.00280599 Jy ! Clean residual min=-0.012721 max=0.014616 Jy/beam ! Clean residual mean=-0.000085 rms=0.003710 Jy/beam ! Combined flux in latest and established models = 1.59314 Jy selfcal ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 107 components and 1.59314 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.129420Jy sigma=0.363450 ! Fit after self-cal, rms=0.129292Jy sigma=0.363326 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.000264326 Jy ! Component: 100 - total flux cleaned = 0.000228442 Jy ! Component: 150 - total flux cleaned = 0.000224953 Jy ! Component: 200 - total flux cleaned = 0.00035658 Jy ! Total flux subtracted in 200 components = 0.00035658 Jy ! Clean residual min=-0.012654 max=0.015069 Jy/beam ! Clean residual mean=-0.000088 rms=0.003649 Jy/beam ! Combined flux in latest and established models = 1.59349 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 118 components and 1.59349 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.128936Jy sigma=0.362253 ! Fit after self-cal, rms=0.128835Jy sigma=0.362164 ! Inverting map ! Added new window around map position (16, 15.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0141249 Jy ! Component: 100 - total flux cleaned = 0.0183524 Jy ! Component: 150 - total flux cleaned = 0.0200458 Jy ! Component: 200 - total flux cleaned = 0.0206369 Jy ! Total flux subtracted in 200 components = 0.0206369 Jy ! Clean residual min=-0.012298 max=0.010930 Jy/beam ! Clean residual mean=-0.000244 rms=0.003097 Jy/beam ! Combined flux in latest and established models = 1.61413 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 153 components and 1.61413 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126770Jy sigma=0.356692 ! Fit after self-cal, rms=0.126521Jy sigma=0.356147 ! 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.126521Jy sigma=0.356147 ! Fit after self-cal, rms=0.123390Jy sigma=0.345294 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.806 mas, bmaj=17.06 mas, bpa=-5.502 degrees ! Estimated noise=8.5917 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00577916 Jy ! Component: 100 - total flux cleaned = 0.00741908 Jy ! Component: 150 - total flux cleaned = 0.00796106 Jy ! Component: 200 - total flux cleaned = 0.00812454 Jy ! Total flux subtracted in 200 components = 0.00812454 Jy ! Clean residual min=-0.009476 max=0.009371 Jy/beam ! Clean residual mean=-0.000038 rms=0.002705 Jy/beam ! Combined flux in latest and established models = 1.62225 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 170 components and 1.62225 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122549Jy sigma=0.343499 ! Fit after self-cal, rms=0.122449Jy sigma=0.343401 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.0015082 Jy ! Component: 100 - total flux cleaned = 0.00151003 Jy ! Component: 150 - total flux cleaned = 0.00151071 Jy ! Component: 200 - total flux cleaned = 0.00164982 Jy ! Total flux subtracted in 200 components = 0.00164982 Jy ! Clean residual min=-0.009313 max=0.009221 Jy/beam ! Clean residual mean=-0.000049 rms=0.002660 Jy/beam ! Combined flux in latest and established models = 1.6239 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 184 components and 1.6239 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121967Jy sigma=0.342558 ! Fit after self-cal, rms=0.121927Jy sigma=0.342524 ! 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.121927Jy sigma=0.342524 ! Fit after self-cal, rms=0.121750Jy sigma=0.342258 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.112 mas, bmaj=22.37 mas, bpa=4.484 degrees ! Estimated noise=10.329 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00456234 Jy ! Component: 100 - total flux cleaned = 0.00481909 Jy ! Component: 150 - total flux cleaned = 0.00493194 Jy ! Component: 200 - total flux cleaned = 0.00502686 Jy ! Total flux subtracted in 200 components = 0.00502686 Jy ! Clean residual min=-0.009747 max=0.009275 Jy/beam ! Clean residual mean=-0.000019 rms=0.002901 Jy/beam ! Combined flux in latest and established models = 1.62893 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 205 components and 1.62893 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121448Jy sigma=0.341991 ! Fit after self-cal, rms=0.121403Jy sigma=0.341898 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.00102199 Jy ! Component: 100 - total flux cleaned = 0.00081807 Jy ! Component: 150 - total flux cleaned = 0.000816712 Jy ! Component: 200 - total flux cleaned = 0.000815152 Jy ! Total flux subtracted in 200 components = 0.000815152 Jy ! Clean residual min=-0.009656 max=0.009114 Jy/beam ! Clean residual mean=-0.000029 rms=0.002867 Jy/beam ! Combined flux in latest and established models = 1.62975 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 223 components and 1.62975 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.121211Jy sigma=0.341771 ! Fit after self-cal, rms=0.121188Jy sigma=0.341720 ! 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:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA 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:MK,true ! 3:MK 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: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 7 telescope corrections were flagged in sub-array 3. ! A total of 3 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 3. ! A total of 3 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 3. ! A total of 3 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 3. ! A total of 3 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.449117Jy sigma=1.502803 ! Fit after self-cal, rms=0.295708Jy sigma=0.787496 selfant "",false ! 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: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:KK self-cal status=correctable weight=1 ! 3:MK 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: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 14 telescope corrections were flagged in sub-array 3. ! A total of 23 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 14 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! A total of 26 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 14 telescope corrections were flagged in sub-array 3. ! A total of 23 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 14 telescope corrections were flagged in sub-array 3. ! A total of 29 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.295708Jy sigma=0.787496 ! Fit after self-cal, rms=0.127878Jy sigma=0.355407 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.863 mas, bmaj=12.76 mas, bpa=-7.696 degrees ! Estimated noise=10.4456 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.00964359 Jy ! Component: 100 - total flux cleaned = -0.0123181 Jy ! Component: 150 - total flux cleaned = -0.0149629 Jy ! Component: 200 - total flux cleaned = -0.0168481 Jy ! Total flux subtracted in 200 components = -0.0168481 Jy ! Clean residual min=-0.020958 max=0.013487 Jy/beam ! Clean residual mean=-0.000061 rms=0.003830 Jy/beam ! Combined flux in latest and established models = 1.6129 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 252 components and 1.6129 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125104Jy sigma=0.347668 ! Fit after self-cal, rms=0.124048Jy sigma=0.345609 ! 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.35 mas, bmaj=15.97 mas, bpa=-7.598 degrees ! Estimated noise=7.87081 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.000103762 Jy ! Component: 100 - total flux cleaned = 6.99193e-05 Jy ! Component: 150 - total flux cleaned = -9.21466e-05 Jy ! Component: 200 - total flux cleaned = -0.000241541 Jy ! Component: 250 - total flux cleaned = -0.000379662 Jy ! Component: 300 - total flux cleaned = -0.000381493 Jy ! Component: 350 - total flux cleaned = -0.000379741 Jy ! Component: 400 - total flux cleaned = -0.000382205 Jy ! Total flux subtracted in 400 components = -0.000382205 Jy ! Clean residual min=-0.013243 max=0.009948 Jy/beam ! Clean residual mean=-0.000037 rms=0.003211 Jy/beam ! Combined flux in latest and established models = 1.61252 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 285 components and 1.61252 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.123180Jy sigma=0.343891 ! Fit after self-cal, rms=0.123160Jy sigma=0.343751 ! 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.000789565 Jy ! Component: 100 - total flux cleaned = 0.000907387 Jy ! Component: 150 - total flux cleaned = 0.000681387 Jy ! Component: 200 - total flux cleaned = 0.000681648 Jy ! Component: 250 - total flux cleaned = 0.000683107 Jy ! Component: 300 - total flux cleaned = 0.000682784 Jy ! Component: 350 - total flux cleaned = 0.000682077 Jy ! Component: 400 - total flux cleaned = 0.000684572 Jy ! Total flux subtracted in 400 components = 0.000684572 Jy ! Clean residual min=-0.013190 max=0.009700 Jy/beam ! Clean residual mean=-0.000040 rms=0.003170 Jy/beam ! Combined flux in latest and established models = 1.6132 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 305 components and 1.6132 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.122699Jy sigma=0.342963 ! Fit after self-cal, rms=0.122692Jy sigma=0.342899 ! 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 21 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 5. ! A total of 8 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 21 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 36 telescope corrections were flagged in sub-array 5. ! A total of 8 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 21 telescope corrections were flagged in sub-array 3. ! A total of 26 telescope corrections were flagged in sub-array 5. ! A total of 8 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 21 telescope corrections were flagged in sub-array 3. ! A total of 32 telescope corrections were flagged in sub-array 5. ! A total of 8 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.445568Jy sigma=1.487222 ! Fit after self-cal, rms=0.121811Jy sigma=0.339571 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.368 mas, bmaj=16.01 mas, bpa=-7.533 degrees ! Estimated noise=7.90581 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00703718 Jy ! Component: 100 - total flux cleaned = -0.00738501 Jy ! Component: 150 - total flux cleaned = -0.00752776 Jy ! Component: 200 - total flux cleaned = -0.0075179 Jy ! Component: 250 - total flux cleaned = -0.00741018 Jy ! Component: 300 - total flux cleaned = -0.0075073 Jy ! Component: 350 - total flux cleaned = -0.00759676 Jy ! Component: 400 - total flux cleaned = -0.00759692 Jy ! Total flux subtracted in 400 components = -0.00759692 Jy ! Clean residual min=-0.011723 max=0.008853 Jy/beam ! Clean residual mean=-0.000031 rms=0.002692 Jy/beam ! Combined flux in latest and established models = 1.6056 Jy selfcal ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 331 components and 1.6056 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120989Jy sigma=0.336847 ! Fit after self-cal, rms=0.120959Jy sigma=0.336771 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.000378132 Jy ! Component: 100 - total flux cleaned = 0.000288377 Jy ! Component: 150 - total flux cleaned = 0.000116042 Jy ! Component: 200 - total flux cleaned = 0.000115367 Jy ! Component: 250 - total flux cleaned = 0.000114389 Jy ! Component: 300 - total flux cleaned = 3.44153e-05 Jy ! Component: 350 - total flux cleaned = 3.69404e-05 Jy ! Component: 400 - total flux cleaned = -3.77064e-05 Jy ! Total flux subtracted in 400 components = -3.77064e-05 Jy ! Clean residual min=-0.011586 max=0.008603 Jy/beam ! Clean residual mean=-0.000030 rms=0.002667 Jy/beam ! Combined flux in latest and established models = 1.60557 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 346 components and 1.60557 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.120698Jy sigma=0.336271 ! Fit after self-cal, rms=0.120687Jy sigma=0.336241 ! 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.120687Jy sigma=0.336241 ! Fit after self-cal, rms=0.120219Jy sigma=0.334391 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.378 mas, bmaj=16.03 mas, bpa=-7.5 degrees ! Estimated noise=7.92404 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00458143 Jy ! Component: 100 - total flux cleaned = -0.00432248 Jy ! Component: 150 - total flux cleaned = -0.00410355 Jy ! Component: 200 - total flux cleaned = -0.0040127 Jy ! Component: 250 - total flux cleaned = -0.00400949 Jy ! Component: 300 - total flux cleaned = -0.00409014 Jy ! Component: 350 - total flux cleaned = -0.00417019 Jy ! Component: 400 - total flux cleaned = -0.00424798 Jy ! Total flux subtracted in 400 components = -0.00424798 Jy ! Clean residual min=-0.010652 max=0.008168 Jy/beam ! Clean residual mean=-0.000033 rms=0.002458 Jy/beam ! Combined flux in latest and established models = 1.60132 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 373 components and 1.60132 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119808Jy sigma=0.333027 ! Fit after self-cal, rms=0.119778Jy sigma=0.332983 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.000409065 Jy ! Component: 100 - total flux cleaned = 0.000409418 Jy ! Component: 150 - total flux cleaned = 0.000336063 Jy ! Component: 200 - total flux cleaned = 0.000336173 Jy ! Component: 250 - total flux cleaned = 0.000335945 Jy ! Component: 300 - total flux cleaned = 0.000334971 Jy ! Component: 350 - total flux cleaned = 0.000268673 Jy ! Component: 400 - total flux cleaned = 0.000270116 Jy ! Total flux subtracted in 400 components = 0.000270116 Jy ! Clean residual min=-0.010503 max=0.008098 Jy/beam ! Clean residual mean=-0.000034 rms=0.002436 Jy/beam ! Combined flux in latest and established models = 1.60159 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 388 components and 1.60159 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119599Jy sigma=0.332614 ! Fit after self-cal, rms=0.119585Jy sigma=0.332591 ! 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.119585Jy sigma=0.332591 ! Fit after self-cal, rms=0.117925Jy sigma=0.325975 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=9.057 mas, bmaj=22.28 mas, bpa=2.479 degrees ! Estimated noise=9.90654 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00223561 Jy ! Component: 100 - total flux cleaned = -0.00190502 Jy ! Component: 150 - total flux cleaned = -0.00176768 Jy ! Component: 200 - total flux cleaned = -0.00176391 Jy ! Component: 250 - total flux cleaned = -0.00176162 Jy ! Component: 300 - total flux cleaned = -0.00176217 Jy ! Component: 350 - total flux cleaned = -0.00170783 Jy ! Component: 400 - total flux cleaned = -0.00170978 Jy ! Total flux subtracted in 400 components = -0.00170978 Jy ! Clean residual min=-0.011051 max=0.009425 Jy/beam ! Clean residual mean=0.000052 rms=0.002983 Jy/beam ! Combined flux in latest and established models = 1.59988 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 403 components and 1.59988 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117775Jy sigma=0.325602 ! Fit after self-cal, rms=0.117739Jy sigma=0.325547 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.00026939 Jy ! Component: 100 - total flux cleaned = 0.000330915 Jy ! Component: 150 - total flux cleaned = 0.000331922 Jy ! Component: 200 - total flux cleaned = 0.000331032 Jy ! Component: 250 - total flux cleaned = 0.000283295 Jy ! Component: 300 - total flux cleaned = 0.00033056 Jy ! Component: 350 - total flux cleaned = 0.00033144 Jy ! Component: 400 - total flux cleaned = 0.00033296 Jy ! Total flux subtracted in 400 components = 0.00033296 Jy ! Clean residual min=-0.011063 max=0.009323 Jy/beam ! Clean residual mean=0.000048 rms=0.002975 Jy/beam ! Combined flux in latest and established models = 1.60021 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 421 components and 1.60021 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.117697Jy sigma=0.325418 ! Fit after self-cal, rms=0.117668Jy sigma=0.325390 ! 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.384 mas, bmaj=16.05 mas, bpa=-7.47 degrees ! Estimated noise=7.93692 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.384 x 16.05 at -7.47 degrees (North through East) ! Clean map min=-0.012689 max=1.0937 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.117668Jy sigma=0.325390 ! Fit after self-cal, rms=0.109781Jy sigma=0.298944 ! Inverting map and beam ! Estimated beam: bmin=3.386 mas, bmaj=16.02 mas, bpa=-7.482 degrees ! Estimated noise=7.93782 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00335641 Jy ! Component: 100 - total flux cleaned = -0.00386998 Jy ! Component: 150 - total flux cleaned = -0.00413189 Jy ! Component: 200 - total flux cleaned = -0.00412901 Jy ! Component: 250 - total flux cleaned = -0.00420285 Jy ! Component: 300 - total flux cleaned = -0.00420421 Jy ! Component: 350 - total flux cleaned = -0.00427342 Jy ! Component: 400 - total flux cleaned = -0.00427117 Jy ! Total flux subtracted in 400 components = -0.00427117 Jy ! Clean residual min=-0.009226 max=0.008191 Jy/beam ! Clean residual mean=0.000058 rms=0.002278 Jy/beam ! Combined flux in latest and established models = 1.59594 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 453 components and 1.59594 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109526Jy sigma=0.298139 ! Fit after self-cal, rms=0.109513Jy sigma=0.298112 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -9.2626e-05 Jy ! Component: 100 - total flux cleaned = -0.000153232 Jy ! Component: 150 - total flux cleaned = -9.47336e-05 Jy ! Component: 200 - total flux cleaned = -0.00015159 Jy ! Component: 250 - total flux cleaned = -0.000208957 Jy ! Component: 300 - total flux cleaned = -0.000261936 Jy ! Component: 350 - total flux cleaned = -0.000313415 Jy ! Component: 400 - total flux cleaned = -0.000313139 Jy ! Total flux subtracted in 400 components = -0.000313139 Jy ! Clean residual min=-0.009135 max=0.008136 Jy/beam ! Clean residual mean=0.000061 rms=0.002257 Jy/beam ! Combined flux in latest and established models = 1.59563 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 463 components and 1.59563 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109390Jy sigma=0.297802 ! Fit after self-cal, rms=0.109382Jy sigma=0.297782 ! 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.109382Jy sigma=0.297782 ! Fit after self-cal, rms=0.109563Jy sigma=0.297577 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.879 mas, bmaj=12.76 mas, bpa=-7.55 degrees ! Estimated noise=10.6606 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.860627 Jy ! Component: 100 - total flux cleaned = 1.12703 Jy ! Component: 150 - total flux cleaned = 1.27766 Jy ! Component: 200 - total flux cleaned = 1.36748 Jy ! Component: 250 - total flux cleaned = 1.42422 Jy ! Component: 300 - total flux cleaned = 1.46191 Jy ! Component: 350 - total flux cleaned = 1.48932 Jy ! Component: 400 - total flux cleaned = 1.51048 Jy ! Total flux subtracted in 400 components = 1.51048 Jy ! Clean residual min=-0.018761 max=0.018367 Jy/beam ! Clean residual mean=0.000141 rms=0.004058 Jy/beam ! Combined flux in latest and established models = 1.51048 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 95 model components to the UV plane model. ! The established model now contains 95 components and 1.51048 Jy ! Inverting map and beam ! Estimated beam: bmin=3.39 mas, bmaj=16.03 mas, bpa=-7.461 degrees ! Estimated noise=7.94335 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 (8, 10.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0292439 Jy ! Component: 100 - total flux cleaned = 0.0462475 Jy ! Component: 150 - total flux cleaned = 0.057971 Jy ! Component: 200 - total flux cleaned = 0.0667743 Jy ! Component: 250 - total flux cleaned = 0.0731005 Jy ! Component: 300 - total flux cleaned = 0.0747262 Jy ! Component: 350 - total flux cleaned = 0.0757973 Jy ! Component: 400 - total flux cleaned = 0.0767872 Jy ! Total flux subtracted in 400 components = 0.0767872 Jy ! Clean residual min=-0.010135 max=0.008361 Jy/beam ! Clean residual mean=0.000112 rms=0.002361 Jy/beam ! Combined flux in latest and established models = 1.58727 Jy ! Adding 77 model components to the UV plane model. ! The established model now contains 170 components and 1.58727 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=9.081 mas, bmaj=22.35 mas, bpa=2.574 degrees ! Estimated noise=9.88586 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00299497 Jy ! Component: 100 - total flux cleaned = 0.00448955 Jy ! Component: 150 - total flux cleaned = 0.00533794 Jy ! Component: 200 - total flux cleaned = 0.00584953 Jy ! Component: 250 - total flux cleaned = 0.00631593 Jy ! Component: 300 - total flux cleaned = 0.00674285 Jy ! Component: 350 - total flux cleaned = 0.00704197 Jy ! Component: 400 - total flux cleaned = 0.00732672 Jy ! Total flux subtracted in 400 components = 0.00732672 Jy ! Clean residual min=-0.010838 max=0.008481 Jy/beam ! Clean residual mean=0.000112 rms=0.002691 Jy/beam ! Combined flux in latest and established models = 1.5946 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.000176178 Jy ! Component: 100 - total flux cleaned = 0.000261858 Jy ! Component: 150 - total flux cleaned = 0.00034355 Jy ! Component: 200 - total flux cleaned = 0.000423798 Jy ! Component: 250 - total flux cleaned = 0.000500797 Jy ! Component: 300 - total flux cleaned = 0.000502874 Jy ! Component: 350 - total flux cleaned = 0.000507801 Jy ! Component: 400 - total flux cleaned = 0.00057964 Jy ! Total flux subtracted in 400 components = 0.00057964 Jy ! Clean residual min=-0.010927 max=0.008380 Jy/beam ! Clean residual mean=0.000109 rms=0.002666 Jy/beam ! Combined flux in latest and established models = 1.59518 Jy ! Adding 52 model components to the UV plane model. ! The established model now contains 217 components and 1.59518 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.39 mas, bmaj=16.03 mas, bpa=-7.461 degrees ! Estimated noise=7.94335 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000222725 Jy ! Component: 100 - total flux cleaned = -0.000674128 Jy ! Component: 150 - total flux cleaned = -0.000666847 Jy ! Component: 200 - total flux cleaned = -0.000542284 Jy ! Component: 250 - total flux cleaned = -0.000426823 Jy ! Component: 300 - total flux cleaned = -0.000425533 Jy ! Component: 350 - total flux cleaned = -0.000424531 Jy ! Component: 400 - total flux cleaned = -0.000328882 Jy ! Total flux subtracted in 400 components = -0.000328882 Jy ! Clean residual min=-0.009212 max=0.008411 Jy/beam ! Clean residual mean=0.000062 rms=0.002221 Jy/beam ! Combined flux in latest and established models = 1.59485 Jy selfcal ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 262 components and 1.59485 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.109839Jy sigma=0.298142 ! Fit after self-cal, rms=0.109759Jy sigma=0.297665 wmodel /scr/pima/rdv73_uvs/J1257-3155_S_map.mod ! Writing 262 model components to file: /scr/pima/rdv73_uvs/J1257-3155_S_map.mod wobs /scr/pima/rdv73_uvs/J1257-3155_S_uvs.fits ! Applying 886 buffered edits. ! Writing UV FITS file: /scr/pima/rdv73_uvs/J1257-3155_S_uvs.fits wwins /scr/pima/rdv73_uvs/J1257-3155_S_map.win ! wwins: Wrote 4 windows to /scr/pima/rdv73_uvs/J1257-3155_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.00152452 Jy ! Component: 100 - total flux cleaned = -0.00211736 Jy ! Component: 150 - total flux cleaned = -0.00171222 Jy ! Component: 200 - total flux cleaned = -0.00106287 Jy ! Component: 250 - total flux cleaned = -0.000564605 Jy ! Component: 300 - total flux cleaned = -0.000683652 Jy ! Component: 350 - total flux cleaned = -0.000798578 Jy ! Component: 400 - total flux cleaned = -0.00124726 Jy ! Component: 450 - total flux cleaned = -0.00157308 Jy ! Component: 500 - total flux cleaned = -0.00220716 Jy ! Component: 550 - total flux cleaned = -0.00251481 Jy ! Component: 600 - total flux cleaned = -0.00301447 Jy ! Component: 650 - total flux cleaned = -0.00359949 Jy ! Component: 700 - total flux cleaned = -0.00398192 Jy ! Component: 750 - total flux cleaned = -0.00454358 Jy ! Component: 800 - total flux cleaned = -0.00500112 Jy ! Component: 850 - total flux cleaned = -0.00553951 Jy ! Component: 900 - total flux cleaned = -0.00580324 Jy ! Component: 950 - total flux cleaned = -0.0058903 Jy ! Component: 1000 - total flux cleaned = -0.00597413 Jy ! Component: 1050 - total flux cleaned = -0.00572246 Jy ! Component: 1100 - total flux cleaned = -0.00547531 Jy ! Component: 1150 - total flux cleaned = -0.00523262 Jy ! Component: 1200 - total flux cleaned = -0.00491261 Jy ! Component: 1250 - total flux cleaned = -0.00483423 Jy ! Component: 1300 - total flux cleaned = -0.00428809 Jy ! Component: 1350 - total flux cleaned = -0.004365 Jy ! Component: 1400 - total flux cleaned = -0.00375765 Jy ! Component: 1450 - total flux cleaned = -0.003309 Jy ! Component: 1500 - total flux cleaned = -0.00301333 Jy ! Component: 1550 - total flux cleaned = -0.00250317 Jy ! Component: 1600 - total flux cleaned = -0.0022152 Jy ! Component: 1650 - total flux cleaned = -0.00185874 Jy ! Component: 1700 - total flux cleaned = -0.00150632 Jy ! Component: 1750 - total flux cleaned = -0.00108773 Jy ! Component: 1800 - total flux cleaned = -0.000949573 Jy ! Component: 1850 - total flux cleaned = -0.000607043 Jy ! Component: 1900 - total flux cleaned = -0.000539216 Jy ! Component: 1950 - total flux cleaned = -0.000270932 Jy ! Component: 2000 - total flux cleaned = 6.16308e-05 Jy ! Component: 2050 - total flux cleaned = 0.000193422 Jy ! Component: 2100 - total flux cleaned = 0.000520215 Jy ! Component: 2150 - total flux cleaned = 0.000713693 Jy ! Component: 2200 - total flux cleaned = 0.000842652 Jy ! Component: 2250 - total flux cleaned = 0.00122532 Jy ! Component: 2300 - total flux cleaned = 0.00128794 Jy ! Component: 2350 - total flux cleaned = 0.00160153 Jy ! Component: 2400 - total flux cleaned = 0.00184997 Jy ! Component: 2450 - total flux cleaned = 0.00209662 Jy ! Component: 2500 - total flux cleaned = 0.00215828 Jy ! Component: 2550 - total flux cleaned = 0.00240126 Jy ! Component: 2600 - total flux cleaned = 0.00264245 Jy ! Component: 2650 - total flux cleaned = 0.00294089 Jy ! Component: 2700 - total flux cleaned = 0.00323717 Jy ! Component: 2750 - total flux cleaned = 0.00329595 Jy ! Component: 2800 - total flux cleaned = 0.00364686 Jy ! Component: 2850 - total flux cleaned = 0.00376262 Jy ! Component: 2900 - total flux cleaned = 0.00410806 Jy ! Component: 2950 - total flux cleaned = 0.00422241 Jy ! Component: 3000 - total flux cleaned = 0.00439245 Jy ! Component: 3050 - total flux cleaned = 0.00461805 Jy ! Component: 3100 - total flux cleaned = 0.0050086 Jy ! Component: 3150 - total flux cleaned = 0.00511944 Jy ! Component: 3200 - total flux cleaned = 0.00517442 Jy ! Component: 3250 - total flux cleaned = 0.00555717 Jy ! Component: 3300 - total flux cleaned = 0.00577361 Jy ! Component: 3350 - total flux cleaned = 0.00588115 Jy ! Component: 3400 - total flux cleaned = 0.00598776 Jy ! Component: 3450 - total flux cleaned = 0.00614709 Jy ! Component: 3500 - total flux cleaned = 0.0063059 Jy ! Component: 3550 - total flux cleaned = 0.00646369 Jy ! Component: 3600 - total flux cleaned = 0.00662029 Jy ! Component: 3650 - total flux cleaned = 0.00682771 Jy ! Component: 3700 - total flux cleaned = 0.00713647 Jy ! Component: 3750 - total flux cleaned = 0.0072897 Jy ! Component: 3800 - total flux cleaned = 0.0073912 Jy ! Component: 3850 - total flux cleaned = 0.00754267 Jy ! Component: 3900 - total flux cleaned = 0.00759271 Jy ! Component: 3950 - total flux cleaned = 0.00774241 Jy ! Component: 4000 - total flux cleaned = 0.00779209 Jy ! Component: 4050 - total flux cleaned = 0.00779251 Jy ! Total flux subtracted in 4096 components = 0.00789058 Jy ! Clean residual min=-0.002626 max=0.002657 Jy/beam ! Clean residual mean=0.000066 rms=0.001020 Jy/beam ! Combined flux in latest and established models = 1.60274 Jy keep ! Adding 322 model components to the UV plane model. ! The established model now contains 584 components and 1.60274 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv73_uvs/J1257-3155_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.39 x 16.03 at -7.461 degrees (North through East) ! Clean map min=-0.010256 max=1.1106 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv73_uvs/J1257-3155_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.1105 rms= 0.00102 Jy quit ! Quitting program ! Log file /scr/pima/rdv73_uvs/J1257-3155_S_uvs.log closed on Sat Jan 9 17:49:56 2016