! Started logfile: J1427-4206_X_1999_03_08_pus_uvs.log on Thu Jan 3 18:36:42 2008 obs J1427-4206_X_1999_03_08_pus_uva.fits ! Reading UV FITS file: J1427-4206_X_1999_03_08_pus_uva.fits ! AN table 1: 36 integrations on 190 of 190 possible baselines. ! AN table 2: 9 integrations on 190 of 190 possible baselines. ! AN table 3: Unused. ! AN table 4: 25 integrations on 190 of 190 possible baselines. ! AN table 5: Unused. ! AN table 6: Unused. ! Apparent sampling: 0.0506767 visibilities/baseline/integration-bin. ! *** This seems a bit low - see "help observe" on the binwid argument. ! Found source: J1427-4206 ! ! 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 2585 lines of history. ! ! Reading 2696 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 J1427-4206_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 720 telescope corrections were flagged in sub-array 1. ! A total of 180 telescope corrections were flagged in sub-array 2. ! A total of 319 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 2. ! A total of 720 telescope corrections were flagged in sub-array 1. ! A total of 180 telescope corrections were flagged in sub-array 2. ! A total of 319 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 3. ! A total of 720 telescope corrections were flagged in sub-array 1. ! A total of 180 telescope corrections were flagged in sub-array 2. ! A total of 319 telescope corrections were flagged in sub-array 3. ! ! Correcting IF 4. ! A total of 720 telescope corrections were flagged in sub-array 1. ! A total of 180 telescope corrections were flagged in sub-array 2. ! A total of 319 telescope corrections were flagged in sub-array 3. ! ! Fit before self-cal, rms=1.469661Jy sigma=0.013708 ! Fit after self-cal, rms=0.865752Jy sigma=0.009038 ! 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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.5776 mas, bmaj=1.91 mas, bpa=-10 degrees ! Estimated noise=2329.27 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.27173 Jy ! Component: 100 - total flux cleaned = 1.59615 Jy ! Total flux subtracted in 100 components = 1.59615 Jy ! Clean residual min=-0.086741 max=0.175636 Jy/beam ! Clean residual mean=0.001402 rms=0.038485 Jy/beam ! Combined flux in latest and established models = 1.59615 Jy ! Performing phase self-cal ! Adding 6 model components to the UV plane model. ! The established model now contains 6 components and 1.59615 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.447225Jy sigma=0.004652 ! Fit after self-cal, rms=0.447404Jy sigma=0.004648 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.795 mas, bmaj=2.584 mas, bpa=-14.43 degrees ! Estimated noise=1580.51 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.270404 Jy ! Component: 100 - total flux cleaned = 0.389552 Jy ! Component: 150 - total flux cleaned = 0.445008 Jy ! Component: 200 - total flux cleaned = 0.468033 Jy ! Total flux subtracted in 200 components = 0.468033 Jy ! Clean residual min=-0.062337 max=0.147500 Jy/beam ! Clean residual mean=0.002122 rms=0.035195 Jy/beam ! Combined flux in latest and established models = 2.06418 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 15 components and 2.06418 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.339204Jy sigma=0.002664 ! Fit after self-cal, rms=0.339006Jy sigma=0.002663 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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* ! ! 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* 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* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 0.94 GC 1.00* GG 1.00* ! HN 1.00* KP 0.96 LA 0.98 NL 1.11 ! OV 0.99 PT 0.94 SC 0.94 WF 1.00* ! GN 1.37 HH 1.00* KK 1.00* MC 1.00* ! MK 1.01 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* ! ! 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* 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* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 0.94 GC 1.00* GG 1.00* ! HN 1.00* KP 0.95 LA 0.98 NL 1.11 ! OV 0.99 PT 0.94 SC 0.97 WF 1.00* ! GN 1.34 HH 1.00* KK 1.00* MC 1.00* ! MK 1.02 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* ! ! 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* 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* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 0.93 GC 1.00* GG 1.00* ! HN 1.00* KP 0.96 LA 0.98 NL 1.06 ! OV 1.00 PT 0.95 SC 1.00 WF 1.00* ! GN 1.33 HH 1.00* KK 1.00* MC 1.00* ! MK 0.98 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* ! ! 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* 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* ! ! Telescope amplitude corrections in sub-array 3: ! BR 1.00* FD 0.93 GC 1.00* GG 1.00* ! HN 1.00* KP 0.96 LA 0.98 NL 1.20 ! OV 0.99 PT 0.93 SC 1.00 WF 1.00* ! GN 1.32 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.339006Jy sigma=0.002663 ! Fit after self-cal, rms=0.283672Jy sigma=0.001632 dynam = dynam - 0.5 flux_cutoff = imstat(rms) * dynam ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.809 mas, bmaj=2.56 mas, bpa=-14.61 degrees ! Estimated noise=1596.67 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 (1.6, 1.8). ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0426357 Jy ! Component: 100 - total flux cleaned = 0.0434556 Jy ! Component: 150 - total flux cleaned = 0.0434453 Jy ! Component: 200 - total flux cleaned = 0.0439953 Jy ! Total flux subtracted in 200 components = 0.0439953 Jy ! Clean residual min=-0.018414 max=0.021839 Jy/beam ! Clean residual mean=-0.000790 rms=0.005636 Jy/beam ! Combined flux in latest and established models = 2.10818 Jy ! Performing phase self-cal ! Adding 10 model components to the UV plane model. ! The established model now contains 24 components and 2.10818 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.272385Jy sigma=0.001529 ! Fit after self-cal, rms=0.271857Jy sigma=0.001519 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.271857Jy sigma=0.001519 ! Fit after self-cal, rms=0.263005Jy sigma=0.001504 dynam = dynam - 0.75 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.8034 mas, bmaj=2.546 mas, bpa=-14.61 degrees ! Estimated noise=1587.92 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.020047 Jy ! Component: 100 - total flux cleaned = 0.0226559 Jy ! Component: 150 - total flux cleaned = 0.0220205 Jy ! Component: 200 - total flux cleaned = 0.0220361 Jy ! Total flux subtracted in 200 components = 0.0220361 Jy ! Clean residual min=-0.016537 max=0.012425 Jy/beam ! Clean residual mean=-0.000697 rms=0.004306 Jy/beam ! Combined flux in latest and established models = 2.13021 Jy selfcal ! Performing phase self-cal ! Adding 21 model components to the UV plane model. ! The established model now contains 39 components and 2.13021 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.260432Jy sigma=0.001483 ! Fit after self-cal, rms=0.260036Jy sigma=0.001480 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.00298739 Jy ! Component: 100 - total flux cleaned = 0.00326475 Jy ! Component: 150 - total flux cleaned = 0.00326573 Jy ! Component: 200 - total flux cleaned = 0.00283269 Jy ! Total flux subtracted in 200 components = 0.00283269 Jy ! Clean residual min=-0.016018 max=0.010390 Jy/beam ! Clean residual mean=-0.000738 rms=0.004144 Jy/beam ! Combined flux in latest and established models = 2.13305 Jy ! Performing phase self-cal ! Adding 15 model components to the UV plane model. ! The established model now contains 45 components and 2.13305 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.259311Jy sigma=0.001474 ! Fit after self-cal, rms=0.259224Jy sigma=0.001473 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.259224Jy sigma=0.001473 ! Fit after self-cal, rms=0.255747Jy sigma=0.001469 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.8009 mas, bmaj=2.539 mas, bpa=-14.58 degrees ! Estimated noise=1584.77 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00957794 Jy ! Component: 100 - total flux cleaned = 0.0115902 Jy ! Component: 150 - total flux cleaned = 0.0112337 Jy ! Component: 200 - total flux cleaned = 0.0109085 Jy ! Total flux subtracted in 200 components = 0.0109085 Jy ! Clean residual min=-0.015241 max=0.009500 Jy/beam ! Clean residual mean=-0.000590 rms=0.003837 Jy/beam ! Combined flux in latest and established models = 2.14396 Jy selfcal ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 54 components and 2.14396 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.254962Jy sigma=0.001465 ! Fit after self-cal, rms=0.254808Jy sigma=0.001464 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.00197379 Jy ! Component: 100 - total flux cleaned = 0.00197724 Jy ! Component: 150 - total flux cleaned = 0.0019737 Jy ! Component: 200 - total flux cleaned = 0.00185339 Jy ! Total flux subtracted in 200 components = 0.00185339 Jy ! Clean residual min=-0.014842 max=0.009262 Jy/beam ! Clean residual mean=-0.000618 rms=0.003781 Jy/beam ! Combined flux in latest and established models = 2.14581 Jy ! Performing phase self-cal ! Adding 12 model components to the UV plane model. ! The established model now contains 57 components and 2.14581 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.254549Jy sigma=0.001462 ! Fit after self-cal, rms=0.254509Jy sigma=0.001461 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.254509Jy sigma=0.001461 ! Fit after self-cal, rms=0.252911Jy sigma=0.001460 clean clean_niter,clean_gain ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=2.343 mas, bmaj=9.72 mas, bpa=-3.157 degrees ! Estimated noise=1939.34 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00551533 Jy ! Component: 100 - total flux cleaned = 0.00564988 Jy ! Component: 150 - total flux cleaned = 0.0057664 Jy ! Component: 200 - total flux cleaned = 0.00576809 Jy ! Total flux subtracted in 200 components = 0.00576809 Jy ! Clean residual min=-0.014665 max=0.012407 Jy/beam ! Clean residual mean=-0.001024 rms=0.006061 Jy/beam ! Combined flux in latest and established models = 2.15158 Jy selfcal ! Performing phase self-cal ! Adding 13 model components to the UV plane model. ! The established model now contains 69 components and 2.15158 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.252605Jy sigma=0.001460 ! Fit after self-cal, rms=0.252530Jy sigma=0.001459 flux_cutoff = imstat(rms) * dynam ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.000966592 Jy ! Component: 100 - total flux cleaned = 0.00108126 Jy ! Component: 150 - total flux cleaned = 0.00116537 Jy ! Component: 200 - total flux cleaned = 0.00124239 Jy ! Total flux subtracted in 200 components = 0.00124239 Jy ! Clean residual min=-0.014591 max=0.012481 Jy/beam ! Clean residual mean=-0.001061 rms=0.006118 Jy/beam ! Combined flux in latest and established models = 2.15282 Jy ! Performing phase self-cal ! Adding 14 model components to the UV plane model. ! The established model now contains 73 components and 2.15282 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.252423Jy sigma=0.001459 ! Fit after self-cal, rms=0.252429Jy sigma=0.001459 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.7997 mas, bmaj=2.536 mas, bpa=-14.57 degrees ! Estimated noise=1584.38 mJy/beam. ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7997 x 2.536 at -14.57 degrees (North through East) ! Clean map min=-0.013525 max=1.7423 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.252429Jy sigma=0.001459 ! Fit after self-cal, rms=0.226360Jy sigma=0.001273 ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.7984 mas, bmaj=2.533 mas, bpa=-14.6 degrees ! Estimated noise=1580.96 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.00562223 Jy ! Component: 100 - total flux cleaned = 0.0059119 Jy ! Component: 150 - total flux cleaned = 0.00591482 Jy ! Component: 200 - total flux cleaned = 0.00591635 Jy ! Total flux subtracted in 200 components = 0.00591635 Jy ! Clean residual min=-0.012786 max=0.009278 Jy/beam ! Clean residual mean=-0.000534 rms=0.003459 Jy/beam ! Combined flux in latest and established models = 2.15874 Jy ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 87 components and 2.15874 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225993Jy sigma=0.001271 ! Fit after self-cal, rms=0.225927Jy sigma=0.001271 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0017742 Jy ! Component: 100 - total flux cleaned = 0.00177284 Jy ! Component: 150 - total flux cleaned = 0.00177364 Jy ! Component: 200 - total flux cleaned = 0.00177605 Jy ! Total flux subtracted in 200 components = 0.00177605 Jy ! Clean residual min=-0.012762 max=0.009378 Jy/beam ! Clean residual mean=-0.000560 rms=0.003457 Jy/beam ! Combined flux in latest and established models = 2.16051 Jy ! Performing phase self-cal ! Adding 16 model components to the UV plane model. ! The established model now contains 89 components and 2.16051 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225792Jy sigma=0.001270 ! Fit after self-cal, rms=0.225770Jy sigma=0.001270 ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Performing amp+phase self-cal ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225770Jy sigma=0.001270 ! Fit after self-cal, rms=0.224854Jy sigma=0.001269 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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.5724 mas, bmaj=1.87 mas, bpa=-9.611 degrees ! Estimated noise=2326.48 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 1.27022 Jy ! Component: 100 - total flux cleaned = 1.60593 Jy ! Component: 150 - total flux cleaned = 1.7728 Jy ! Component: 200 - total flux cleaned = 1.87425 Jy ! Total flux subtracted in 200 components = 1.87425 Jy ! Clean residual min=-0.032484 max=0.055664 Jy/beam ! Clean residual mean=0.000427 rms=0.015322 Jy/beam ! Combined flux in latest and established models = 1.87425 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 26 model components to the UV plane model. ! The established model now contains 26 components and 1.87425 Jy ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.7979 mas, bmaj=2.532 mas, bpa=-14.59 degrees ! Estimated noise=1581.41 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.121505 Jy ! Component: 100 - total flux cleaned = 0.186768 Jy ! Component: 150 - total flux cleaned = 0.224817 Jy ! Component: 200 - total flux cleaned = 0.247725 Jy ! Total flux subtracted in 200 components = 0.247725 Jy ! Clean residual min=-0.014742 max=0.014887 Jy/beam ! Clean residual mean=0.000090 rms=0.005430 Jy/beam ! Combined flux in latest and established models = 2.12198 Jy ! Adding 43 model components to the UV plane model. ! The established model now contains 59 components and 2.12198 Jy ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=2.336 mas, bmaj=9.67 mas, bpa=-3.126 degrees ! Estimated noise=1943.28 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 0.0218067 Jy ! Component: 100 - total flux cleaned = 0.0315762 Jy ! Component: 150 - total flux cleaned = 0.0363682 Jy ! Component: 200 - total flux cleaned = 0.0370948 Jy ! Total flux subtracted in 200 components = 0.0370948 Jy ! Clean residual min=-0.014804 max=0.010553 Jy/beam ! Clean residual mean=-0.000811 rms=0.005827 Jy/beam ! Combined flux in latest and established models = 2.15907 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.000316732 Jy ! Component: 100 - total flux cleaned = 0.000611942 Jy ! Component: 150 - total flux cleaned = 0.000885193 Jy ! Component: 200 - total flux cleaned = 0.00120756 Jy ! Total flux subtracted in 200 components = 0.00120756 Jy ! Clean residual min=-0.014560 max=0.010419 Jy/beam ! Clean residual mean=-0.000847 rms=0.005770 Jy/beam ! Combined flux in latest and established models = 2.16028 Jy ! Adding 22 model components to the UV plane model. ! The established model now contains 76 components and 2.16028 Jy ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec uvtaper 0,0 ! No UV-taper is currently set. clean ! Inverting map and beam ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! Estimated beam: bmin=0.7979 mas, bmaj=2.532 mas, bpa=-14.59 degrees ! Estimated noise=1581.41 mJy/beam. ! clean: niter=200 gain=0.03 cutoff=0 ! Component: 050 - total flux cleaned = 8.93827e-05 Jy ! Component: 100 - total flux cleaned = 0.000313274 Jy ! Component: 150 - total flux cleaned = 0.000497502 Jy ! Component: 200 - total flux cleaned = 0.000922554 Jy ! Total flux subtracted in 200 components = 0.000922554 Jy ! Clean residual min=-0.012682 max=0.009094 Jy/beam ! Clean residual mean=-0.000434 rms=0.003411 Jy/beam ! Combined flux in latest and established models = 2.1612 Jy selfcal ! Performing phase self-cal ! Adding 20 model components to the UV plane model. ! The established model now contains 94 components and 2.1612 Jy ! ! Correcting IF 1. ! ! Correcting IF 2. ! ! Correcting IF 3. ! ! Correcting IF 4. ! ! Fit before self-cal, rms=0.225116Jy sigma=0.001269 ! Fit after self-cal, rms=0.224765Jy sigma=0.001267 wmodel J1427-4206_X_1999_03_08_pus_map.mod ! Writing 94 model components to file: J1427-4206_X_1999_03_08_pus_map.mod wobs J1427-4206_X_1999_03_08_pus_uvs.fits ! Writing UV FITS file: J1427-4206_X_1999_03_08_pus_uvs.fits wwins J1427-4206_X_1999_03_08_pus_map.win ! wwins: Wrote 2 windows to J1427-4206_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.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! clean: niter=1024 gain=0.01 cutoff=0 ! Component: 050 - total flux cleaned = -0.00425334 Jy ! Component: 100 - total flux cleaned = -0.00472617 Jy ! Component: 150 - total flux cleaned = -0.00576027 Jy ! Component: 200 - total flux cleaned = -0.00669828 Jy ! Component: 250 - total flux cleaned = -0.00780826 Jy ! Component: 300 - total flux cleaned = -0.0087169 Jy ! Component: 350 - total flux cleaned = -0.00978119 Jy ! Component: 400 - total flux cleaned = -0.0111987 Jy ! Component: 450 - total flux cleaned = -0.0125561 Jy ! Component: 500 - total flux cleaned = -0.0137699 Jy ! Component: 550 - total flux cleaned = -0.0147636 Jy ! Component: 600 - total flux cleaned = -0.0158144 Jy ! Component: 650 - total flux cleaned = -0.0164944 Jy ! Component: 700 - total flux cleaned = -0.0171563 Jy ! Component: 750 - total flux cleaned = -0.0179592 Jy ! Component: 800 - total flux cleaned = -0.0183508 Jy ! Component: 850 - total flux cleaned = -0.0191154 Jy ! Component: 900 - total flux cleaned = -0.0197111 Jy ! Component: 950 - total flux cleaned = -0.0202935 Jy ! Component: 1000 - total flux cleaned = -0.0209339 Jy ! Total flux subtracted in 1024 components = -0.0212841 Jy ! Clean residual min=-0.003853 max=0.003615 Jy/beam ! Clean residual mean=-0.000315 rms=0.001540 Jy/beam ! Combined flux in latest and established models = 2.13992 Jy keep ! Adding 63 model components to the UV plane model. ! The established model now contains 157 components and 2.13992 Jy ![Exited script file: muppet] wmap "J1427-4206_X_1999_03_08_pus\_map.fits" ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec ! restore: Substituting estimate of restoring beam from last 'invert'. ! Restoring with beam: 0.7979 x 2.532 at -14.59 degrees (North through East) ! Clean map min=-0.0085768 max=1.7375 Jy/beam ! Writing clean map to FITS file: J1427-4206_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 1.73754 invert ! Inverting map ! Your choice of large map pixels excluded 0.439% of the data. ! The x-axis pixel size should ideally be below 0.1989 milli-arcsec float image_rms image_rms = imstat(rms) print "The final image rms is",image_rms ! The final image rms is 0.0015388 device "J1427-4206_X_1999_03_08_pus\_map.eps/VPS" ! Attempting to open device: 'J1427-4206_X_1999_03_08_pus_map.eps/VPS' loglevs 100*3*image_rms/peak_flux ! The new contour levels are: ! -0.265685 0.265685 0.531371 1.06274 2.12548 4.25097 8.50193 17.0039 34.0077 68.0155 print "clev=",int(3*image_rms*10000)/10000 ! clev= 0.0046 print "peak=",int(1000*peak_flux)/1000 ! peak= 1.737 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.7979 x 2.532 at -14.59 degrees (North through East) ! Clean map min=-0.0085768 max=1.7375 Jy/beam mapcolor color, 0.9 ! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5. device "J1427-4206_X_1999_03_08_pus\_map.gif/vgif" ! Attempting to open device: 'J1427-4206_X_1999_03_08_pus_map.gif/vgif' mapplot cln quit ! Quitting program ! Log file J1427-4206_X_1999_03_08_pus_uvs.log closed on Thu Jan 3 18:36:45 2008