! Started logfile: /scr/pima/rdv62_uvs/J1407+2827_X_uvs.log on Fri Nov 22 22:38:51 2013 obs /scr/pima/rdv62_uvs/J1407+2827_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J1407+2827_X_uva.fits ! AN table 1: 2 integrations on 10 of 10 possible baselines. ! AN table 2: 13 integrations on 36 of 36 possible baselines. ! AN table 3: 9 integrations on 45 of 45 possible baselines. ! AN table 4: 4 integrations on 28 of 28 possible baselines. ! AN table 5: 2 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.605479 visibilities/baseline/integration-bin. ! Found source: J1407+2827 ! ! 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 8e+06 1 8e+06 ! 02 2 8.47599e+09 8e+06 1 8e+06 ! 03 3 8.79099e+09 8e+06 1 8e+06 ! 04 4 8.89599e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2652 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/rdv62_uvs/J1407+2827_X] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN 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:NL 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 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 41 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 41 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.457254Jy sigma=19.802323 ! Fit after self-cal, rms=0.460540Jy sigma=4.584705 ! 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 2:NY ! Flagging baselines of antenna 2:NY in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 4:KK ! Flagging baselines of antenna 4: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=1.303 mas, bmaj=2.028 mas, bpa=31.96 degrees ! Estimated noise=2.17955 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.801711 Jy ! Component: 100 - total flux cleaned = 1.05217 Jy ! Total flux subtracted in 100 components = 1.05217 Jy ! Clean residual min=-0.048368 max=0.111630 Jy/beam ! Clean residual mean=0.000398 rms=0.013049 Jy/beam ! Combined flux in latest and established models = 1.05218 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 22 components and 1.05218 Jy ! ! Correcting IF 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 9 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.297258Jy sigma=4.570016 ! Fit after self-cal, rms=0.296936Jy sigma=4.568412 ! Inverting map and beam ! Estimated beam: bmin=1.303 mas, bmaj=2.028 mas, bpa=31.96 degrees ! Estimated noise=2.17955 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.130399 Jy ! Component: 100 - total flux cleaned = 0.209042 Jy ! Total flux subtracted in 100 components = 0.209042 Jy ! Clean residual min=-0.030183 max=0.061884 Jy/beam ! Clean residual mean=0.000178 rms=0.007907 Jy/beam ! Combined flux in latest and established models = 1.26122 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 40 components and 1.26122 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.199943Jy sigma=2.917570 ! Fit after self-cal, rms=0.199356Jy sigma=2.914958 ! Inverting map ! Added new window around map position (-5.2, -3.8). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0617077 Jy ! Component: 100 - total flux cleaned = 0.104235 Jy ! Total flux subtracted in 100 components = 0.104235 Jy ! Clean residual min=-0.023901 max=0.030253 Jy/beam ! Clean residual mean=0.000090 rms=0.006235 Jy/beam ! Combined flux in latest and established models = 1.36545 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 55 components and 1.36545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168144Jy sigma=2.380524 ! Fit after self-cal, rms=0.166794Jy sigma=2.367542 ! Inverting map print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.819 mas, bmaj=2.527 mas, bpa=47.8 degrees ! Estimated noise=1.23837 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.0697758 Jy ! Component: 100 - total flux cleaned = 0.112866 Jy ! Component: 150 - total flux cleaned = 0.14167 Jy ! Component: 200 - total flux cleaned = 0.161386 Jy ! Total flux subtracted in 200 components = 0.161386 Jy ! Clean residual min=-0.022971 max=0.029694 Jy/beam ! Clean residual mean=0.000032 rms=0.006404 Jy/beam ! Combined flux in latest and established models = 1.52684 Jy ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 86 components and 1.52684 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.138378Jy sigma=1.858558 ! Fit after self-cal, rms=0.136984Jy sigma=1.843185 ! Inverting map print "=== Finished natural weighting clean for VLBI sites ===" ! === Finished natural weighting clean for VLBI sites === gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.04 LA 0.98 NL 1.02 ! PT 1.02 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.95 FD 0.98 HN 0.98 KP 0.96 ! LA 0.94 NL 0.97 NY 1.00* PT 0.96 ! SC 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.96 FD 0.96 HN 0.97 KP 0.96 ! LA 0.92 NL 0.97 NY 1.00* OV 1.00 ! PT 0.95 SC 1.01 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.95 FD 0.97 KK 1.00* KP 1.01 ! LA 0.94 MK 0.94 OV 0.97 PT 0.95 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.98 FD 1.08 HN 1.01 KP 1.04 ! LA 0.97 MK 1.09 NL 0.99 OV 1.00 ! PT 0.97 SC 1.16 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.06 FD 1.11 LA 0.94 NL 1.01 ! PT 1.08 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 0.96 FD 0.98 HN 1.02 KP 0.91 ! LA 0.95 NL 0.96 NY 1.00* PT 1.01 ! SC 0.97 ! ! Telescope amplitude corrections in sub-array 3: ! BR 0.96 FD 0.96 HN 0.97 KP 0.92 ! LA 0.93 NL 0.94 NY 1.00* OV 0.98 ! PT 1.01 SC 1.02 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 0.94 FD 0.98 KK 1.00* KP 0.94 ! LA 0.94 MK 0.99 OV 0.96 PT 1.02 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 0.96 FD 1.08 HN 1.08 KP 1.00 ! LA 0.98 MK 1.09 NL 0.98 OV 0.99 ! PT 1.03 SC 1.09 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.10 FD 1.14 LA 1.09 NL 1.08 ! PT 1.14 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.05 HN 1.02 KP 0.87 ! LA 1.03 NL 1.02 NY 1.00* PT 1.08 ! SC 1.05 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.07 FD 1.04 HN 1.01 KP 0.89 ! LA 1.02 NL 1.00 NY 1.00* OV 1.09 ! PT 1.05 SC 0.98 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.01 FD 1.06 KK 1.00* KP 0.92 ! LA 1.03 MK 0.99 OV 1.05 PT 1.06 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.05 FD 1.17 HN 1.08 KP 0.94 ! LA 1.07 MK 1.15 NL 1.04 OV 1.08 ! PT 1.10 SC 1.05 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.16 LA 1.05 NL 1.10 ! PT 1.20 ! ! A total of 6 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.01 FD 1.07 HN 1.09 KP 0.94 ! LA 1.02 NL 1.06 NY 1.00* PT 1.07 ! SC 1.04 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.03 FD 1.07 HN 1.08 KP 0.95 ! LA 1.00 NL 1.04 NY 1.00* OV 1.06 ! PT 1.05 SC 1.03 ! ! A total of 3 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.02 FD 1.08 KK 1.00* KP 0.98 ! LA 1.01 MK 1.05 OV 1.04 PT 1.07 ! ! A total of 3 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 FD 1.20 HN 1.17 KP 1.02 ! LA 1.05 MK 1.17 NL 1.06 OV 1.07 ! PT 1.07 SC 1.15 ! ! ! Fit before self-cal, rms=0.136984Jy sigma=1.843185 ! Fit after self-cal, rms=0.090862Jy sigma=1.044802 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.839 mas, bmaj=2.568 mas, bpa=48.81 degrees ! Estimated noise=1.24777 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 (-7, -3.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0154187 Jy ! Component: 100 - total flux cleaned = 0.0217326 Jy ! Component: 150 - total flux cleaned = 0.0229178 Jy ! Component: 200 - total flux cleaned = 0.0236069 Jy ! Total flux subtracted in 200 components = 0.0236069 Jy ! Clean residual min=-0.008922 max=0.008445 Jy/beam ! Clean residual mean=-0.000030 rms=0.001680 Jy/beam ! Combined flux in latest and established models = 1.55044 Jy ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 108 components and 1.55044 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.088056Jy sigma=1.004557 ! Fit after self-cal, rms=0.087839Jy sigma=1.000030 ! 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.087839Jy sigma=1.000030 ! Fit after self-cal, rms=0.087676Jy sigma=0.997089 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.84 mas, bmaj=2.569 mas, bpa=48.78 degrees ! Estimated noise=1.2473 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00486906 Jy ! Component: 100 - total flux cleaned = 0.00587088 Jy ! Component: 150 - total flux cleaned = 0.00601265 Jy ! Component: 200 - total flux cleaned = 0.00614179 Jy ! Total flux subtracted in 200 components = 0.00614179 Jy ! Clean residual min=-0.008627 max=0.007597 Jy/beam ! Clean residual mean=-0.000017 rms=0.001449 Jy/beam ! Combined flux in latest and established models = 1.55659 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 121 components and 1.55659 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.086768Jy sigma=0.989163 ! Fit after self-cal, rms=0.086685Jy sigma=0.988189 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.00165232 Jy ! Component: 100 - total flux cleaned = 0.00230546 Jy ! Component: 150 - total flux cleaned = 0.00242105 Jy ! Component: 200 - total flux cleaned = 0.00253229 Jy ! Total flux subtracted in 200 components = 0.00253229 Jy ! Clean residual min=-0.008225 max=0.006847 Jy/beam ! Clean residual mean=-0.000018 rms=0.001389 Jy/beam ! Combined flux in latest and established models = 1.55912 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 128 components and 1.55912 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.086024Jy sigma=0.985042 ! Fit after self-cal, rms=0.085976Jy sigma=0.984701 ! 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.085976Jy sigma=0.984701 ! Fit after self-cal, rms=0.085947Jy sigma=0.983960 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.841 mas, bmaj=2.569 mas, bpa=48.78 degrees ! Estimated noise=1.24685 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000723168 Jy ! Component: 100 - total flux cleaned = 0.000961797 Jy ! Component: 150 - total flux cleaned = 0.00106869 Jy ! Component: 200 - total flux cleaned = 0.0011731 Jy ! Total flux subtracted in 200 components = 0.0011731 Jy ! Clean residual min=-0.007961 max=0.006713 Jy/beam ! Clean residual mean=-0.000006 rms=0.001344 Jy/beam ! Combined flux in latest and established models = 1.56029 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 132 components and 1.56029 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.085364Jy sigma=0.981270 ! Fit after self-cal, rms=0.085319Jy sigma=0.980990 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.000379897 Jy ! Component: 100 - total flux cleaned = 0.000704081 Jy ! Component: 150 - total flux cleaned = 0.000902302 Jy ! Component: 200 - total flux cleaned = 0.000998602 Jy ! Total flux subtracted in 200 components = 0.000998602 Jy ! Clean residual min=-0.007634 max=0.006389 Jy/beam ! Clean residual mean=-0.000006 rms=0.001318 Jy/beam ! Combined flux in latest and established models = 1.56129 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 138 components and 1.56129 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084836Jy sigma=0.978990 ! Fit after self-cal, rms=0.084799Jy sigma=0.978783 ! 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.084799Jy sigma=0.978783 ! Fit after self-cal, rms=0.084772Jy sigma=0.978618 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.157 mas, bmaj=3.446 mas, bpa=0.6839 degrees ! Estimated noise=1.21755 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000355694 Jy ! Component: 100 - total flux cleaned = 0.0010699 Jy ! Component: 150 - total flux cleaned = 0.00165742 Jy ! Component: 200 - total flux cleaned = 0.00196763 Jy ! Total flux subtracted in 200 components = 0.00196763 Jy ! Clean residual min=-0.005253 max=0.004530 Jy/beam ! Clean residual mean=-0.000002 rms=0.001301 Jy/beam ! Combined flux in latest and established models = 1.56326 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 147 components and 1.56326 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084785Jy sigma=0.976324 ! Fit after self-cal, rms=0.084612Jy sigma=0.975809 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 = 5.52537e-05 Jy ! Component: 100 - total flux cleaned = 0.000384129 Jy ! Component: 150 - total flux cleaned = 0.000553956 Jy ! Component: 200 - total flux cleaned = 0.00070636 Jy ! Total flux subtracted in 200 components = 0.00070636 Jy ! Clean residual min=-0.005161 max=0.004637 Jy/beam ! Clean residual mean=-0.000002 rms=0.001281 Jy/beam ! Combined flux in latest and established models = 1.56396 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 152 components and 1.56396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.084640Jy sigma=0.974401 ! Fit after self-cal, rms=0.084496Jy sigma=0.974006 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1: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 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 selfant 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3: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:BR,true ! 4:BR self-cal status=fixed weight=1 selfant 4:FD,true ! 4:FD self-cal status=fixed weight=1 selfant 4:KP,true ! 4:KP self-cal status=fixed weight=1 selfant 4:LA,true ! 4:LA self-cal status=fixed weight=1 selfant 4:MK,true ! 4:MK self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 4:PT,true ! 4:PT 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:HN,true ! 5:HN 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:NL,true ! 5:NL 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 ![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 4 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 2. ! A total of 4 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 3. ! A total of 4 telescope corrections were flagged in sub-array 2. ! ! Correcting IF 4. ! A total of 4 telescope corrections were flagged in sub-array 2. ! ! Fit before self-cal, rms=0.304219Jy sigma=1.153572 ! Fit after self-cal, rms=0.226336Jy sigma=0.900655 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD 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 ! 2:BR self-cal status=correctable weight=1 ! 2:FD self-cal status=correctable weight=1 ! 2:HN self-cal status=correctable weight=1 ! 2:KP self-cal status=correctable weight=1 ! 2:LA self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:NY self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 3:PT self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 4:BR self-cal status=correctable weight=1 ! 4:FD self-cal status=correctable weight=1 ! 4:KK self-cal status=correctable weight=1 ! 4:KP self-cal status=correctable weight=1 ! 4:LA self-cal status=correctable weight=1 ! 4:MK self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 4:PT self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN 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:NL 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 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 40 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 40 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 40 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 40 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.226336Jy sigma=0.900655 ! Fit after self-cal, rms=0.227102Jy sigma=0.895795 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=1.299 mas, bmaj=1.831 mas, bpa=25.45 degrees ! Estimated noise=2.36277 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.00120513 Jy ! Component: 100 - total flux cleaned = 0.00201031 Jy ! Component: 150 - total flux cleaned = 0.00275059 Jy ! Component: 200 - total flux cleaned = 0.00321555 Jy ! Total flux subtracted in 200 components = 0.00321555 Jy ! Clean residual min=-0.013335 max=0.010237 Jy/beam ! Clean residual mean=-0.000003 rms=0.001803 Jy/beam ! Combined flux in latest and established models = 1.56718 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 162 components and 1.56718 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225327Jy sigma=0.895309 ! Fit after self-cal, rms=0.224578Jy sigma=0.893337 ! 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=1.79 mas, bmaj=2.108 mas, bpa=32.04 degrees ! Estimated noise=1.35279 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.000193117 Jy ! Component: 100 - total flux cleaned = -2.5941e-05 Jy ! Component: 150 - total flux cleaned = -0.00017258 Jy ! Component: 200 - total flux cleaned = -0.000439856 Jy ! Component: 250 - total flux cleaned = -0.000570577 Jy ! Component: 300 - total flux cleaned = -0.000694182 Jy ! Component: 350 - total flux cleaned = -0.000814433 Jy ! Component: 400 - total flux cleaned = -0.000820309 Jy ! Total flux subtracted in 400 components = -0.000820309 Jy ! Clean residual min=-0.007784 max=0.006122 Jy/beam ! Clean residual mean=-0.000001 rms=0.001285 Jy/beam ! Combined flux in latest and established models = 1.56636 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 173 components and 1.56636 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.221705Jy sigma=0.891107 ! Fit after self-cal, rms=0.221413Jy sigma=0.890591 ! 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.000362728 Jy ! Component: 100 - total flux cleaned = -0.000233231 Jy ! Component: 150 - total flux cleaned = -0.000351108 Jy ! Component: 200 - total flux cleaned = -0.00023431 Jy ! Component: 250 - total flux cleaned = -0.000234902 Jy ! Component: 300 - total flux cleaned = -0.00023154 Jy ! Component: 350 - total flux cleaned = -0.000231584 Jy ! Component: 400 - total flux cleaned = -0.000238713 Jy ! Total flux subtracted in 400 components = -0.000238713 Jy ! Clean residual min=-0.007238 max=0.006006 Jy/beam ! Clean residual mean=-0.000000 rms=0.001269 Jy/beam ! Combined flux in latest and established models = 1.56612 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 175 components and 1.56612 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.219194Jy sigma=0.889353 ! Fit after self-cal, rms=0.218972Jy sigma=0.888974 ! 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 44 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 44 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 44 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 44 telescope corrections were flagged in sub-array 2. ! A total of 16 telescope corrections were flagged in sub-array 3. ! A total of 4 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.225003Jy sigma=0.893939 ! Fit after self-cal, rms=0.213107Jy sigma=0.888650 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.786 mas, bmaj=2.098 mas, bpa=31.4 degrees ! Estimated noise=1.35209 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000577343 Jy ! Component: 100 - total flux cleaned = 0.000583501 Jy ! Component: 150 - total flux cleaned = 0.000567632 Jy ! Component: 200 - total flux cleaned = 0.000445427 Jy ! Component: 250 - total flux cleaned = 0.000341208 Jy ! Component: 300 - total flux cleaned = 0.000237669 Jy ! Component: 350 - total flux cleaned = 0.000231584 Jy ! Component: 400 - total flux cleaned = 0.000231064 Jy ! Total flux subtracted in 400 components = 0.000231064 Jy ! Clean residual min=-0.006948 max=0.005930 Jy/beam ! Clean residual mean=0.000004 rms=0.001254 Jy/beam ! Combined flux in latest and established models = 1.56635 Jy selfcal ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 177 components and 1.56635 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.210711Jy sigma=0.887777 ! Fit after self-cal, rms=0.210489Jy sigma=0.887478 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.00015727 Jy ! Component: 100 - total flux cleaned = -0.000165804 Jy ! Component: 150 - total flux cleaned = -0.000156056 Jy ! Component: 200 - total flux cleaned = -5.43045e-05 Jy ! Component: 250 - total flux cleaned = 3.80293e-05 Jy ! Component: 300 - total flux cleaned = 3.29705e-05 Jy ! Component: 350 - total flux cleaned = 2.83255e-05 Jy ! Component: 400 - total flux cleaned = 0.000111916 Jy ! Total flux subtracted in 400 components = 0.000111916 Jy ! Clean residual min=-0.006642 max=0.005686 Jy/beam ! Clean residual mean=0.000003 rms=0.001240 Jy/beam ! Combined flux in latest and established models = 1.56646 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 179 components and 1.56646 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.208705Jy sigma=0.886589 ! Fit after self-cal, rms=0.208559Jy sigma=0.886351 ! 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.208559Jy sigma=0.886351 ! Fit after self-cal, rms=0.204478Jy sigma=0.886213 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.783 mas, bmaj=2.09 mas, bpa=30.94 degrees ! Estimated noise=1.35144 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000310716 Jy ! Component: 100 - total flux cleaned = 0.000424611 Jy ! Component: 150 - total flux cleaned = 0.000411653 Jy ! Component: 200 - total flux cleaned = 0.00031238 Jy ! Component: 250 - total flux cleaned = 0.000308652 Jy ! Component: 300 - total flux cleaned = 0.000307966 Jy ! Component: 350 - total flux cleaned = 0.000311434 Jy ! Component: 400 - total flux cleaned = 0.000220524 Jy ! Total flux subtracted in 400 components = 0.000220524 Jy ! Clean residual min=-0.006421 max=0.005576 Jy/beam ! Clean residual mean=0.000004 rms=0.001230 Jy/beam ! Combined flux in latest and established models = 1.56668 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 181 components and 1.56668 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.202575Jy sigma=0.885531 ! Fit after self-cal, rms=0.202425Jy sigma=0.885316 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.00013043 Jy ! Component: 100 - total flux cleaned = -3.01368e-05 Jy ! Component: 150 - total flux cleaned = -3.13987e-05 Jy ! Component: 200 - total flux cleaned = -2.80382e-05 Jy ! Component: 250 - total flux cleaned = 0.000136189 Jy ! Component: 300 - total flux cleaned = 5.42718e-05 Jy ! Component: 350 - total flux cleaned = 0.000138412 Jy ! Component: 400 - total flux cleaned = 0.000130895 Jy ! Total flux subtracted in 400 components = 0.000130895 Jy ! Clean residual min=-0.006176 max=0.005300 Jy/beam ! Clean residual mean=0.000004 rms=0.001218 Jy/beam ! Combined flux in latest and established models = 1.56681 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 185 components and 1.56681 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.200945Jy sigma=0.884524 ! Fit after self-cal, rms=0.200840Jy sigma=0.884333 ! 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.200840Jy sigma=0.884333 ! Fit after self-cal, rms=0.197661Jy sigma=0.863682 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.066 mas, bmaj=3.443 mas, bpa=5.11 degrees ! Estimated noise=1.2108 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00130275 Jy ! Component: 100 - total flux cleaned = 0.00141202 Jy ! Component: 150 - total flux cleaned = 0.00141482 Jy ! Component: 200 - total flux cleaned = 0.00141552 Jy ! Component: 250 - total flux cleaned = 0.00133694 Jy ! Component: 300 - total flux cleaned = 0.00130303 Jy ! Component: 350 - total flux cleaned = 0.00130545 Jy ! Component: 400 - total flux cleaned = 0.00133818 Jy ! Total flux subtracted in 400 components = 0.00133818 Jy ! Clean residual min=-0.004800 max=0.004137 Jy/beam ! Clean residual mean=0.000004 rms=0.001248 Jy/beam ! Combined flux in latest and established models = 1.56815 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 195 components and 1.56815 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.197952Jy sigma=0.862459 ! Fit after self-cal, rms=0.197884Jy sigma=0.862230 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 = -4.16471e-05 Jy ! Component: 100 - total flux cleaned = 0.000163441 Jy ! Component: 150 - total flux cleaned = 0.000163443 Jy ! Component: 200 - total flux cleaned = 0.00016355 Jy ! Component: 250 - total flux cleaned = 0.000194019 Jy ! Component: 300 - total flux cleaned = 0.000194219 Jy ! Component: 350 - total flux cleaned = 0.000248086 Jy ! Component: 400 - total flux cleaned = 0.000248638 Jy ! Total flux subtracted in 400 components = 0.000248638 Jy ! Clean residual min=-0.004761 max=0.004198 Jy/beam ! Clean residual mean=0.000004 rms=0.001238 Jy/beam ! Combined flux in latest and established models = 1.5684 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 200 components and 1.5684 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.198032Jy sigma=0.861502 ! Fit after self-cal, rms=0.197957Jy sigma=0.861307 ! 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=1.781 mas, bmaj=2.084 mas, bpa=30.6 degrees ! Estimated noise=1.35194 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.781 x 2.084 at 30.6 degrees (North through East) ! Clean map min=-0.0048786 max=1.0685 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.197957Jy sigma=0.861307 ! Fit after self-cal, rms=0.208936Jy sigma=0.776725 ! Inverting map and beam ! Estimated beam: bmin=1.782 mas, bmaj=2.087 mas, bpa=30.69 degrees ! Estimated noise=1.35287 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -4.53644e-05 Jy ! Component: 100 - total flux cleaned = 5.32892e-05 Jy ! Component: 150 - total flux cleaned = 4.69948e-05 Jy ! Component: 200 - total flux cleaned = 5.13146e-05 Jy ! Component: 250 - total flux cleaned = 5.1157e-05 Jy ! Component: 300 - total flux cleaned = 4.92954e-05 Jy ! Component: 350 - total flux cleaned = 0.000123957 Jy ! Component: 400 - total flux cleaned = 0.00012301 Jy ! Total flux subtracted in 400 components = 0.00012301 Jy ! Clean residual min=-0.005618 max=0.004888 Jy/beam ! Clean residual mean=0.000012 rms=0.001184 Jy/beam ! Combined flux in latest and established models = 1.56852 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 201 components and 1.56852 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.207389Jy sigma=0.776119 ! Fit after self-cal, rms=0.207325Jy sigma=0.775958 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000110036 Jy ! Component: 100 - total flux cleaned = -2.49103e-05 Jy ! Component: 150 - total flux cleaned = 0.000137742 Jy ! Component: 200 - total flux cleaned = 0.000296166 Jy ! Component: 250 - total flux cleaned = 0.000369766 Jy ! Component: 300 - total flux cleaned = 0.000441616 Jy ! Component: 350 - total flux cleaned = 0.000438902 Jy ! Component: 400 - total flux cleaned = 0.000435272 Jy ! Total flux subtracted in 400 components = 0.000435272 Jy ! Clean residual min=-0.005542 max=0.004757 Jy/beam ! Clean residual mean=0.000012 rms=0.001178 Jy/beam ! Combined flux in latest and established models = 1.56896 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 203 components and 1.56896 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.206108Jy sigma=0.775414 ! Fit after self-cal, rms=0.206046Jy sigma=0.775278 ! 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.206046Jy sigma=0.775278 ! Fit after self-cal, rms=0.203135Jy sigma=0.775076 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=1.294 mas, bmaj=1.816 mas, bpa=24.93 degrees ! Estimated noise=2.3598 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.824241 Jy ! Component: 100 - total flux cleaned = 1.08128 Jy ! Component: 150 - total flux cleaned = 1.2108 Jy ! Component: 200 - total flux cleaned = 1.29416 Jy ! Component: 250 - total flux cleaned = 1.35257 Jy ! Component: 300 - total flux cleaned = 1.39697 Jy ! Component: 350 - total flux cleaned = 1.43149 Jy ! Component: 400 - total flux cleaned = 1.4587 Jy ! Total flux subtracted in 400 components = 1.4587 Jy ! Clean residual min=-0.012331 max=0.018522 Jy/beam ! Clean residual mean=0.000084 rms=0.002838 Jy/beam ! Combined flux in latest and established models = 1.45871 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 57 model components to the UV plane model. ! The established model now contains 57 components and 1.45871 Jy ! Inverting map and beam ! Estimated beam: bmin=1.781 mas, bmaj=2.082 mas, bpa=30.38 degrees ! Estimated noise=1.35281 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.032678 Jy ! Component: 100 - total flux cleaned = 0.0541318 Jy ! Component: 150 - total flux cleaned = 0.0693076 Jy ! Component: 200 - total flux cleaned = 0.0804119 Jy ! Component: 250 - total flux cleaned = 0.0886371 Jy ! Component: 300 - total flux cleaned = 0.0937741 Jy ! Component: 350 - total flux cleaned = 0.0957715 Jy ! Component: 400 - total flux cleaned = 0.0973936 Jy ! Total flux subtracted in 400 components = 0.0973936 Jy ! Clean residual min=-0.007206 max=0.008234 Jy/beam ! Clean residual mean=0.000016 rms=0.001400 Jy/beam ! Combined flux in latest and established models = 1.5561 Jy ! Adding 41 model components to the UV plane model. ! The established model now contains 93 components and 1.5561 Jy ! Inverting map ! Added new window around map position (-0.6, 3). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00576031 Jy ! Component: 100 - total flux cleaned = 0.00671076 Jy ! Component: 150 - total flux cleaned = 0.00706713 Jy ! Component: 200 - total flux cleaned = 0.0072421 Jy ! Component: 250 - total flux cleaned = 0.0075765 Jy ! Component: 300 - total flux cleaned = 0.00773214 Jy ! Component: 350 - total flux cleaned = 0.00773007 Jy ! Component: 400 - total flux cleaned = 0.00803629 Jy ! Total flux subtracted in 400 components = 0.00803629 Jy ! Clean residual min=-0.006205 max=0.005152 Jy/beam ! Clean residual mean=0.000011 rms=0.001215 Jy/beam ! Combined flux in latest and established models = 1.56413 Jy ! Adding 36 model components to the UV plane model. ! The established model now contains 116 components and 1.56413 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=3.067 mas, bmaj=3.446 mas, bpa=5.129 degrees ! Estimated noise=1.21129 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0019751 Jy ! Component: 100 - total flux cleaned = 0.00252008 Jy ! Component: 150 - total flux cleaned = 0.00290453 Jy ! Component: 200 - total flux cleaned = 0.00304441 Jy ! Component: 250 - total flux cleaned = 0.0032338 Jy ! Component: 300 - total flux cleaned = 0.00340989 Jy ! Component: 350 - total flux cleaned = 0.00346419 Jy ! Component: 400 - total flux cleaned = 0.00346342 Jy ! Total flux subtracted in 400 components = 0.00346342 Jy ! Clean residual min=-0.004845 max=0.004076 Jy/beam ! Clean residual mean=0.000014 rms=0.001180 Jy/beam ! Combined flux in latest and established models = 1.5676 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.000146042 Jy ! Component: 100 - total flux cleaned = 0.000145616 Jy ! Component: 150 - total flux cleaned = 0.000145644 Jy ! Component: 200 - total flux cleaned = 0.00023419 Jy ! Component: 250 - total flux cleaned = 0.000276573 Jy ! Component: 300 - total flux cleaned = 0.000357865 Jy ! Component: 350 - total flux cleaned = 0.000357031 Jy ! Component: 400 - total flux cleaned = 0.00043372 Jy ! Total flux subtracted in 400 components = 0.00043372 Jy ! Clean residual min=-0.004824 max=0.004036 Jy/beam ! Clean residual mean=0.000013 rms=0.001165 Jy/beam ! Combined flux in latest and established models = 1.56803 Jy ! Adding 49 model components to the UV plane model. ! The established model now contains 161 components and 1.56803 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=1.781 mas, bmaj=2.082 mas, bpa=30.38 degrees ! Estimated noise=1.35281 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -6.54609e-06 Jy ! Component: 100 - total flux cleaned = -6.08274e-06 Jy ! Component: 150 - total flux cleaned = -0.000139254 Jy ! Component: 200 - total flux cleaned = -0.000267875 Jy ! Component: 250 - total flux cleaned = -0.000269488 Jy ! Component: 300 - total flux cleaned = -0.000518204 Jy ! Component: 350 - total flux cleaned = -0.00051281 Jy ! Component: 400 - total flux cleaned = -0.000509135 Jy ! Total flux subtracted in 400 components = -0.000509135 Jy ! Clean residual min=-0.005368 max=0.004537 Jy/beam ! Clean residual mean=0.000008 rms=0.001133 Jy/beam ! Combined flux in latest and established models = 1.56752 Jy selfcal ! Performing phase self-cal ! Adding 26 model components to the UV plane model. ! The established model now contains 174 components and 1.56752 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.203894Jy sigma=0.776335 ! Fit after self-cal, rms=0.203281Jy sigma=0.774699 wmodel /scr/pima/rdv62_uvs/J1407+2827_X_map.mod ! Writing 174 model components to file: /scr/pima/rdv62_uvs/J1407+2827_X_map.mod wobs /scr/pima/rdv62_uvs/J1407+2827_X_uvs.fits ! Applying 2216 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J1407+2827_X_uvs.fits wwins /scr/pima/rdv62_uvs/J1407+2827_X_map.win ! wwins: Wrote 4 windows to /scr/pima/rdv62_uvs/J1407+2827_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.00242583 Jy ! Component: 100 - total flux cleaned = -0.00395632 Jy ! Component: 150 - total flux cleaned = -0.00483496 Jy ! Component: 200 - total flux cleaned = -0.00536925 Jy ! Component: 250 - total flux cleaned = -0.00595707 Jy ! Component: 300 - total flux cleaned = -0.00631122 Jy ! Component: 350 - total flux cleaned = -0.00644998 Jy ! Component: 400 - total flux cleaned = -0.00651706 Jy ! Component: 450 - total flux cleaned = -0.00638716 Jy ! Component: 500 - total flux cleaned = -0.00613363 Jy ! Component: 550 - total flux cleaned = -0.005885 Jy ! Component: 600 - total flux cleaned = -0.00557828 Jy ! Component: 650 - total flux cleaned = -0.00521481 Jy ! Component: 700 - total flux cleaned = -0.00479829 Jy ! Component: 750 - total flux cleaned = -0.00438811 Jy ! Component: 800 - total flux cleaned = -0.00415686 Jy ! Component: 850 - total flux cleaned = -0.00375705 Jy ! Component: 900 - total flux cleaned = -0.00336225 Jy ! Component: 950 - total flux cleaned = -0.00313902 Jy ! Component: 1000 - total flux cleaned = -0.0028094 Jy ! Component: 1050 - total flux cleaned = -0.00264657 Jy ! Component: 1100 - total flux cleaned = -0.00221593 Jy ! Component: 1150 - total flux cleaned = -0.00189571 Jy ! Component: 1200 - total flux cleaned = -0.00142166 Jy ! Component: 1250 - total flux cleaned = -0.00105674 Jy ! Component: 1300 - total flux cleaned = -0.00069546 Jy ! Component: 1350 - total flux cleaned = -0.000235136 Jy ! Component: 1400 - total flux cleaned = 0.00022219 Jy ! Component: 1450 - total flux cleaned = 0.000474718 Jy ! Component: 1500 - total flux cleaned = 0.00097483 Jy ! Component: 1550 - total flux cleaned = 0.00142054 Jy ! Component: 1600 - total flux cleaned = 0.00191238 Jy ! Component: 1650 - total flux cleaned = 0.00220577 Jy ! Component: 1700 - total flux cleaned = 0.00278756 Jy ! Component: 1750 - total flux cleaned = 0.00307626 Jy ! Component: 1800 - total flux cleaned = 0.00341111 Jy ! Component: 1850 - total flux cleaned = 0.00374355 Jy ! Component: 1900 - total flux cleaned = 0.0039801 Jy ! Component: 1950 - total flux cleaned = 0.00416776 Jy ! Component: 2000 - total flux cleaned = 0.00454067 Jy ! Component: 2050 - total flux cleaned = 0.00491106 Jy ! Component: 2100 - total flux cleaned = 0.00514128 Jy ! Component: 2150 - total flux cleaned = 0.00541584 Jy ! Component: 2200 - total flux cleaned = 0.00573452 Jy ! Component: 2250 - total flux cleaned = 0.00591538 Jy ! Component: 2300 - total flux cleaned = 0.00596029 Jy ! Component: 2350 - total flux cleaned = 0.00640754 Jy ! Component: 2400 - total flux cleaned = 0.00654068 Jy ! Component: 2450 - total flux cleaned = 0.00667343 Jy ! Component: 2500 - total flux cleaned = 0.00693769 Jy ! Component: 2550 - total flux cleaned = 0.00724415 Jy ! Component: 2600 - total flux cleaned = 0.00737473 Jy ! Component: 2650 - total flux cleaned = 0.00763488 Jy ! Component: 2700 - total flux cleaned = 0.00780768 Jy ! Component: 2750 - total flux cleaned = 0.00806514 Jy ! Component: 2800 - total flux cleaned = 0.00819314 Jy ! Component: 2850 - total flux cleaned = 0.00836336 Jy ! Component: 2900 - total flux cleaned = 0.0084056 Jy ! Component: 2950 - total flux cleaned = 0.00882657 Jy ! Component: 3000 - total flux cleaned = 0.00907783 Jy ! Component: 3050 - total flux cleaned = 0.00920262 Jy ! Component: 3100 - total flux cleaned = 0.00953444 Jy ! Component: 3150 - total flux cleaned = 0.00969965 Jy ! Component: 3200 - total flux cleaned = 0.00969959 Jy ! Component: 3250 - total flux cleaned = 0.0101093 Jy ! Component: 3300 - total flux cleaned = 0.0100685 Jy ! Component: 3350 - total flux cleaned = 0.0103125 Jy ! Component: 3400 - total flux cleaned = 0.0104743 Jy ! Component: 3450 - total flux cleaned = 0.0106358 Jy ! Component: 3500 - total flux cleaned = 0.0106356 Jy ! Component: 3550 - total flux cleaned = 0.0110361 Jy ! Component: 3600 - total flux cleaned = 0.0110363 Jy ! Component: 3650 - total flux cleaned = 0.0111557 Jy ! Component: 3700 - total flux cleaned = 0.0114328 Jy ! Component: 3750 - total flux cleaned = 0.0116693 Jy ! Component: 3800 - total flux cleaned = 0.0117088 Jy ! Component: 3850 - total flux cleaned = 0.0118653 Jy ! Component: 3900 - total flux cleaned = 0.0120214 Jy ! Component: 3950 - total flux cleaned = 0.0122159 Jy ! Component: 4000 - total flux cleaned = 0.0123325 Jy ! Component: 4050 - total flux cleaned = 0.0124098 Jy ! Total flux subtracted in 4096 components = 0.0125254 Jy ! Clean residual min=-0.001921 max=0.001951 Jy/beam ! Clean residual mean=0.000013 rms=0.000646 Jy/beam ! Combined flux in latest and established models = 1.58005 Jy keep ! Adding 239 model components to the UV plane model. ! The established model now contains 406 components and 1.58005 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J1407+2827_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 1.781 x 2.082 at 30.38 degrees (North through East) ! Clean map min=-0.004025 max=1.0746 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J1407+2827_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= 1.0745 rms= 0.00064 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J1407+2827_X_uvs.log closed on Fri Nov 22 22:38:57 2013