! Started logfile: /scr/pima/rdv94_uvs/J0856-1105_S_uvs.log on Tue Oct 8 12:16:20 2013 obs /scr/pima/rdv94_uvs/J0856-1105_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv94_uvs/J0856-1105_S_uva.fits ! AN table 1: 6 integrations on 45 of 45 possible baselines. ! AN table 2: 2 integrations on 15 of 15 possible baselines. ! AN table 3: 5 integrations on 1 of 1 possible baselines. ! AN table 4: 80 integrations on 3 of 3 possible baselines. ! AN table 5: 10 integrations on 10 of 10 possible baselines. ! AN table 6: 10 integrations on 15 of 15 possible baselines. ! AN table 7: 20 integrations on 6 of 6 possible baselines. ! AN table 8: 12 integrations on 45 of 45 possible baselines. ! AN table 9: 3 integrations on 3 of 3 possible baselines. ! Apparent sampling: 0.678962 visibilities/baseline/integration-bin. ! Found source: J0856-1105 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 1.6e+07 1 1.6e+07 ! 02 2 2.26299e+09 1.6e+07 1 1.6e+07 ! 03 3 2.35299e+09 1.6e+07 1 1.6e+07 ! 04 4 2.37299e+09 1.6e+07 1 1.6e+07 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 3976 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.5 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv94_uvs/J0856-1105_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:MK 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:WF 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:OV self-cal status=correctable weight=1 ! 3:HH self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:ON self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:ON self-cal status=correctable weight=1 ! 5:WF self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:ON self-cal status=correctable weight=1 ! 6:WF self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:HH self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WF self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HH self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:ON self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:WF self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:OV self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 90 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 11 telescope corrections were flagged in sub-array 7. ! A total of 57 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 90 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 11 telescope corrections were flagged in sub-array 7. ! A total of 57 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 90 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 11 telescope corrections were flagged in sub-array 7. ! A total of 57 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 90 telescope corrections were flagged in sub-array 4. ! A total of 10 telescope corrections were flagged in sub-array 5. ! A total of 17 telescope corrections were flagged in sub-array 6. ! A total of 11 telescope corrections were flagged in sub-array 7. ! A total of 57 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=1.192294Jy sigma=6.068635 ! Fit after self-cal, rms=0.625744Jy sigma=3.708746 ! 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:WF ! Flagging baselines of antenna 1:WF in the currently selected channels. flag 3:HH ! Flagging baselines of antenna 3:HH in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 4:HH ! Flagging baselines of antenna 4:HH in the currently selected channels. flag 4:ON ! Flagging baselines of antenna 4:ON in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 5:HH ! Flagging baselines of antenna 5:HH in the currently selected channels. flag 5:ON ! Flagging baselines of antenna 5:ON in the currently selected channels. flag 5:WF ! Flagging baselines of antenna 5:WF in the currently selected channels. flag 5:WZ ! Flagging baselines of antenna 5:WZ in the currently selected channels. flag 6:HH ! Flagging baselines of antenna 6:HH in the currently selected channels. flag 6:ON ! Flagging baselines of antenna 6:ON in the currently selected channels. flag 6:WF ! Flagging baselines of antenna 6:WF in the currently selected channels. flag 6:WZ ! Flagging baselines of antenna 6:WZ in the currently selected channels. flag 7:HH ! Flagging baselines of antenna 7:HH in the currently selected channels. flag 7:ON ! Flagging baselines of antenna 7:ON in the currently selected channels. flag 7:WF ! Flagging baselines of antenna 7:WF in the currently selected channels. flag 7:WZ ! Flagging baselines of antenna 7:WZ in the currently selected channels. flag 8:HH ! Flagging baselines of antenna 8:HH in the currently selected channels. flag 8:ON ! Flagging baselines of antenna 8:ON in the currently selected channels. flag 8:WF ! Flagging baselines of antenna 8:WF in the currently selected channels. flag 8:WZ ! Flagging baselines of antenna 8: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=3.077 mas, bmaj=12.84 mas, bpa=-17.3 degrees ! Estimated noise=2.94144 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.257254 Jy ! Component: 100 - total flux cleaned = 0.313421 Jy ! Total flux subtracted in 100 components = 0.313421 Jy ! Clean residual min=-0.009045 max=0.016393 Jy/beam ! Clean residual mean=-0.000014 rms=0.002880 Jy/beam ! Combined flux in latest and established models = 0.313421 Jy ! Performing phase self-cal ! Adding 4 model components to the UV plane model. ! The established model now contains 4 components and 0.313421 Jy ! ! Correcting IF 1. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 150 telescope corrections were flagged in sub-array 4. ! A total of 40 telescope corrections were flagged in sub-array 5. ! A total of 43 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 31 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 150 telescope corrections were flagged in sub-array 4. ! A total of 40 telescope corrections were flagged in sub-array 5. ! A total of 43 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 31 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 150 telescope corrections were flagged in sub-array 4. ! A total of 40 telescope corrections were flagged in sub-array 5. ! A total of 43 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 31 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 6 telescope corrections were flagged in sub-array 1. ! A total of 150 telescope corrections were flagged in sub-array 4. ! A total of 40 telescope corrections were flagged in sub-array 5. ! A total of 43 telescope corrections were flagged in sub-array 6. ! A total of 69 telescope corrections were flagged in sub-array 7. ! A total of 31 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.080358Jy sigma=0.671077 ! Fit after self-cal, rms=0.079191Jy sigma=0.672575 ! Inverting map and beam ! Estimated beam: bmin=3.075 mas, bmaj=12.83 mas, bpa=-17.3 degrees ! Estimated noise=2.94598 mJy/beam. print "=== Finished uniform weighting clean ===" ! === Finished uniform weighting clean === uvw 0,-1 ! Uniform weighting is not currently selected. ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. win_mult = win_mult * 1.6 clean_niter = clean_niter * 2 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.95 mas, bmaj=14.18 mas, bpa=-16.62 degrees ! Estimated noise=2.32592 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.0214449 Jy ! Component: 100 - total flux cleaned = 0.0312334 Jy ! Component: 150 - total flux cleaned = 0.0334052 Jy ! Component: 200 - total flux cleaned = 0.0343445 Jy ! Total flux subtracted in 200 components = 0.0343445 Jy ! Clean residual min=-0.009895 max=0.008651 Jy/beam ! Clean residual mean=-0.000103 rms=0.002428 Jy/beam ! Combined flux in latest and established models = 0.347765 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 38 components and 0.347765 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.073959Jy sigma=0.621384 ! Fit after self-cal, rms=0.073936Jy sigma=0.621178 ! 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.94 FD 1.07 HN 1.07 KP 1.15 ! LA 1.22 MK 0.91 NL 0.92 OV 0.95 ! PT 1.02 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.90 FD 1.16 KP 1.14 LA 1.29 ! MK 0.89 OV 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* HN 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! HH 0.00* ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.10 HH 1.00* HN 1.14 KP 1.10 ! LA 1.16 NL 0.88 ON 1.00* OV 0.97 ! WF 1.00* WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! KP 0.00* MK 0.00* OV 0.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.93 FD 1.14 HN 1.17 KP 1.13 ! LA 1.12 MK 0.98 NL 0.97 OV 0.99 ! PT 1.02 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.94 FD 1.13 KP 1.11 LA 1.21 ! MK 0.87 OV 1.03 ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* HN 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! HH 0.00* ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.07 HH 1.00* HN 1.17 KP 1.16 ! LA 1.09 NL 0.94 ON 1.00* OV 1.02 ! WF 1.00* WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! KP 0.00* MK 0.00* OV 0.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.99 FD 1.15 HN 1.14 KP 1.00 ! LA 0.89 MK 0.99 NL 0.98 OV 0.98 ! PT 1.04 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.98 FD 1.17 KP 1.01 LA 0.87 ! MK 0.93 OV 1.06 ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* HN 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! HH 0.00* ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.09 HH 1.00* HN 1.11 KP 0.99 ! LA 0.87 NL 0.99 ON 1.00* OV 0.99 ! WF 1.00* WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! KP 0.00* MK 0.00* OV 0.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.97 FD 1.10 HN 0.98 KP 0.94 ! LA 0.73 MK 1.08 NL 0.96 OV 0.84 ! PT 0.91 WF 1.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.86 FD 1.17 KP 0.96 LA 0.74 ! MK 0.99 OV 0.85 ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* ON 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 5: ! HH 0.00* HN 0.00* ON 0.00* WF 0.00* ! WZ 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! HH 0.00* HN 0.00* NL 0.00* ON 0.00* ! WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 7: ! HH 0.00* ON 0.00* WF 0.00* WZ 0.00* ! ! A total of 3 telescope corrections were flagged in sub-array 8. ! Telescope amplitude corrections in sub-array 8: ! FD 1.06 HH 1.00* HN 0.99 KP 0.90 ! LA 0.73 NL 0.89 ON 1.00* OV 0.84 ! WF 1.00* WZ 1.00* ! ! A total of 9 telescope corrections were flagged in sub-array 9. ! Telescope amplitude corrections in sub-array 9: ! KP 0.00* MK 0.00* OV 0.00* ! ! ! Fit before self-cal, rms=0.073936Jy sigma=0.621178 ! Fit after self-cal, rms=0.044458Jy sigma=0.371719 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.984 mas, bmaj=14.11 mas, bpa=-16.26 degrees ! Estimated noise=2.3556 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.000478589 Jy ! Component: 100 - total flux cleaned = 0.000486797 Jy ! Component: 150 - total flux cleaned = 0.000448602 Jy ! Component: 200 - total flux cleaned = 0.000342689 Jy ! Total flux subtracted in 200 components = 0.000342689 Jy ! Clean residual min=-0.004484 max=0.004709 Jy/beam ! Clean residual mean=-0.000005 rms=0.001011 Jy/beam ! Combined flux in latest and established models = 0.348108 Jy ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 58 components and 0.348108 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044307Jy sigma=0.370964 ! Fit after self-cal, rms=0.044368Jy sigma=0.370613 ! 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.044368Jy sigma=0.370613 ! Fit after self-cal, rms=0.044312Jy sigma=0.370472 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.987 mas, bmaj=14.1 mas, bpa=-16.25 degrees ! Estimated noise=2.35161 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000147632 Jy ! Component: 100 - total flux cleaned = -0.000147098 Jy ! Component: 150 - total flux cleaned = -0.000210213 Jy ! Component: 200 - total flux cleaned = -0.000263233 Jy ! Total flux subtracted in 200 components = -0.000263233 Jy ! Clean residual min=-0.004559 max=0.004719 Jy/beam ! Clean residual mean=-0.000005 rms=0.001002 Jy/beam ! Combined flux in latest and established models = 0.347845 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 69 components and 0.347845 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.044240Jy sigma=0.370185 ! Fit after self-cal, rms=0.044240Jy sigma=0.370151 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) ! Added new window around map position (1, -73). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00176201 Jy ! Component: 100 - total flux cleaned = 0.00198674 Jy ! Component: 150 - total flux cleaned = 0.00230076 Jy ! Component: 200 - total flux cleaned = 0.00257953 Jy ! Total flux subtracted in 200 components = 0.00257953 Jy ! Clean residual min=-0.004189 max=0.003851 Jy/beam ! Clean residual mean=-0.000007 rms=0.000935 Jy/beam ! Combined flux in latest and established models = 0.350424 Jy ! Performing phase self-cal ! Adding 35 model components to the UV plane model. ! The established model now contains 103 components and 0.350424 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.043731Jy sigma=0.365698 ! Fit after self-cal, rms=0.043682Jy sigma=0.365333 ! 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.043682Jy sigma=0.365333 ! Fit after self-cal, rms=0.042255Jy sigma=0.346011 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.991 mas, bmaj=14.1 mas, bpa=-16.26 degrees ! Estimated noise=2.35121 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000251114 Jy ! Component: 100 - total flux cleaned = 0.000303445 Jy ! Component: 150 - total flux cleaned = 0.000395816 Jy ! Component: 200 - total flux cleaned = 0.000595014 Jy ! Total flux subtracted in 200 components = 0.000595014 Jy ! Clean residual min=-0.002684 max=0.002608 Jy/beam ! Clean residual mean=-0.000003 rms=0.000685 Jy/beam ! Combined flux in latest and established models = 0.351019 Jy selfcal ! Performing phase self-cal ! Adding 19 model components to the UV plane model. ! The established model now contains 120 components and 0.351019 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.042141Jy sigma=0.345186 ! Fit after self-cal, rms=0.042146Jy sigma=0.345154 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.39444e-05 Jy ! Component: 100 - total flux cleaned = 0.000191046 Jy ! Component: 150 - total flux cleaned = 0.000285695 Jy ! Component: 200 - total flux cleaned = 0.000373632 Jy ! Total flux subtracted in 200 components = 0.000373632 Jy ! Clean residual min=-0.002679 max=0.002535 Jy/beam ! Clean residual mean=-0.000003 rms=0.000680 Jy/beam ! Combined flux in latest and established models = 0.351393 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 132 components and 0.351393 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.042097Jy sigma=0.344813 ! Fit after self-cal, rms=0.042098Jy sigma=0.344802 ! 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.042098Jy sigma=0.344802 ! Fit after self-cal, rms=0.042058Jy sigma=0.344730 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.054 mas, bmaj=15.82 mas, bpa=-7.659 degrees ! Estimated noise=2.63289 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000283118 Jy ! Component: 100 - total flux cleaned = 0.000517235 Jy ! Component: 150 - total flux cleaned = 0.000644337 Jy ! Component: 200 - total flux cleaned = 0.000798964 Jy ! Total flux subtracted in 200 components = 0.000798964 Jy ! Clean residual min=-0.002830 max=0.002442 Jy/beam ! Clean residual mean=-0.000006 rms=0.000719 Jy/beam ! Combined flux in latest and established models = 0.352192 Jy selfcal ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 159 components and 0.352192 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.042038Jy sigma=0.344546 ! Fit after self-cal, rms=0.042040Jy sigma=0.344530 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.57068e-05 Jy ! Component: 100 - total flux cleaned = 7.51857e-05 Jy ! Component: 150 - total flux cleaned = 0.000147086 Jy ! Component: 200 - total flux cleaned = 0.000196858 Jy ! Total flux subtracted in 200 components = 0.000196858 Jy ! Clean residual min=-0.002873 max=0.002421 Jy/beam ! Clean residual mean=-0.000007 rms=0.000717 Jy/beam ! Combined flux in latest and established models = 0.352389 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 182 components and 0.352389 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.042033Jy sigma=0.344480 ! Fit after self-cal, rms=0.042035Jy sigma=0.344473 ! 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:MK,true ! 1:MK 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:OV,true ! 2:OV self-cal status=fixed weight=1 selfant 5:HN,true ! 5:HN self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:NL,true ! 6:NL self-cal status=fixed weight=1 selfant 8:FD,true ! 8:FD self-cal status=fixed weight=1 selfant 8:HN,true ! 8:HN self-cal status=fixed weight=1 selfant 8:KP,true ! 8:KP self-cal status=fixed weight=1 selfant 8:LA,true ! 8:LA self-cal status=fixed weight=1 selfant 8:NL,true ! 8:NL self-cal status=fixed weight=1 selfant 8:OV,true ! 8:OV self-cal status=fixed weight=1 selfant 9:KP,true ! 9:KP self-cal status=fixed weight=1 selfant 9:MK,true ! 9:MK self-cal status=fixed weight=1 selfant 9:OV,true ! 9:OV self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 20 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 20 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 20 telescope corrections were flagged in sub-array 8. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 20 telescope corrections were flagged in sub-array 5. ! A total of 16 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 20 telescope corrections were flagged in sub-array 8. ! ! Fit before self-cal, rms=0.395067Jy sigma=0.947458 ! Fit after self-cal, rms=0.097697Jy sigma=0.473652 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:MK 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:WF 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:OV self-cal status=correctable weight=1 ! 3:HH self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:ON self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 5:HH self-cal status=correctable weight=1 ! 5:HN self-cal status=correctable weight=1 ! 5:ON self-cal status=correctable weight=1 ! 5:WF self-cal status=correctable weight=1 ! 5:WZ self-cal status=correctable weight=1 ! 6:HH self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:ON self-cal status=correctable weight=1 ! 6:WF self-cal status=correctable weight=1 ! 6:WZ self-cal status=correctable weight=1 ! 7:HH self-cal status=correctable weight=1 ! 7:ON self-cal status=correctable weight=1 ! 7:WF self-cal status=correctable weight=1 ! 7:WZ self-cal status=correctable weight=1 ! 8:FD self-cal status=correctable weight=1 ! 8:HH self-cal status=correctable weight=1 ! 8:HN self-cal status=correctable weight=1 ! 8:KP self-cal status=correctable weight=1 ! 8:LA self-cal status=correctable weight=1 ! 8:NL self-cal status=correctable weight=1 ! 8:ON self-cal status=correctable weight=1 ! 8:OV self-cal status=correctable weight=1 ! 8:WF self-cal status=correctable weight=1 ! 8:WZ self-cal status=correctable weight=1 ! 9:KP self-cal status=correctable weight=1 ! 9:MK self-cal status=correctable weight=1 ! 9:OV 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 5 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 5 telescope corrections were flagged in sub-array 5. ! A total of 10 telescope corrections were flagged in sub-array 6. ! A total of 40 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.097697Jy sigma=0.473652 ! Fit after self-cal, rms=0.099698Jy sigma=0.464885 uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=2.123 mas, bmaj=3.585 mas, bpa=9.431 degrees ! Estimated noise=2.71504 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 (92, -26). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00325604 Jy ! Component: 100 - total flux cleaned = 0.00342412 Jy ! Component: 150 - total flux cleaned = 0.00402064 Jy ! Component: 200 - total flux cleaned = 0.00428486 Jy ! Total flux subtracted in 200 components = 0.00428486 Jy ! Clean residual min=-0.003842 max=0.004422 Jy/beam ! Clean residual mean=0.000006 rms=0.000975 Jy/beam ! Combined flux in latest and established models = 0.356674 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 212 components and 0.356674 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099034Jy sigma=0.462755 ! Fit after self-cal, rms=0.100203Jy sigma=0.456790 ! Inverting map ! Added new window around map position (-65, -75.5). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00086099 Jy ! Component: 100 - total flux cleaned = 0.000728774 Jy ! Component: 150 - total flux cleaned = 0.000483674 Jy ! Component: 200 - total flux cleaned = 0.000369677 Jy ! Total flux subtracted in 200 components = 0.000369677 Jy ! Clean residual min=-0.003580 max=0.004137 Jy/beam ! Clean residual mean=0.000005 rms=0.000915 Jy/beam ! Combined flux in latest and established models = 0.357043 Jy ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 230 components and 0.357043 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099840Jy sigma=0.455423 ! Fit after self-cal, rms=0.099835Jy sigma=0.455361 ! Inverting map ! Added new window around map position (-1.5, 55). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00251323 Jy ! Component: 100 - total flux cleaned = 0.00229591 Jy ! Component: 150 - total flux cleaned = 0.00199665 Jy ! Component: 200 - total flux cleaned = 0.00190367 Jy ! Total flux subtracted in 200 components = 0.00190367 Jy ! Clean residual min=-0.003476 max=0.003535 Jy/beam ! Clean residual mean=0.000003 rms=0.000878 Jy/beam ! Combined flux in latest and established models = 0.358947 Jy ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 248 components and 0.358947 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099597Jy sigma=0.454448 ! Fit after self-cal, rms=0.099585Jy sigma=0.454410 ! 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=2.41 mas, bmaj=3.903 mas, bpa=11.24 degrees ! Estimated noise=2.12535 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.000772435 Jy ! Component: 100 - total flux cleaned = 0.00132721 Jy ! Component: 150 - total flux cleaned = 0.00165546 Jy ! Component: 200 - total flux cleaned = 0.00188236 Jy ! Component: 250 - total flux cleaned = 0.00209281 Jy ! Component: 300 - total flux cleaned = 0.00235847 Jy ! Component: 350 - total flux cleaned = 0.00242089 Jy ! Component: 400 - total flux cleaned = 0.00248197 Jy ! Total flux subtracted in 400 components = 0.00248197 Jy ! Clean residual min=-0.003234 max=0.003096 Jy/beam ! Clean residual mean=-0.000002 rms=0.000751 Jy/beam ! Combined flux in latest and established models = 0.361429 Jy ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 293 components and 0.361429 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099232Jy sigma=0.453303 ! Fit after self-cal, rms=0.099248Jy sigma=0.453239 ! Inverting map dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000337486 Jy ! Component: 100 - total flux cleaned = 0.000518511 Jy ! Component: 150 - total flux cleaned = 0.000632475 Jy ! Component: 200 - total flux cleaned = 0.000688341 Jy ! Component: 250 - total flux cleaned = 0.000795893 Jy ! Component: 300 - total flux cleaned = 0.000899748 Jy ! Component: 350 - total flux cleaned = 0.000999795 Jy ! Component: 400 - total flux cleaned = 0.00109737 Jy ! Total flux subtracted in 400 components = 0.00109737 Jy ! Clean residual min=-0.003082 max=0.003039 Jy/beam ! Clean residual mean=-0.000015 rms=0.000732 Jy/beam ! Combined flux in latest and established models = 0.362526 Jy ! Performing phase self-cal ! Adding 67 model components to the UV plane model. ! The established model now contains 317 components and 0.362526 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.099069Jy sigma=0.452691 ! Fit after self-cal, rms=0.099073Jy sigma=0.452673 ! 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. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 25 telescope corrections were flagged in sub-array 5. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 60 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 25 telescope corrections were flagged in sub-array 5. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 60 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 25 telescope corrections were flagged in sub-array 5. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 60 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 3. ! A total of 240 telescope corrections were flagged in sub-array 4. ! A total of 25 telescope corrections were flagged in sub-array 5. ! A total of 26 telescope corrections were flagged in sub-array 6. ! A total of 32 telescope corrections were flagged in sub-array 7. ! A total of 60 telescope corrections were flagged in sub-array 8. ! A total of 9 telescope corrections were flagged in sub-array 9. ! ! Fit before self-cal, rms=0.295380Jy sigma=0.802189 ! Fit after self-cal, rms=0.098221Jy sigma=0.452042 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.404 mas, bmaj=3.896 mas, bpa=11.1 degrees ! Estimated noise=2.12156 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000682227 Jy ! Component: 100 - total flux cleaned = 0.000867907 Jy ! Component: 150 - total flux cleaned = 0.00103761 Jy ! Component: 200 - total flux cleaned = 0.0010903 Jy ! Component: 250 - total flux cleaned = 0.00109013 Jy ! Component: 300 - total flux cleaned = 0.00113913 Jy ! Component: 350 - total flux cleaned = 0.00104718 Jy ! Component: 400 - total flux cleaned = 0.00109313 Jy ! Total flux subtracted in 400 components = 0.00109313 Jy ! Clean residual min=-0.002913 max=0.002866 Jy/beam ! Clean residual mean=-0.000036 rms=0.000694 Jy/beam ! Combined flux in latest and established models = 0.363619 Jy selfcal ! Performing phase self-cal ! Adding 53 model components to the UV plane model. ! The established model now contains 335 components and 0.363619 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098044Jy sigma=0.451494 ! Fit after self-cal, rms=0.098053Jy sigma=0.451473 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.000141704 Jy ! Component: 100 - total flux cleaned = 0.00014163 Jy ! Component: 150 - total flux cleaned = 0.000141924 Jy ! Component: 200 - total flux cleaned = 0.000183257 Jy ! Component: 250 - total flux cleaned = 0.000222643 Jy ! Component: 300 - total flux cleaned = 0.000261923 Jy ! Component: 350 - total flux cleaned = 0.000300718 Jy ! Component: 400 - total flux cleaned = 0.000337414 Jy ! Total flux subtracted in 400 components = 0.000337414 Jy ! Clean residual min=-0.002888 max=0.002856 Jy/beam ! Clean residual mean=-0.000039 rms=0.000685 Jy/beam ! Combined flux in latest and established models = 0.363957 Jy ! Performing phase self-cal ! Adding 61 model components to the UV plane model. ! The established model now contains 341 components and 0.363957 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097956Jy sigma=0.451149 ! Fit after self-cal, rms=0.097962Jy sigma=0.451141 ! 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.097962Jy sigma=0.451141 ! Fit after self-cal, rms=0.097551Jy sigma=0.447841 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.401 mas, bmaj=3.892 mas, bpa=11.07 degrees ! Estimated noise=2.12025 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000410279 Jy ! Component: 100 - total flux cleaned = 0.000412272 Jy ! Component: 150 - total flux cleaned = 0.000455556 Jy ! Component: 200 - total flux cleaned = 0.000454984 Jy ! Component: 250 - total flux cleaned = 0.000417228 Jy ! Component: 300 - total flux cleaned = 0.000416296 Jy ! Component: 350 - total flux cleaned = 0.000417062 Jy ! Component: 400 - total flux cleaned = 0.00038268 Jy ! Total flux subtracted in 400 components = 0.00038268 Jy ! Clean residual min=-0.002831 max=0.002864 Jy/beam ! Clean residual mean=-0.000044 rms=0.000664 Jy/beam ! Combined flux in latest and established models = 0.364339 Jy selfcal ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 369 components and 0.364339 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097452Jy sigma=0.447522 ! Fit after self-cal, rms=0.097484Jy sigma=0.447490 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) ! Added new window around map position (126.5, -120.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00157555 Jy ! Component: 100 - total flux cleaned = 0.00171618 Jy ! Component: 150 - total flux cleaned = 0.00171953 Jy ! Component: 200 - total flux cleaned = 0.00172125 Jy ! Component: 250 - total flux cleaned = 0.00155192 Jy ! Component: 300 - total flux cleaned = 0.00147386 Jy ! Component: 350 - total flux cleaned = 0.00151078 Jy ! Component: 400 - total flux cleaned = 0.00158082 Jy ! Total flux subtracted in 400 components = 0.00158082 Jy ! Clean residual min=-0.002828 max=0.002691 Jy/beam ! Clean residual mean=-0.000049 rms=0.000645 Jy/beam ! Combined flux in latest and established models = 0.36592 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 392 components and 0.36592 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097288Jy sigma=0.446744 ! Fit after self-cal, rms=0.097281Jy sigma=0.446694 ! 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.097281Jy sigma=0.446694 ! Fit after self-cal, rms=0.097948Jy sigma=0.431303 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=8.01 mas, bmaj=16.67 mas, bpa=-7.037 degrees ! Estimated noise=2.50728 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000914609 Jy ! Component: 100 - total flux cleaned = 0.00160528 Jy ! Component: 150 - total flux cleaned = 0.00212606 Jy ! Component: 200 - total flux cleaned = 0.00256478 Jy ! Component: 250 - total flux cleaned = 0.00283195 Jy ! Component: 300 - total flux cleaned = 0.00304417 Jy ! Component: 350 - total flux cleaned = 0.00327668 Jy ! Component: 400 - total flux cleaned = 0.00346479 Jy ! Total flux subtracted in 400 components = 0.00346479 Jy ! Clean residual min=-0.002489 max=0.002075 Jy/beam ! Clean residual mean=-0.000008 rms=0.000648 Jy/beam ! Combined flux in latest and established models = 0.369385 Jy selfcal ! Performing phase self-cal ! Adding 55 model components to the UV plane model. ! The established model now contains 437 components and 0.369385 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097952Jy sigma=0.431155 ! Fit after self-cal, rms=0.098024Jy sigma=0.431047 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.000275038 Jy ! Component: 100 - total flux cleaned = 0.000405509 Jy ! Component: 150 - total flux cleaned = 0.000557818 Jy ! Component: 200 - total flux cleaned = 0.000666634 Jy ! Component: 250 - total flux cleaned = 0.000739167 Jy ! Component: 300 - total flux cleaned = 0.000780708 Jy ! Component: 350 - total flux cleaned = 0.000833267 Jy ! Component: 400 - total flux cleaned = 0.000871169 Jy ! Total flux subtracted in 400 components = 0.000871169 Jy ! Clean residual min=-0.002546 max=0.002059 Jy/beam ! Clean residual mean=-0.000029 rms=0.000645 Jy/beam ! Combined flux in latest and established models = 0.370256 Jy ! Performing phase self-cal ! Adding 49 model components to the UV plane model. ! The established model now contains 458 components and 0.370256 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.098039Jy sigma=0.431099 ! Fit after self-cal, rms=0.098037Jy sigma=0.431094 ! 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=2.403 mas, bmaj=3.899 mas, bpa=10.87 degrees ! Estimated noise=2.12979 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.403 x 3.899 at 10.87 degrees (North through East) ! Clean map min=-0.0046584 max=0.32174 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 print "=== CLearing model and starting over ===" ! === CLearing model and starting over === clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=2.119 mas, bmaj=3.584 mas, bpa=9.045 degrees ! Estimated noise=2.71846 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.250104 Jy ! Component: 100 - total flux cleaned = 0.305257 Jy ! Component: 150 - total flux cleaned = 0.323575 Jy ! Component: 200 - total flux cleaned = 0.332828 Jy ! Component: 250 - total flux cleaned = 0.337378 Jy ! Component: 300 - total flux cleaned = 0.339848 Jy ! Component: 350 - total flux cleaned = 0.342035 Jy ! Component: 400 - total flux cleaned = 0.343989 Jy ! Total flux subtracted in 400 components = 0.343989 Jy ! Clean residual min=-0.003788 max=0.004267 Jy/beam ! Clean residual mean=0.000020 rms=0.001011 Jy/beam ! Combined flux in latest and established models = 0.343989 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 38 model components to the UV plane model. ! The established model now contains 38 components and 0.343989 Jy ! Inverting map and beam ! Estimated beam: bmin=2.403 mas, bmaj=3.899 mas, bpa=10.87 degrees ! Estimated noise=2.12979 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.00657629 Jy ! Component: 100 - total flux cleaned = 0.0102464 Jy ! Component: 150 - total flux cleaned = 0.0128965 Jy ! Component: 200 - total flux cleaned = 0.0145701 Jy ! Component: 250 - total flux cleaned = 0.0157231 Jy ! Component: 300 - total flux cleaned = 0.0166698 Jy ! Component: 350 - total flux cleaned = 0.0176587 Jy ! Component: 400 - total flux cleaned = 0.0183796 Jy ! Total flux subtracted in 400 components = 0.0183796 Jy ! Clean residual min=-0.002889 max=0.003359 Jy/beam ! Clean residual mean=0.000076 rms=0.000768 Jy/beam ! Combined flux in latest and established models = 0.362369 Jy ! Adding 55 model components to the UV plane model. ! The established model now contains 82 components and 0.362369 Jy ! Inverting map ! Added new window around map position (97, -41). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00154068 Jy ! Component: 100 - total flux cleaned = 0.00182044 Jy ! Component: 150 - total flux cleaned = 0.0020871 Jy ! Component: 200 - total flux cleaned = 0.0023413 Jy ! Component: 250 - total flux cleaned = 0.00258437 Jy ! Component: 300 - total flux cleaned = 0.0027433 Jy ! Component: 350 - total flux cleaned = 0.00297439 Jy ! Component: 400 - total flux cleaned = 0.00326806 Jy ! Total flux subtracted in 400 components = 0.00326806 Jy ! Clean residual min=-0.002844 max=0.002906 Jy/beam ! Clean residual mean=0.000045 rms=0.000713 Jy/beam ! Combined flux in latest and established models = 0.365637 Jy ! Adding 58 model components to the UV plane model. ! The established model now contains 111 components and 0.365637 Jy ! Inverting map ! Added new window around map position (19.5, 109.5). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00105452 Jy ! Component: 100 - total flux cleaned = 0.000792903 Jy ! Component: 150 - total flux cleaned = 0.000939508 Jy ! Component: 200 - total flux cleaned = 0.00114903 Jy ! Component: 250 - total flux cleaned = 0.00128422 Jy ! Component: 300 - total flux cleaned = 0.00147937 Jy ! Component: 350 - total flux cleaned = 0.00160685 Jy ! Component: 400 - total flux cleaned = 0.00179256 Jy ! Total flux subtracted in 400 components = 0.00179256 Jy ! Clean residual min=-0.002825 max=0.002632 Jy/beam ! Clean residual mean=0.000024 rms=0.000678 Jy/beam ! Combined flux in latest and established models = 0.36743 Jy ! Adding 75 model components to the UV plane model. ! The established model now contains 151 components and 0.36743 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=8.01 mas, bmaj=16.67 mas, bpa=-7.037 degrees ! Estimated noise=2.50728 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0015696 Jy ! Component: 100 - total flux cleaned = 0.00186737 Jy ! Component: 150 - total flux cleaned = 0.00212071 Jy ! Component: 200 - total flux cleaned = 0.0023832 Jy ! Component: 250 - total flux cleaned = 0.00248445 Jy ! Component: 300 - total flux cleaned = 0.00261042 Jy ! Component: 350 - total flux cleaned = 0.00272805 Jy ! Component: 400 - total flux cleaned = 0.00283863 Jy ! Total flux subtracted in 400 components = 0.00283863 Jy ! Clean residual min=-0.002757 max=0.001946 Jy/beam ! Clean residual mean=-0.000026 rms=0.000659 Jy/beam ! Combined flux in latest and established models = 0.370268 Jy flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam keep new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.000103715 Jy ! Component: 100 - total flux cleaned = 0.000201696 Jy ! Component: 150 - total flux cleaned = 0.000271428 Jy ! Component: 200 - total flux cleaned = 0.00033849 Jy ! Component: 250 - total flux cleaned = 0.000381729 Jy ! Component: 300 - total flux cleaned = 0.000444106 Jy ! Component: 350 - total flux cleaned = 0.000483866 Jy ! Component: 400 - total flux cleaned = 0.00050309 Jy ! Total flux subtracted in 400 components = 0.00050309 Jy ! Clean residual min=-0.002717 max=0.001913 Jy/beam ! Clean residual mean=-0.000041 rms=0.000649 Jy/beam ! Combined flux in latest and established models = 0.370771 Jy ! Adding 87 model components to the UV plane model. ! The established model now contains 237 components and 0.370771 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=2.403 mas, bmaj=3.899 mas, bpa=10.87 degrees ! Estimated noise=2.12979 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000223122 Jy ! Component: 100 - total flux cleaned = -0.000355215 Jy ! Component: 150 - total flux cleaned = -0.000476554 Jy ! Component: 200 - total flux cleaned = -0.000534342 Jy ! Component: 250 - total flux cleaned = -0.000591563 Jy ! Component: 300 - total flux cleaned = -0.000647175 Jy ! Component: 350 - total flux cleaned = -0.000751474 Jy ! Component: 400 - total flux cleaned = -0.000750235 Jy ! Total flux subtracted in 400 components = -0.000750235 Jy ! Clean residual min=-0.002771 max=0.002578 Jy/beam ! Clean residual mean=-0.000008 rms=0.000650 Jy/beam ! Combined flux in latest and established models = 0.370021 Jy selfcal ! Performing phase self-cal ! Adding 77 model components to the UV plane model. ! The established model now contains 265 components and 0.370021 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.097877Jy sigma=0.430701 ! Fit after self-cal, rms=0.097792Jy sigma=0.430553 wmodel /scr/pima/rdv94_uvs/J0856-1105_S_map.mod ! Writing 265 model components to file: /scr/pima/rdv94_uvs/J0856-1105_S_map.mod wobs /scr/pima/rdv94_uvs/J0856-1105_S_uvs.fits ! Applying 832 buffered edits. ! Writing UV FITS file: /scr/pima/rdv94_uvs/J0856-1105_S_uvs.fits wwins /scr/pima/rdv94_uvs/J0856-1105_S_map.win ! wwins: Wrote 8 windows to /scr/pima/rdv94_uvs/J0856-1105_S_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! clean: niter=4096 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.000404939 Jy ! Component: 100 - total flux cleaned = -0.000592036 Jy ! Component: 150 - total flux cleaned = -0.00063619 Jy ! Component: 200 - total flux cleaned = -0.00059301 Jy ! Component: 250 - total flux cleaned = -0.000550331 Jy ! Component: 300 - total flux cleaned = -0.000508286 Jy ! Component: 350 - total flux cleaned = -0.00042593 Jy ! Component: 400 - total flux cleaned = -0.000305175 Jy ! Component: 450 - total flux cleaned = -0.000186845 Jy ! Component: 500 - total flux cleaned = -0.00014798 Jy ! Component: 550 - total flux cleaned = -0.000109357 Jy ! Component: 600 - total flux cleaned = -0.000146862 Jy ! Component: 650 - total flux cleaned = -0.000108983 Jy ! Component: 700 - total flux cleaned = -3.52736e-05 Jy ! Component: 750 - total flux cleaned = -0.000108827 Jy ! Component: 800 - total flux cleaned = -7.24527e-05 Jy ! Component: 850 - total flux cleaned = -0.00010839 Jy ! Component: 900 - total flux cleaned = -3.77685e-05 Jy ! Component: 950 - total flux cleaned = 3.20863e-05 Jy ! Component: 1000 - total flux cleaned = -0.000107041 Jy ! Component: 1050 - total flux cleaned = -3.8267e-05 Jy ! Component: 1100 - total flux cleaned = -0.000106442 Jy ! Component: 1150 - total flux cleaned = -7.25541e-05 Jy ! Component: 1200 - total flux cleaned = -3.91242e-05 Jy ! Component: 1250 - total flux cleaned = -7.21199e-05 Jy ! Component: 1300 - total flux cleaned = -6.02163e-06 Jy ! Component: 1350 - total flux cleaned = -3.82748e-05 Jy ! Component: 1400 - total flux cleaned = 2.65894e-05 Jy ! Component: 1450 - total flux cleaned = -3.78187e-05 Jy ! Component: 1500 - total flux cleaned = -5.72778e-06 Jy ! Component: 1550 - total flux cleaned = -5.65918e-06 Jy ! Component: 1600 - total flux cleaned = -3.7267e-05 Jy ! Component: 1650 - total flux cleaned = -9.97883e-05 Jy ! Component: 1700 - total flux cleaned = -0.000130721 Jy ! Component: 1750 - total flux cleaned = -0.000130689 Jy ! Component: 1800 - total flux cleaned = -0.000191737 Jy ! Component: 1850 - total flux cleaned = -0.000252594 Jy ! Component: 1900 - total flux cleaned = -0.000222281 Jy ! Component: 1950 - total flux cleaned = -0.000342394 Jy ! Component: 2000 - total flux cleaned = -0.000282744 Jy ! Component: 2050 - total flux cleaned = -0.000401516 Jy ! Component: 2100 - total flux cleaned = -0.000372072 Jy ! Component: 2150 - total flux cleaned = -0.000489138 Jy ! Component: 2200 - total flux cleaned = -0.000459917 Jy ! Component: 2250 - total flux cleaned = -0.000546558 Jy ! Component: 2300 - total flux cleaned = -0.000546509 Jy ! Component: 2350 - total flux cleaned = -0.000574993 Jy ! Component: 2400 - total flux cleaned = -0.000574972 Jy ! Component: 2450 - total flux cleaned = -0.000659613 Jy ! Component: 2500 - total flux cleaned = -0.00065952 Jy ! Component: 2550 - total flux cleaned = -0.00068739 Jy ! Component: 2600 - total flux cleaned = -0.000715164 Jy ! Component: 2650 - total flux cleaned = -0.000715346 Jy ! Component: 2700 - total flux cleaned = -0.000797837 Jy ! Component: 2750 - total flux cleaned = -0.000825127 Jy ! Component: 2800 - total flux cleaned = -0.000852141 Jy ! Component: 2850 - total flux cleaned = -0.00087927 Jy ! Component: 2900 - total flux cleaned = -0.000852235 Jy ! Component: 2950 - total flux cleaned = -0.000959467 Jy ! Component: 3000 - total flux cleaned = -0.00090594 Jy ! Component: 3050 - total flux cleaned = -0.000932611 Jy ! Component: 3100 - total flux cleaned = -0.000985489 Jy ! Component: 3150 - total flux cleaned = -0.000985551 Jy ! Component: 3200 - total flux cleaned = -0.00101174 Jy ! Component: 3250 - total flux cleaned = -0.000985795 Jy ! Component: 3300 - total flux cleaned = -0.00108934 Jy ! Component: 3350 - total flux cleaned = -0.00108932 Jy ! Component: 3400 - total flux cleaned = -0.00106363 Jy ! Component: 3450 - total flux cleaned = -0.00108919 Jy ! Component: 3500 - total flux cleaned = -0.00106374 Jy ! Component: 3550 - total flux cleaned = -0.00113959 Jy ! Component: 3600 - total flux cleaned = -0.0011647 Jy ! Component: 3650 - total flux cleaned = -0.00116459 Jy ! Component: 3700 - total flux cleaned = -0.00113959 Jy ! Component: 3750 - total flux cleaned = -0.0011892 Jy ! Component: 3800 - total flux cleaned = -0.00118913 Jy ! Component: 3850 - total flux cleaned = -0.00118905 Jy ! Component: 3900 - total flux cleaned = -0.00121365 Jy ! Component: 3950 - total flux cleaned = -0.00123807 Jy ! Component: 4000 - total flux cleaned = -0.00123801 Jy ! Component: 4050 - total flux cleaned = -0.00121388 Jy ! Total flux subtracted in 4096 components = -0.00126217 Jy ! Clean residual min=-0.001668 max=0.001538 Jy/beam ! Clean residual mean=-0.000016 rms=0.000434 Jy/beam ! Combined flux in latest and established models = 0.368759 Jy keep ! Adding 307 model components to the UV plane model. ! The established model now contains 572 components and 0.368759 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv94_uvs/J0856-1105_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 2.403 x 3.899 at 10.87 degrees (North through East) ! Clean map min=-0.0036559 max=0.31657 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv94_uvs/J0856-1105_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.3165 rms= 0.00043 Jy quit ! Quitting program ! Log file /scr/pima/rdv94_uvs/J0856-1105_S_uvs.log closed on Tue Oct 8 12:16:26 2013