! Started logfile: /scr/pima/rdv94_uvs/J1521+4336_X_uvs.log on Tue Oct 8 12:21:40 2013 obs /scr/pima/rdv94_uvs/J1521+4336_X_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J1521+4336_X_uva.fits ! AN table 1: 12 integrations on 36 of 36 possible baselines. ! AN table 2: 7 integrations on 28 of 28 possible baselines. ! AN table 3: 6 integrations on 36 of 36 possible baselines. ! AN table 4: 9 integrations on 21 of 21 possible baselines. ! AN table 5: 7 integrations on 21 of 21 possible baselines. ! Apparent sampling: 0.957627 visibilities/baseline/integration-bin. ! Found source: J1521+4336 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40599e+09 1.6e+07 1 1.6e+07 ! 02 2 8.47599e+09 1.6e+07 1 1.6e+07 ! 03 3 8.79099e+09 1.6e+07 1 1.6e+07 ! 04 4 8.89599e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 4520 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/rdv94_uvs/J1521+4336_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: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:ON self-cal status=correctable weight=1 ! 1:OV 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: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:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR 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:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:WZ 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 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=1.168623Jy sigma=11.680580 ! Fit after self-cal, rms=0.412349Jy sigma=3.762130 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 1:ON ! Flagging baselines of antenna 1:ON in the currently selected channels. 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:NY ! Flagging baselines of antenna 4:NY in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ 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.149 mas, bmaj=1.418 mas, bpa=-18.79 degrees ! Estimated noise=1.86708 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.487154 Jy ! Component: 100 - total flux cleaned = 0.596726 Jy ! Total flux subtracted in 100 components = 0.596726 Jy ! Clean residual min=-0.020470 max=0.037693 Jy/beam ! Clean residual mean=-0.000166 rms=0.004849 Jy/beam ! Combined flux in latest and established models = 0.596726 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 0.596726 Jy ! ! Correcting IF 1. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 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 12 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 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 12 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 12 telescope corrections were flagged in sub-array 1. ! A total of 7 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.141109Jy sigma=1.324649 ! Fit after self-cal, rms=0.140928Jy sigma=1.323021 ! Inverting map and beam ! Estimated beam: bmin=1.149 mas, bmaj=1.418 mas, bpa=-18.79 degrees ! Estimated noise=1.86708 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0372458 Jy ! Component: 100 - total flux cleaned = 0.054501 Jy ! Total flux subtracted in 100 components = 0.054501 Jy ! Clean residual min=-0.017589 max=0.017222 Jy/beam ! Clean residual mean=-0.000207 rms=0.003938 Jy/beam ! Combined flux in latest and established models = 0.651227 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 22 components and 0.651227 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.133040Jy sigma=1.146976 ! Fit after self-cal, rms=0.133007Jy sigma=1.146456 ! 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.528 mas, bmaj=1.712 mas, bpa=-17.67 degrees ! Estimated noise=1.18941 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.0166715 Jy ! Component: 100 - total flux cleaned = 0.0242757 Jy ! Component: 150 - total flux cleaned = 0.028176 Jy ! Component: 200 - total flux cleaned = 0.0288886 Jy ! Total flux subtracted in 200 components = 0.0288886 Jy ! Clean residual min=-0.016698 max=0.015989 Jy/beam ! Clean residual mean=-0.000567 rms=0.004117 Jy/beam ! Combined flux in latest and established models = 0.680115 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 34 components and 0.680115 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.132625Jy sigma=1.113500 ! Fit after self-cal, rms=0.132618Jy sigma=1.113334 ! 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 0.95 FD 0.89 HN 0.91 KP 0.96 ! LA 1.04 NL 1.00 ON 1.00* OV 0.97 ! PT 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.97 FD 0.90 KP 0.95 LA 1.06 ! MK 1.04 NL 0.98 NY 1.00* OV 0.96 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00 FD 0.89 HN 0.90 KP 1.80 ! LA 1.08 MK 0.96 NL 1.00 NY 1.00* ! OV 0.97 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.00 KP 0.96 LA 1.12 MK 0.91 ! NL 1.03 NY 1.00* OV 0.99 ! ! Telescope amplitude corrections in sub-array 5: ! BR 0.97 FD 0.93 HN 0.91 LA 1.14 ! NL 1.05 PT 1.26 WZ 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 1.35 FD 0.89 HN 0.94 KP 0.91 ! LA 1.01 NL 1.02 ON 1.00* OV 0.96 ! PT 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.38 FD 0.90 KP 0.92 LA 1.02 ! MK 1.09 NL 0.99 NY 1.00* OV 0.94 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.41 FD 0.89 HN 0.88 KP 1.78 ! LA 1.07 MK 0.99 NL 0.99 NY 1.00* ! OV 0.95 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.39 KP 0.93 LA 1.04 MK 0.93 ! NL 1.08 NY 1.00* OV 0.98 ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.16 FD 0.94 HN 0.95 LA 1.10 ! NL 1.03 PT 1.38 WZ 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 1.03 FD 0.89 HN 0.95 KP 0.81 ! LA 1.08 NL 1.05 ON 1.00* OV 1.02 ! PT 0.99 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.03 FD 0.92 KP 0.83 LA 1.06 ! MK 1.09 NL 1.02 NY 1.00* OV 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.07 FD 0.89 HN 0.92 KP 1.60 ! LA 1.16 MK 0.98 NL 1.01 NY 1.00* ! OV 1.01 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.07 KP 0.82 LA 1.12 MK 0.92 ! NL 1.11 NY 1.00* OV 1.05 ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.03 FD 0.97 HN 0.92 LA 1.17 ! NL 1.09 PT 1.36 WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 1.00 FD 0.92 HN 0.96 KP 0.86 ! LA 1.04 NL 1.08 ON 1.00* OV 0.99 ! PT 1.04 ! ! Telescope amplitude corrections in sub-array 2: ! BR 1.04 FD 1.00 KP 0.86 LA 1.06 ! MK 1.08 NL 1.04 NY 1.00* OV 1.00 ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.04 FD 1.00 HN 0.91 KP 2.16 ! LA 1.15 MK 1.05 NL 1.03 NY 1.00* ! OV 0.97 ! ! Telescope amplitude corrections in sub-array 4: ! BR 1.12 KP 0.86 LA 1.12 MK 0.94 ! NL 1.13 NY 1.00* OV 1.01 ! ! Telescope amplitude corrections in sub-array 5: ! BR 1.02 FD 1.04 HN 0.97 LA 1.12 ! NL 1.09 PT 1.41 WZ 1.00* ! ! ! Fit before self-cal, rms=0.132618Jy sigma=1.113334 ! Fit after self-cal, rms=0.100421Jy sigma=0.704203 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=1.525 mas, bmaj=1.707 mas, bpa=-18.74 degrees ! Estimated noise=1.18707 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.00180339 Jy ! Component: 100 - total flux cleaned = -0.00172694 Jy ! Component: 150 - total flux cleaned = -0.00172544 Jy ! Component: 200 - total flux cleaned = -0.00184015 Jy ! Total flux subtracted in 200 components = -0.00184015 Jy ! Clean residual min=-0.002990 max=0.003293 Jy/beam ! Clean residual mean=-0.000022 rms=0.000774 Jy/beam ! Combined flux in latest and established models = 0.678275 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 43 components and 0.678275 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100306Jy sigma=0.702532 ! Fit after self-cal, rms=0.100451Jy sigma=0.701991 ! 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.100451Jy sigma=0.701991 ! Fit after self-cal, rms=0.100450Jy sigma=0.701764 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.527 mas, bmaj=1.707 mas, bpa=-18.37 degrees ! Estimated noise=1.1863 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00138117 Jy ! Component: 100 - total flux cleaned = -0.00138237 Jy ! Component: 150 - total flux cleaned = -0.00138486 Jy ! Component: 200 - total flux cleaned = -0.00153204 Jy ! Total flux subtracted in 200 components = -0.00153204 Jy ! Clean residual min=-0.002846 max=0.003063 Jy/beam ! Clean residual mean=-0.000022 rms=0.000750 Jy/beam ! Combined flux in latest and established models = 0.676743 Jy selfcal ! Performing phase self-cal ! Adding 9 model components to the UV plane model. ! The established model now contains 43 components and 0.676743 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100375Jy sigma=0.700672 ! Fit after self-cal, rms=0.100370Jy sigma=0.700585 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.87364e-05 Jy ! Component: 100 - total flux cleaned = -3.17265e-05 Jy ! Component: 150 - total flux cleaned = -0.000156601 Jy ! Component: 200 - total flux cleaned = -0.00023368 Jy ! Total flux subtracted in 200 components = -0.00023368 Jy ! Clean residual min=-0.002795 max=0.003029 Jy/beam ! Clean residual mean=-0.000021 rms=0.000742 Jy/beam ! Combined flux in latest and established models = 0.67651 Jy ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 43 components and 0.67651 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100330Jy sigma=0.700070 ! Fit after self-cal, rms=0.100330Jy sigma=0.700024 ! 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.100330Jy sigma=0.700024 ! Fit after self-cal, rms=0.100399Jy sigma=0.699803 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=1.529 mas, bmaj=1.707 mas, bpa=-17.99 degrees ! Estimated noise=1.18615 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.0010323 Jy ! Component: 100 - total flux cleaned = -0.00102914 Jy ! Component: 150 - total flux cleaned = -0.00102731 Jy ! Component: 200 - total flux cleaned = -0.00106311 Jy ! Total flux subtracted in 200 components = -0.00106311 Jy ! Clean residual min=-0.002790 max=0.002904 Jy/beam ! Clean residual mean=-0.000020 rms=0.000726 Jy/beam ! Combined flux in latest and established models = 0.675446 Jy selfcal ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 45 components and 0.675446 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100363Jy sigma=0.699239 ! Fit after self-cal, rms=0.100363Jy sigma=0.699206 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 = 1.0692e-06 Jy ! Component: 100 - total flux cleaned = -9.63425e-05 Jy ! Component: 150 - total flux cleaned = -0.000186254 Jy ! Component: 200 - total flux cleaned = -0.000271879 Jy ! Total flux subtracted in 200 components = -0.000271879 Jy ! Clean residual min=-0.002770 max=0.002901 Jy/beam ! Clean residual mean=-0.000020 rms=0.000722 Jy/beam ! Combined flux in latest and established models = 0.675175 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 47 components and 0.675174 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100343Jy sigma=0.698947 ! Fit after self-cal, rms=0.100344Jy sigma=0.698928 ! 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.100344Jy sigma=0.698928 ! Fit after self-cal, rms=0.100375Jy sigma=0.698799 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.557 mas, bmaj=2.791 mas, bpa=47.46 degrees ! Estimated noise=1.327 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000743735 Jy ! Component: 100 - total flux cleaned = -0.000758307 Jy ! Component: 150 - total flux cleaned = -0.000770346 Jy ! Component: 200 - total flux cleaned = -0.000781087 Jy ! Total flux subtracted in 200 components = -0.000781087 Jy ! Clean residual min=-0.002975 max=0.002753 Jy/beam ! Clean residual mean=-0.000035 rms=0.000735 Jy/beam ! Combined flux in latest and established models = 0.674393 Jy selfcal ! Performing phase self-cal ! Adding 8 model components to the UV plane model. ! The established model now contains 51 components and 0.674393 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100363Jy sigma=0.698620 ! Fit after self-cal, rms=0.100362Jy sigma=0.698591 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 = -2.46053e-05 Jy ! Component: 100 - total flux cleaned = -5.24031e-05 Jy ! Component: 150 - total flux cleaned = -6.33124e-05 Jy ! Component: 200 - total flux cleaned = -8.19549e-05 Jy ! Total flux subtracted in 200 components = -8.19549e-05 Jy ! Clean residual min=-0.002955 max=0.002758 Jy/beam ! Clean residual mean=-0.000035 rms=0.000734 Jy/beam ! Combined flux in latest and established models = 0.674311 Jy ! Performing phase self-cal ! Adding 7 model components to the UV plane model. ! The established model now contains 51 components and 0.674311 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100354Jy sigma=0.698528 ! Fit after self-cal, rms=0.100353Jy sigma=0.698509 ! 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 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: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 3:BR,true ! 3:BR self-cal status=fixed weight=1 selfant 3:FD,true ! 3:FD self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:KP,true ! 3:KP self-cal status=fixed weight=1 selfant 3:LA,true ! 3:LA self-cal status=fixed weight=1 selfant 3:MK,true ! 3:MK self-cal status=fixed weight=1 selfant 3:NL,true ! 3:NL self-cal status=fixed weight=1 selfant 3:OV,true ! 3:OV self-cal status=fixed weight=1 selfant 4:BR,true ! 4:BR 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:NL,true ! 4:NL self-cal status=fixed weight=1 selfant 4:OV,true ! 4:OV self-cal status=fixed weight=1 selfant 5:BR,true ! 5:BR self-cal status=fixed weight=1 selfant 5:FD,true ! 5:FD self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:NL,true ! 5:NL self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT 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. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 7 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 9 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=0.110994Jy sigma=0.745042 ! Fit after self-cal, rms=0.102093Jy sigma=0.686508 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:ON self-cal status=correctable weight=1 ! 1:OV 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: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:NY self-cal status=correctable weight=1 ! 2:OV self-cal status=correctable weight=1 ! 3:BR self-cal status=correctable weight=1 ! 3:FD self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:KP self-cal status=correctable weight=1 ! 3:LA self-cal status=correctable weight=1 ! 3:MK self-cal status=correctable weight=1 ! 3:NL self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:OV self-cal status=correctable weight=1 ! 4:BR 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:NL self-cal status=correctable weight=1 ! 4:NY self-cal status=correctable weight=1 ! 4:OV self-cal status=correctable weight=1 ! 5:BR self-cal status=correctable weight=1 ! 5:FD self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:NL self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 5:WZ 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 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.102093Jy sigma=0.686508 ! Fit after self-cal, rms=0.102007Jy sigma=0.685854 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.7458 mas, bmaj=1.109 mas, bpa=-1.493 degrees ! Estimated noise=1.82968 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 (-33, 7.2). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0026927 Jy ! Component: 100 - total flux cleaned = 0.00229353 Jy ! Component: 150 - total flux cleaned = 0.00206078 Jy ! Component: 200 - total flux cleaned = 0.00205474 Jy ! Total flux subtracted in 200 components = 0.00205474 Jy ! Clean residual min=-0.005311 max=0.005012 Jy/beam ! Clean residual mean=-0.000006 rms=0.001254 Jy/beam ! Combined flux in latest and established models = 0.676366 Jy ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 67 components and 0.676366 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101969Jy sigma=0.685562 ! Fit after self-cal, rms=0.101568Jy sigma=0.683800 ! Inverting map uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=0.9479 mas, bmaj=1.419 mas, bpa=-0.6171 degrees ! Estimated noise=1.14128 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.00116575 Jy ! Component: 100 - total flux cleaned = -0.0015877 Jy ! Component: 150 - total flux cleaned = -0.00182204 Jy ! Component: 200 - total flux cleaned = -0.00199602 Jy ! Component: 250 - total flux cleaned = -0.00211841 Jy ! Component: 300 - total flux cleaned = -0.00223521 Jy ! Component: 350 - total flux cleaned = -0.00230961 Jy ! Component: 400 - total flux cleaned = -0.00238137 Jy ! Total flux subtracted in 400 components = -0.00238137 Jy ! Clean residual min=-0.002846 max=0.002859 Jy/beam ! Clean residual mean=-0.000021 rms=0.000686 Jy/beam ! Combined flux in latest and established models = 0.673985 Jy ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 83 components and 0.673985 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101468Jy sigma=0.682800 ! Fit after self-cal, rms=0.101471Jy sigma=0.682730 ! Inverting map ! Added new window around map position (48.4, 41.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00231246 Jy ! Component: 100 - total flux cleaned = 0.00217735 Jy ! Component: 150 - total flux cleaned = 0.00223558 Jy ! Component: 200 - total flux cleaned = 0.00218046 Jy ! Component: 250 - total flux cleaned = 0.00218198 Jy ! Component: 300 - total flux cleaned = 0.00223102 Jy ! Component: 350 - total flux cleaned = 0.00222976 Jy ! Component: 400 - total flux cleaned = 0.00222889 Jy ! Total flux subtracted in 400 components = 0.00222889 Jy ! Clean residual min=-0.002640 max=0.002819 Jy/beam ! Clean residual mean=-0.000021 rms=0.000661 Jy/beam ! Combined flux in latest and established models = 0.676214 Jy ! Performing phase self-cal ! Adding 28 model components to the UV plane model. ! The established model now contains 108 components and 0.676214 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101291Jy sigma=0.681096 ! Fit after self-cal, rms=0.101289Jy sigma=0.680978 ! Inverting map ! Added new window around map position (-17.8, -37). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000887734 Jy ! Component: 100 - total flux cleaned = 0.00103889 Jy ! Component: 150 - total flux cleaned = 0.00115669 Jy ! Component: 200 - total flux cleaned = 0.00126697 Jy ! Component: 250 - total flux cleaned = 0.00137029 Jy ! Component: 300 - total flux cleaned = 0.00127663 Jy ! Component: 350 - total flux cleaned = 0.00123294 Jy ! Component: 400 - total flux cleaned = 0.00123341 Jy ! Total flux subtracted in 400 components = 0.00123341 Jy ! Clean residual min=-0.002380 max=0.002759 Jy/beam ! Clean residual mean=-0.000023 rms=0.000629 Jy/beam ! Combined flux in latest and established models = 0.677447 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 138 components and 0.677447 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.101127Jy sigma=0.679170 ! Fit after self-cal, rms=0.101126Jy sigma=0.679083 ! Inverting map ! Added new window around map position (-25.4, 21.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000307209 Jy ! Component: 100 - total flux cleaned = -0.000310782 Jy ! Component: 150 - total flux cleaned = -0.000495678 Jy ! Component: 200 - total flux cleaned = -0.000608333 Jy ! Component: 250 - total flux cleaned = -0.000557734 Jy ! Component: 300 - total flux cleaned = -0.000605177 Jy ! Component: 350 - total flux cleaned = -0.000651473 Jy ! Component: 400 - total flux cleaned = -0.000694305 Jy ! Total flux subtracted in 400 components = -0.000694305 Jy ! Clean residual min=-0.002418 max=0.002494 Jy/beam ! Clean residual mean=-0.000022 rms=0.000591 Jy/beam ! Combined flux in latest and established models = 0.676753 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 167 components and 0.676753 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100962Jy sigma=0.677245 ! Fit after self-cal, rms=0.100967Jy sigma=0.677165 ! Inverting map ! Added new window around map position (-50.8, -26.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00111034 Jy ! Component: 100 - total flux cleaned = 0.000903146 Jy ! Component: 150 - total flux cleaned = 0.000903073 Jy ! Component: 200 - total flux cleaned = 0.000902149 Jy ! Component: 250 - total flux cleaned = 0.000944164 Jy ! Component: 300 - total flux cleaned = 0.000900217 Jy ! Component: 350 - total flux cleaned = 0.000938798 Jy ! Component: 400 - total flux cleaned = 0.000818498 Jy ! Total flux subtracted in 400 components = 0.000818498 Jy ! Clean residual min=-0.002425 max=0.002354 Jy/beam ! Clean residual mean=-0.000024 rms=0.000571 Jy/beam ! Combined flux in latest and established models = 0.677571 Jy ! Performing phase self-cal ! Adding 64 model components to the UV plane model. ! The established model now contains 195 components and 0.677571 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100862Jy sigma=0.676003 ! Fit after self-cal, rms=0.100856Jy sigma=0.675965 ! 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 = 6.11979e-07 Jy ! Component: 100 - total flux cleaned = 9.79771e-07 Jy ! Component: 150 - total flux cleaned = -7.66959e-05 Jy ! Component: 200 - total flux cleaned = -7.61184e-05 Jy ! Component: 250 - total flux cleaned = -0.000113886 Jy ! Component: 300 - total flux cleaned = -0.000150363 Jy ! Component: 350 - total flux cleaned = -4.22421e-05 Jy ! Component: 400 - total flux cleaned = -0.000114767 Jy ! Total flux subtracted in 400 components = -0.000114767 Jy ! Clean residual min=-0.002371 max=0.002283 Jy/beam ! Clean residual mean=-0.000023 rms=0.000563 Jy/beam ! Combined flux in latest and established models = 0.677456 Jy ! Performing phase self-cal ! Adding 80 model components to the UV plane model. ! The established model now contains 211 components and 0.677456 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100795Jy sigma=0.675392 ! Fit after self-cal, rms=0.100795Jy sigma=0.675386 ! 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 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 2. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 1 telescope corrections were flagged in sub-array 2. ! A total of 3 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 3. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Correcting IF 4. ! A total of 2 telescope corrections were flagged in sub-array 1. ! A total of 8 telescope corrections were flagged in sub-array 2. ! A total of 6 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 2 telescope corrections were flagged in sub-array 5. ! ! Fit before self-cal, rms=0.100795Jy sigma=0.675386 ! Fit after self-cal, rms=0.100883Jy sigma=0.674849 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9486 mas, bmaj=1.42 mas, bpa=-0.5615 degrees ! Estimated noise=1.14082 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000368857 Jy ! Component: 100 - total flux cleaned = 0.000578022 Jy ! Component: 150 - total flux cleaned = 0.000658095 Jy ! Component: 200 - total flux cleaned = 0.00073546 Jy ! Component: 250 - total flux cleaned = 0.0008101 Jy ! Component: 300 - total flux cleaned = 0.000846327 Jy ! Component: 350 - total flux cleaned = 0.000846235 Jy ! Component: 400 - total flux cleaned = 0.000882063 Jy ! Total flux subtracted in 400 components = 0.000882063 Jy ! Clean residual min=-0.002295 max=0.002291 Jy/beam ! Clean residual mean=-0.000026 rms=0.000544 Jy/beam ! Combined flux in latest and established models = 0.678338 Jy selfcal ! Performing phase self-cal ! Adding 66 model components to the UV plane model. ! The established model now contains 225 components and 0.678338 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100821Jy sigma=0.674222 ! Fit after self-cal, rms=0.100822Jy sigma=0.674207 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 = -3.9663e-05 Jy ! Component: 100 - total flux cleaned = 0.000102317 Jy ! Component: 150 - total flux cleaned = 0.000102223 Jy ! Component: 200 - total flux cleaned = 0.000169594 Jy ! Component: 250 - total flux cleaned = 0.00013572 Jy ! Component: 300 - total flux cleaned = 0.000168947 Jy ! Component: 350 - total flux cleaned = 0.00023395 Jy ! Component: 400 - total flux cleaned = 0.000299108 Jy ! Total flux subtracted in 400 components = 0.000299108 Jy ! Clean residual min=-0.002278 max=0.002238 Jy/beam ! Clean residual mean=-0.000026 rms=0.000538 Jy/beam ! Combined flux in latest and established models = 0.678637 Jy ! Performing phase self-cal ! Adding 73 model components to the UV plane model. ! The established model now contains 227 components and 0.678637 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100775Jy sigma=0.673755 ! Fit after self-cal, rms=0.100777Jy sigma=0.673750 ! 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.100777Jy sigma=0.673750 ! Fit after self-cal, rms=0.100804Jy sigma=0.673700 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.9487 mas, bmaj=1.42 mas, bpa=-0.5727 degrees ! Estimated noise=1.14056 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000148049 Jy ! Component: 100 - total flux cleaned = 0.000352739 Jy ! Component: 150 - total flux cleaned = 0.000452005 Jy ! Component: 200 - total flux cleaned = 0.000484622 Jy ! Component: 250 - total flux cleaned = 0.000580619 Jy ! Component: 300 - total flux cleaned = 0.000643963 Jy ! Component: 350 - total flux cleaned = 0.000706339 Jy ! Component: 400 - total flux cleaned = 0.000738069 Jy ! Total flux subtracted in 400 components = 0.000738069 Jy ! Clean residual min=-0.002261 max=0.002230 Jy/beam ! Clean residual mean=-0.000026 rms=0.000533 Jy/beam ! Combined flux in latest and established models = 0.679375 Jy selfcal ! Performing phase self-cal ! Adding 84 model components to the UV plane model. ! The established model now contains 236 components and 0.679375 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100763Jy sigma=0.673275 ! Fit after self-cal, rms=0.100764Jy sigma=0.673270 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 = 5.89854e-05 Jy ! Component: 100 - total flux cleaned = 5.96828e-05 Jy ! Component: 150 - total flux cleaned = 0.000120394 Jy ! Component: 200 - total flux cleaned = 0.000180267 Jy ! Component: 250 - total flux cleaned = 0.000209947 Jy ! Component: 300 - total flux cleaned = 0.000239497 Jy ! Component: 350 - total flux cleaned = 0.000269005 Jy ! Component: 400 - total flux cleaned = 0.000268675 Jy ! Total flux subtracted in 400 components = 0.000268675 Jy ! Clean residual min=-0.002239 max=0.002216 Jy/beam ! Clean residual mean=-0.000026 rms=0.000528 Jy/beam ! Combined flux in latest and established models = 0.679644 Jy ! Performing phase self-cal ! Adding 93 model components to the UV plane model. ! The established model now contains 243 components and 0.679644 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.100726Jy sigma=0.672925 ! Fit after self-cal, rms=0.100728Jy sigma=0.672922 ! 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.100728Jy sigma=0.672922 ! Fit after self-cal, rms=0.098047Jy sigma=0.657327 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.552 mas, bmaj=2.79 mas, bpa=46.8 degrees ! Estimated noise=1.3233 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000256939 Jy ! Component: 100 - total flux cleaned = -0.00057874 Jy ! Component: 150 - total flux cleaned = -0.00081333 Jy ! Component: 200 - total flux cleaned = -0.000983764 Jy ! Component: 250 - total flux cleaned = -0.001102 Jy ! Component: 300 - total flux cleaned = -0.00113932 Jy ! Component: 350 - total flux cleaned = -0.00117593 Jy ! Component: 400 - total flux cleaned = -0.00120984 Jy ! Total flux subtracted in 400 components = -0.00120984 Jy ! Clean residual min=-0.002417 max=0.002027 Jy/beam ! Clean residual mean=-0.000051 rms=0.000553 Jy/beam ! Combined flux in latest and established models = 0.678434 Jy selfcal ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 270 components and 0.678434 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098040Jy sigma=0.657182 ! Fit after self-cal, rms=0.098038Jy sigma=0.657162 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 = -5.28217e-06 Jy ! Component: 100 - total flux cleaned = 1.17044e-05 Jy ! Component: 150 - total flux cleaned = 4.35449e-05 Jy ! Component: 200 - total flux cleaned = 2.75799e-05 Jy ! Component: 250 - total flux cleaned = 8.90026e-05 Jy ! Component: 300 - total flux cleaned = 0.000118104 Jy ! Component: 350 - total flux cleaned = 0.000131835 Jy ! Component: 400 - total flux cleaned = 0.00017329 Jy ! Total flux subtracted in 400 components = 0.00017329 Jy ! Clean residual min=-0.002381 max=0.002059 Jy/beam ! Clean residual mean=-0.000050 rms=0.000550 Jy/beam ! Combined flux in latest and established models = 0.678607 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 283 components and 0.678607 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098036Jy sigma=0.657092 ! Fit after self-cal, rms=0.098037Jy sigma=0.657088 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "=== Finished amplitude self-cal for all stations ===" ! === Finished amplitude self-cal for all stations === restore ! Inverting map and beam ! Estimated beam: bmin=0.9489 mas, bmaj=1.42 mas, bpa=-0.5837 degrees ! Estimated noise=1.14079 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.9489 x 1.42 at -0.5837 degrees (North through East) ! Clean map min=-0.003076 max=0.61848 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.098037Jy sigma=0.657088 ! Fit after self-cal, rms=0.093389Jy sigma=0.608215 ! Inverting map and beam ! Estimated beam: bmin=0.9489 mas, bmaj=1.421 mas, bpa=-0.5764 degrees ! Estimated noise=1.14004 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -5.83031e-05 Jy ! Component: 100 - total flux cleaned = -9.49284e-05 Jy ! Component: 150 - total flux cleaned = -0.000196496 Jy ! Component: 200 - total flux cleaned = -0.000197351 Jy ! Component: 250 - total flux cleaned = -0.00025945 Jy ! Component: 300 - total flux cleaned = -0.000229571 Jy ! Component: 350 - total flux cleaned = -0.000288862 Jy ! Component: 400 - total flux cleaned = -0.000260144 Jy ! Total flux subtracted in 400 components = -0.000260144 Jy ! Clean residual min=-0.002252 max=0.002257 Jy/beam ! Clean residual mean=-0.000025 rms=0.000512 Jy/beam ! Combined flux in latest and established models = 0.678347 Jy ! Performing phase self-cal ! Adding 54 model components to the UV plane model. ! The established model now contains 305 components and 0.678347 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093347Jy sigma=0.607759 ! Fit after self-cal, rms=0.093349Jy sigma=0.607752 ! Inverting map ! Added new window around map position (11, -36.8). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00131004 Jy ! Component: 100 - total flux cleaned = 0.0014334 Jy ! Component: 150 - total flux cleaned = 0.00143432 Jy ! Component: 200 - total flux cleaned = 0.00148305 Jy ! Component: 250 - total flux cleaned = 0.00152669 Jy ! Component: 300 - total flux cleaned = 0.00144064 Jy ! Component: 350 - total flux cleaned = 0.001401 Jy ! Component: 400 - total flux cleaned = 0.00140126 Jy ! Total flux subtracted in 400 components = 0.00140126 Jy ! Clean residual min=-0.002106 max=0.002103 Jy/beam ! Clean residual mean=-0.000026 rms=0.000493 Jy/beam ! Combined flux in latest and established models = 0.679748 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 335 components and 0.679748 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093170Jy sigma=0.606553 ! Fit after self-cal, rms=0.093173Jy sigma=0.606515 ! 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.093173Jy sigma=0.606515 ! Fit after self-cal, rms=0.093110Jy sigma=0.606452 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=0.7465 mas, bmaj=1.112 mas, bpa=-1.445 degrees ! Estimated noise=1.83035 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.479232 Jy ! Component: 100 - total flux cleaned = 0.588154 Jy ! Component: 150 - total flux cleaned = 0.626845 Jy ! Component: 200 - total flux cleaned = 0.644784 Jy ! Component: 250 - total flux cleaned = 0.655108 Jy ! Component: 300 - total flux cleaned = 0.661839 Jy ! Component: 350 - total flux cleaned = 0.665186 Jy ! Component: 400 - total flux cleaned = 0.66829 Jy ! Total flux subtracted in 400 components = 0.66829 Jy ! Clean residual min=-0.005409 max=0.004991 Jy/beam ! Clean residual mean=-0.000004 rms=0.001216 Jy/beam ! Combined flux in latest and established models = 0.66829 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 45 model components to the UV plane model. ! The established model now contains 45 components and 0.66829 Jy ! Inverting map and beam ! Estimated beam: bmin=0.949 mas, bmaj=1.421 mas, bpa=-0.5638 degrees ! Estimated noise=1.1396 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.00513955 Jy ! Component: 100 - total flux cleaned = 0.00760969 Jy ! Component: 150 - total flux cleaned = 0.00925619 Jy ! Component: 200 - total flux cleaned = 0.0103061 Jy ! Component: 250 - total flux cleaned = 0.011188 Jy ! Component: 300 - total flux cleaned = 0.0120193 Jy ! Component: 350 - total flux cleaned = 0.0123683 Jy ! Component: 400 - total flux cleaned = 0.0124503 Jy ! Total flux subtracted in 400 components = 0.0124503 Jy ! Clean residual min=-0.002662 max=0.002788 Jy/beam ! Clean residual mean=-0.000028 rms=0.000630 Jy/beam ! Combined flux in latest and established models = 0.680741 Jy ! Adding 48 model components to the UV plane model. ! The established model now contains 88 components and 0.680741 Jy ! Inverting map ! Added new window around map position (-4.4, 39.2). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00126326 Jy ! Component: 100 - total flux cleaned = 0.00158128 Jy ! Component: 150 - total flux cleaned = 0.00165551 Jy ! Component: 200 - total flux cleaned = 0.00180502 Jy ! Component: 250 - total flux cleaned = 0.00194975 Jy ! Component: 300 - total flux cleaned = 0.00215857 Jy ! Component: 350 - total flux cleaned = 0.00202348 Jy ! Component: 400 - total flux cleaned = 0.00215427 Jy ! Total flux subtracted in 400 components = 0.00215427 Jy ! Clean residual min=-0.002500 max=0.002483 Jy/beam ! Clean residual mean=-0.000029 rms=0.000585 Jy/beam ! Combined flux in latest and established models = 0.682895 Jy ! Adding 66 model components to the UV plane model. ! The established model now contains 113 components and 0.682895 Jy ! Inverting map uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.552 mas, bmaj=2.789 mas, bpa=46.82 degrees ! Estimated noise=1.32152 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000244765 Jy ! Component: 100 - total flux cleaned = -0.000627582 Jy ! Component: 150 - total flux cleaned = -0.000876358 Jy ! Component: 200 - total flux cleaned = -0.00100898 Jy ! Component: 250 - total flux cleaned = -0.00121297 Jy ! Component: 300 - total flux cleaned = -0.00125063 Jy ! Component: 350 - total flux cleaned = -0.00128714 Jy ! Component: 400 - total flux cleaned = -0.00125189 Jy ! Total flux subtracted in 400 components = -0.00125189 Jy ! Clean residual min=-0.002137 max=0.001913 Jy/beam ! Clean residual mean=-0.000057 rms=0.000552 Jy/beam ! Combined flux in latest and established models = 0.681643 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 = 1.04277e-06 Jy ! Component: 100 - total flux cleaned = 3.33709e-05 Jy ! Component: 150 - total flux cleaned = 3.32898e-05 Jy ! Component: 200 - total flux cleaned = 9.27188e-05 Jy ! Component: 250 - total flux cleaned = 6.28743e-05 Jy ! Component: 300 - total flux cleaned = 6.24117e-05 Jy ! Component: 350 - total flux cleaned = -2.10406e-05 Jy ! Component: 400 - total flux cleaned = -2.00375e-05 Jy ! Total flux subtracted in 400 components = -2.00375e-05 Jy ! Clean residual min=-0.002138 max=0.001921 Jy/beam ! Clean residual mean=-0.000057 rms=0.000538 Jy/beam ! Combined flux in latest and established models = 0.681623 Jy ! Adding 68 model components to the UV plane model. ! The established model now contains 176 components and 0.681623 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=0.949 mas, bmaj=1.421 mas, bpa=-0.5638 degrees ! Estimated noise=1.1396 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.8227e-06 Jy ! Component: 100 - total flux cleaned = 3.64577e-06 Jy ! Component: 150 - total flux cleaned = 6.36266e-05 Jy ! Component: 200 - total flux cleaned = -4.981e-05 Jy ! Component: 250 - total flux cleaned = 0.000173004 Jy ! Component: 300 - total flux cleaned = 0.000224934 Jy ! Component: 350 - total flux cleaned = 0.000278477 Jy ! Component: 400 - total flux cleaned = 0.000330622 Jy ! Total flux subtracted in 400 components = 0.000330622 Jy ! Clean residual min=-0.002200 max=0.002304 Jy/beam ! Clean residual mean=-0.000028 rms=0.000534 Jy/beam ! Combined flux in latest and established models = 0.681954 Jy selfcal ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 213 components and 0.681954 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.093566Jy sigma=0.610312 ! Fit after self-cal, rms=0.093526Jy sigma=0.609994 wmodel /scr/pima/rdv94_uvs/J1521+4336_X_map.mod ! Writing 213 model components to file: /scr/pima/rdv94_uvs/J1521+4336_X_map.mod wobs /scr/pima/rdv94_uvs/J1521+4336_X_uvs.fits ! Applying 2401 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J1521+4336_X_uvs.fits wwins /scr/pima/rdv94_uvs/J1521+4336_X_map.win ! wwins: Wrote 8 windows to /scr/pima/rdv94_uvs/J1521+4336_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 = -7.26651e-05 Jy ! Component: 100 - total flux cleaned = -0.000231968 Jy ! Component: 150 - total flux cleaned = -0.000385851 Jy ! Component: 200 - total flux cleaned = -0.000499149 Jy ! Component: 250 - total flux cleaned = -0.000824872 Jy ! Component: 300 - total flux cleaned = -0.00103649 Jy ! Component: 350 - total flux cleaned = -0.00120966 Jy ! Component: 400 - total flux cleaned = -0.00134619 Jy ! Component: 450 - total flux cleaned = -0.00148082 Jy ! Component: 500 - total flux cleaned = -0.00158066 Jy ! Component: 550 - total flux cleaned = -0.00171175 Jy ! Component: 600 - total flux cleaned = -0.00180856 Jy ! Component: 650 - total flux cleaned = -0.00193698 Jy ! Component: 700 - total flux cleaned = -0.00193663 Jy ! Component: 750 - total flux cleaned = -0.00199962 Jy ! Component: 800 - total flux cleaned = -0.00212365 Jy ! Component: 850 - total flux cleaned = -0.00212357 Jy ! Component: 900 - total flux cleaned = -0.00221511 Jy ! Component: 950 - total flux cleaned = -0.0022451 Jy ! Component: 1000 - total flux cleaned = -0.00233469 Jy ! Component: 1050 - total flux cleaned = -0.00242369 Jy ! Component: 1100 - total flux cleaned = -0.00257103 Jy ! Component: 1150 - total flux cleaned = -0.00257095 Jy ! Component: 1200 - total flux cleaned = -0.00268637 Jy ! Component: 1250 - total flux cleaned = -0.00277239 Jy ! Component: 1300 - total flux cleaned = -0.002687 Jy ! Component: 1350 - total flux cleaned = -0.00285667 Jy ! Component: 1400 - total flux cleaned = -0.00285675 Jy ! Component: 1450 - total flux cleaned = -0.00291224 Jy ! Component: 1500 - total flux cleaned = -0.00299492 Jy ! Component: 1550 - total flux cleaned = -0.00304969 Jy ! Component: 1600 - total flux cleaned = -0.00307688 Jy ! Component: 1650 - total flux cleaned = -0.0030769 Jy ! Component: 1700 - total flux cleaned = -0.00313051 Jy ! Component: 1750 - total flux cleaned = -0.00305053 Jy ! Component: 1800 - total flux cleaned = -0.0031302 Jy ! Component: 1850 - total flux cleaned = -0.00307743 Jy ! Component: 1900 - total flux cleaned = -0.00307734 Jy ! Component: 1950 - total flux cleaned = -0.00312904 Jy ! Component: 2000 - total flux cleaned = -0.00307707 Jy ! Component: 2050 - total flux cleaned = -0.00302602 Jy ! Component: 2100 - total flux cleaned = -0.0031022 Jy ! Component: 2150 - total flux cleaned = -0.00300061 Jy ! Component: 2200 - total flux cleaned = -0.00305101 Jy ! Component: 2250 - total flux cleaned = -0.00300121 Jy ! Component: 2300 - total flux cleaned = -0.00302605 Jy ! Component: 2350 - total flux cleaned = -0.00305094 Jy ! Component: 2400 - total flux cleaned = -0.00317395 Jy ! Component: 2450 - total flux cleaned = -0.0031251 Jy ! Component: 2500 - total flux cleaned = -0.00314958 Jy ! Component: 2550 - total flux cleaned = -0.00322209 Jy ! Component: 2600 - total flux cleaned = -0.00322228 Jy ! Component: 2650 - total flux cleaned = -0.00324633 Jy ! Component: 2700 - total flux cleaned = -0.00324641 Jy ! Component: 2750 - total flux cleaned = -0.00329383 Jy ! Component: 2800 - total flux cleaned = -0.00324661 Jy ! Component: 2850 - total flux cleaned = -0.0033171 Jy ! Component: 2900 - total flux cleaned = -0.00345706 Jy ! Component: 2950 - total flux cleaned = -0.00343391 Jy ! Component: 3000 - total flux cleaned = -0.0034107 Jy ! Component: 3050 - total flux cleaned = -0.00343359 Jy ! Component: 3100 - total flux cleaned = -0.00345658 Jy ! Component: 3150 - total flux cleaned = -0.00359319 Jy ! Component: 3200 - total flux cleaned = -0.00347986 Jy ! Component: 3250 - total flux cleaned = -0.00343471 Jy ! Component: 3300 - total flux cleaned = -0.00350208 Jy ! Component: 3350 - total flux cleaned = -0.00354691 Jy ! Component: 3400 - total flux cleaned = -0.00370266 Jy ! Component: 3450 - total flux cleaned = -0.00365847 Jy ! Component: 3500 - total flux cleaned = -0.00363635 Jy ! Component: 3550 - total flux cleaned = -0.00381217 Jy ! Component: 3600 - total flux cleaned = -0.00379017 Jy ! Component: 3650 - total flux cleaned = -0.00383382 Jy ! Component: 3700 - total flux cleaned = -0.00398549 Jy ! Component: 3750 - total flux cleaned = -0.00402855 Jy ! Component: 3800 - total flux cleaned = -0.00405007 Jy ! Component: 3850 - total flux cleaned = -0.00413571 Jy ! Component: 3900 - total flux cleaned = -0.00413571 Jy ! Component: 3950 - total flux cleaned = -0.00424184 Jy ! Component: 4000 - total flux cleaned = -0.00432628 Jy ! Component: 4050 - total flux cleaned = -0.00434732 Jy ! Total flux subtracted in 4096 components = -0.00445223 Jy ! Clean residual min=-0.001328 max=0.001150 Jy/beam ! Clean residual mean=-0.000022 rms=0.000340 Jy/beam ! Combined flux in latest and established models = 0.677501 Jy keep ! Adding 281 model components to the UV plane model. ! The established model now contains 494 components and 0.677501 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J1521+4336_X\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.949 x 1.421 at -0.5638 degrees (North through East) ! Clean map min=-0.0025337 max=0.61728 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J1521+4336_X_map.fits float peak_flux peak_flux = peak(flux,max) invert ! Inverting map float image_peak image_peak = int(10000*peak_flux)/10000 float image_rms image_rms = int(100000*imstat(rms))/100000 print "=========================================" ! ========================================= print "The image peak= ", image_peak, " rms= ", image_rms, " Jy" ! The image peak= 0.6172 rms= 0.00033 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J1521+4336_X_uvs.log closed on Tue Oct 8 12:21:47 2013