! Started logfile: J0538-4405_S_2003_07_09_pus_uvs.log on Sun Apr 7 19:39:10 2013 obs J0538-4405_S_2003_07_09_pus_uva.fits ! Reading UV FITS file: J0538-4405_S_2003_07_09_pus_uva.fits ! AN table 1: 64 integrations on 171 of 171 possible baselines. ! AN table 2: 33 integrations on 171 of 171 possible baselines. ! AN table 3: 45 integrations on 171 of 171 possible baselines. ! Apparent sampling: 0.0192324 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0538-4405 ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! Warning: Unknown AN-table MNTSTA value (3). ! Will assume that it is ground-based. ! ! 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.23697e+09 8e+06 1 8e+06 ! 02 2 2.26697e+09 8e+06 1 8e+06 ! 03 3 2.35697e+09 8e+06 1 8e+06 ! 04 4 2.37697e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 1516 lines of history. ! ! Reading 1868 visibilities. select I ! Polarization I is unavailable. ! 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 = 512 print "field_cell =",field_cell ! field_cell = 0.5 print "taper_size =",taper_size ! taper_size = 20 ![@muppet J0538-4405_S_2003_07_09_pus] integer clean_niter clean_niter = 100 float clean_gain clean_gain = 0.03 float dynam dynam = 6.0 float soltime1 soltime1 = 120 float soltime2 soltime2 = 20 float thresh thresh = 0.5 float win_mult win_mult = 1.8 float time_av time_av = 16 float old_peak float new_peak float flux_cutoff mapsize field_size, field_cell ! Map grid = 512x512 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 1216 telescope corrections were flagged in sub-array 1. ! A total of 627 telescope corrections were flagged in sub-array 2. ! A total of 700 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 1216 telescope corrections were flagged in sub-array 1. ! A total of 627 telescope corrections were flagged in sub-array 2. ! A total of 700 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 1216 telescope corrections were flagged in sub-array 1. ! A total of 627 telescope corrections were flagged in sub-array 2. ! A total of 700 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 1216 telescope corrections were flagged in sub-array 1. ! A total of 627 telescope corrections were flagged in sub-array 2. ! A total of 700 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=3.135017Jy sigma=0.020097 ! Fit after self-cal, rms=3.291465Jy sigma=0.021958 ! clrmod: Cleared the established, tentative and continuum models. ! Redundant starting model cleared. 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.707 mas, bmaj=4.123 mas, bpa=24.12 degrees ! Estimated noise=6041.89 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=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.13417 Jy ! Component: 100 - total flux cleaned = 3.93975 Jy ! Total flux subtracted in 100 components = 3.93975 Jy ! Clean residual min=-0.245401 max=0.295756 Jy/beam ! Clean residual mean=0.001813 rms=0.092078 Jy/beam ! Combined flux in latest and established models = 3.93975 Jy ! Performing phase self-cal ! Adding 18 model components to the UV plane model. ! The established model now contains 18 components and 3.93975 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=1.041845Jy sigma=0.006263 ! Fit after self-cal, rms=1.041577Jy sigma=0.006262 ! 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.973 mas, bmaj=6.808 mas, bpa=10.9 degrees ! Estimated noise=3215.32 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.460355 Jy ! Component: 100 - total flux cleaned = 0.591157 Jy ! Component: 150 - total flux cleaned = 0.597512 Jy ! Component: 200 - total flux cleaned = 0.588353 Jy ! Total flux subtracted in 200 components = 0.588353 Jy ! Clean residual min=-0.071487 max=0.070970 Jy/beam ! Clean residual mean=0.002267 rms=0.029562 Jy/beam ! Combined flux in latest and established models = 4.5281 Jy ! Performing phase self-cal ! Adding 42 model components to the UV plane model. ! The established model now contains 58 components and 4.5281 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.851413Jy sigma=0.004891 ! Fit after self-cal, rms=0.851135Jy sigma=0.004891 ! Inverting map print "************** FINISHED NATURAL WEIGHTING CLEAN **************" ! ************** FINISHED NATURAL WEIGHTING CLEAN ************** ! gscale true ! Performing overall amplitude self-cal ! ! Correcting IF 1. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.15 GC 1.00* GG 1.00* ! HN 1.00* KP 0.93 LA 1.20 MA 1.00* ! MK 0.95 NL 1.00* NY 1.00* ON 1.00* ! OV 0.90 PT 1.02 SC 1.00* TC 1.03 ! TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 2. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.06 GC 1.00* GG 1.00* ! HN 1.00* KP 1.02 LA 1.14 MA 1.00* ! MK 1.02 NL 1.00* NY 1.00* ON 1.00* ! OV 1.08 PT 1.00 SC 1.00* TC 1.01 ! TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 3. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 1.01 GC 1.00* GG 1.00* ! HN 1.00* KP 1.01 LA 0.95 MA 1.00* ! MK 1.06 NL 1.00* NY 1.00* ON 1.00* ! OV 1.03 PT 1.01 SC 1.00* TC 0.99 ! TS 1.00* WF 1.00* WZ 1.00* ! ! ! Correcting IF 4. ! Telescope amplitude corrections in sub-array 1: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! Telescope amplitude corrections in sub-array 2: ! BR 0.00* FD 0.00* GC 0.00* GG 0.00* ! HN 0.00* KP 0.00* LA 0.00* MA 0.00* ! MK 0.00* NL 0.00* NY 0.00* ON 0.00* ! OV 0.00* PT 0.00* SC 0.00* TC 0.00* ! TS 0.00* WF 0.00* WZ 0.00* ! ! A total of 15 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 0.91 GC 1.00* GG 1.00* ! HN 1.00* KP 1.00 LA 0.77 MA 1.00* ! MK 1.03 NL 1.00* NY 1.00* ON 1.00* ! OV 0.84 PT 1.00 SC 1.00* TC 1.09 ! TS 1.00* WF 1.00* WZ 1.00* ! ! ! Fit before self-cal, rms=0.851135Jy sigma=0.004891 ! Fit after self-cal, rms=0.518317Jy sigma=0.001755 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Estimated beam: bmin=4.078 mas, bmaj=7.171 mas, bpa=8.947 degrees ! Estimated noise=3202.56 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.00254702 Jy ! Component: 100 - total flux cleaned = -0.00373888 Jy ! Component: 150 - total flux cleaned = -0.00785475 Jy ! Component: 200 - total flux cleaned = -0.00876311 Jy ! Total flux subtracted in 200 components = -0.00876311 Jy ! Clean residual min=-0.014473 max=0.014717 Jy/beam ! Clean residual mean=0.000172 rms=0.005476 Jy/beam ! Combined flux in latest and established models = 4.51934 Jy ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 88 components and 4.51934 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.512563Jy sigma=0.001691 ! Fit after self-cal, rms=0.513194Jy sigma=0.001689 ! Inverting map selfcal true, true, soltime1 ! Performing amp+phase self-cal over 120 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.513194Jy sigma=0.001689 ! Fit after self-cal, rms=0.511641Jy sigma=0.001682 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.077 mas, bmaj=7.163 mas, bpa=8.975 degrees ! Estimated noise=3203.92 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00995615 Jy ! Component: 100 - total flux cleaned = -0.0139876 Jy ! Component: 150 - total flux cleaned = -0.0139799 Jy ! Component: 200 - total flux cleaned = -0.013979 Jy ! Total flux subtracted in 200 components = -0.013979 Jy ! Clean residual min=-0.011274 max=0.011165 Jy/beam ! Clean residual mean=-0.000660 rms=0.003867 Jy/beam ! Combined flux in latest and established models = 4.50536 Jy selfcal ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 117 components and 4.50536 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.509299Jy sigma=0.001664 ! Fit after self-cal, rms=0.509240Jy sigma=0.001664 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.00251707 Jy ! Component: 100 - total flux cleaned = -0.00693424 Jy ! Component: 150 - total flux cleaned = -0.00890881 Jy ! Component: 200 - total flux cleaned = -0.0103014 Jy ! Total flux subtracted in 200 components = -0.0103014 Jy ! Clean residual min=-0.007508 max=0.007496 Jy/beam ! Clean residual mean=-0.000050 rms=0.002603 Jy/beam ! Combined flux in latest and established models = 4.49506 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 149 components and 4.49506 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.507602Jy sigma=0.001654 ! Fit after self-cal, rms=0.507630Jy sigma=0.001654 ! Inverting map selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.507630Jy sigma=0.001654 ! Fit after self-cal, rms=0.507480Jy sigma=0.001649 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=4.078 mas, bmaj=7.166 mas, bpa=8.957 degrees ! Estimated noise=3203.61 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0103999 Jy ! Component: 100 - total flux cleaned = 0.012717 Jy ! Component: 150 - total flux cleaned = 0.0151912 Jy ! Component: 200 - total flux cleaned = 0.0167528 Jy ! Total flux subtracted in 200 components = 0.0167528 Jy ! Clean residual min=-0.006346 max=0.006373 Jy/beam ! Clean residual mean=-0.000222 rms=0.002285 Jy/beam ! Combined flux in latest and established models = 4.51181 Jy selfcal ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 177 components and 4.51181 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.506291Jy sigma=0.001641 ! Fit after self-cal, rms=0.506286Jy sigma=0.001641 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.0034047 Jy ! Component: 100 - total flux cleaned = -0.0027202 Jy ! Component: 150 - total flux cleaned = -0.00303713 Jy ! Component: 200 - total flux cleaned = -0.00302765 Jy ! Total flux subtracted in 200 components = -0.00302765 Jy ! Clean residual min=-0.004960 max=0.004940 Jy/beam ! Clean residual mean=0.000048 rms=0.001820 Jy/beam ! Combined flux in latest and established models = 4.50878 Jy ! Performing phase self-cal ! Adding 31 model components to the UV plane model. ! The established model now contains 194 components and 4.50878 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.505356Jy sigma=0.001637 ! Fit after self-cal, rms=0.505348Jy sigma=0.001637 ! 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 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.505348Jy sigma=0.001637 ! Fit after self-cal, rms=0.504954Jy sigma=0.001635 clean clean_niter,clean_gain ! Inverting map and beam ! Estimated beam: bmin=11.5 mas, bmaj=46.99 mas, bpa=12.63 degrees ! Estimated noise=3437.26 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00599074 Jy ! Component: 100 - total flux cleaned = 0.00671575 Jy ! Component: 150 - total flux cleaned = 0.00693271 Jy ! Component: 200 - total flux cleaned = 0.00672164 Jy ! Total flux subtracted in 200 components = 0.00672164 Jy ! Clean residual min=-0.003314 max=0.003338 Jy/beam ! Clean residual mean=0.000085 rms=0.001055 Jy/beam ! Combined flux in latest and established models = 4.5155 Jy selfcal ! Performing phase self-cal ! Adding 39 model components to the UV plane model. ! The established model now contains 227 components and 4.5155 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.504803Jy sigma=0.001633 ! Fit after self-cal, rms=0.504794Jy sigma=0.001633 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.000695994 Jy ! Component: 100 - total flux cleaned = 0.000104923 Jy ! Component: 150 - total flux cleaned = 0.00123994 Jy ! Component: 200 - total flux cleaned = 0.00254836 Jy ! Total flux subtracted in 200 components = 0.00254836 Jy ! Clean residual min=-0.001565 max=0.001550 Jy/beam ! Clean residual mean=0.000098 rms=0.000674 Jy/beam ! Combined flux in latest and established models = 4.51805 Jy ! Performing phase self-cal ! Adding 36 model components to the UV plane model. ! The established model now contains 259 components and 4.51805 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.504784Jy sigma=0.001632 ! Fit after self-cal, rms=0.504787Jy sigma=0.001632 ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. win_mult = win_mult / 1.2 print "************** FINISHED AMPLITUDE SELF-CAL **************" ! ************** FINISHED AMPLITUDE SELF-CAL ************** restore ! Inverting map and beam ! Estimated beam: bmin=4.078 mas, bmaj=7.165 mas, bpa=8.964 degrees ! Estimated noise=3203.5 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.078 x 7.165 at 8.964 degrees (North through East) ! Clean map min=-0.036159 max=3.4292 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.504787Jy sigma=0.001632 ! Fit after self-cal, rms=0.466639Jy sigma=0.001119 ! Inverting map and beam ! Estimated beam: bmin=4.075 mas, bmaj=7.167 mas, bpa=9.008 degrees ! Estimated noise=3202.71 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00432918 Jy ! Component: 100 - total flux cleaned = -0.00573816 Jy ! Component: 150 - total flux cleaned = -0.00796705 Jy ! Component: 200 - total flux cleaned = -0.0094588 Jy ! Total flux subtracted in 200 components = -0.0094588 Jy ! Clean residual min=-0.004819 max=0.004821 Jy/beam ! Clean residual mean=-0.000146 rms=0.001704 Jy/beam ! Combined flux in latest and established models = 4.50859 Jy ! Performing phase self-cal ! Adding 29 model components to the UV plane model. ! The established model now contains 281 components and 4.50859 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.465409Jy sigma=0.001114 ! Fit after self-cal, rms=0.465429Jy sigma=0.001114 ! Inverting map ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.000657877 Jy ! Component: 100 - total flux cleaned = 0.000120896 Jy ! Component: 150 - total flux cleaned = -0.000122437 Jy ! Component: 200 - total flux cleaned = -0.000602602 Jy ! Total flux subtracted in 200 components = -0.000602602 Jy ! Clean residual min=-0.004036 max=0.003954 Jy/beam ! Clean residual mean=-0.000072 rms=0.001455 Jy/beam ! Combined flux in latest and established models = 4.50799 Jy ! Performing phase self-cal ! Adding 27 model components to the UV plane model. ! The established model now contains 298 components and 4.50799 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.464424Jy sigma=0.001112 ! Fit after self-cal, rms=0.464413Jy sigma=0.001112 ! 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.464413Jy sigma=0.001112 ! Fit after self-cal, rms=0.464579Jy sigma=0.001111 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.729 mas, bmaj=4.162 mas, bpa=23.36 degrees ! Estimated noise=6315.57 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 3.10246 Jy ! Component: 100 - total flux cleaned = 3.88942 Jy ! Component: 150 - total flux cleaned = 4.19734 Jy ! Component: 200 - total flux cleaned = 4.35856 Jy ! Total flux subtracted in 200 components = 4.35856 Jy ! Clean residual min=-0.063058 max=0.076773 Jy/beam ! Clean residual mean=0.000410 rms=0.024647 Jy/beam ! Combined flux in latest and established models = 4.35856 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 43 model components to the UV plane model. ! The established model now contains 43 components and 4.35856 Jy ! Inverting map and beam ! Estimated beam: bmin=4.075 mas, bmaj=7.169 mas, bpa=9 degrees ! Estimated noise=3203.07 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=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.132622 Jy ! Component: 100 - total flux cleaned = 0.196052 Jy ! Component: 150 - total flux cleaned = 0.217533 Jy ! Component: 200 - total flux cleaned = 0.215089 Jy ! Total flux subtracted in 200 components = 0.215089 Jy ! Clean residual min=-0.019552 max=0.019575 Jy/beam ! Clean residual mean=-0.001945 rms=0.008673 Jy/beam ! Combined flux in latest and established models = 4.57365 Jy ! Adding 45 model components to the UV plane model. ! The established model now contains 85 components and 4.57365 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=11.49 mas, bmaj=47 mas, bpa=12.65 degrees ! Estimated noise=3437.94 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00478767 Jy ! Component: 100 - total flux cleaned = -0.0149374 Jy ! Component: 150 - total flux cleaned = -0.0240948 Jy ! Component: 200 - total flux cleaned = -0.0322989 Jy ! Total flux subtracted in 200 components = -0.0322989 Jy ! Clean residual min=-0.005955 max=0.005976 Jy/beam ! Clean residual mean=-0.001800 rms=0.002688 Jy/beam ! Combined flux in latest and established models = 4.54135 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=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00661417 Jy ! Component: 100 - total flux cleaned = -0.0127754 Jy ! Component: 150 - total flux cleaned = -0.0175314 Jy ! Component: 200 - total flux cleaned = -0.021433 Jy ! Total flux subtracted in 200 components = -0.021433 Jy ! Clean residual min=-0.003917 max=0.003907 Jy/beam ! Clean residual mean=-0.001173 rms=0.001612 Jy/beam ! Combined flux in latest and established models = 4.51991 Jy ! Adding 103 model components to the UV plane model. ! The established model now contains 187 components and 4.51992 Jy ! Inverting map uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Estimated beam: bmin=4.075 mas, bmaj=7.169 mas, bpa=9 degrees ! Estimated noise=3203.07 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = -0.00654624 Jy ! Component: 100 - total flux cleaned = -0.0101271 Jy ! Component: 150 - total flux cleaned = -0.0150341 Jy ! Component: 200 - total flux cleaned = -0.0195281 Jy ! Total flux subtracted in 200 components = -0.0195281 Jy ! Clean residual min=-0.008904 max=0.008908 Jy/beam ! Clean residual mean=-0.000164 rms=0.003360 Jy/beam ! Combined flux in latest and established models = 4.50039 Jy selfcal ! Performing phase self-cal ! Adding 24 model components to the UV plane model. ! The established model now contains 210 components and 4.50039 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.472914Jy sigma=0.001140 ! Fit after self-cal, rms=0.472263Jy sigma=0.001139 wmodel J0538-4405_S_2003_07_09_pus_map.mod ! Writing 210 model components to file: J0538-4405_S_2003_07_09_pus_map.mod wobs J0538-4405_S_2003_07_09_pus_uvs.fits ! Writing UV FITS file: J0538-4405_S_2003_07_09_pus_uvs.fits wwins J0538-4405_S_2003_07_09_pus_map.win ! wwins: No CLEAN windows to write. ! Error occured in command: wwins ![Exited script file: muppet] ! wmap "J0538-4405_S_2003_07_09_pus\_map.fits" ! Inverting map ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.075 x 7.169 at 9 degrees (North through East) ! Clean map min=-0.018962 max=3.6008 Jy/beam ! Writing clean map to FITS file: J0538-4405_S_2003_07_09_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 3.60077 invert ! Inverting map float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00329282 device "J0538-4405_S_2003_07_09_pus\_map.eps/VPS" ! Attempting to open device: 'J0538-4405_S_2003_07_09_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.274343 0.274343 0.548686 1.09737 2.19475 4.38949 8.77898 17.558 35.1159 70.2318 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0098 print "peak=",int(1000*peak_flux)/1000 ! peak= 3.6 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 4.075 x 7.169 at 9 degrees (North through East) ! Clean map min=-0.018962 max=3.6008 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. quit ! Quitting program ! Log file J0538-4405_S_2003_07_09_pus_uvs.log closed on Sun Apr 7 19:39:14 2013