! Started logfile: /scr/pima/rdv95_uvs/J1543-0757_S_uvs.log on Wed Nov 6 20:47:38 2013 obs /scr/pima/rdv95_uvs/J1543-0757_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv95_uvs/J1543-0757_S_uva.fits ! AN table 1: 22 integrations on 36 of 36 possible baselines. ! AN table 2: 22 integrations on 45 of 45 possible baselines. ! Apparent sampling: 1 visibilities/baseline/integration-bin. ! Found source: J1543-0757 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 7128 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv95_uvs/J1543-0757_S] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 800000 float soltime2 soltime2 = 15 float soltime3 soltime3 = 1 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float old_peak float new_peak float flux_cutoff selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 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:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.206755Jy sigma=12.925842 ! Fit after self-cal, rms=0.611292Jy sigma=3.408824 ! 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 \] ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.6 mas, bmaj=9.487 mas, bpa=12.02 degrees ! Estimated noise=2.96428 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.4). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.579369 Jy ! Component: 100 - total flux cleaned = 0.790126 Jy ! Total flux subtracted in 100 components = 0.790126 Jy ! Clean residual min=-0.081375 max=0.232643 Jy/beam ! Clean residual mean=0.002043 rms=0.034283 Jy/beam ! Combined flux in latest and established models = 0.790126 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 33 components and 0.790126 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.295045Jy sigma=3.622586 ! Fit after self-cal, rms=0.295092Jy sigma=3.616246 ! Inverting map ! Added new window around map position (-7.6, -44). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.19667 Jy ! Component: 100 - total flux cleaned = 0.295556 Jy ! Total flux subtracted in 100 components = 0.295556 Jy ! Clean residual min=-0.055243 max=0.059300 Jy/beam ! Clean residual mean=0.001167 rms=0.017986 Jy/beam ! Combined flux in latest and established models = 1.08568 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 75 components and 1.08568 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.180058Jy sigma=2.150065 ! Fit after self-cal, rms=0.166558Jy sigma=2.055944 ! 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=6.285 mas, bmaj=11.75 mas, bpa=10.24 degrees ! Estimated noise=1.14189 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.0961646 Jy ! Component: 100 - total flux cleaned = 0.156867 Jy ! Component: 150 - total flux cleaned = 0.19917 Jy ! Component: 200 - total flux cleaned = 0.229913 Jy ! Total flux subtracted in 200 components = 0.229913 Jy ! Clean residual min=-0.023675 max=0.028234 Jy/beam ! Clean residual mean=0.000528 rms=0.008311 Jy/beam ! Combined flux in latest and established models = 1.31559 Jy ! Performing phase self-cal ! Adding 105 model components to the UV plane model. ! The established model now contains 176 components and 1.3156 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.124912Jy sigma=1.351355 ! Fit after self-cal, rms=0.121040Jy sigma=1.321041 ! 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 0.95 HN 1.07 KP 0.94 ! LA 1.14 NL 1.02 OV 0.94 PT 0.99 ! SC 0.95 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.99 FD 0.97 HN 1.08 KP 0.91 ! LA 1.13 MK 1.18 NL 0.96 OV 0.95 ! PT 0.99 SC 1.01 ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.01 FD 1.02 HN 1.09 KP 0.99 ! LA 1.04 NL 1.08 OV 0.99 PT 0.98 ! SC 1.01 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.02 FD 1.02 HN 1.13 KP 0.98 ! LA 1.04 MK 1.09 NL 1.00 OV 0.99 ! PT 0.99 SC 1.08 ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.11 FD 1.05 HN 1.06 KP 1.00 ! LA 0.96 NL 1.15 OV 1.04 PT 1.00 ! SC 0.98 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.07 FD 1.06 HN 1.11 KP 1.00 ! LA 0.97 MK 1.04 NL 1.05 OV 1.01 ! PT 1.00 SC 1.06 ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.04 FD 1.02 HN 0.92 KP 0.93 ! LA 0.82 NL 1.09 OV 0.86 PT 0.89 ! SC 1.05 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 1.10 HN 0.98 KP 0.93 ! LA 0.80 MK 1.17 NL 0.99 OV 0.86 ! PT 0.89 SC 1.09 ! ! ! Fit before self-cal, rms=0.121040Jy sigma=1.321041 ! Fit after self-cal, rms=0.087741Jy sigma=0.622198 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=6.422 mas, bmaj=12.05 mas, bpa=10.54 degrees ! Estimated noise=1.12745 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 (-14.8, -35.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0180913 Jy ! Component: 100 - total flux cleaned = 0.0238513 Jy ! Component: 150 - total flux cleaned = 0.028333 Jy ! Component: 200 - total flux cleaned = 0.0298584 Jy ! Total flux subtracted in 200 components = 0.0298584 Jy ! Clean residual min=-0.008434 max=0.010361 Jy/beam ! Clean residual mean=0.000107 rms=0.002691 Jy/beam ! Combined flux in latest and established models = 1.34545 Jy ! Performing phase self-cal ! Adding 44 model components to the UV plane model. ! The established model now contains 218 components and 1.34545 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.078612Jy sigma=0.542776 ! Fit after self-cal, rms=0.077893Jy sigma=0.528664 ! 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.077893Jy sigma=0.528664 ! Fit after self-cal, rms=0.076883Jy sigma=0.523673 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.39 mas, bmaj=12 mas, bpa=10.74 degrees ! Estimated noise=1.1231 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00571095 Jy ! Component: 100 - total flux cleaned = 0.00823201 Jy ! Component: 150 - total flux cleaned = 0.0102758 Jy ! Component: 200 - total flux cleaned = 0.011734 Jy ! Total flux subtracted in 200 components = 0.011734 Jy ! Clean residual min=-0.006294 max=0.007231 Jy/beam ! Clean residual mean=0.000069 rms=0.001833 Jy/beam ! Combined flux in latest and established models = 1.35719 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 242 components and 1.35719 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.074371Jy sigma=0.507191 ! Fit after self-cal, rms=0.074312Jy sigma=0.503751 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.00364059 Jy ! Component: 100 - total flux cleaned = 0.00453392 Jy ! Component: 150 - total flux cleaned = 0.00533969 Jy ! Component: 200 - total flux cleaned = 0.00608814 Jy ! Total flux subtracted in 200 components = 0.00608814 Jy ! Clean residual min=-0.005281 max=0.005945 Jy/beam ! Clean residual mean=0.000046 rms=0.001574 Jy/beam ! Combined flux in latest and established models = 1.36328 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 267 components and 1.36328 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.072934Jy sigma=0.497963 ! Fit after self-cal, rms=0.072967Jy sigma=0.497091 ! 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.072967Jy sigma=0.497091 ! Fit after self-cal, rms=0.072523Jy sigma=0.494671 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.385 mas, bmaj=11.99 mas, bpa=10.92 degrees ! Estimated noise=1.121 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0017689 Jy ! Component: 100 - total flux cleaned = 0.00256243 Jy ! Component: 150 - total flux cleaned = 0.00343999 Jy ! Component: 200 - total flux cleaned = 0.00438342 Jy ! Total flux subtracted in 200 components = 0.00438342 Jy ! Clean residual min=-0.004592 max=0.005239 Jy/beam ! Clean residual mean=0.000040 rms=0.001335 Jy/beam ! Combined flux in latest and established models = 1.36766 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 285 components and 1.36766 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.071432Jy sigma=0.490008 ! Fit after self-cal, rms=0.071493Jy sigma=0.489529 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.00139517 Jy ! Component: 100 - total flux cleaned = 0.00202728 Jy ! Component: 150 - total flux cleaned = 0.00284088 Jy ! Component: 200 - total flux cleaned = 0.00338427 Jy ! Total flux subtracted in 200 components = 0.00338427 Jy ! Clean residual min=-0.004355 max=0.004791 Jy/beam ! Clean residual mean=0.000029 rms=0.001240 Jy/beam ! Combined flux in latest and established models = 1.37104 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 293 components and 1.37104 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070742Jy sigma=0.487057 ! Fit after self-cal, rms=0.070801Jy sigma=0.486848 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 15 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070801Jy sigma=0.486848 ! Fit after self-cal, rms=0.070648Jy sigma=0.485820 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.17 mas, bmaj=15.78 mas, bpa=11.84 degrees ! Estimated noise=0.95556 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00340166 Jy ! Component: 100 - total flux cleaned = 0.00474745 Jy ! Component: 150 - total flux cleaned = 0.00515723 Jy ! Component: 200 - total flux cleaned = 0.00545904 Jy ! Total flux subtracted in 200 components = 0.00545904 Jy ! Clean residual min=-0.003513 max=0.004264 Jy/beam ! Clean residual mean=0.000031 rms=0.001070 Jy/beam ! Combined flux in latest and established models = 1.3765 Jy selfcal ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 312 components and 1.3765 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070478Jy sigma=0.483414 ! Fit after self-cal, rms=0.070465Jy sigma=0.483222 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.000966864 Jy ! Component: 100 - total flux cleaned = 0.00116894 Jy ! Component: 150 - total flux cleaned = 0.00129494 Jy ! Component: 200 - total flux cleaned = 0.00141232 Jy ! Total flux subtracted in 200 components = 0.00141232 Jy ! Clean residual min=-0.003310 max=0.004050 Jy/beam ! Clean residual mean=0.000026 rms=0.001014 Jy/beam ! Combined flux in latest and established models = 1.37791 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 326 components and 1.37791 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.070375Jy sigma=0.482023 ! Fit after self-cal, rms=0.070380Jy sigma=0.481952 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for VLBA stations ===" ! === Finished amplitude self-cal for VLBA stations === unflag * ! Unflagging all baselines in the currently selected channels. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. integer lun lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_selfant_process.csh /tmp/ff.ff" ![@/tmp/ff.ff] selfant 1:BR,true ! 1:BR self-cal status=fixed weight=1 selfant 1:FD,true ! 1:FD self-cal status=fixed weight=1 selfant 1:HN,true ! 1:HN self-cal status=fixed weight=1 selfant 1:KP,true ! 1:KP self-cal status=fixed weight=1 selfant 1:LA,true ! 1:LA self-cal status=fixed weight=1 selfant 1:NL,true ! 1:NL self-cal status=fixed weight=1 selfant 1:OV,true ! 1:OV self-cal status=fixed weight=1 selfant 1:PT,true ! 1:PT self-cal status=fixed weight=1 selfant 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 2:BR,true ! 2:BR self-cal status=fixed weight=1 selfant 2:FD,true ! 2:FD self-cal status=fixed weight=1 selfant 2:HN,true ! 2:HN self-cal status=fixed weight=1 selfant 2:KP,true ! 2:KP self-cal status=fixed weight=1 selfant 2:LA,true ! 2:LA self-cal status=fixed weight=1 selfant 2:MK,true ! 2:MK self-cal status=fixed weight=1 selfant 2:NL,true ! 2:NL self-cal status=fixed weight=1 selfant 2:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 2:PT,true ! 2:PT self-cal status=fixed weight=1 selfant 2:SC,true ! 2:SC self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] 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.070380Jy sigma=0.481952 ! Fit after self-cal, rms=0.070380Jy sigma=0.481952 selfant "",false ! 1:BR self-cal status=correctable weight=1 ! 1:FD self-cal status=correctable weight=1 ! 1:HN self-cal status=correctable weight=1 ! 1:KP self-cal status=correctable weight=1 ! 1:LA self-cal status=correctable weight=1 ! 1:NL self-cal status=correctable weight=1 ! 1:OV self-cal status=correctable weight=1 ! 1:PT self-cal status=correctable weight=1 ! 1:SC self-cal status=correctable weight=1 ! 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:MK self-cal status=correctable weight=1 ! 2:NL self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 2:PT self-cal status=correctable weight=1 ! 2:SC self-cal status=correctable weight=1 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.070380Jy sigma=0.481952 ! Fit after self-cal, rms=0.070199Jy sigma=0.481443 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.672 mas, bmaj=9.721 mas, bpa=13.09 degrees ! Estimated noise=2.9139 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 (14.8, 69.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000822564 Jy ! Component: 100 - total flux cleaned = 0.00124791 Jy ! Component: 150 - total flux cleaned = 0.00183086 Jy ! Component: 200 - total flux cleaned = 0.00218416 Jy ! Total flux subtracted in 200 components = 0.00218416 Jy ! Clean residual min=-0.005466 max=0.004880 Jy/beam ! Clean residual mean=0.000015 rms=0.001255 Jy/beam ! Combined flux in latest and established models = 1.3801 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 344 components and 1.3801 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068804Jy sigma=0.479881 ! Fit after self-cal, rms=0.068929Jy sigma=0.479489 ! 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=6.383 mas, bmaj=11.98 mas, bpa=11.1 degrees ! Estimated noise=1.11967 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 (54.8, -102.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00186985 Jy ! Component: 100 - total flux cleaned = 0.00297034 Jy ! Component: 150 - total flux cleaned = 0.00325166 Jy ! Component: 200 - total flux cleaned = 0.00359521 Jy ! Component: 250 - total flux cleaned = 0.00399611 Jy ! Component: 300 - total flux cleaned = 0.00429658 Jy ! Component: 350 - total flux cleaned = 0.00451121 Jy ! Component: 400 - total flux cleaned = 0.0044427 Jy ! Total flux subtracted in 400 components = 0.0044427 Jy ! Clean residual min=-0.002933 max=0.003089 Jy/beam ! Clean residual mean=0.000021 rms=0.000799 Jy/beam ! Combined flux in latest and established models = 1.38454 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 384 components and 1.38454 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.068175Jy sigma=0.475761 ! Fit after self-cal, rms=0.068205Jy sigma=0.475404 ! 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) ! Added new window around map position (-49.2, -30.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00304493 Jy ! Component: 100 - total flux cleaned = 0.00366668 Jy ! Component: 150 - total flux cleaned = 0.00366607 Jy ! Component: 200 - total flux cleaned = 0.00359137 Jy ! Component: 250 - total flux cleaned = 0.00345074 Jy ! Component: 300 - total flux cleaned = 0.00351746 Jy ! Component: 350 - total flux cleaned = 0.00332695 Jy ! Component: 400 - total flux cleaned = 0.00320491 Jy ! Total flux subtracted in 400 components = 0.00320491 Jy ! Clean residual min=-0.002541 max=0.002456 Jy/beam ! Clean residual mean=0.000021 rms=0.000632 Jy/beam ! Combined flux in latest and established models = 1.38775 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 418 components and 1.38775 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.067562Jy sigma=0.473136 ! Fit after self-cal, rms=0.067586Jy sigma=0.472922 ! Inverting map print "=== Finished natural weighting clean for al sites ===" ! === Finished natural weighting clean for al sites === unflag * ! Unflagging all baselines in the currently selected channels. ! Applying 2560 buffered edits. uncalib false, false, true ! uncal: All telescope correction flags have been un-done. selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.067586Jy sigma=0.472922 ! Fit after self-cal, rms=0.067522Jy sigma=0.472681 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.384 mas, bmaj=11.98 mas, bpa=11.07 degrees ! Estimated noise=1.11967 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0010718 Jy ! Component: 100 - total flux cleaned = 0.00126128 Jy ! Component: 150 - total flux cleaned = 0.00126041 Jy ! Component: 200 - total flux cleaned = 0.00143192 Jy ! Component: 250 - total flux cleaned = 0.00137675 Jy ! Component: 300 - total flux cleaned = 0.00132304 Jy ! Component: 350 - total flux cleaned = 0.00142668 Jy ! Component: 400 - total flux cleaned = 0.00137447 Jy ! Total flux subtracted in 400 components = 0.00137447 Jy ! Clean residual min=-0.002243 max=0.002340 Jy/beam ! Clean residual mean=0.000016 rms=0.000546 Jy/beam ! Combined flux in latest and established models = 1.38912 Jy selfcal ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 451 components and 1.38912 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.067122Jy sigma=0.471719 ! Fit after self-cal, rms=0.067156Jy sigma=0.471639 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.000789294 Jy ! Component: 100 - total flux cleaned = 0.000942832 Jy ! Component: 150 - total flux cleaned = 0.000943205 Jy ! Component: 200 - total flux cleaned = 0.000943957 Jy ! Component: 250 - total flux cleaned = 0.00103616 Jy ! Component: 300 - total flux cleaned = 0.000991195 Jy ! Component: 350 - total flux cleaned = 0.00107962 Jy ! Component: 400 - total flux cleaned = 0.00112332 Jy ! Total flux subtracted in 400 components = 0.00112332 Jy ! Clean residual min=-0.002182 max=0.002127 Jy/beam ! Clean residual mean=0.000013 rms=0.000502 Jy/beam ! Combined flux in latest and established models = 1.39024 Jy ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 477 components and 1.39024 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066847Jy sigma=0.471072 ! Fit after self-cal, rms=0.066872Jy sigma=0.471022 ! 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.066872Jy sigma=0.471022 ! Fit after self-cal, rms=0.066865Jy sigma=0.470820 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=6.386 mas, bmaj=11.98 mas, bpa=11.07 degrees ! Estimated noise=1.11962 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000509307 Jy ! Component: 100 - total flux cleaned = 0.000686924 Jy ! Component: 150 - total flux cleaned = 0.000771176 Jy ! Component: 200 - total flux cleaned = 0.000771824 Jy ! Component: 250 - total flux cleaned = 0.000854574 Jy ! Component: 300 - total flux cleaned = 0.000935255 Jy ! Component: 350 - total flux cleaned = 0.000974897 Jy ! Component: 400 - total flux cleaned = 0.00101342 Jy ! Total flux subtracted in 400 components = 0.00101342 Jy ! Clean residual min=-0.002155 max=0.001979 Jy/beam ! Clean residual mean=0.000010 rms=0.000472 Jy/beam ! Combined flux in latest and established models = 1.39126 Jy selfcal ! Performing phase self-cal ! Adding 68 model components to the UV plane model. ! The established model now contains 497 components and 1.39126 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066602Jy sigma=0.470416 ! Fit after self-cal, rms=0.066624Jy sigma=0.470386 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.000428034 Jy ! Component: 100 - total flux cleaned = 0.000429037 Jy ! Component: 150 - total flux cleaned = 0.000429539 Jy ! Component: 200 - total flux cleaned = 0.000430483 Jy ! Component: 250 - total flux cleaned = 0.000467062 Jy ! Component: 300 - total flux cleaned = 0.000466786 Jy ! Component: 350 - total flux cleaned = 0.000466447 Jy ! Component: 400 - total flux cleaned = 0.000467451 Jy ! Total flux subtracted in 400 components = 0.000467451 Jy ! Clean residual min=-0.002097 max=0.001911 Jy/beam ! Clean residual mean=0.000009 rms=0.000449 Jy/beam ! Combined flux in latest and established models = 1.39172 Jy ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 513 components and 1.39172 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066405Jy sigma=0.470094 ! Fit after self-cal, rms=0.066421Jy sigma=0.470076 ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime3 ! Performing amp+phase self-cal over 1 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.066421Jy sigma=0.470076 ! Fit after self-cal, rms=0.065467Jy sigma=0.450034 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.16 mas, bmaj=15.76 mas, bpa=11.78 degrees ! Estimated noise=0.955938 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000118328 Jy ! Component: 100 - total flux cleaned = -0.000259389 Jy ! Component: 150 - total flux cleaned = -0.000285401 Jy ! Component: 200 - total flux cleaned = -0.000369428 Jy ! Component: 250 - total flux cleaned = -0.000368431 Jy ! Component: 300 - total flux cleaned = -0.000368133 Jy ! Component: 350 - total flux cleaned = -0.000403467 Jy ! Component: 400 - total flux cleaned = -0.00038479 Jy ! Total flux subtracted in 400 components = -0.00038479 Jy ! Clean residual min=-0.001577 max=0.001431 Jy/beam ! Clean residual mean=0.000011 rms=0.000353 Jy/beam ! Combined flux in latest and established models = 1.39134 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 555 components and 1.39134 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065452Jy sigma=0.449810 ! Fit after self-cal, rms=0.065461Jy sigma=0.449799 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 = 2.46264e-05 Jy ! Component: 100 - total flux cleaned = -9.30459e-06 Jy ! Component: 150 - total flux cleaned = -4.16358e-05 Jy ! Component: 200 - total flux cleaned = -7.31292e-05 Jy ! Component: 250 - total flux cleaned = -0.000103829 Jy ! Component: 300 - total flux cleaned = -0.000103641 Jy ! Component: 350 - total flux cleaned = -0.000132536 Jy ! Component: 400 - total flux cleaned = -0.0001038 Jy ! Total flux subtracted in 400 components = -0.0001038 Jy ! Clean residual min=-0.001550 max=0.001377 Jy/beam ! Clean residual mean=0.000010 rms=0.000340 Jy/beam ! Combined flux in latest and established models = 1.39123 Jy ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 593 components and 1.39123 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.065434Jy sigma=0.449691 ! Fit after self-cal, rms=0.065437Jy sigma=0.449687 ! 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=6.391 mas, bmaj=11.99 mas, bpa=11.12 degrees ! Estimated noise=1.11994 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 6.391 x 11.99 at 11.12 degrees (North through East) ! Clean map min=-0.0024219 max=0.70292 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.065437Jy sigma=0.449687 ! Fit after self-cal, rms=0.063263Jy sigma=0.391234 ! Inverting map and beam ! Estimated beam: bmin=6.386 mas, bmaj=11.98 mas, bpa=11.06 degrees ! Estimated noise=1.11952 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000166718 Jy ! Component: 100 - total flux cleaned = 0.000348307 Jy ! Component: 150 - total flux cleaned = 0.000486721 Jy ! Component: 200 - total flux cleaned = 0.00062108 Jy ! Component: 250 - total flux cleaned = 0.000751859 Jy ! Component: 300 - total flux cleaned = 0.000815933 Jy ! Component: 350 - total flux cleaned = 0.00087858 Jy ! Component: 400 - total flux cleaned = 0.000908904 Jy ! Total flux subtracted in 400 components = 0.000908904 Jy ! Clean residual min=-0.001827 max=0.001757 Jy/beam ! Clean residual mean=0.000008 rms=0.000390 Jy/beam ! Combined flux in latest and established models = 1.39214 Jy ! Performing phase self-cal ! Adding 50 model components to the UV plane model. ! The established model now contains 624 components and 1.39214 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.063038Jy sigma=0.391011 ! Fit after self-cal, rms=0.063052Jy sigma=0.390998 ! Inverting map ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000292963 Jy ! Component: 100 - total flux cleaned = 0.000322586 Jy ! Component: 150 - total flux cleaned = 0.000381829 Jy ! Component: 200 - total flux cleaned = 0.00038139 Jy ! Component: 250 - total flux cleaned = 0.00043843 Jy ! Component: 300 - total flux cleaned = 0.000467096 Jy ! Component: 350 - total flux cleaned = 0.000410918 Jy ! Component: 400 - total flux cleaned = 0.000410687 Jy ! Total flux subtracted in 400 components = 0.000410687 Jy ! Clean residual min=-0.001782 max=0.001754 Jy/beam ! Clean residual mean=0.000007 rms=0.000382 Jy/beam ! Combined flux in latest and established models = 1.39255 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 650 components and 1.39256 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.062883Jy sigma=0.390841 ! Fit after self-cal, rms=0.062896Jy sigma=0.390832 ! 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.062896Jy sigma=0.390832 ! Fit after self-cal, rms=0.062964Jy sigma=0.390818 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=4.675 mas, bmaj=9.727 mas, bpa=13.03 degrees ! Estimated noise=2.91371 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.570347 Jy ! Component: 100 - total flux cleaned = 0.820492 Jy ! Component: 150 - total flux cleaned = 0.969234 Jy ! Component: 200 - total flux cleaned = 1.06674 Jy ! Component: 250 - total flux cleaned = 1.13519 Jy ! Component: 300 - total flux cleaned = 1.18573 Jy ! Component: 350 - total flux cleaned = 1.2241 Jy ! Component: 400 - total flux cleaned = 1.25391 Jy ! Total flux subtracted in 400 components = 1.25391 Jy ! Clean residual min=-0.015391 max=0.019588 Jy/beam ! Clean residual mean=0.000269 rms=0.005182 Jy/beam ! Combined flux in latest and established models = 1.25391 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 145 model components to the UV plane model. ! The established model now contains 145 components and 1.25391 Jy ! Inverting map and beam ! Estimated beam: bmin=6.39 mas, bmaj=11.99 mas, bpa=11.09 degrees ! Estimated noise=1.11965 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.0323011 Jy ! Component: 100 - total flux cleaned = 0.0549684 Jy ! Component: 150 - total flux cleaned = 0.0718162 Jy ! Component: 200 - total flux cleaned = 0.0845663 Jy ! Component: 250 - total flux cleaned = 0.0943421 Jy ! Component: 300 - total flux cleaned = 0.10199 Jy ! Component: 350 - total flux cleaned = 0.108315 Jy ! Component: 400 - total flux cleaned = 0.113295 Jy ! Total flux subtracted in 400 components = 0.113295 Jy ! Clean residual min=-0.004566 max=0.006019 Jy/beam ! Clean residual mean=0.000054 rms=0.001397 Jy/beam ! Combined flux in latest and established models = 1.3672 Jy ! Adding 151 model components to the UV plane model. ! The established model now contains 280 components and 1.3672 Jy ! Inverting map ! Added new window around map position (6.8, -3.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00420426 Jy ! Component: 100 - total flux cleaned = 0.00683599 Jy ! Component: 150 - total flux cleaned = 0.00826456 Jy ! Component: 200 - total flux cleaned = 0.00937833 Jy ! Component: 250 - total flux cleaned = 0.010652 Jy ! Component: 300 - total flux cleaned = 0.0121802 Jy ! Component: 350 - total flux cleaned = 0.0136001 Jy ! Component: 400 - total flux cleaned = 0.0145282 Jy ! Total flux subtracted in 400 components = 0.0145282 Jy ! Clean residual min=-0.002435 max=0.002750 Jy/beam ! Clean residual mean=0.000026 rms=0.000753 Jy/beam ! Combined flux in latest and established models = 1.38173 Jy ! Adding 82 model components to the UV plane model. ! The established model now contains 357 components and 1.38173 Jy ! Inverting map ! Added new window around map position (-26.8, -21.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00280876 Jy ! Component: 100 - total flux cleaned = 0.00414387 Jy ! Component: 150 - total flux cleaned = 0.0046792 Jy ! Component: 200 - total flux cleaned = 0.00543505 Jy ! Component: 250 - total flux cleaned = 0.00631042 Jy ! Component: 300 - total flux cleaned = 0.00668546 Jy ! Component: 350 - total flux cleaned = 0.00719427 Jy ! Component: 400 - total flux cleaned = 0.007541 Jy ! Total flux subtracted in 400 components = 0.007541 Jy ! Clean residual min=-0.001957 max=0.002271 Jy/beam ! Clean residual mean=0.000018 rms=0.000559 Jy/beam ! Combined flux in latest and established models = 1.38927 Jy ! Adding 80 model components to the UV plane model. ! The established model now contains 415 components and 1.38927 Jy ! Inverting map ! Added new window around map position (-32.8, 84). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00121893 Jy ! Component: 100 - total flux cleaned = 0.00143358 Jy ! Component: 150 - total flux cleaned = 0.00163548 Jy ! Component: 200 - total flux cleaned = 0.00150511 Jy ! Component: 250 - total flux cleaned = 0.0016268 Jy ! Component: 300 - total flux cleaned = 0.00156767 Jy ! Component: 350 - total flux cleaned = 0.00174013 Jy ! Component: 400 - total flux cleaned = 0.00179468 Jy ! Total flux subtracted in 400 components = 0.00179468 Jy ! Clean residual min=-0.001730 max=0.001685 Jy/beam ! Clean residual mean=0.000012 rms=0.000451 Jy/beam ! Combined flux in latest and established models = 1.39107 Jy ! Adding 71 model components to the UV plane model. ! The established model now contains 466 components and 1.39107 Jy ! Inverting map ! Added new window around map position (17.2, -93.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000747213 Jy ! Component: 100 - total flux cleaned = 0.00114389 Jy ! Component: 150 - total flux cleaned = 0.00147031 Jy ! Component: 200 - total flux cleaned = 0.00168057 Jy ! Component: 250 - total flux cleaned = 0.0018324 Jy ! Component: 300 - total flux cleaned = 0.00198039 Jy ! Component: 350 - total flux cleaned = 0.00202895 Jy ! Component: 400 - total flux cleaned = 0.00202949 Jy ! Total flux subtracted in 400 components = 0.00202949 Jy ! Clean residual min=-0.001591 max=0.001481 Jy/beam ! Clean residual mean=0.000010 rms=0.000393 Jy/beam ! Combined flux in latest and established models = 1.3931 Jy ! Adding 69 model components to the UV plane model. ! The established model now contains 509 components and 1.3931 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=10.15 mas, bmaj=15.76 mas, bpa=11.77 degrees ! Estimated noise=0.955827 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000278087 Jy ! Component: 100 - total flux cleaned = -0.000311248 Jy ! Component: 150 - total flux cleaned = -0.00025228 Jy ! Component: 200 - total flux cleaned = -0.000167524 Jy ! Component: 250 - total flux cleaned = -8.76326e-05 Jy ! Component: 300 - total flux cleaned = -1.18338e-05 Jy ! Component: 350 - total flux cleaned = 1.23499e-05 Jy ! Component: 400 - total flux cleaned = 5.78573e-05 Jy ! Total flux subtracted in 400 components = 5.78573e-05 Jy ! Clean residual min=-0.001089 max=0.001203 Jy/beam ! Clean residual mean=0.000011 rms=0.000285 Jy/beam ! Combined flux in latest and established models = 1.39315 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) ! Added new window around map position (-61.6, -94). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000177775 Jy ! Component: 100 - total flux cleaned = 0.000300277 Jy ! Component: 150 - total flux cleaned = 0.000325466 Jy ! Component: 200 - total flux cleaned = 0.000349374 Jy ! Component: 250 - total flux cleaned = 0.000395698 Jy ! Component: 300 - total flux cleaned = 0.000460641 Jy ! Component: 350 - total flux cleaned = 0.000608652 Jy ! Component: 400 - total flux cleaned = 0.000691819 Jy ! Total flux subtracted in 400 components = 0.000691819 Jy ! Clean residual min=-0.001051 max=0.000844 Jy/beam ! Clean residual mean=0.000010 rms=0.000239 Jy/beam ! Combined flux in latest and established models = 1.39385 Jy ! Adding 124 model components to the UV plane model. ! The established model now contains 631 components and 1.39385 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=6.39 mas, bmaj=11.99 mas, bpa=11.09 degrees ! Estimated noise=1.11965 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000182126 Jy ! Component: 100 - total flux cleaned = 9.38937e-05 Jy ! Component: 150 - total flux cleaned = 5.33101e-05 Jy ! Component: 200 - total flux cleaned = 9.53263e-05 Jy ! Component: 250 - total flux cleaned = 0.000252366 Jy ! Component: 300 - total flux cleaned = 0.000289381 Jy ! Component: 350 - total flux cleaned = 0.000287328 Jy ! Component: 400 - total flux cleaned = 0.000250295 Jy ! Total flux subtracted in 400 components = 0.000250295 Jy ! Clean residual min=-0.001273 max=0.001052 Jy/beam ! Clean residual mean=0.000008 rms=0.000295 Jy/beam ! Combined flux in latest and established models = 1.3941 Jy selfcal ! Performing phase self-cal ! Adding 62 model components to the UV plane model. ! The established model now contains 675 components and 1.3941 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.062704Jy sigma=0.389300 ! Fit after self-cal, rms=0.062522Jy sigma=0.389094 wmodel /scr/pima/rdv95_uvs/J1543-0757_S_map.mod ! Writing 675 model components to file: /scr/pima/rdv95_uvs/J1543-0757_S_map.mod wobs /scr/pima/rdv95_uvs/J1543-0757_S_uvs.fits ! Applying 1004 buffered edits. ! Writing UV FITS file: /scr/pima/rdv95_uvs/J1543-0757_S_uvs.fits wwins /scr/pima/rdv95_uvs/J1543-0757_S_map.win ! wwins: Wrote 11 windows to /scr/pima/rdv95_uvs/J1543-0757_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 2.25124e-05 Jy ! Component: 100 - total flux cleaned = 6.16476e-05 Jy ! Component: 150 - total flux cleaned = 4.20745e-05 Jy ! Component: 200 - total flux cleaned = -1.2827e-05 Jy ! Component: 250 - total flux cleaned = -8.30741e-05 Jy ! Component: 300 - total flux cleaned = -0.000117057 Jy ! Component: 350 - total flux cleaned = -0.000100669 Jy ! Component: 400 - total flux cleaned = -0.000100293 Jy ! Component: 450 - total flux cleaned = -5.23328e-05 Jy ! Component: 500 - total flux cleaned = -2.1102e-05 Jy ! Component: 550 - total flux cleaned = -5.74732e-06 Jy ! Component: 600 - total flux cleaned = -2.08385e-05 Jy ! Component: 650 - total flux cleaned = -2.06467e-05 Jy ! Component: 700 - total flux cleaned = -2.02788e-05 Jy ! Component: 750 - total flux cleaned = -2.00692e-05 Jy ! Component: 800 - total flux cleaned = -5.80842e-06 Jy ! Component: 850 - total flux cleaned = -2.00371e-05 Jy ! Component: 900 - total flux cleaned = 2.1958e-05 Jy ! Component: 950 - total flux cleaned = 4.95565e-05 Jy ! Component: 1000 - total flux cleaned = 0.000132031 Jy ! Component: 1050 - total flux cleaned = 0.000172829 Jy ! Component: 1100 - total flux cleaned = 0.00018631 Jy ! Component: 1150 - total flux cleaned = 0.000199833 Jy ! Component: 1200 - total flux cleaned = 0.000186484 Jy ! Component: 1250 - total flux cleaned = 0.000186686 Jy ! Component: 1300 - total flux cleaned = 0.000186706 Jy ! Component: 1350 - total flux cleaned = 0.000225653 Jy ! Component: 1400 - total flux cleaned = 0.000212716 Jy ! Component: 1450 - total flux cleaned = 0.00023834 Jy ! Component: 1500 - total flux cleaned = 0.000200024 Jy ! Component: 1550 - total flux cleaned = 0.000250651 Jy ! Component: 1600 - total flux cleaned = 0.000212986 Jy ! Component: 1650 - total flux cleaned = 0.000225592 Jy ! Component: 1700 - total flux cleaned = 0.000238025 Jy ! Component: 1750 - total flux cleaned = 0.00022574 Jy ! Component: 1800 - total flux cleaned = 0.000274737 Jy ! Component: 1850 - total flux cleaned = 0.000262492 Jy ! Component: 1900 - total flux cleaned = 0.000262481 Jy ! Component: 1950 - total flux cleaned = 0.000250686 Jy ! Component: 2000 - total flux cleaned = 0.000310326 Jy ! Component: 2050 - total flux cleaned = 0.000298597 Jy ! Component: 2100 - total flux cleaned = 0.000310417 Jy ! Component: 2150 - total flux cleaned = 0.000333941 Jy ! Component: 2200 - total flux cleaned = 0.000345549 Jy ! Component: 2250 - total flux cleaned = 0.00032231 Jy ! Component: 2300 - total flux cleaned = 0.000322251 Jy ! Component: 2350 - total flux cleaned = 0.000310738 Jy ! Component: 2400 - total flux cleaned = 0.00033352 Jy ! Component: 2450 - total flux cleaned = 0.000322026 Jy ! Component: 2500 - total flux cleaned = 0.000321987 Jy ! Component: 2550 - total flux cleaned = 0.000321994 Jy ! Component: 2600 - total flux cleaned = 0.000310973 Jy ! Component: 2650 - total flux cleaned = 0.000333216 Jy ! Component: 2700 - total flux cleaned = 0.000355319 Jy ! Component: 2750 - total flux cleaned = 0.000322336 Jy ! Component: 2800 - total flux cleaned = 0.000355218 Jy ! Component: 2850 - total flux cleaned = 0.000322539 Jy ! Component: 2900 - total flux cleaned = 0.000311717 Jy ! Component: 2950 - total flux cleaned = 0.000333248 Jy ! Component: 3000 - total flux cleaned = 0.000311725 Jy ! Component: 3050 - total flux cleaned = 0.000311727 Jy ! Component: 3100 - total flux cleaned = 0.000343604 Jy ! Component: 3150 - total flux cleaned = 0.000311832 Jy ! Component: 3200 - total flux cleaned = 0.000311827 Jy ! Component: 3250 - total flux cleaned = 0.000311795 Jy ! Component: 3300 - total flux cleaned = 0.000332643 Jy ! Component: 3350 - total flux cleaned = 0.000332563 Jy ! Component: 3400 - total flux cleaned = 0.000291124 Jy ! Component: 3450 - total flux cleaned = 0.000322113 Jy ! Component: 3500 - total flux cleaned = 0.000352821 Jy ! Component: 3550 - total flux cleaned = 0.00033246 Jy ! Component: 3600 - total flux cleaned = 0.000342646 Jy ! Component: 3650 - total flux cleaned = 0.000342594 Jy ! Component: 3700 - total flux cleaned = 0.000352641 Jy ! Component: 3750 - total flux cleaned = 0.000372611 Jy ! Component: 3800 - total flux cleaned = 0.000372567 Jy ! Component: 3850 - total flux cleaned = 0.000362553 Jy ! Component: 3900 - total flux cleaned = 0.000352616 Jy ! Component: 3950 - total flux cleaned = 0.000342761 Jy ! Component: 4000 - total flux cleaned = 0.000382046 Jy ! Component: 4050 - total flux cleaned = 0.000352616 Jy ! Total flux subtracted in 4096 components = 0.000372112 Jy ! Clean residual min=-0.000564 max=0.000542 Jy/beam ! Clean residual mean=0.000006 rms=0.000188 Jy/beam ! Combined flux in latest and established models = 1.39447 Jy keep ! Adding 199 model components to the UV plane model. ! The established model now contains 860 components and 1.39447 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv95_uvs/J1543-0757_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 6.39 x 11.99 at 11.09 degrees (North through East) ! Clean map min=-0.0027824 max=0.70546 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv95_uvs/J1543-0757_S_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.7054 rms= 0.00018 Jy quit ! Quitting program ! Log file /scr/pima/rdv95_uvs/J1543-0757_S_uvs.log closed on Wed Nov 6 20:47:46 2013