! Started logfile: /scr/pima/rdv93_uvs/J0922-3959_X_uvs.log on Tue Oct 8 15:45:34 2013 obs /scr/pima/rdv93_uvs/J0922-3959_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv93_uvs/J0922-3959_X_uva.fits ! AN table 1: 9 integrations on 21 of 21 possible baselines. ! AN table 2: 9 integrations on 15 of 15 possible baselines. ! AN table 3: 11 integrations on 21 of 21 possible baselines. ! AN table 4: 10 integrations on 28 of 28 possible baselines. ! AN table 5: 10 integrations on 15 of 15 possible baselines. ! AN table 6: 17 integrations on 1 of 1 possible baselines. ! Apparent sampling: 0.845309 visibilities/baseline/integration-bin. ! Found source: J0922-3959 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3388 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv93_uvs/J0922-3959_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:FD self-cal status=correctable weight=1 ! 1:FT self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:FT 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 ! 5:SC self-cal status=correctable weight=1 ! 6:FT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 26 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 17 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 26 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 17 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 26 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 17 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 26 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 17 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=1.445013Jy sigma=24.274906 ! Fit after self-cal, rms=0.317760Jy sigma=4.467950 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:FT ! Flagging baselines of antenna 1:FT in the currently selected channels. flag 2:FT ! Flagging baselines of antenna 2:FT in the currently selected channels. flag 3:FT ! Flagging baselines of antenna 3:FT in the currently selected channels. flag 4:FT ! Flagging baselines of antenna 4:FT in the currently selected channels. flag 5:FT ! Flagging baselines of antenna 5:FT in the currently selected channels. flag 6:FT ! Flagging baselines of antenna 6:FT in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.6832 mas, bmaj=2.514 mas, bpa=4.851 degrees ! Estimated noise=1.69071 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.772693 Jy ! Component: 100 - total flux cleaned = 1.00404 Jy ! Total flux subtracted in 100 components = 1.00404 Jy ! Clean residual min=-0.059764 max=0.103437 Jy/beam ! Clean residual mean=0.000070 rms=0.014895 Jy/beam ! Combined flux in latest and established models = 1.00404 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 19 components and 1.00404 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.245127Jy sigma=4.189650 ! Fit after self-cal, rms=0.244719Jy sigma=4.186509 ! Inverting map ! Added new window around map position (-0.2, -2.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.104882 Jy ! Component: 100 - total flux cleaned = 0.167348 Jy ! Total flux subtracted in 100 components = 0.167348 Jy ! Clean residual min=-0.036951 max=0.061763 Jy/beam ! Clean residual mean=0.000054 rms=0.008958 Jy/beam ! Combined flux in latest and established models = 1.17139 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 39 components and 1.17139 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.197653Jy sigma=2.668997 ! Fit after self-cal, rms=0.197005Jy sigma=2.656462 ! Inverting map ! Added new window around map position (-1.8, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0535491 Jy ! Component: 100 - total flux cleaned = 0.0857356 Jy ! Total flux subtracted in 100 components = 0.0857356 Jy ! Clean residual min=-0.029901 max=0.051045 Jy/beam ! Clean residual mean=0.000053 rms=0.007151 Jy/beam ! Combined flux in latest and established models = 1.25713 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 48 components and 1.25713 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.176739Jy sigma=2.201674 ! Fit after self-cal, rms=0.175949Jy sigma=2.161892 ! 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=0.885 mas, bmaj=3.202 mas, bpa=4.322 degrees ! Estimated noise=1.0699 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.0417584 Jy ! Component: 100 - total flux cleaned = 0.0636355 Jy ! Component: 150 - total flux cleaned = 0.0723589 Jy ! Component: 200 - total flux cleaned = 0.0744832 Jy ! Total flux subtracted in 200 components = 0.0744832 Jy ! Clean residual min=-0.031804 max=0.033355 Jy/beam ! Clean residual mean=0.000056 rms=0.006664 Jy/beam ! Combined flux in latest and established models = 1.33161 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 59 components and 1.33161 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166485Jy sigma=1.931980 ! Fit after self-cal, rms=0.166299Jy sigma=1.912283 ! 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.98 FT 1.00* KP 0.97 LA 1.00 ! NL 1.35 PT 1.05 SC 0.96 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.94 FT 1.00* LA 1.01 NL 1.12 ! PT 0.97 SC 1.02 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 0.97 FT 1.00* LA 1.00 NL 1.05 ! OV 0.96 PT 0.98 SC 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.99 FT 1.00* LA 0.98 MK 1.02 ! NL 1.30 OV 0.94 PT 0.99 SC 1.01 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.94 FT 1.00* MK 1.06 OV 0.95 ! PT 0.99 SC 0.99 ! ! Telescope amplitude corrections in sub-array 6: ! FT 0.00* SC 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! FD 1.03 FT 1.00* KP 0.87 LA 0.99 ! NL 1.34 PT 1.00 SC 0.99 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.97 FT 1.00* LA 0.99 NL 1.14 ! PT 0.97 SC 1.04 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 0.98 FT 1.00* LA 0.98 NL 1.06 ! OV 0.95 PT 0.98 SC 1.04 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.99 FT 1.00* LA 0.95 MK 1.02 ! NL 1.34 OV 0.93 PT 1.00 SC 1.05 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.96 FT 1.00* MK 1.07 OV 0.93 ! PT 1.00 SC 1.02 ! ! Telescope amplitude corrections in sub-array 6: ! FT 0.00* SC 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! FD 1.00 FT 1.00* KP 0.83 LA 1.00 ! NL 1.36 PT 1.01 SC 0.99 ! ! Telescope amplitude corrections in sub-array 2: ! FD 0.97 FT 1.00* LA 1.08 NL 1.16 ! PT 0.98 SC 1.01 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 0.96 FT 1.00* LA 1.04 NL 1.13 ! OV 1.00 PT 0.99 SC 1.00 ! ! Telescope amplitude corrections in sub-array 4: ! FD 0.98 FT 1.00* LA 1.02 MK 1.01 ! NL 1.45 OV 0.96 PT 1.00 SC 1.03 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 0.94 FT 1.00* MK 1.05 OV 1.00 ! PT 1.00 SC 1.00 ! ! Telescope amplitude corrections in sub-array 6: ! FT 0.00* SC 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! FD 1.09 FT 1.00* KP 0.84 LA 1.05 ! NL 1.36 PT 1.06 SC 1.01 ! ! Telescope amplitude corrections in sub-array 2: ! FD 1.03 FT 1.00* LA 1.05 NL 1.21 ! PT 1.00 SC 1.03 ! ! A total of 3 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! FD 1.03 FT 1.00* LA 1.00 NL 1.14 ! OV 0.98 PT 1.02 SC 1.03 ! ! Telescope amplitude corrections in sub-array 4: ! FD 1.07 FT 1.00* LA 0.99 MK 1.09 ! NL 1.51 OV 0.95 PT 1.03 SC 1.00 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! FD 1.03 FT 1.00* MK 1.05 OV 0.96 ! PT 1.05 SC 1.02 ! ! Telescope amplitude corrections in sub-array 6: ! FT 0.00* SC 0.00* ! ! ! Fit before self-cal, rms=0.166299Jy sigma=1.912283 ! Fit after self-cal, rms=0.129563Jy sigma=1.329205 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.8884 mas, bmaj=3.253 mas, bpa=4.381 degrees ! Estimated noise=1.09049 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.00928479 Jy ! Component: 100 - total flux cleaned = -0.00975176 Jy ! Component: 150 - total flux cleaned = -0.0103464 Jy ! Component: 200 - total flux cleaned = -0.010705 Jy ! Total flux subtracted in 200 components = -0.010705 Jy ! Clean residual min=-0.009541 max=0.009174 Jy/beam ! Clean residual mean=0.000098 rms=0.002060 Jy/beam ! Combined flux in latest and established models = 1.3209 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 73 components and 1.3209 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127551Jy sigma=1.314448 ! Fit after self-cal, rms=0.127496Jy sigma=1.311845 ! 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.127496Jy sigma=1.311845 ! Fit after self-cal, rms=0.128007Jy sigma=1.308055 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8894 mas, bmaj=3.258 mas, bpa=4.354 degrees ! Estimated noise=1.09225 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00532238 Jy ! Component: 100 - total flux cleaned = -0.00585964 Jy ! Component: 150 - total flux cleaned = -0.00617743 Jy ! Component: 200 - total flux cleaned = -0.00646931 Jy ! Total flux subtracted in 200 components = -0.00646931 Jy ! Clean residual min=-0.008979 max=0.008163 Jy/beam ! Clean residual mean=0.000110 rms=0.001912 Jy/beam ! Combined flux in latest and established models = 1.31443 Jy selfcal ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 79 components and 1.31444 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.127027Jy sigma=1.301957 ! Fit after self-cal, rms=0.126878Jy sigma=1.300875 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.00131659 Jy ! Component: 100 - total flux cleaned = -0.00157562 Jy ! Component: 150 - total flux cleaned = -0.00192678 Jy ! Component: 200 - total flux cleaned = -0.0020331 Jy ! Total flux subtracted in 200 components = -0.0020331 Jy ! Clean residual min=-0.008804 max=0.008072 Jy/beam ! Clean residual mean=0.000108 rms=0.001894 Jy/beam ! Combined flux in latest and established models = 1.3124 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 85 components and 1.3124 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126429Jy sigma=1.298525 ! Fit after self-cal, rms=0.126392Jy sigma=1.298185 ! 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.126392Jy sigma=1.298185 ! Fit after self-cal, rms=0.127118Jy sigma=1.295848 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8899 mas, bmaj=3.261 mas, bpa=4.352 degrees ! Estimated noise=1.09336 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00333147 Jy ! Component: 100 - total flux cleaned = -0.00401871 Jy ! Component: 150 - total flux cleaned = -0.0043765 Jy ! Component: 200 - total flux cleaned = -0.00459512 Jy ! Total flux subtracted in 200 components = -0.00459512 Jy ! Clean residual min=-0.008508 max=0.007703 Jy/beam ! Clean residual mean=0.000112 rms=0.001816 Jy/beam ! Combined flux in latest and established models = 1.30781 Jy selfcal ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 92 components and 1.30781 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126531Jy sigma=1.293409 ! Fit after self-cal, rms=0.126487Jy sigma=1.292902 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.000817155 Jy ! Component: 100 - total flux cleaned = -0.00101387 Jy ! Component: 150 - total flux cleaned = -0.00128298 Jy ! Component: 200 - total flux cleaned = -0.00144454 Jy ! Total flux subtracted in 200 components = -0.00144454 Jy ! Clean residual min=-0.008391 max=0.007659 Jy/beam ! Clean residual mean=0.000111 rms=0.001804 Jy/beam ! Combined flux in latest and established models = 1.30636 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 97 components and 1.30636 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126209Jy sigma=1.291728 ! Fit after self-cal, rms=0.126204Jy sigma=1.291506 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126204Jy sigma=1.291506 ! Fit after self-cal, rms=0.126824Jy sigma=1.290031 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.604 mas, bmaj=9.834 mas, bpa=-1.519 degrees ! Estimated noise=1.21048 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00291815 Jy ! Component: 100 - total flux cleaned = -0.00325116 Jy ! Component: 150 - total flux cleaned = -0.00352571 Jy ! Component: 200 - total flux cleaned = -0.00371157 Jy ! Total flux subtracted in 200 components = -0.00371157 Jy ! Clean residual min=-0.009067 max=0.008951 Jy/beam ! Clean residual mean=0.000205 rms=0.002226 Jy/beam ! Combined flux in latest and established models = 1.30265 Jy selfcal ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 101 components and 1.30265 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126593Jy sigma=1.289435 ! Fit after self-cal, rms=0.126536Jy sigma=1.289044 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -3.89887e-05 Jy ! Component: 100 - total flux cleaned = -0.000241992 Jy ! Component: 150 - total flux cleaned = -0.000413905 Jy ! Component: 200 - total flux cleaned = -0.00056799 Jy ! Total flux subtracted in 200 components = -0.00056799 Jy ! Clean residual min=-0.009046 max=0.009056 Jy/beam ! Clean residual mean=0.000204 rms=0.002231 Jy/beam ! Combined flux in latest and established models = 1.30208 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 103 components and 1.30208 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126403Jy sigma=1.288858 ! Fit after self-cal, rms=0.126358Jy sigma=1.288577 ! 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:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1: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:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 3:PT,true ! 3:PT self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD 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 5:SC,true ! 5:SC self-cal status=fixed weight=1 selfant 6:SC,true ! 6:SC 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 9 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 10 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 10 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 10 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 11 telescope corrections were flagged in sub-array 3. ! A total of 10 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.193715Jy sigma=3.344347 ! Fit after self-cal, rms=0.193715Jy sigma=3.344347 selfant "",false ! 1:FD self-cal status=correctable weight=1 ! 1:FT self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:FT self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:FT self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:FT self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:NL self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:FT 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 ! 5:SC self-cal status=correctable weight=1 ! 6:FT self-cal status=correctable weight=1 ! 6:SC 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 1 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 7 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 7 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 7 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 1 telescope corrections were flagged in sub-array 1. ! A total of 5 telescope corrections were flagged in sub-array 2. ! A total of 18 telescope corrections were flagged in sub-array 3. ! A total of 7 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.193715Jy sigma=3.344347 ! Fit after self-cal, rms=0.126391Jy sigma=1.287670 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.684 mas, bmaj=2.535 mas, bpa=5.071 degrees ! Estimated noise=1.74837 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.00121741 Jy ! Component: 100 - total flux cleaned = -0.000808058 Jy ! Component: 150 - total flux cleaned = -0.000810303 Jy ! Component: 200 - total flux cleaned = -0.000959029 Jy ! Total flux subtracted in 200 components = -0.000959029 Jy ! Clean residual min=-0.010580 max=0.009679 Jy/beam ! Clean residual mean=0.000045 rms=0.002367 Jy/beam ! Combined flux in latest and established models = 1.30112 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 115 components and 1.30112 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.126062Jy sigma=1.287763 ! Fit after self-cal, rms=0.126086Jy sigma=1.285884 ! 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=0.8903 mas, bmaj=3.266 mas, bpa=4.357 degrees ! Estimated noise=1.09313 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.00266148 Jy ! Component: 100 - total flux cleaned = -0.00299949 Jy ! Component: 150 - total flux cleaned = -0.00321623 Jy ! Component: 200 - total flux cleaned = -0.00328444 Jy ! Component: 250 - total flux cleaned = -0.00334618 Jy ! Component: 300 - total flux cleaned = -0.00334804 Jy ! Component: 350 - total flux cleaned = -0.00334983 Jy ! Component: 400 - total flux cleaned = -0.00340459 Jy ! Total flux subtracted in 400 components = -0.00340459 Jy ! Clean residual min=-0.007608 max=0.007178 Jy/beam ! Clean residual mean=0.000118 rms=0.001702 Jy/beam ! Combined flux in latest and established models = 1.29772 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 122 components and 1.29772 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125708Jy sigma=1.284664 ! Fit after self-cal, rms=0.125721Jy sigma=1.284431 ! 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 = -3.63114e-05 Jy ! Component: 100 - total flux cleaned = -9.93295e-05 Jy ! Component: 150 - total flux cleaned = -4.46602e-05 Jy ! Component: 200 - total flux cleaned = -4.53547e-05 Jy ! Component: 250 - total flux cleaned = -9.61498e-05 Jy ! Component: 300 - total flux cleaned = -9.51931e-05 Jy ! Component: 350 - total flux cleaned = -9.53147e-05 Jy ! Component: 400 - total flux cleaned = -9.5629e-05 Jy ! Total flux subtracted in 400 components = -9.5629e-05 Jy ! Clean residual min=-0.007475 max=0.007245 Jy/beam ! Clean residual mean=0.000118 rms=0.001688 Jy/beam ! Combined flux in latest and established models = 1.29762 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 125 components and 1.29762 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125556Jy sigma=1.283517 ! Fit after self-cal, rms=0.125582Jy sigma=1.283301 ! 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 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 29 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 29 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 29 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 14 telescope corrections were flagged in sub-array 2. ! A total of 29 telescope corrections were flagged in sub-array 3. ! A total of 17 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 34 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.193334Jy sigma=3.342078 ! Fit after self-cal, rms=0.126353Jy sigma=1.281972 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8906 mas, bmaj=3.269 mas, bpa=4.367 degrees ! Estimated noise=1.09362 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00203084 Jy ! Component: 100 - total flux cleaned = -0.00253196 Jy ! Component: 150 - total flux cleaned = -0.00282718 Jy ! Component: 200 - total flux cleaned = -0.00309318 Jy ! Component: 250 - total flux cleaned = -0.00327704 Jy ! Component: 300 - total flux cleaned = -0.00344857 Jy ! Component: 350 - total flux cleaned = -0.00350255 Jy ! Component: 400 - total flux cleaned = -0.00355407 Jy ! Total flux subtracted in 400 components = -0.00355407 Jy ! Clean residual min=-0.007329 max=0.007141 Jy/beam ! Clean residual mean=0.000112 rms=0.001642 Jy/beam ! Combined flux in latest and established models = 1.29407 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 127 components and 1.29407 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125973Jy sigma=1.281125 ! Fit after self-cal, rms=0.126005Jy sigma=1.280868 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.000123264 Jy ! Component: 100 - total flux cleaned = -0.000237347 Jy ! Component: 150 - total flux cleaned = -0.000342799 Jy ! Component: 200 - total flux cleaned = -0.000340835 Jy ! Component: 250 - total flux cleaned = -0.000340944 Jy ! Component: 300 - total flux cleaned = -0.000386775 Jy ! Component: 350 - total flux cleaned = -0.000430418 Jy ! Component: 400 - total flux cleaned = -0.000388453 Jy ! Total flux subtracted in 400 components = -0.000388453 Jy ! Clean residual min=-0.007242 max=0.007279 Jy/beam ! Clean residual mean=0.000112 rms=0.001627 Jy/beam ! Combined flux in latest and established models = 1.29368 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 128 components and 1.29368 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125844Jy sigma=1.280169 ! Fit after self-cal, rms=0.125869Jy sigma=1.279992 ! Inverting map ! Added new window around map position (-22.2, 34.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00676372 Jy ! Component: 100 - total flux cleaned = 0.00853592 Jy ! Component: 150 - total flux cleaned = 0.00928843 Jy ! Component: 200 - total flux cleaned = 0.00964691 Jy ! Component: 250 - total flux cleaned = 0.00989495 Jy ! Component: 300 - total flux cleaned = 0.0101266 Jy ! Component: 350 - total flux cleaned = 0.0102698 Jy ! Component: 400 - total flux cleaned = 0.0104052 Jy ! Total flux subtracted in 400 components = 0.0104052 Jy ! Clean residual min=-0.006493 max=0.006117 Jy/beam ! Clean residual mean=0.000125 rms=0.001470 Jy/beam ! Combined flux in latest and established models = 1.30409 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 166 components and 1.30409 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125209Jy sigma=1.272848 ! Fit after self-cal, rms=0.125156Jy sigma=1.271507 ! 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.125156Jy sigma=1.271507 ! Fit after self-cal, rms=0.125399Jy sigma=1.269314 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.8911 mas, bmaj=3.272 mas, bpa=4.375 degrees ! Estimated noise=1.09377 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00327076 Jy ! Component: 100 - total flux cleaned = 0.00397989 Jy ! Component: 150 - total flux cleaned = 0.00405647 Jy ! Component: 200 - total flux cleaned = 0.00412662 Jy ! Component: 250 - total flux cleaned = 0.00419114 Jy ! Component: 300 - total flux cleaned = 0.00425373 Jy ! Component: 350 - total flux cleaned = 0.00425306 Jy ! Component: 400 - total flux cleaned = 0.00425515 Jy ! Total flux subtracted in 400 components = 0.00425515 Jy ! Clean residual min=-0.005979 max=0.005485 Jy/beam ! Clean residual mean=0.000111 rms=0.001370 Jy/beam ! Combined flux in latest and established models = 1.30834 Jy selfcal ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 188 components and 1.30834 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.125030Jy sigma=1.266953 ! Fit after self-cal, rms=0.125027Jy sigma=1.266337 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.00137529 Jy ! Component: 100 - total flux cleaned = 0.00168268 Jy ! Component: 150 - total flux cleaned = 0.00173615 Jy ! Component: 200 - total flux cleaned = 0.0017367 Jy ! Component: 250 - total flux cleaned = 0.001788 Jy ! Component: 300 - total flux cleaned = 0.0017875 Jy ! Component: 350 - total flux cleaned = 0.00178816 Jy ! Component: 400 - total flux cleaned = 0.00178777 Jy ! Total flux subtracted in 400 components = 0.00178777 Jy ! Clean residual min=-0.005852 max=0.005475 Jy/beam ! Clean residual mean=0.000114 rms=0.001367 Jy/beam ! Combined flux in latest and established models = 1.31013 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 197 components and 1.31013 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124839Jy sigma=1.265615 ! Fit after self-cal, rms=0.124856Jy sigma=1.265470 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124856Jy sigma=1.265470 ! Fit after self-cal, rms=0.124309Jy sigma=1.207901 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.616 mas, bmaj=9.903 mas, bpa=-1.446 degrees ! Estimated noise=1.20761 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000521743 Jy ! Component: 100 - total flux cleaned = 0.00067555 Jy ! Component: 150 - total flux cleaned = 0.000849604 Jy ! Component: 200 - total flux cleaned = 0.000948021 Jy ! Component: 250 - total flux cleaned = 0.00107524 Jy ! Component: 300 - total flux cleaned = 0.00115119 Jy ! Component: 350 - total flux cleaned = 0.00122275 Jy ! Component: 400 - total flux cleaned = 0.00125558 Jy ! Total flux subtracted in 400 components = 0.00125558 Jy ! Clean residual min=-0.006692 max=0.005205 Jy/beam ! Clean residual mean=0.000204 rms=0.001793 Jy/beam ! Combined flux in latest and established models = 1.31139 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 208 components and 1.31139 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124115Jy sigma=1.208372 ! Fit after self-cal, rms=0.124163Jy sigma=1.207911 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.000784949 Jy ! Component: 100 - total flux cleaned = 0.00118903 Jy ! Component: 150 - total flux cleaned = 0.00135146 Jy ! Component: 200 - total flux cleaned = 0.00145654 Jy ! Component: 250 - total flux cleaned = 0.00148789 Jy ! Component: 300 - total flux cleaned = 0.00151793 Jy ! Component: 350 - total flux cleaned = 0.00151771 Jy ! Component: 400 - total flux cleaned = 0.00154364 Jy ! Total flux subtracted in 400 components = 0.00154364 Jy ! Clean residual min=-0.006495 max=0.005064 Jy/beam ! Clean residual mean=0.000209 rms=0.001814 Jy/beam ! Combined flux in latest and established models = 1.31293 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 213 components and 1.31293 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124072Jy sigma=1.208354 ! Fit after self-cal, rms=0.124082Jy sigma=1.208134 ! 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=0.8915 mas, bmaj=3.275 mas, bpa=4.386 degrees ! Estimated noise=1.09365 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8915 x 3.275 at 4.386 degrees (North through East) ! Clean map min=-0.014527 max=0.94864 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.124082Jy sigma=1.208134 ! Fit after self-cal, rms=0.119696Jy sigma=1.013163 ! Inverting map and beam ! Estimated beam: bmin=0.8917 mas, bmaj=3.279 mas, bpa=4.379 degrees ! Estimated noise=1.09316 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000651727 Jy ! Component: 100 - total flux cleaned = 0.000665349 Jy ! Component: 150 - total flux cleaned = 0.000605329 Jy ! Component: 200 - total flux cleaned = 0.000550338 Jy ! Component: 250 - total flux cleaned = 0.000548249 Jy ! Component: 300 - total flux cleaned = 0.000646199 Jy ! Component: 350 - total flux cleaned = 0.000739738 Jy ! Component: 400 - total flux cleaned = 0.00078551 Jy ! Total flux subtracted in 400 components = 0.00078551 Jy ! Clean residual min=-0.005098 max=0.004730 Jy/beam ! Clean residual mean=0.000101 rms=0.001321 Jy/beam ! Combined flux in latest and established models = 1.31371 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 227 components and 1.31371 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119445Jy sigma=1.012080 ! Fit after self-cal, rms=0.119481Jy sigma=1.011760 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000479307 Jy ! Component: 100 - total flux cleaned = 0.000625795 Jy ! Component: 150 - total flux cleaned = 0.000717101 Jy ! Component: 200 - total flux cleaned = 0.000759637 Jy ! Component: 250 - total flux cleaned = 0.000845062 Jy ! Component: 300 - total flux cleaned = 0.000843873 Jy ! Component: 350 - total flux cleaned = 0.000924496 Jy ! Component: 400 - total flux cleaned = 0.00100365 Jy ! Total flux subtracted in 400 components = 0.00100365 Jy ! Clean residual min=-0.004971 max=0.004728 Jy/beam ! Clean residual mean=0.000103 rms=0.001318 Jy/beam ! Combined flux in latest and established models = 1.31472 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 232 components and 1.31472 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119332Jy sigma=1.011163 ! Fit after self-cal, rms=0.119362Jy sigma=1.011047 ! 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.119362Jy sigma=1.011047 ! Fit after self-cal, rms=0.119614Jy sigma=1.010522 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=0.684 mas, bmaj=2.542 mas, bpa=5.09 degrees ! Estimated noise=1.74732 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.793349 Jy ! Component: 100 - total flux cleaned = 1.02612 Jy ! Component: 150 - total flux cleaned = 1.12485 Jy ! Component: 200 - total flux cleaned = 1.1782 Jy ! Component: 250 - total flux cleaned = 1.21225 Jy ! Component: 300 - total flux cleaned = 1.23508 Jy ! Component: 350 - total flux cleaned = 1.25069 Jy ! Component: 400 - total flux cleaned = 1.26169 Jy ! Total flux subtracted in 400 components = 1.26169 Jy ! Clean residual min=-0.011429 max=0.012749 Jy/beam ! Clean residual mean=0.000020 rms=0.002542 Jy/beam ! Combined flux in latest and established models = 1.26169 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 77 model components to the UV plane model. ! The established model now contains 77 components and 1.26169 Jy ! Inverting map and beam ! Estimated beam: bmin=0.8919 mas, bmaj=3.279 mas, bpa=4.37 degrees ! Estimated noise=1.09277 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 (0.4, 4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0153437 Jy ! Component: 100 - total flux cleaned = 0.0260789 Jy ! Component: 150 - total flux cleaned = 0.0342722 Jy ! Component: 200 - total flux cleaned = 0.0382368 Jy ! Component: 250 - total flux cleaned = 0.0411436 Jy ! Component: 300 - total flux cleaned = 0.0425354 Jy ! Component: 350 - total flux cleaned = 0.0430862 Jy ! Component: 400 - total flux cleaned = 0.0437724 Jy ! Total flux subtracted in 400 components = 0.0437724 Jy ! Clean residual min=-0.005683 max=0.005727 Jy/beam ! Clean residual mean=0.000100 rms=0.001519 Jy/beam ! Combined flux in latest and established models = 1.30546 Jy ! Adding 63 model components to the UV plane model. ! The established model now contains 134 components and 1.30546 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.619 mas, bmaj=9.921 mas, bpa=-1.438 degrees ! Estimated noise=1.2058 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00355034 Jy ! Component: 100 - total flux cleaned = 0.00508268 Jy ! Component: 150 - total flux cleaned = 0.00566776 Jy ! Component: 200 - total flux cleaned = 0.00615108 Jy ! Component: 250 - total flux cleaned = 0.00640128 Jy ! Component: 300 - total flux cleaned = 0.00669711 Jy ! Component: 350 - total flux cleaned = 0.00695472 Jy ! Component: 400 - total flux cleaned = 0.00718518 Jy ! Total flux subtracted in 400 components = 0.00718518 Jy ! Clean residual min=-0.005726 max=0.005146 Jy/beam ! Clean residual mean=0.000194 rms=0.001809 Jy/beam ! Combined flux in latest and established models = 1.31265 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.0002157 Jy ! Component: 100 - total flux cleaned = 0.000416042 Jy ! Component: 150 - total flux cleaned = 0.000559438 Jy ! Component: 200 - total flux cleaned = 0.000652843 Jy ! Component: 250 - total flux cleaned = 0.000785424 Jy ! Component: 300 - total flux cleaned = 0.000871382 Jy ! Component: 350 - total flux cleaned = 0.000954794 Jy ! Component: 400 - total flux cleaned = 0.000997309 Jy ! Total flux subtracted in 400 components = 0.000997309 Jy ! Clean residual min=-0.005707 max=0.005171 Jy/beam ! Clean residual mean=0.000194 rms=0.001791 Jy/beam ! Combined flux in latest and established models = 1.31364 Jy ! Adding 55 model components to the UV plane model. ! The established model now contains 183 components and 1.31364 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.8919 mas, bmaj=3.279 mas, bpa=4.37 degrees ! Estimated noise=1.09277 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000329103 Jy ! Component: 100 - total flux cleaned = -0.000462423 Jy ! Component: 150 - total flux cleaned = -0.000343224 Jy ! Component: 200 - total flux cleaned = -0.000343931 Jy ! Component: 250 - total flux cleaned = -0.0003417 Jy ! Component: 300 - total flux cleaned = -0.000339829 Jy ! Component: 350 - total flux cleaned = -0.000242219 Jy ! Component: 400 - total flux cleaned = -0.000144521 Jy ! Total flux subtracted in 400 components = -0.000144521 Jy ! Clean residual min=-0.005063 max=0.005414 Jy/beam ! Clean residual mean=0.000106 rms=0.001313 Jy/beam ! Combined flux in latest and established models = 1.3135 Jy selfcal ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 208 components and 1.3135 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119458Jy sigma=1.013376 ! Fit after self-cal, rms=0.119452Jy sigma=1.011128 wmodel /scr/pima/rdv93_uvs/J0922-3959_X_map.mod ! Writing 208 model components to file: /scr/pima/rdv93_uvs/J0922-3959_X_map.mod wobs /scr/pima/rdv93_uvs/J0922-3959_X_uvs.fits ! Applying 2070 buffered edits. ! Writing UV FITS file: /scr/pima/rdv93_uvs/J0922-3959_X_uvs.fits wwins /scr/pima/rdv93_uvs/J0922-3959_X_map.win ! wwins: Wrote 5 windows to /scr/pima/rdv93_uvs/J0922-3959_X_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000387352 Jy ! Component: 100 - total flux cleaned = 0.000556668 Jy ! Component: 150 - total flux cleaned = 0.000958101 Jy ! Component: 200 - total flux cleaned = 0.00141954 Jy ! Component: 250 - total flux cleaned = 0.0020107 Jy ! Component: 300 - total flux cleaned = 0.00279366 Jy ! Component: 350 - total flux cleaned = 0.00320717 Jy ! Component: 400 - total flux cleaned = 0.00387874 Jy ! Component: 450 - total flux cleaned = 0.00433559 Jy ! Component: 500 - total flux cleaned = 0.00503743 Jy ! Component: 550 - total flux cleaned = 0.00566016 Jy ! Component: 600 - total flux cleaned = 0.00645061 Jy ! Component: 650 - total flux cleaned = 0.0070429 Jy ! Component: 700 - total flux cleaned = 0.00739139 Jy ! Component: 750 - total flux cleaned = 0.00784871 Jy ! Component: 800 - total flux cleaned = 0.00824121 Jy ! Component: 850 - total flux cleaned = 0.00857195 Jy ! Component: 900 - total flux cleaned = 0.00889727 Jy ! Component: 950 - total flux cleaned = 0.0091105 Jy ! Component: 1000 - total flux cleaned = 0.00947685 Jy ! Component: 1050 - total flux cleaned = 0.00957976 Jy ! Component: 1100 - total flux cleaned = 0.00993552 Jy ! Component: 1150 - total flux cleaned = 0.010286 Jy ! Component: 1200 - total flux cleaned = 0.0104834 Jy ! Component: 1250 - total flux cleaned = 0.0108246 Jy ! Component: 1300 - total flux cleaned = 0.0111595 Jy ! Component: 1350 - total flux cleaned = 0.0114902 Jy ! Component: 1400 - total flux cleaned = 0.0116765 Jy ! Component: 1450 - total flux cleaned = 0.0118143 Jy ! Component: 1500 - total flux cleaned = 0.0120416 Jy ! Component: 1550 - total flux cleaned = 0.0122661 Jy ! Component: 1600 - total flux cleaned = 0.0124875 Jy ! Component: 1650 - total flux cleaned = 0.0125315 Jy ! Component: 1700 - total flux cleaned = 0.0127056 Jy ! Component: 1750 - total flux cleaned = 0.0128774 Jy ! Component: 1800 - total flux cleaned = 0.0131748 Jy ! Component: 1850 - total flux cleaned = 0.0134266 Jy ! Component: 1900 - total flux cleaned = 0.0135929 Jy ! Component: 1950 - total flux cleaned = 0.0138393 Jy ! Component: 2000 - total flux cleaned = 0.0140429 Jy ! Component: 2050 - total flux cleaned = 0.0142844 Jy ! Component: 2100 - total flux cleaned = 0.0144438 Jy ! Component: 2150 - total flux cleaned = 0.0146017 Jy ! Component: 2200 - total flux cleaned = 0.0147972 Jy ! Component: 2250 - total flux cleaned = 0.0148364 Jy ! Component: 2300 - total flux cleaned = 0.0149909 Jy ! Component: 2350 - total flux cleaned = 0.0152202 Jy ! Component: 2400 - total flux cleaned = 0.0153338 Jy ! Component: 2450 - total flux cleaned = 0.0153339 Jy ! Component: 2500 - total flux cleaned = 0.015446 Jy ! Component: 2550 - total flux cleaned = 0.0156679 Jy ! Component: 2600 - total flux cleaned = 0.0157417 Jy ! Component: 2650 - total flux cleaned = 0.0159598 Jy ! Component: 2700 - total flux cleaned = 0.015996 Jy ! Component: 2750 - total flux cleaned = 0.0162109 Jy ! Component: 2800 - total flux cleaned = 0.016353 Jy ! Component: 2850 - total flux cleaned = 0.0164935 Jy ! Component: 2900 - total flux cleaned = 0.0165984 Jy ! Component: 2950 - total flux cleaned = 0.0166333 Jy ! Component: 3000 - total flux cleaned = 0.016771 Jy ! Component: 3050 - total flux cleaned = 0.0168738 Jy ! Component: 3100 - total flux cleaned = 0.0169079 Jy ! Component: 3150 - total flux cleaned = 0.0170434 Jy ! Component: 3200 - total flux cleaned = 0.0171775 Jy ! Component: 3250 - total flux cleaned = 0.0172438 Jy ! Component: 3300 - total flux cleaned = 0.0173099 Jy ! Component: 3350 - total flux cleaned = 0.0174089 Jy ! Component: 3400 - total flux cleaned = 0.017507 Jy ! Component: 3450 - total flux cleaned = 0.0176043 Jy ! Component: 3500 - total flux cleaned = 0.0175721 Jy ! Component: 3550 - total flux cleaned = 0.0176685 Jy ! Component: 3600 - total flux cleaned = 0.0177639 Jy ! Component: 3650 - total flux cleaned = 0.0178586 Jy ! Component: 3700 - total flux cleaned = 0.0179213 Jy ! Component: 3750 - total flux cleaned = 0.0178899 Jy ! Component: 3800 - total flux cleaned = 0.0179832 Jy ! Component: 3850 - total flux cleaned = 0.0179829 Jy ! Component: 3900 - total flux cleaned = 0.0180749 Jy ! Component: 3950 - total flux cleaned = 0.0180138 Jy ! Component: 4000 - total flux cleaned = 0.0179832 Jy ! Component: 4050 - total flux cleaned = 0.0180134 Jy ! Total flux subtracted in 4096 components = 0.0180736 Jy ! Clean residual min=-0.001619 max=0.001700 Jy/beam ! Clean residual mean=0.000030 rms=0.000530 Jy/beam ! Combined flux in latest and established models = 1.33157 Jy keep ! Adding 251 model components to the UV plane model. ! The established model now contains 459 components and 1.33157 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv93_uvs/J0922-3959_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.8919 x 3.279 at 4.37 degrees (North through East) ! Clean map min=-0.0061331 max=0.96756 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv93_uvs/J0922-3959_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.9675 rms= 0.00051 Jy quit ! Quitting program ! Log file /scr/pima/rdv93_uvs/J0922-3959_X_uvs.log closed on Tue Oct 8 15:45:40 2013