! Started logfile: /scr/pima/rdv62_uvs/J2225-0457_S_uvs.log on Fri Nov 22 22:42:03 2013 obs /scr/pima/rdv62_uvs/J2225-0457_S_uva.fits ! Reading UV FITS file: /scr/pima/rdv62_uvs/J2225-0457_S_uva.fits ! AN table 1: 7 integrations on 45 of 45 possible baselines. ! AN table 2: 3 integrations on 3 of 3 possible baselines. ! AN table 3: 13 integrations on 15 of 15 possible baselines. ! AN table 4: 9 integrations on 10 of 10 possible baselines. ! AN table 5: 3 integrations on 3 of 3 possible baselines. ! AN table 6: 21 integrations on 45 of 45 possible baselines. ! Apparent sampling: 0.372361 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J2225-0457 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 2.23299e+09 8e+06 1 8e+06 ! 02 2 2.26299e+09 8e+06 1 8e+06 ! 03 3 2.35299e+09 8e+06 1 8e+06 ! 04 4 2.37299e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 0 lines of history. ! ! Reading 2328 visibilities. select RR ! Selecting polarization: RR, channels: 1..4 ! Reading IF 1 channels: 1..1 ! Reading IF 2 channels: 2..2 ! Reading IF 3 channels: 3..3 ! Reading IF 4 channels: 4..4 print "field_size =",field_size ! field_size = 1024 print "field_cell =",field_cell ! field_cell = 0.4 print "taper_size =",taper_size ! taper_size = 20 ![@/home/lpetrov/bin/rdv_muppet /scr/pima/rdv62_uvs/J2225-0457_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:SC self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:HH self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 system "rm /tmp/ff.ff" mapsize field_size, field_cell ! Map grid = 1024x1024 pixels with 0.400x0.400 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 51 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 51 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 51 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 10 telescope corrections were flagged in sub-array 1. ! A total of 51 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=2.402510Jy sigma=54.686071 ! Fit after self-cal, rms=1.279251Jy sigma=31.983063 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. print "=== Flag out IVS stations" ! === Flag out IVS stations integer lun string cfg_file lun = outfile ( "/tmp/ff.ff" ) do i=1,nsub do j=1,ntel(i) fprintf lun, "%s\n", telname(j,i) end do end do close ( lun ) system "/home/lpetrov/bin/dimap_flag_process.csh /tmp/ff.ff" ![@/tmp/ff.ff \] flag 2:HH ! Flagging baselines of antenna 2:HH in the currently selected channels. flag 2:WZ ! Flagging baselines of antenna 2:WZ in the currently selected channels. flag 2:ZC ! Flagging baselines of antenna 2:ZC in the currently selected channels. flag 3:HH ! Flagging baselines of antenna 3:HH in the currently selected channels. flag 3:NY ! Flagging baselines of antenna 3:NY in the currently selected channels. flag 3:WZ ! Flagging baselines of antenna 3:WZ in the currently selected channels. flag 3:ZC ! Flagging baselines of antenna 3:ZC in the currently selected channels. flag 4:HH ! Flagging baselines of antenna 4:HH in the currently selected channels. flag 4:WZ ! Flagging baselines of antenna 4:WZ in the currently selected channels. flag 4:ZC ! Flagging baselines of antenna 4:ZC in the currently selected channels. flag 6:KK ! Flagging baselines of antenna 6:KK in the currently selected channels. ![Exited script file: /tmp/ff.ff] uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.106 mas, bmaj=7.151 mas, bpa=-5.665 degrees ! Estimated noise=1.0624 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 = 1.73212 Jy ! Component: 100 - total flux cleaned = 2.23804 Jy ! Total flux subtracted in 100 components = 2.23804 Jy ! Clean residual min=-0.071638 max=0.218695 Jy/beam ! Clean residual mean=-0.000350 rms=0.026085 Jy/beam ! Combined flux in latest and established models = 2.23804 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 10 components and 2.23804 Jy ! ! Correcting IF 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 27 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 27 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 27 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 27 telescope corrections were flagged in sub-array 3. ! A total of 27 telescope corrections were flagged in sub-array 4. ! A total of 3 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.531019Jy sigma=14.268057 ! Fit after self-cal, rms=0.534310Jy sigma=14.425583 ! Inverting map and beam ! Estimated beam: bmin=3.092 mas, bmaj=7.216 mas, bpa=-5.463 degrees ! Estimated noise=1.07051 mJy/beam. ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.241238 Jy ! Component: 100 - total flux cleaned = 0.360801 Jy ! Total flux subtracted in 100 components = 0.360801 Jy ! Clean residual min=-0.059347 max=0.100015 Jy/beam ! Clean residual mean=-0.000074 rms=0.017517 Jy/beam ! Combined flux in latest and established models = 2.59884 Jy ! Performing phase self-cal ! Adding 22 model components to the UV plane model. ! The established model now contains 29 components and 2.59884 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.419192Jy sigma=11.063271 ! Fit after self-cal, rms=0.418177Jy sigma=11.038723 ! 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=3.538 mas, bmaj=8.112 mas, bpa=-5.678 degrees ! Estimated noise=0.799885 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.107274 Jy ! Component: 100 - total flux cleaned = 0.161624 Jy ! Component: 150 - total flux cleaned = 0.192487 Jy ! Component: 200 - total flux cleaned = 0.211234 Jy ! Total flux subtracted in 200 components = 0.211234 Jy ! Clean residual min=-0.072015 max=0.143834 Jy/beam ! Clean residual mean=0.000294 rms=0.024014 Jy/beam ! Combined flux in latest and established models = 2.81007 Jy ! Performing phase self-cal ! Adding 45 model components to the UV plane model. ! The established model now contains 68 components and 2.81008 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.399103Jy sigma=10.356941 ! Fit after self-cal, rms=0.398750Jy sigma=10.352821 ! Inverting map ! Added new window around map position (49.2, -18). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.134345 Jy ! Component: 100 - total flux cleaned = 0.204976 Jy ! Component: 150 - total flux cleaned = 0.249671 Jy ! Component: 200 - total flux cleaned = 0.27973 Jy ! Total flux subtracted in 200 components = 0.27973 Jy ! Clean residual min=-0.054325 max=0.074342 Jy/beam ! Clean residual mean=0.000090 rms=0.016181 Jy/beam ! Combined flux in latest and established models = 3.08981 Jy ! Performing phase self-cal ! Adding 75 model components to the UV plane model. ! The established model now contains 143 components and 3.0898 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.365579Jy sigma=9.145630 ! Fit after self-cal, rms=0.361605Jy sigma=9.076739 ! 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. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 0.99 HN 0.98 KP 0.99 ! LA 1.12 MK 1.01 NL 1.03 OV 1.01 ! PT 1.04 SC 0.93 ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* SC 0.00* WZ 0.00* ! ZC 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.96 FD 1.00 HN 0.94 KK 1.00* ! LA 1.16 MK 1.01 NL 1.10 OV 0.96 ! PT 1.05 SC 0.92 ! ! ! Correcting IF 2. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.02 FD 1.06 HN 1.11 KP 1.01 ! LA 1.03 MK 1.02 NL 1.06 OV 1.04 ! PT 1.04 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* SC 0.00* WZ 0.00* ! ZC 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 0.99 FD 1.06 HN 1.07 KK 1.00* ! LA 1.04 MK 1.06 NL 1.11 OV 0.99 ! PT 1.03 SC 0.98 ! ! ! Correcting IF 3. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 1.07 FD 1.06 HN 1.04 KP 0.99 ! LA 0.83 MK 1.03 NL 1.13 OV 1.07 ! PT 1.05 SC 1.00 ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* SC 0.00* WZ 0.00* ! ZC 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.00 FD 1.06 HN 0.98 KK 1.00* ! LA 0.86 MK 1.04 NL 1.15 OV 1.03 ! PT 1.04 SC 0.99 ! ! ! Correcting IF 4. ! A total of 3 telescope corrections were flagged in sub-array 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.98 FD 1.02 HN 0.99 KP 0.93 ! LA 0.71 MK 0.98 NL 1.03 OV 0.87 ! PT 1.00 SC 1.10 ! ! Telescope amplitude corrections in sub-array 2: ! HH 0.00* WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 3: ! HH 0.00* HN 0.00* NY 0.00* SC 0.00* ! WZ 0.00* ZC 0.00* ! ! Telescope amplitude corrections in sub-array 4: ! HH 0.00* HN 0.00* SC 0.00* WZ 0.00* ! ZC 0.00* ! ! A total of 9 telescope corrections were flagged in sub-array 5. ! Telescope amplitude corrections in sub-array 5: ! KP 0.00* LA 0.00* PT 0.00* ! ! Telescope amplitude corrections in sub-array 6: ! BR 1.01 FD 1.03 HN 0.95 KK 1.00* ! LA 0.75 MK 1.00 NL 1.01 OV 0.86 ! PT 0.98 SC 1.07 ! ! ! Fit before self-cal, rms=0.361605Jy sigma=9.076739 ! Fit after self-cal, rms=0.180736Jy sigma=4.385692 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=3.549 mas, bmaj=8.098 mas, bpa=-5.492 degrees ! Estimated noise=0.805175 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.0323608 Jy ! Component: 100 - total flux cleaned = 0.0471798 Jy ! Component: 150 - total flux cleaned = 0.0557492 Jy ! Component: 200 - total flux cleaned = 0.0596182 Jy ! Total flux subtracted in 200 components = 0.0596182 Jy ! Clean residual min=-0.027208 max=0.040286 Jy/beam ! Clean residual mean=0.000080 rms=0.008928 Jy/beam ! Combined flux in latest and established models = 3.14942 Jy ! Performing phase self-cal ! Adding 32 model components to the UV plane model. ! The established model now contains 173 components and 3.14942 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.175478Jy sigma=4.197488 ! Fit after self-cal, rms=0.175176Jy sigma=4.186481 ! 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.175176Jy sigma=4.186481 ! Fit after self-cal, rms=0.175067Jy sigma=4.180304 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.55 mas, bmaj=8.101 mas, bpa=-5.471 degrees ! Estimated noise=0.80458 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000350027 Jy ! Component: 100 - total flux cleaned = 0.00413452 Jy ! Component: 150 - total flux cleaned = 0.00627749 Jy ! Component: 200 - total flux cleaned = 0.00825107 Jy ! Total flux subtracted in 200 components = 0.00825107 Jy ! Clean residual min=-0.027504 max=0.038021 Jy/beam ! Clean residual mean=0.000058 rms=0.008559 Jy/beam ! Combined flux in latest and established models = 3.15767 Jy selfcal ! Performing phase self-cal ! Adding 41 model components to the UV plane model. ! The established model now contains 193 components and 3.15767 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.173537Jy sigma=4.141208 ! Fit after self-cal, rms=0.173402Jy sigma=4.138466 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 (58.4, -19.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0343716 Jy ! Component: 100 - total flux cleaned = 0.0439356 Jy ! Component: 150 - total flux cleaned = 0.0458769 Jy ! Component: 200 - total flux cleaned = 0.0490501 Jy ! Total flux subtracted in 200 components = 0.0490501 Jy ! Clean residual min=-0.023689 max=0.036143 Jy/beam ! Clean residual mean=0.000011 rms=0.007542 Jy/beam ! Combined flux in latest and established models = 3.20672 Jy ! Performing phase self-cal ! Adding 48 model components to the UV plane model. ! The established model now contains 240 components and 3.20673 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.167858Jy sigma=3.946385 ! Fit after self-cal, rms=0.167315Jy sigma=3.934663 ! Inverting map ! Added new window around map position (5.6, -2.4). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0328022 Jy ! Component: 100 - total flux cleaned = 0.0343149 Jy ! Component: 150 - total flux cleaned = 0.036519 Jy ! Component: 200 - total flux cleaned = 0.0394861 Jy ! Total flux subtracted in 200 components = 0.0394861 Jy ! Clean residual min=-0.021801 max=0.024402 Jy/beam ! Clean residual mean=0.000032 rms=0.006866 Jy/beam ! Combined flux in latest and established models = 3.24621 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 263 components and 3.24621 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.160603Jy sigma=3.763167 ! Fit after self-cal, rms=0.158948Jy sigma=3.738633 ! 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.158948Jy sigma=3.738633 ! Fit after self-cal, rms=0.150580Jy sigma=3.524807 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.549 mas, bmaj=8.105 mas, bpa=-5.582 degrees ! Estimated noise=0.80363 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0181558 Jy ! Component: 100 - total flux cleaned = 0.0302844 Jy ! Component: 150 - total flux cleaned = 0.0338435 Jy ! Component: 200 - total flux cleaned = 0.0358313 Jy ! Total flux subtracted in 200 components = 0.0358313 Jy ! Clean residual min=-0.018701 max=0.019733 Jy/beam ! Clean residual mean=0.000040 rms=0.005490 Jy/beam ! Combined flux in latest and established models = 3.28204 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 289 components and 3.28204 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147604Jy sigma=3.445806 ! Fit after self-cal, rms=0.147049Jy sigma=3.439337 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00398641 Jy ! Component: 100 - total flux cleaned = 0.00624089 Jy ! Component: 150 - total flux cleaned = 0.00793941 Jy ! Component: 200 - total flux cleaned = 0.00887491 Jy ! Total flux subtracted in 200 components = 0.00887491 Jy ! Clean residual min=-0.017953 max=0.018465 Jy/beam ! Clean residual mean=0.000039 rms=0.005283 Jy/beam ! Combined flux in latest and established models = 3.29092 Jy ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 309 components and 3.29092 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.145576Jy sigma=3.402900 ! Fit after self-cal, rms=0.145265Jy sigma=3.400012 ! 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.145265Jy sigma=3.400012 ! Fit after self-cal, rms=0.144544Jy sigma=3.389263 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.786 mas, bmaj=12.84 mas, bpa=-4.448 degrees ! Estimated noise=0.857225 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.014974 Jy ! Component: 100 - total flux cleaned = 0.0225241 Jy ! Component: 150 - total flux cleaned = 0.0275228 Jy ! Component: 200 - total flux cleaned = 0.0301066 Jy ! Total flux subtracted in 200 components = 0.0301066 Jy ! Clean residual min=-0.022860 max=0.026271 Jy/beam ! Clean residual mean=0.000072 rms=0.008114 Jy/beam ! Combined flux in latest and established models = 3.32102 Jy selfcal ! Performing phase self-cal ! Adding 46 model components to the UV plane model. ! The established model now contains 341 components and 3.32103 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143913Jy sigma=3.349641 ! Fit after self-cal, rms=0.143556Jy sigma=3.343735 flux_cutoff = imstat(rms) * dynam ! Inverting map repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00268745 Jy ! Component: 100 - total flux cleaned = 0.00388914 Jy ! Component: 150 - total flux cleaned = 0.00495701 Jy ! Component: 200 - total flux cleaned = 0.00591518 Jy ! Total flux subtracted in 200 components = 0.00591518 Jy ! Clean residual min=-0.022304 max=0.026101 Jy/beam ! Clean residual mean=0.000069 rms=0.007963 Jy/beam ! Combined flux in latest and established models = 3.32694 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 358 components and 3.32694 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.143263Jy sigma=3.327239 ! Fit after self-cal, rms=0.143050Jy sigma=3.324609 ! 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 1:SC,true ! 1:SC self-cal status=fixed weight=1 selfant 3:HN,true ! 3:HN self-cal status=fixed weight=1 selfant 3:SC,true ! 3:SC self-cal status=fixed weight=1 selfant 4:HN,true ! 4:HN self-cal status=fixed weight=1 selfant 4:SC,true ! 4:SC self-cal status=fixed weight=1 selfant 5:KP,true ! 5:KP self-cal status=fixed weight=1 selfant 5:LA,true ! 5:LA self-cal status=fixed weight=1 selfant 5:PT,true ! 5:PT self-cal status=fixed weight=1 selfant 6:BR,true ! 6:BR self-cal status=fixed weight=1 selfant 6:FD,true ! 6:FD self-cal status=fixed weight=1 selfant 6:HN,true ! 6:HN self-cal status=fixed weight=1 selfant 6:LA,true ! 6:LA self-cal status=fixed weight=1 selfant 6:MK,true ! 6:MK self-cal status=fixed weight=1 selfant 6:NL,true ! 6:NL self-cal status=fixed weight=1 selfant 6:OV,true ! 6:OV self-cal status=fixed weight=1 selfant 6:PT,true ! 6:PT self-cal status=fixed weight=1 selfant 6:SC,true ! 6:SC self-cal status=fixed weight=1 ![Exited script file: /tmp/ff.ff] selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 40 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 40 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 40 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 40 telescope corrections were flagged in sub-array 3. ! A total of 18 telescope corrections were flagged in sub-array 4. ! A total of 18 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.717718Jy sigma=11.820127 ! Fit after self-cal, rms=0.668438Jy sigma=11.808824 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:SC self-cal status=correctable weight=1 ! 2:HH self-cal status=correctable weight=1 ! 2:WZ self-cal status=correctable weight=1 ! 2:ZC self-cal status=correctable weight=1 ! 3:HH self-cal status=correctable weight=1 ! 3:HN self-cal status=correctable weight=1 ! 3:NY self-cal status=correctable weight=1 ! 3:SC self-cal status=correctable weight=1 ! 3:WZ self-cal status=correctable weight=1 ! 3:ZC self-cal status=correctable weight=1 ! 4:HH self-cal status=correctable weight=1 ! 4:HN self-cal status=correctable weight=1 ! 4:SC self-cal status=correctable weight=1 ! 4:WZ self-cal status=correctable weight=1 ! 4:ZC self-cal status=correctable weight=1 ! 5:KP self-cal status=correctable weight=1 ! 5:LA self-cal status=correctable weight=1 ! 5:PT self-cal status=correctable weight=1 ! 6:BR self-cal status=correctable weight=1 ! 6:FD self-cal status=correctable weight=1 ! 6:HN self-cal status=correctable weight=1 ! 6:KK self-cal status=correctable weight=1 ! 6:LA self-cal status=correctable weight=1 ! 6:MK self-cal status=correctable weight=1 ! 6:NL self-cal status=correctable weight=1 ! 6:OV self-cal status=correctable weight=1 ! 6:PT self-cal status=correctable weight=1 ! 6:SC self-cal status=correctable weight=1 selfcal true, true, soltime1 ! Performing amp+phase self-cal over 800000 minute time intervals ! ! Correcting IF 1. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 135 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 135 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 135 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 20 telescope corrections were flagged in sub-array 3. ! A total of 12 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 135 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.668438Jy sigma=11.808824 ! Fit after self-cal, rms=0.181491Jy sigma=3.147428 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.769 mas, bmaj=5.262 mas, bpa=-4.502 degrees ! Estimated noise=1.14381 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.0154121 Jy ! Component: 100 - total flux cleaned = -0.020357 Jy ! Component: 150 - total flux cleaned = -0.0236759 Jy ! Component: 200 - total flux cleaned = -0.0282374 Jy ! Total flux subtracted in 200 components = -0.0282374 Jy ! Clean residual min=-0.016097 max=0.018684 Jy/beam ! Clean residual mean=0.000410 rms=0.004736 Jy/beam ! Combined flux in latest and established models = 3.2987 Jy ! Performing phase self-cal ! Adding 23 model components to the UV plane model. ! The established model now contains 375 components and 3.2987 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.177692Jy sigma=3.129850 ! Fit after self-cal, rms=0.177972Jy sigma=3.126556 ! 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=3.222 mas, bmaj=6.449 mas, bpa=-6.096 degrees ! Estimated noise=0.814613 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.000281766 Jy ! Component: 100 - total flux cleaned = 0.00279404 Jy ! Component: 150 - total flux cleaned = 0.00432344 Jy ! Component: 200 - total flux cleaned = 0.00543975 Jy ! Component: 250 - total flux cleaned = 0.0069974 Jy ! Component: 300 - total flux cleaned = 0.00772508 Jy ! Component: 350 - total flux cleaned = 0.00864695 Jy ! Component: 400 - total flux cleaned = 0.00908469 Jy ! Total flux subtracted in 400 components = 0.00908469 Jy ! Clean residual min=-0.015635 max=0.013950 Jy/beam ! Clean residual mean=0.000356 rms=0.004523 Jy/beam ! Combined flux in latest and established models = 3.30779 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 407 components and 3.30779 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.174655Jy sigma=3.097740 ! Fit after self-cal, rms=0.174385Jy sigma=3.095111 ! 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.000285968 Jy ! Component: 100 - total flux cleaned = 0.00184772 Jy ! Component: 150 - total flux cleaned = 0.00225524 Jy ! Component: 200 - total flux cleaned = 0.00205738 Jy ! Component: 250 - total flux cleaned = 0.00225004 Jy ! Component: 300 - total flux cleaned = 0.00225111 Jy ! Component: 350 - total flux cleaned = 0.00225195 Jy ! Component: 400 - total flux cleaned = 0.0022555 Jy ! Total flux subtracted in 400 components = 0.0022555 Jy ! Clean residual min=-0.015318 max=0.014032 Jy/beam ! Clean residual mean=0.000356 rms=0.004467 Jy/beam ! Combined flux in latest and established models = 3.31004 Jy ! Performing phase self-cal ! Adding 56 model components to the UV plane model. ! The established model now contains 426 components and 3.31004 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.172466Jy sigma=3.084103 ! Fit after self-cal, rms=0.172360Jy sigma=3.083070 ! 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 13 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 60 telescope corrections were flagged in sub-array 3. ! A total of 30 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 2. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 60 telescope corrections were flagged in sub-array 3. ! A total of 30 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 3. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 60 telescope corrections were flagged in sub-array 3. ! A total of 30 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Correcting IF 4. ! A total of 13 telescope corrections were flagged in sub-array 1. ! A total of 9 telescope corrections were flagged in sub-array 2. ! A total of 60 telescope corrections were flagged in sub-array 3. ! A total of 30 telescope corrections were flagged in sub-array 4. ! A total of 9 telescope corrections were flagged in sub-array 5. ! A total of 153 telescope corrections were flagged in sub-array 6. ! ! Fit before self-cal, rms=0.676688Jy sigma=11.466224 ! Fit after self-cal, rms=0.171689Jy sigma=3.073904 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.228 mas, bmaj=6.462 mas, bpa=-6.05 degrees ! Estimated noise=0.815428 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00618919 Jy ! Component: 100 - total flux cleaned = -0.00619503 Jy ! Component: 150 - total flux cleaned = -0.00681652 Jy ! Component: 200 - total flux cleaned = -0.00681403 Jy ! Component: 250 - total flux cleaned = -0.00720174 Jy ! Component: 300 - total flux cleaned = -0.00756344 Jy ! Component: 350 - total flux cleaned = -0.00791122 Jy ! Component: 400 - total flux cleaned = -0.00791391 Jy ! Total flux subtracted in 400 components = -0.00791391 Jy ! Clean residual min=-0.015233 max=0.014858 Jy/beam ! Clean residual mean=0.000323 rms=0.004377 Jy/beam ! Combined flux in latest and established models = 3.30213 Jy selfcal ! Performing phase self-cal ! Adding 52 model components to the UV plane model. ! The established model now contains 439 components and 3.30213 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.169685Jy sigma=3.060365 ! Fit after self-cal, rms=0.169640Jy sigma=3.059465 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.000861477 Jy ! Component: 100 - total flux cleaned = 0.00103201 Jy ! Component: 150 - total flux cleaned = 0.00119553 Jy ! Component: 200 - total flux cleaned = 0.00104033 Jy ! Component: 250 - total flux cleaned = 0.000889455 Jy ! Component: 300 - total flux cleaned = 0.00103799 Jy ! Component: 350 - total flux cleaned = 0.000890367 Jy ! Component: 400 - total flux cleaned = 0.000891892 Jy ! Total flux subtracted in 400 components = 0.000891892 Jy ! Clean residual min=-0.015226 max=0.014841 Jy/beam ! Clean residual mean=0.000323 rms=0.004355 Jy/beam ! Combined flux in latest and established models = 3.30302 Jy ! Performing phase self-cal ! Adding 60 model components to the UV plane model. ! The established model now contains 455 components and 3.30302 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.168337Jy sigma=3.053148 ! Fit after self-cal, rms=0.168294Jy sigma=3.052742 ! 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.168294Jy sigma=3.052742 ! Fit after self-cal, rms=0.167766Jy sigma=3.044277 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=3.234 mas, bmaj=6.47 mas, bpa=-6.004 degrees ! Estimated noise=0.81584 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00451575 Jy ! Component: 100 - total flux cleaned = -0.00549987 Jy ! Component: 150 - total flux cleaned = -0.0060168 Jy ! Component: 200 - total flux cleaned = -0.00617575 Jy ! Component: 250 - total flux cleaned = -0.00647741 Jy ! Component: 300 - total flux cleaned = -0.00647832 Jy ! Component: 350 - total flux cleaned = -0.00633445 Jy ! Component: 400 - total flux cleaned = -0.00647438 Jy ! Total flux subtracted in 400 components = -0.00647438 Jy ! Clean residual min=-0.014972 max=0.015063 Jy/beam ! Clean residual mean=0.000297 rms=0.004294 Jy/beam ! Combined flux in latest and established models = 3.29655 Jy selfcal ! Performing phase self-cal ! Adding 59 model components to the UV plane model. ! The established model now contains 468 components and 3.29655 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.166353Jy sigma=3.035204 ! Fit after self-cal, rms=0.166302Jy sigma=3.034614 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 (-66, -66.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0163307 Jy ! Component: 100 - total flux cleaned = 0.0278364 Jy ! Component: 150 - total flux cleaned = 0.0367891 Jy ! Component: 200 - total flux cleaned = 0.0418097 Jy ! Component: 250 - total flux cleaned = 0.045109 Jy ! Component: 300 - total flux cleaned = 0.0468964 Jy ! Component: 350 - total flux cleaned = 0.0489047 Jy ! Component: 400 - total flux cleaned = 0.051116 Jy ! Total flux subtracted in 400 components = 0.051116 Jy ! Clean residual min=-0.013489 max=0.013843 Jy/beam ! Clean residual mean=0.000287 rms=0.003926 Jy/beam ! Combined flux in latest and established models = 3.34766 Jy ! Performing phase self-cal ! Adding 63 model components to the UV plane model. ! The established model now contains 523 components and 3.34766 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.163354Jy sigma=2.962557 ! Fit after self-cal, rms=0.163284Jy sigma=2.954059 ! 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.163284Jy sigma=2.954059 ! Fit after self-cal, rms=0.151958Jy sigma=2.566056 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=7.672 mas, bmaj=12.6 mas, bpa=-5.167 degrees ! Estimated noise=0.850673 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0138241 Jy ! Component: 100 - total flux cleaned = 0.0183109 Jy ! Component: 150 - total flux cleaned = 0.0221997 Jy ! Component: 200 - total flux cleaned = 0.0257191 Jy ! Component: 250 - total flux cleaned = 0.0289889 Jy ! Component: 300 - total flux cleaned = 0.0316638 Jy ! Component: 350 - total flux cleaned = 0.0326935 Jy ! Component: 400 - total flux cleaned = 0.0334897 Jy ! Total flux subtracted in 400 components = 0.0334897 Jy ! Clean residual min=-0.017067 max=0.022571 Jy/beam ! Clean residual mean=0.000496 rms=0.006492 Jy/beam ! Combined flux in latest and established models = 3.38115 Jy selfcal ! Performing phase self-cal ! Adding 91 model components to the UV plane model. ! The established model now contains 604 components and 3.38115 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.151863Jy sigma=2.534154 ! Fit after self-cal, rms=0.151676Jy sigma=2.528611 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 (91.6, -54.4). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0193215 Jy ! Component: 100 - total flux cleaned = 0.0266582 Jy ! Component: 150 - total flux cleaned = 0.0302532 Jy ! Component: 200 - total flux cleaned = 0.033868 Jy ! Component: 250 - total flux cleaned = 0.0364604 Jy ! Component: 300 - total flux cleaned = 0.039549 Jy ! Component: 350 - total flux cleaned = 0.0411647 Jy ! Component: 400 - total flux cleaned = 0.0430093 Jy ! Total flux subtracted in 400 components = 0.0430093 Jy ! Clean residual min=-0.016769 max=0.020106 Jy/beam ! Clean residual mean=0.000466 rms=0.005625 Jy/beam ! Combined flux in latest and established models = 3.42416 Jy ! Performing phase self-cal ! Adding 154 model components to the UV plane model. ! The established model now contains 747 components and 3.42416 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.149783Jy sigma=2.441798 ! Fit after self-cal, rms=0.149636Jy sigma=2.435536 ! Inverting map ! Added new window around map position (-29.2, -88). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0195502 Jy ! Component: 100 - total flux cleaned = 0.0286632 Jy ! Component: 150 - total flux cleaned = 0.0345981 Jy ! Component: 200 - total flux cleaned = 0.0383087 Jy ! Component: 250 - total flux cleaned = 0.0408439 Jy ! Component: 300 - total flux cleaned = 0.0426671 Jy ! Component: 350 - total flux cleaned = 0.0447 Jy ! Component: 400 - total flux cleaned = 0.0455336 Jy ! Total flux subtracted in 400 components = 0.0455336 Jy ! Clean residual min=-0.012464 max=0.017737 Jy/beam ! Clean residual mean=0.000405 rms=0.004727 Jy/beam ! Combined flux in latest and established models = 3.4697 Jy ! Performing phase self-cal ! Adding 163 model components to the UV plane model. ! The established model now contains 893 components and 3.4697 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.148130Jy sigma=2.357591 ! Fit after self-cal, rms=0.148057Jy sigma=2.354002 ! Inverting map ! Added new window around map position (-80.8, -85.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0168034 Jy ! Component: 100 - total flux cleaned = 0.0259195 Jy ! Component: 150 - total flux cleaned = 0.034875 Jy ! Component: 200 - total flux cleaned = 0.0430804 Jy ! Component: 250 - total flux cleaned = 0.0507732 Jy ! Component: 300 - total flux cleaned = 0.0570597 Jy ! Component: 350 - total flux cleaned = 0.0627211 Jy ! Component: 400 - total flux cleaned = 0.0669365 Jy ! Total flux subtracted in 400 components = 0.0669365 Jy ! Clean residual min=-0.013634 max=0.014058 Jy/beam ! Clean residual mean=0.000222 rms=0.004083 Jy/beam ! Combined flux in latest and established models = 3.53663 Jy ! Performing phase self-cal ! Adding 163 model components to the UV plane model. ! The established model now contains 1043 components and 3.53663 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.147091Jy sigma=2.281736 ! Fit after self-cal, rms=0.147002Jy sigma=2.279490 ! Inverting map ! Added new window around map position (97.6, 15.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.01152 Jy ! Component: 100 - total flux cleaned = 0.0166259 Jy ! Component: 150 - total flux cleaned = 0.0217499 Jy ! Component: 200 - total flux cleaned = 0.0273823 Jy ! Component: 250 - total flux cleaned = 0.0317605 Jy ! Component: 300 - total flux cleaned = 0.0359341 Jy ! Component: 350 - total flux cleaned = 0.0402159 Jy ! Component: 400 - total flux cleaned = 0.0448671 Jy ! Total flux subtracted in 400 components = 0.0448671 Jy ! Clean residual min=-0.010881 max=0.010866 Jy/beam ! Clean residual mean=0.000272 rms=0.003399 Jy/beam ! Combined flux in latest and established models = 3.5815 Jy ! Performing phase self-cal ! Adding 139 model components to the UV plane model. ! The established model now contains 1170 components and 3.5815 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.146066Jy sigma=2.220750 ! Fit after self-cal, rms=0.145871Jy sigma=2.216985 ! 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=3.237 mas, bmaj=6.47 mas, bpa=-5.964 degrees ! Estimated noise=0.815767 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.237 x 6.47 at -5.964 degrees (North through East) ! Clean map min=-0.015863 max=2.2079 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.145871Jy sigma=2.216985 ! Fit after self-cal, rms=0.123126Jy sigma=1.835304 ! Inverting map and beam ! Estimated beam: bmin=3.236 mas, bmaj=6.475 mas, bpa=-5.97 degrees ! Estimated noise=0.814999 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00894659 Jy ! Component: 100 - total flux cleaned = 0.0128717 Jy ! Component: 150 - total flux cleaned = 0.015587 Jy ! Component: 200 - total flux cleaned = 0.018994 Jy ! Component: 250 - total flux cleaned = 0.0225044 Jy ! Component: 300 - total flux cleaned = 0.0256098 Jy ! Component: 350 - total flux cleaned = 0.0288561 Jy ! Component: 400 - total flux cleaned = 0.0315179 Jy ! Total flux subtracted in 400 components = 0.0315179 Jy ! Clean residual min=-0.007651 max=0.008368 Jy/beam ! Clean residual mean=0.000133 rms=0.002098 Jy/beam ! Combined flux in latest and established models = 3.61302 Jy ! Performing phase self-cal ! Adding 85 model components to the UV plane model. ! The established model now contains 1249 components and 3.61302 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.119514Jy sigma=1.779316 ! Fit after self-cal, rms=0.119445Jy sigma=1.777497 ! Inverting map ! Added new window around map position (-34.4, 51.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00700508 Jy ! Component: 100 - total flux cleaned = 0.00894892 Jy ! Component: 150 - total flux cleaned = 0.0110439 Jy ! Component: 200 - total flux cleaned = 0.0135216 Jy ! Component: 250 - total flux cleaned = 0.0157126 Jy ! Component: 300 - total flux cleaned = 0.0180755 Jy ! Component: 350 - total flux cleaned = 0.0206057 Jy ! Component: 400 - total flux cleaned = 0.0228727 Jy ! Total flux subtracted in 400 components = 0.0228727 Jy ! Clean residual min=-0.007179 max=0.008284 Jy/beam ! Clean residual mean=0.000125 rms=0.001931 Jy/beam ! Combined flux in latest and established models = 3.63589 Jy ! Performing phase self-cal ! Adding 99 model components to the UV plane model. ! The established model now contains 1300 components and 3.63589 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.116918Jy sigma=1.741781 ! Fit after self-cal, rms=0.116864Jy sigma=1.740519 ! Inverting map ! Added new window around map position (-58.8, 18). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00640417 Jy ! Component: 100 - total flux cleaned = 0.0090584 Jy ! Component: 150 - total flux cleaned = 0.012828 Jy ! Component: 200 - total flux cleaned = 0.0166732 Jy ! Component: 250 - total flux cleaned = 0.01904 Jy ! Component: 300 - total flux cleaned = 0.0219315 Jy ! Component: 350 - total flux cleaned = 0.0245858 Jy ! Component: 400 - total flux cleaned = 0.0275608 Jy ! Total flux subtracted in 400 components = 0.0275608 Jy ! Clean residual min=-0.006704 max=0.007172 Jy/beam ! Clean residual mean=0.000111 rms=0.001753 Jy/beam ! Combined flux in latest and established models = 3.66345 Jy ! Performing phase self-cal ! Adding 103 model components to the UV plane model. ! The established model now contains 1348 components and 3.66345 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.114720Jy sigma=1.709449 ! Fit after self-cal, rms=0.114667Jy sigma=1.708561 ! Inverting map ! Added new window around map position (14.8, 17.6). ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00219113 Jy ! Component: 100 - total flux cleaned = 0.00323529 Jy ! Component: 150 - total flux cleaned = 0.00546144 Jy ! Component: 200 - total flux cleaned = 0.00818356 Jy ! Component: 250 - total flux cleaned = 0.0108222 Jy ! Component: 300 - total flux cleaned = 0.0130278 Jy ! Component: 350 - total flux cleaned = 0.015551 Jy ! Component: 400 - total flux cleaned = 0.0180298 Jy ! Total flux subtracted in 400 components = 0.0180298 Jy ! Clean residual min=-0.006485 max=0.006277 Jy/beam ! Clean residual mean=0.000112 rms=0.001604 Jy/beam ! Combined flux in latest and established models = 3.68148 Jy ! Performing phase self-cal ! Adding 81 model components to the UV plane model. ! The established model now contains 1388 components and 3.68148 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.112489Jy sigma=1.678712 ! Fit after self-cal, rms=0.112461Jy sigma=1.678032 ! 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.112461Jy sigma=1.678032 ! Fit after self-cal, rms=0.112174Jy sigma=1.671280 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.782 mas, bmaj=5.295 mas, bpa=-4.414 degrees ! Estimated noise=1.14534 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.71019 Jy ! Component: 100 - total flux cleaned = 2.22658 Jy ! Component: 150 - total flux cleaned = 2.45939 Jy ! Component: 200 - total flux cleaned = 2.58865 Jy ! Component: 250 - total flux cleaned = 2.6795 Jy ! Component: 300 - total flux cleaned = 2.7506 Jy ! Component: 350 - total flux cleaned = 2.80808 Jy ! Component: 400 - total flux cleaned = 2.85592 Jy ! Total flux subtracted in 400 components = 2.85592 Jy ! Clean residual min=-0.029121 max=0.029377 Jy/beam ! Clean residual mean=0.000208 rms=0.007953 Jy/beam ! Combined flux in latest and established models = 2.85591 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 84 model components to the UV plane model. ! The established model now contains 84 components and 2.85591 Jy ! Inverting map and beam ! Estimated beam: bmin=3.236 mas, bmaj=6.479 mas, bpa=-5.98 degrees ! Estimated noise=0.815017 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.0725306 Jy ! Component: 100 - total flux cleaned = 0.127238 Jy ! Component: 150 - total flux cleaned = 0.172769 Jy ! Component: 200 - total flux cleaned = 0.205423 Jy ! Component: 250 - total flux cleaned = 0.233766 Jy ! Component: 300 - total flux cleaned = 0.260038 Jy ! Component: 350 - total flux cleaned = 0.281627 Jy ! Component: 400 - total flux cleaned = 0.297549 Jy ! Total flux subtracted in 400 components = 0.297549 Jy ! Clean residual min=-0.027514 max=0.023750 Jy/beam ! Clean residual mean=0.000294 rms=0.007066 Jy/beam ! Combined flux in latest and established models = 3.15346 Jy ! Adding 187 model components to the UV plane model. ! The established model now contains 255 components and 3.15346 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=7.674 mas, bmaj=12.61 mas, bpa=-5.155 degrees ! Estimated noise=0.849846 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0406562 Jy ! Component: 100 - total flux cleaned = 0.0713305 Jy ! Component: 150 - total flux cleaned = 0.0970339 Jy ! Component: 200 - total flux cleaned = 0.118944 Jy ! Component: 250 - total flux cleaned = 0.140801 Jy ! Component: 300 - total flux cleaned = 0.160827 Jy ! Component: 350 - total flux cleaned = 0.180213 Jy ! Component: 400 - total flux cleaned = 0.19586 Jy ! Total flux subtracted in 400 components = 0.19586 Jy ! Clean residual min=-0.022829 max=0.020608 Jy/beam ! Clean residual mean=0.000478 rms=0.007164 Jy/beam ! Combined flux in latest and established models = 3.34932 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.0161613 Jy ! Component: 100 - total flux cleaned = 0.0314317 Jy ! Component: 150 - total flux cleaned = 0.0459216 Jy ! Component: 200 - total flux cleaned = 0.0590644 Jy ! Component: 250 - total flux cleaned = 0.0721681 Jy ! Component: 300 - total flux cleaned = 0.0840912 Jy ! Component: 350 - total flux cleaned = 0.0960149 Jy ! Component: 400 - total flux cleaned = 0.107408 Jy ! Total flux subtracted in 400 components = 0.107408 Jy ! Clean residual min=-0.016971 max=0.015515 Jy/beam ! Clean residual mean=0.000350 rms=0.005233 Jy/beam ! Combined flux in latest and established models = 3.45673 Jy ! Adding 388 model components to the UV plane model. ! The established model now contains 630 components and 3.45673 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=3.236 mas, bmaj=6.479 mas, bpa=-5.98 degrees ! Estimated noise=0.815017 mJy/beam. ! clean: niter=400 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00169087 Jy ! Component: 100 - total flux cleaned = 0.00783497 Jy ! Component: 150 - total flux cleaned = 0.0152508 Jy ! Component: 200 - total flux cleaned = 0.0222375 Jy ! Component: 250 - total flux cleaned = 0.0292979 Jy ! Component: 300 - total flux cleaned = 0.0360497 Jy ! Component: 350 - total flux cleaned = 0.0425438 Jy ! Component: 400 - total flux cleaned = 0.0491504 Jy ! Total flux subtracted in 400 components = 0.0491504 Jy ! Clean residual min=-0.012491 max=0.011117 Jy/beam ! Clean residual mean=0.000166 rms=0.002835 Jy/beam ! Combined flux in latest and established models = 3.50588 Jy selfcal ! Performing phase self-cal ! Adding 88 model components to the UV plane model. ! The established model now contains 716 components and 3.50588 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.115555Jy sigma=1.920846 ! Fit after self-cal, rms=0.113679Jy sigma=1.894557 wmodel /scr/pima/rdv62_uvs/J2225-0457_S_map.mod ! Writing 716 model components to file: /scr/pima/rdv62_uvs/J2225-0457_S_map.mod wobs /scr/pima/rdv62_uvs/J2225-0457_S_uvs.fits ! Applying 675 buffered edits. ! Writing UV FITS file: /scr/pima/rdv62_uvs/J2225-0457_S_uvs.fits wwins /scr/pima/rdv62_uvs/J2225-0457_S_map.win ! wwins: Wrote 12 windows to /scr/pima/rdv62_uvs/J2225-0457_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.000250169 Jy ! Component: 100 - total flux cleaned = 0.000292622 Jy ! Component: 150 - total flux cleaned = -0.000203648 Jy ! Component: 200 - total flux cleaned = -0.000359672 Jy ! Component: 250 - total flux cleaned = -0.000962237 Jy ! Component: 300 - total flux cleaned = -0.00183693 Jy ! Component: 350 - total flux cleaned = -0.00254522 Jy ! Component: 400 - total flux cleaned = -0.00350713 Jy ! Component: 450 - total flux cleaned = -0.00430893 Jy ! Component: 500 - total flux cleaned = -0.00509494 Jy ! Component: 550 - total flux cleaned = -0.00599217 Jy ! Component: 600 - total flux cleaned = -0.00712022 Jy ! Component: 650 - total flux cleaned = -0.00760896 Jy ! Component: 700 - total flux cleaned = -0.00760856 Jy ! Component: 750 - total flux cleaned = -0.00784836 Jy ! Component: 800 - total flux cleaned = -0.00831387 Jy ! Component: 850 - total flux cleaned = -0.00796909 Jy ! Component: 900 - total flux cleaned = -0.00774452 Jy ! Component: 950 - total flux cleaned = -0.00785747 Jy ! Component: 1000 - total flux cleaned = -0.00785698 Jy ! Component: 1050 - total flux cleaned = -0.00785661 Jy ! Component: 1100 - total flux cleaned = -0.00731846 Jy ! Component: 1150 - total flux cleaned = -0.00731885 Jy ! Component: 1200 - total flux cleaned = -0.00710895 Jy ! Component: 1250 - total flux cleaned = -0.00669393 Jy ! Component: 1300 - total flux cleaned = -0.00669399 Jy ! Component: 1350 - total flux cleaned = -0.00659157 Jy ! Component: 1400 - total flux cleaned = -0.00598809 Jy ! Component: 1450 - total flux cleaned = -0.00539365 Jy ! Component: 1500 - total flux cleaned = -0.00519747 Jy ! Component: 1550 - total flux cleaned = -0.00471112 Jy ! Component: 1600 - total flux cleaned = -0.00451801 Jy ! Component: 1650 - total flux cleaned = -0.00385157 Jy ! Component: 1700 - total flux cleaned = -0.0031928 Jy ! Component: 1750 - total flux cleaned = -0.00291081 Jy ! Component: 1800 - total flux cleaned = -0.00198591 Jy ! Component: 1850 - total flux cleaned = -0.00161853 Jy ! Component: 1900 - total flux cleaned = -0.000799872 Jy ! Component: 1950 - total flux cleaned = -0.000350245 Jy ! Component: 2000 - total flux cleaned = 0.000186862 Jy ! Component: 2050 - total flux cleaned = 0.000984255 Jy ! Component: 2100 - total flux cleaned = 0.0015118 Jy ! Component: 2150 - total flux cleaned = 0.00203587 Jy ! Component: 2200 - total flux cleaned = 0.0024683 Jy ! Component: 2250 - total flux cleaned = 0.00281172 Jy ! Component: 2300 - total flux cleaned = 0.00332255 Jy ! Component: 2350 - total flux cleaned = 0.003746 Jy ! Component: 2400 - total flux cleaned = 0.00408236 Jy ! Component: 2450 - total flux cleaned = 0.00466561 Jy ! Component: 2500 - total flux cleaned = 0.00516269 Jy ! Component: 2550 - total flux cleaned = 0.00581992 Jy ! Component: 2600 - total flux cleaned = 0.00647379 Jy ! Component: 2650 - total flux cleaned = 0.00704097 Jy ! Component: 2700 - total flux cleaned = 0.00736324 Jy ! Component: 2750 - total flux cleaned = 0.0079256 Jy ! Component: 2800 - total flux cleaned = 0.00872038 Jy ! Component: 2850 - total flux cleaned = 0.00903733 Jy ! Component: 2900 - total flux cleaned = 0.00982248 Jy ! Component: 2950 - total flux cleaned = 0.010213 Jy ! Component: 3000 - total flux cleaned = 0.0107546 Jy ! Component: 3050 - total flux cleaned = 0.0112162 Jy ! Component: 3100 - total flux cleaned = 0.0119042 Jy ! Component: 3150 - total flux cleaned = 0.0125886 Jy ! Component: 3200 - total flux cleaned = 0.013118 Jy ! Component: 3250 - total flux cleaned = 0.0134187 Jy ! Component: 3300 - total flux cleaned = 0.0140911 Jy ! Component: 3350 - total flux cleaned = 0.0147587 Jy ! Component: 3400 - total flux cleaned = 0.0149802 Jy ! Component: 3450 - total flux cleaned = 0.0157144 Jy ! Component: 3500 - total flux cleaned = 0.016225 Jy ! Component: 3550 - total flux cleaned = 0.0169495 Jy ! Component: 3600 - total flux cleaned = 0.0173826 Jy ! Component: 3650 - total flux cleaned = 0.017598 Jy ! Component: 3700 - total flux cleaned = 0.0183108 Jy ! Component: 3750 - total flux cleaned = 0.0189489 Jy ! Component: 3800 - total flux cleaned = 0.019513 Jy ! Component: 3850 - total flux cleaned = 0.019934 Jy ! Component: 3900 - total flux cleaned = 0.0205614 Jy ! Component: 3950 - total flux cleaned = 0.021185 Jy ! Component: 4000 - total flux cleaned = 0.02153 Jy ! Component: 4050 - total flux cleaned = 0.0222159 Jy ! Total flux subtracted in 4096 components = 0.0223526 Jy ! Clean residual min=-0.004661 max=0.003773 Jy/beam ! Clean residual mean=0.000100 rms=0.001457 Jy/beam ! Combined flux in latest and established models = 3.52823 Jy keep ! Adding 330 model components to the UV plane model. ! The established model now contains 1042 components and 3.52823 Jy ![Exited script file: /home/lpetrov/bin/rdv_muppet] wmap "/scr/pima/rdv62_uvs/J2225-0457_S\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 3.236 x 6.479 at -5.98 degrees (North through East) ! Clean map min=-0.010347 max=2.2027 Jy/beam ! Writing clean map to FITS file: /scr/pima/rdv62_uvs/J2225-0457_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= 2.2026 rms= 0.00145 Jy quit ! Quitting program ! Log file /scr/pima/rdv62_uvs/J2225-0457_S_uvs.log closed on Fri Nov 22 22:42:11 2013