! Started logfile: J0241-0815_X_1999_03_08_pus_uvs.log on Thu Jan 3 14:33:50 2008 obs J0241-0815_X_1999_03_08_pus_uva.fits ! Reading UV FITS file: J0241-0815_X_1999_03_08_pus_uva.fits ! AN table 1: 40 integrations on 190 of 190 possible baselines. ! AN table 2: 116 integrations on 190 of 190 possible baselines. ! AN table 3: 159 integrations on 190 of 190 possible baselines. ! AN table 4: 63 integrations on 190 of 190 possible baselines. ! Apparent sampling: 0.0752715 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J0241-0815 ! ! There are 4 IFs, and a total of 4 channels: ! ! IF Channel Frequency Freq offset Number of Overall IF ! origin at origin per channel channels bandwidth ! ------------------------------------------------------------- (Hz) ! 01 1 8.40997e+09 8e+06 1 8e+06 ! 02 2 8.47997e+09 8e+06 1 8e+06 ! 03 3 8.79497e+09 8e+06 1 8e+06 ! 04 4 8.89997e+09 8e+06 1 8e+06 ! ! Polarization(s): RR ! ! Read 2464 lines of history. ! ! Reading 21624 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 = 256 print "field_cell =",field_cell ! field_cell = 0.2 print "taper_size =",taper_size ! taper_size = 70 ![@muppet J0241-0815_X_1999_03_08_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 = 256x256 pixels with 0.200x0.200 milli-arcsec cellsize. startmod ! Applying default point source starting model. ! Performing phase self-cal ! Adding 1 model components to the UV plane model. ! The established model now contains 1 components and 1 Jy ! ! Correcting IF 1. ! A total of 800 telescope corrections were flagged in sub-array 1. ! A total of 1809 telescope corrections were flagged in sub-array 2. ! A total of 2334 telescope corrections were flagged in sub-array 3. ! A total of 954 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 2. ! A total of 800 telescope corrections were flagged in sub-array 1. ! A total of 1809 telescope corrections were flagged in sub-array 2. ! A total of 2334 telescope corrections were flagged in sub-array 3. ! A total of 954 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 3. ! A total of 800 telescope corrections were flagged in sub-array 1. ! A total of 1809 telescope corrections were flagged in sub-array 2. ! A total of 2336 telescope corrections were flagged in sub-array 3. ! A total of 954 telescope corrections were flagged in sub-array 4. ! ! Correcting IF 4. ! A total of 800 telescope corrections were flagged in sub-array 1. ! A total of 1809 telescope corrections were flagged in sub-array 2. ! A total of 2335 telescope corrections were flagged in sub-array 3. ! A total of 954 telescope corrections were flagged in sub-array 4. ! ! Fit before self-cal, rms=1.016932Jy sigma=0.008181 ! Fit after self-cal, rms=0.809057Jy sigma=0.006539 ! 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 ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.6171 mas, bmaj=1.574 mas, bpa=-4.626 degrees ! Estimated noise=1135.25 mJy/beam. repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) ! Added new window around map position (0, 0). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.29119 Jy ! Component: 100 - total flux cleaned = 0.433254 Jy ! Total flux subtracted in 100 components = 0.433254 Jy ! Clean residual min=-0.041703 max=0.154095 Jy/beam ! Clean residual mean=0.007825 rms=0.019490 Jy/beam ! Combined flux in latest and established models = 0.433253 Jy ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 13 components and 0.433253 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.553119Jy sigma=0.004655 ! Fit after self-cal, rms=0.550082Jy sigma=0.004641 ! Inverting map ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Added new window around map position (1, 0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.140761 Jy ! Component: 100 - total flux cleaned = 0.228125 Jy ! Total flux subtracted in 100 components = 0.228125 Jy ! Clean residual min=-0.039143 max=0.118336 Jy/beam ! Clean residual mean=0.005945 rms=0.015236 Jy/beam ! Combined flux in latest and established models = 0.661379 Jy ! Performing phase self-cal ! Adding 17 model components to the UV plane model. ! The established model now contains 30 components and 0.661379 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.475316Jy sigma=0.003735 ! Fit after self-cal, rms=0.471910Jy sigma=0.003704 ! Inverting map ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Added new window around map position (-1.2, -0.2). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.102174 Jy ! Component: 100 - total flux cleaned = 0.167312 Jy ! Total flux subtracted in 100 components = 0.167312 Jy ! Clean residual min=-0.026092 max=0.093994 Jy/beam ! Clean residual mean=0.004487 rms=0.012445 Jy/beam ! Combined flux in latest and established models = 0.828691 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 46 components and 0.828691 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.426625Jy sigma=0.003122 ! Fit after self-cal, rms=0.425738Jy sigma=0.003107 ! Inverting map ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Added new window around map position (2.4, 0.6). ! clean: niter=100 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0867394 Jy ! Component: 100 - total flux cleaned = 0.143033 Jy ! Total flux subtracted in 100 components = 0.143033 Jy ! Clean residual min=-0.026867 max=0.061428 Jy/beam ! Clean residual mean=0.003428 rms=0.010020 Jy/beam ! Combined flux in latest and established models = 0.971724 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 56 components and 0.971724 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.386609Jy sigma=0.002661 ! Fit after self-cal, rms=0.383408Jy sigma=0.002642 ! Inverting map ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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 ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7291 mas, bmaj=2.013 mas, bpa=-5.321 degrees ! Estimated noise=590.574 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.113412 Jy ! Component: 100 - total flux cleaned = 0.188731 Jy ! Component: 150 - total flux cleaned = 0.241615 Jy ! Component: 200 - total flux cleaned = 0.280215 Jy ! Total flux subtracted in 200 components = 0.280215 Jy ! Clean residual min=-0.034797 max=0.095095 Jy/beam ! Clean residual mean=0.005862 rms=0.019811 Jy/beam ! Combined flux in latest and established models = 1.25194 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 83 components and 1.25194 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.339858Jy sigma=0.002155 ! Fit after self-cal, rms=0.336695Jy sigma=0.002131 ! Inverting map ! Your choice of large map pixels excluded 0.302% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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* NL 0.00* ! OV 0.00* PT 0.00* SC 0.00* WF 0.00* ! GN 0.00* HH 0.00* KK 0.00* MC 0.00* ! MK 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 60 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* NL 1.00* ! OV 1.00* PT 1.00* SC 0.87 WF 1.24 ! GN 1.74 HH 1.00* KK 1.00* MC 0.87 ! MK 1.00* NY 1.00* ON 0.89 WZ 0.73 ! ! A total of 51 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.95 FD 0.98 GC 1.10 GG 1.00* ! HN 0.93 KP 1.16 LA 0.88 NL 1.16 ! OV 1.15 PT 0.98 SC 1.13 WF 0.88 ! GN 1.24 HH 1.00* KK 0.90 MC 1.00* ! MK 0.97 NY 1.00* ON 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.17 FD 0.97 GC 1.00* GG 1.00* ! HN 0.92 KP 1.00* LA 1.00* NL 1.22 ! OV 1.00* PT 1.00* SC 0.97 WF 0.86 ! GN 1.08 HH 1.00* KK 1.00* MC 1.00* ! MK 1.00* NY 1.00* ON 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* NL 0.00* ! OV 0.00* PT 0.00* SC 0.00* WF 0.00* ! GN 0.00* HH 0.00* KK 0.00* MC 0.00* ! MK 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 60 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* NL 1.00* ! OV 1.00* PT 1.00* SC 1.13 WF 1.22 ! GN 1.69 HH 1.00* KK 1.00* MC 0.89 ! MK 1.00* NY 1.00* ON 0.88 WZ 0.70 ! ! A total of 51 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.94 FD 0.98 GC 1.09 GG 1.00* ! HN 0.91 KP 1.16 LA 0.88 NL 1.09 ! OV 1.15 PT 0.98 SC 1.13 WF 0.91 ! GN 1.25 HH 1.00* KK 0.87 MC 1.00* ! MK 1.02 NY 1.00* ON 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.14 FD 0.98 GC 1.00* GG 1.00* ! HN 0.87 KP 1.00* LA 1.00* NL 1.20 ! OV 1.00* PT 1.00* SC 0.97 WF 0.89 ! GN 1.19 HH 1.00* KK 1.00* MC 1.00* ! MK 1.00* NY 1.00* ON 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* NL 0.00* ! OV 0.00* PT 0.00* SC 0.00* WF 0.00* ! GN 0.00* HH 0.00* KK 0.00* MC 0.00* ! MK 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 60 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* NL 1.00* ! OV 1.00* PT 1.00* SC 1.44 WF 0.96 ! GN 2.16 HH 1.00* KK 1.00* MC 0.86 ! MK 1.00* NY 1.00* ON 0.90 WZ 0.73 ! ! A total of 51 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.90 FD 0.98 GC 1.04 GG 1.00* ! HN 0.91 KP 1.17 LA 0.89 NL 1.13 ! OV 1.16 PT 0.98 SC 1.18 WF 0.90 ! GN 1.30 HH 1.00* KK 0.85 MC 1.00* ! MK 1.01 NY 1.00* ON 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.15 FD 0.97 GC 1.00* GG 1.00* ! HN 0.87 KP 1.00* LA 1.00* NL 1.13 ! OV 1.00* PT 1.00* SC 1.02 WF 0.90 ! GN 1.14 HH 1.00* KK 1.00* MC 1.00* ! MK 1.00* NY 1.00* ON 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* NL 0.00* ! OV 0.00* PT 0.00* SC 0.00* WF 0.00* ! GN 0.00* HH 0.00* KK 0.00* MC 0.00* ! MK 0.00* NY 0.00* ON 0.00* WZ 0.00* ! ! A total of 60 telescope corrections were flagged in sub-array 2. ! Telescope amplitude corrections in sub-array 2: ! BR 1.00* FD 1.00* GC 1.00* GG 1.00* ! HN 1.00* KP 1.00* LA 1.00* NL 1.00* ! OV 1.00* PT 1.00* SC 1.02 WF 0.77 ! GN 2.58 HH 1.00* KK 1.00* MC 0.85 ! MK 1.00* NY 1.00* ON 0.90 WZ 0.76 ! ! A total of 51 telescope corrections were flagged in sub-array 3. ! Telescope amplitude corrections in sub-array 3: ! BR 0.93 FD 0.98 GC 1.14 GG 1.00* ! HN 0.95 KP 1.18 LA 0.88 NL 1.15 ! OV 1.18 PT 0.97 SC 1.19 WF 0.87 ! GN 1.29 HH 1.00* KK 0.84 MC 1.00* ! MK 1.03 NY 1.00* ON 1.00* WZ 1.00* ! ! A total of 6 telescope corrections were flagged in sub-array 4. ! Telescope amplitude corrections in sub-array 4: ! BR 1.19 FD 0.95 GC 1.00* GG 1.00* ! HN 0.90 KP 1.00* LA 1.00* NL 1.11 ! OV 1.00* PT 1.00* SC 0.99 WF 0.87 ! GN 1.20 HH 1.00* KK 1.00* MC 1.00* ! MK 1.00* NY 1.00* ON 1.00* WZ 1.00* ! ! ! Fit before self-cal, rms=0.336695Jy sigma=0.002131 ! Fit after self-cal, rms=0.302116Jy sigma=0.001674 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7346 mas, bmaj=2.08 mas, bpa=-6.257 degrees ! Estimated noise=599.414 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 (-6.4, -2.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0481305 Jy ! Component: 100 - total flux cleaned = 0.0659523 Jy ! Component: 150 - total flux cleaned = 0.0685312 Jy ! Component: 200 - total flux cleaned = 0.0716726 Jy ! Total flux subtracted in 200 components = 0.0716726 Jy ! Clean residual min=-0.012321 max=0.029298 Jy/beam ! Clean residual mean=0.001432 rms=0.004835 Jy/beam ! Combined flux in latest and established models = 1.32361 Jy ! Performing phase self-cal ! Adding 47 model components to the UV plane model. ! The established model now contains 127 components and 1.32361 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.295584Jy sigma=0.001538 ! Fit after self-cal, rms=0.294386Jy sigma=0.001521 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Added new window around map position (-3.8, -1). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0208433 Jy ! Component: 100 - total flux cleaned = 0.0249354 Jy ! Component: 150 - total flux cleaned = 0.026502 Jy ! Component: 200 - total flux cleaned = 0.027802 Jy ! Total flux subtracted in 200 components = 0.027802 Jy ! Clean residual min=-0.011106 max=0.017303 Jy/beam ! Clean residual mean=0.000874 rms=0.003810 Jy/beam ! Combined flux in latest and established models = 1.35141 Jy ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 144 components and 1.35141 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.290962Jy sigma=0.001473 ! Fit after self-cal, rms=0.290623Jy sigma=0.001468 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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.290623Jy sigma=0.001468 ! Fit after self-cal, rms=0.279501Jy sigma=0.001355 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7356 mas, bmaj=2.066 mas, bpa=-5.932 degrees ! Estimated noise=596.908 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.013729 Jy ! Component: 100 - total flux cleaned = 0.0154556 Jy ! Component: 150 - total flux cleaned = 0.0172442 Jy ! Component: 200 - total flux cleaned = 0.0185037 Jy ! Total flux subtracted in 200 components = 0.0185037 Jy ! Clean residual min=-0.008327 max=0.015982 Jy/beam ! Clean residual mean=0.000604 rms=0.002810 Jy/beam ! Combined flux in latest and established models = 1.36992 Jy selfcal ! Performing phase self-cal ! Adding 30 model components to the UV plane model. ! The established model now contains 158 components and 1.36992 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.278204Jy sigma=0.001336 ! Fit after self-cal, rms=0.277837Jy sigma=0.001334 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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 (3.6, 0.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.013789 Jy ! Component: 100 - total flux cleaned = 0.0165735 Jy ! Component: 150 - total flux cleaned = 0.0168646 Jy ! Component: 200 - total flux cleaned = 0.0165924 Jy ! Total flux subtracted in 200 components = 0.0165924 Jy ! Clean residual min=-0.009338 max=0.009914 Jy/beam ! Clean residual mean=0.000273 rms=0.002534 Jy/beam ! Combined flux in latest and established models = 1.38651 Jy ! Performing phase self-cal ! Adding 34 model components to the UV plane model. ! The established model now contains 177 components and 1.38651 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.276981Jy sigma=0.001317 ! Fit after self-cal, rms=0.276859Jy sigma=0.001316 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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.276859Jy sigma=0.001316 ! Fit after self-cal, rms=0.274313Jy sigma=0.001309 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7362 mas, bmaj=2.067 mas, bpa=-5.897 degrees ! Estimated noise=595.564 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0112939 Jy ! Component: 100 - total flux cleaned = 0.0172418 Jy ! Component: 150 - total flux cleaned = 0.0191525 Jy ! Component: 200 - total flux cleaned = 0.0205801 Jy ! Total flux subtracted in 200 components = 0.0205801 Jy ! Clean residual min=-0.006826 max=0.008526 Jy/beam ! Clean residual mean=0.000380 rms=0.002107 Jy/beam ! Combined flux in latest and established models = 1.40709 Jy selfcal ! Performing phase self-cal ! Adding 43 model components to the UV plane model. ! The established model now contains 193 components and 1.40709 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273733Jy sigma=0.001296 ! Fit after self-cal, rms=0.273656Jy sigma=0.001295 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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.00227932 Jy ! Component: 100 - total flux cleaned = 0.00315573 Jy ! Component: 150 - total flux cleaned = 0.00375823 Jy ! Component: 200 - total flux cleaned = 0.00376613 Jy ! Total flux subtracted in 200 components = 0.00376613 Jy ! Clean residual min=-0.006741 max=0.008078 Jy/beam ! Clean residual mean=0.000318 rms=0.001992 Jy/beam ! Combined flux in latest and established models = 1.41086 Jy ! Performing phase self-cal ! Adding 38 model components to the UV plane model. ! The established model now contains 203 components and 1.41086 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273327Jy sigma=0.001290 ! Fit after self-cal, rms=0.273300Jy sigma=0.001290 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. win_mult = win_mult * 1.2 selfcal true, true, soltime2 ! Performing amp+phase self-cal over 20 minute time intervals ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273300Jy sigma=0.001290 ! Fit after self-cal, rms=0.272898Jy sigma=0.001287 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=2.082 mas, bmaj=3.427 mas, bpa=7.713 degrees ! Estimated noise=796.54 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0122131 Jy ! Component: 100 - total flux cleaned = 0.0141886 Jy ! Component: 150 - total flux cleaned = 0.0152561 Jy ! Component: 200 - total flux cleaned = 0.0159258 Jy ! Total flux subtracted in 200 components = 0.0159258 Jy ! Clean residual min=-0.006764 max=0.010299 Jy/beam ! Clean residual mean=0.000989 rms=0.002820 Jy/beam ! Combined flux in latest and established models = 1.42678 Jy selfcal ! Performing phase self-cal ! Adding 33 model components to the UV plane model. ! The established model now contains 221 components and 1.42678 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272688Jy sigma=0.001280 ! Fit after self-cal, rms=0.272619Jy sigma=0.001279 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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.00165551 Jy ! Component: 100 - total flux cleaned = 0.00210031 Jy ! Component: 150 - total flux cleaned = 0.00210685 Jy ! Component: 200 - total flux cleaned = 0.00211767 Jy ! Total flux subtracted in 200 components = 0.00211767 Jy ! Clean residual min=-0.007025 max=0.009546 Jy/beam ! Clean residual mean=0.000911 rms=0.002533 Jy/beam ! Combined flux in latest and established models = 1.4289 Jy ! Performing phase self-cal ! Adding 37 model components to the UV plane model. ! The established model now contains 237 components and 1.4289 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272534Jy sigma=0.001277 ! Fit after self-cal, rms=0.272547Jy sigma=0.001277 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec 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 ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7364 mas, bmaj=2.066 mas, bpa=-5.84 degrees ! Estimated noise=594.872 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7364 x 2.066 at -5.84 degrees (North through East) ! Clean map min=-0.0053006 max=0.33145 Jy/beam if(peak(flux) > thresh) selfcal true, true clean clean_niter,clean_gain selfcal flux_cutoff = imstat(rms) * dynam repeat if (peak(flux) > flux_cutoff) peakwin win_mult clean clean_niter,clean_gain flux_cutoff = imstat(rms) * dynam selfcal new_peak = peak(flux) until(new_peak<=flux_cutoff) selfcal true, true end if print "************** CLEARING MODEL AND STARTING OVER **********" ! ************** CLEARING MODEL AND STARTING OVER ********** clrmod true ! clrmod: Cleared the established and tentative models. uvw 2,-1 ! Uniform weighting binwidth: 2 (pixels). ! Gridding weights will be scaled by errors raised to the power -1. ! Radial weighting is not currently selected. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.6106 mas, bmaj=1.636 mas, bpa=-5.461 degrees ! Estimated noise=1156.04 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.299188 Jy ! Component: 100 - total flux cleaned = 0.482083 Jy ! Component: 150 - total flux cleaned = 0.620198 Jy ! Component: 200 - total flux cleaned = 0.728753 Jy ! Total flux subtracted in 200 components = 0.728753 Jy ! Clean residual min=-0.024932 max=0.067503 Jy/beam ! Clean residual mean=0.005304 rms=0.012657 Jy/beam ! Combined flux in latest and established models = 0.728753 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 29 model components to the UV plane model. ! The established model now contains 29 components and 0.728753 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7364 mas, bmaj=2.066 mas, bpa=-5.84 degrees ! Estimated noise=594.872 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) ! Added new window around map position (1.6, 0.6). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.167085 Jy ! Component: 100 - total flux cleaned = 0.284364 Jy ! Component: 150 - total flux cleaned = 0.371929 Jy ! Component: 200 - total flux cleaned = 0.439193 Jy ! Total flux subtracted in 200 components = 0.439193 Jy ! Clean residual min=-0.016426 max=0.039652 Jy/beam ! Clean residual mean=0.006106 rms=0.010178 Jy/beam ! Combined flux in latest and established models = 1.16795 Jy ! Adding 65 model components to the UV plane model. ! The established model now contains 80 components and 1.16795 Jy ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec uvtaper 0.3,taper_size ! Gaussian taper: value 0.3 at UV radius = 70 mega-wavelengths. clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=2.082 mas, bmaj=3.427 mas, bpa=7.713 degrees ! Estimated noise=796.54 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.104369 Jy ! Component: 100 - total flux cleaned = 0.165052 Jy ! Component: 150 - total flux cleaned = 0.202135 Jy ! Component: 200 - total flux cleaned = 0.225598 Jy ! Total flux subtracted in 200 components = 0.225598 Jy ! Clean residual min=-0.007520 max=0.014406 Jy/beam ! Clean residual mean=0.002954 rms=0.003897 Jy/beam ! Combined flux in latest and established models = 1.39355 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.0151837 Jy ! Component: 100 - total flux cleaned = 0.0252852 Jy ! Component: 150 - total flux cleaned = 0.030465 Jy ! Component: 200 - total flux cleaned = 0.0318134 Jy ! Total flux subtracted in 200 components = 0.0318134 Jy ! Clean residual min=-0.007394 max=0.009110 Jy/beam ! Clean residual mean=0.001202 rms=0.002524 Jy/beam ! Combined flux in latest and established models = 1.42536 Jy ! Adding 113 model components to the UV plane model. ! The established model now contains 169 components and 1.42536 Jy ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! Estimated beam: bmin=0.7364 mas, bmaj=2.066 mas, bpa=-5.84 degrees ! Estimated noise=594.872 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0043122 Jy ! Component: 100 - total flux cleaned = 0.0049645 Jy ! Component: 150 - total flux cleaned = 0.00498034 Jy ! Component: 200 - total flux cleaned = 0.00497988 Jy ! Total flux subtracted in 200 components = 0.00497988 Jy ! Clean residual min=-0.007214 max=0.007344 Jy/beam ! Clean residual mean=0.000335 rms=0.001849 Jy/beam ! Combined flux in latest and established models = 1.43034 Jy selfcal ! Performing phase self-cal ! Adding 25 model components to the UV plane model. ! The established model now contains 191 components and 1.43034 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.273432Jy sigma=0.001283 ! Fit after self-cal, rms=0.272648Jy sigma=0.001277 wmodel J0241-0815_X_1999_03_08_pus_map.mod ! Writing 191 model components to file: J0241-0815_X_1999_03_08_pus_map.mod wobs J0241-0815_X_1999_03_08_pus_uvs.fits ! Writing UV FITS file: J0241-0815_X_1999_03_08_pus_uvs.fits wwins J0241-0815_X_1999_03_08_pus_map.win ! wwins: Wrote 8 windows to J0241-0815_X_1999_03_08_pus_map.win x = (field_size-8) * field_cell / 4 addwin -x,x,-x,x clean (field_size*4),0.01 ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! clean: niter=1024 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = 0.000696697 Jy ! Component: 100 - total flux cleaned = 0.00154104 Jy ! Component: 150 - total flux cleaned = 0.00224037 Jy ! Component: 200 - total flux cleaned = 0.00281083 Jy ! Component: 250 - total flux cleaned = 0.00345276 Jy ! Component: 300 - total flux cleaned = 0.0042462 Jy ! Component: 350 - total flux cleaned = 0.00492417 Jy ! Component: 400 - total flux cleaned = 0.00549714 Jy ! Component: 450 - total flux cleaned = 0.00597314 Jy ! Component: 500 - total flux cleaned = 0.00659019 Jy ! Component: 550 - total flux cleaned = 0.0071152 Jy ! Component: 600 - total flux cleaned = 0.00770003 Jy ! Component: 650 - total flux cleaned = 0.00812811 Jy ! Component: 700 - total flux cleaned = 0.008686 Jy ! Component: 750 - total flux cleaned = 0.00916374 Jy ! Component: 800 - total flux cleaned = 0.00949926 Jy ! Component: 850 - total flux cleaned = 0.00989443 Jy ! Component: 900 - total flux cleaned = 0.0102185 Jy ! Component: 950 - total flux cleaned = 0.0106658 Jy ! Component: 1000 - total flux cleaned = 0.0109802 Jy ! Total flux subtracted in 1024 components = 0.0111049 Jy ! Clean residual min=-0.003197 max=0.003731 Jy/beam ! Clean residual mean=0.000163 rms=0.001320 Jy/beam ! Combined flux in latest and established models = 1.44144 Jy keep ! Adding 70 model components to the UV plane model. ! The established model now contains 258 components and 1.44144 Jy ![Exited script file: muppet] wmap "J0241-0815_X_1999_03_08_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7364 x 2.066 at -5.84 degrees (North through East) ! Clean map min=-0.0071071 max=0.32696 Jy/beam ! Writing clean map to FITS file: J0241-0815_X_1999_03_08_pus_map.fits float peak_flux peak_flux = peak(flux,max) print "The peak flux is",peak_flux ! The peak flux is 0.32696 invert ! Inverting map ! Your choice of large map pixels excluded 0.309% of the data. ! The x-axis pixel size should ideally be below 0.1919 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.00131686 device "J0241-0815_X_1999_03_08_pus\_map.eps/VPS" ! Attempting to open device: 'J0241-0815_X_1999_03_08_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -1.20828 1.20828 2.41655 4.8331 9.66621 19.3324 38.6648 77.3297 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0039 print "peak=",int(1000*peak_flux)/1000 ! peak= 0.326 mapcolor none ! Mapplot colormap: none, contrast: 1 brightness: 0.5. mapplot cln ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7364 x 2.066 at -5.84 degrees (North through East) ! Clean map min=-0.0071071 max=0.32696 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J0241-0815_X_1999_03_08_pus\_map.gif/vgif" ! Attempting to open device: 'J0241-0815_X_1999_03_08_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J0241-0815_X_1999_03_08_pus_uvs.log closed on Thu Jan 3 14:34:10 2008