! Started logfile: 3C274_U.log on Mon Dec 16 16:57:49 1996
observe 3C274_U.edt
! Reading UV FITS file: 3C274_U.edt
! AN table 1: 54 integrations on 28 of 28 possible baselines.
! Apparent sampling: 0.902116 visibilities/baseline/integration-bin.
! Found source: 1228+126
! 
! Number of spectral-line channels: 1
! 
! IF    Freq (Hz)  Chan BW (Hz)  Full BW (Hz)  Sideband
! 01  1.49105e+10         2e+06         2e+06       USB
! 02  1.49205e+10         2e+06         2e+06       USB
! 03  1.49505e+10         2e+06         2e+06       USB
! 04  1.50105e+10         2e+06         2e+06       USB
! 05  1.51205e+10         2e+06         2e+06       USB
! 06  1.52005e+10         2e+06         2e+06       USB
! 07  1.52505e+10         2e+06         2e+06       USB
! 08  1.52705e+10         2e+06         2e+06       USB
! 
! Polarization(s): RR
! 
! Read 2592 lines of history.
! 
! Reading 10912 visibilities.
! Selecting polarization: RR,  channels: 1..1
print "***********************************"
! *********************************** 
print "*** Using gonzo instruction set ***"
! *** Using gonzo instruction set *** 
print "***********************************"
! *********************************** 
float field_size
field_size = 1024
float field_cell
field_cell = 0.1
integer clean_niter
clean_niter = -100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float win_mult
win_mult = 1.5
float solint
solint = 4.0/6.0
logical doflt
doflt = true
logical doamp
float new_peak
float flux_cutoff
mapsize field_size,field_cell
! Map grid = 1024x1024 pixels with 0.100x0.100 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 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 67 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 5.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 6.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 7.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 8.
!  A total of 64 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.673320Jy  sigma=3.917254
! Fit after  self-cal, rms=0.672787Jy  sigma=3.916694
! 
! Reinstating the default IF.
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
uvweight 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.
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.5399 mas, bmaj=1.13 mas, bpa=30.44 degrees
! Estimated noise=2.29592 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  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.903812 Jy
! Component: 100  -  total flux cleaned = 1.22658 Jy
! Total flux subtracted in 100 components = 1.22658 Jy
! Clean residual min=-0.067299 max=0.196088 Jy/beam
! Clean residual mean=-0.000166 rms=0.020281 Jy/beam
! Combined flux in latest and established models = 1.22658 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 15 model components to the UV plane model.
! The established model now contains 15 components and 1.22658 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.532106Jy  sigma=3.088201
! Fit after  self-cal, rms=0.530862Jy  sigma=3.080623
! 
! Reinstating the default IF.
! Inverting map 
! Added new window around map position (-0.8, 0.3).
! clean: niter=-100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.205442 Jy
! Component: 100  -  total flux cleaned = 0.33546 Jy
! Total flux subtracted in 100 components = 0.33546 Jy
! Clean residual min=-0.045667 max=0.094772 Jy/beam
! Clean residual mean=-0.000111 rms=0.012528 Jy/beam
! Combined flux in latest and established models = 1.56204 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 29 model components to the UV plane model.
! The established model now contains 40 components and 1.56204 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.369776Jy  sigma=2.130052
! Fit after  self-cal, rms=0.366571Jy  sigma=2.110947
! 
! Reinstating the default IF.
! Inverting map 
! clean: niter=-100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0982822 Jy
! Component: 100  -  total flux cleaned = 0.167771 Jy
! Total flux subtracted in 100 components = 0.167771 Jy
! Clean residual min=-0.035195 max=0.057109 Jy/beam
! Clean residual mean=-0.000084 rms=0.009372 Jy/beam
! Combined flux in latest and established models = 1.72981 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 43 model components to the UV plane model.
! The established model now contains 76 components and 1.72981 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.311375Jy  sigma=1.783359
! Fit after  self-cal, rms=0.309911Jy  sigma=1.774872
! 
! Reinstating the default IF.
! Inverting map 
! Added new window around map position (-1.7, 0.4).
! clean: niter=-100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0627675 Jy
! Component: 100  -  total flux cleaned = 0.10877 Jy
! Total flux subtracted in 100 components = 0.10877 Jy
! Clean residual min=-0.030737 max=0.043167 Jy/beam
! Clean residual mean=-0.000063 rms=0.007694 Jy/beam
! Combined flux in latest and established models = 1.83858 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 34 model components to the UV plane model.
! The established model now contains 105 components and 1.83858 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.284033Jy  sigma=1.620824
! Fit after  self-cal, rms=0.282900Jy  sigma=1.614289
! 
! Reinstating the default IF.
! Inverting map 
print "****************************************"
! **************************************** 
print "*** Finished Uniform Weighting CLEAN ***"
! *** Finished Uniform Weighting CLEAN *** 
print "****************************************"
! **************************************** 
uvweight 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.3
clean_niter = clean_niter * 2
dynam = dynam - 1.0
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.6903 mas, bmaj=1.205 mas, bpa=26.96 degrees
! Estimated noise=1.2394 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (-2.6, 0.4).
! clean: niter=-200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0759865 Jy
! Component: 100  -  total flux cleaned = 0.127037 Jy
! Component: 150  -  total flux cleaned = 0.163957 Jy
! Component: 200  -  total flux cleaned = 0.190889 Jy
! Total flux subtracted in 200 components = 0.190889 Jy
! Clean residual min=-0.032148 max=0.039344 Jy/beam
! Clean residual mean=0.000001 rms=0.007914 Jy/beam
! Combined flux in latest and established models = 2.02947 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 68 model components to the UV plane model.
! The established model now contains 159 components and 2.02947 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.256851Jy  sigma=1.458983
! Fit after  self-cal, rms=0.254876Jy  sigma=1.447316
! 
! Reinstating the default IF.
! Inverting map 
! Added new window around map position (-5.7, 2.9).
! clean: niter=-200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0414702 Jy
! Component: 100  -  total flux cleaned = 0.06815 Jy
! Component: 150  -  total flux cleaned = 0.0883389 Jy
! Component: 200  -  total flux cleaned = 0.104255 Jy
! Total flux subtracted in 200 components = 0.104255 Jy
! Clean residual min=-0.029198 max=0.030626 Jy/beam
! Clean residual mean=0.000017 rms=0.006187 Jy/beam
! Combined flux in latest and established models = 2.13373 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 49 model components to the UV plane model.
! The established model now contains 202 components and 2.13373 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.246671Jy  sigma=1.396899
! Fit after  self-cal, rms=0.245839Jy  sigma=1.392007
! 
! Reinstating the default IF.
! Inverting map 
print "****************************************"
! **************************************** 
print "*** Finished Natural Weighting CLEAN ***"
! *** Finished Natural Weighting CLEAN *** 
print "****************************************"
! **************************************** 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.02     FD        0.97     HN        1.00     KP        1.00     
!   LA        1.00     OV        1.00     PT        0.98     SC        0.99     
! 
! 
! Correcting IF 2.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.01     FD        0.97     HN        1.06     KP        0.99     
!   LA        0.98     OV        1.00     PT        1.00     SC        0.99     
! 
! 
! Correcting IF 3.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.04     FD        0.97     HN        1.01     KP        0.97     
!   LA        1.00     OV        1.01     PT        0.99     SC        1.01     
! 
! 
! Correcting IF 4.
!  A total of 15 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.01     FD        0.99     HN        1.03     KP        1.00     
!   LA        1.01     OV        1.02     PT        0.98     SC        1.00     
! 
! 
! Correcting IF 5.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        0.96     HN        1.03     KP        0.98     
!   LA        0.98     OV        0.99     PT        0.97     SC        1.00     
! 
! 
! Correcting IF 6.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        0.99     HN        1.02     KP        0.98     
!   LA        0.98     OV        0.99     PT        0.99     SC        1.00     
! 
! 
! Correcting IF 7.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.02     FD        0.99     HN        1.03     KP        0.99     
!   LA        1.00     OV        1.02     PT        1.00     SC        1.00     
! 
! 
! Correcting IF 8.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.02     FD        0.99     HN        1.03     KP        1.01     
!   LA        1.00     OV        1.04     PT        0.99     SC        0.99     
! 
! 
! Fit before self-cal, rms=0.245839Jy  sigma=1.392007
! Fit after  self-cal, rms=0.241986Jy  sigma=1.374714
! 
! Reinstating the default IF.
clean_niter = -clean_niter
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.6916 mas, bmaj=1.211 mas, bpa=26.76 degrees
! Estimated noise=1.24243 mJy/beam.
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (-1.8, -7.1).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0305236 Jy
! Component: 100  -  total flux cleaned = 0.047555 Jy
! Component: 150  -  total flux cleaned = 0.0601041 Jy
! Component: 200  -  total flux cleaned = 0.0700442 Jy
! Total flux subtracted in 200 components = 0.0700442 Jy
! Clean residual min=-0.024148 max=0.020018 Jy/beam
! Clean residual mean=-0.000009 rms=0.005140 Jy/beam
! Combined flux in latest and established models = 2.20377 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 41 model components to the UV plane model.
! The established model now contains 239 components and 2.20377 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.237827Jy  sigma=1.348845
! Fit after  self-cal, rms=0.237628Jy  sigma=1.347737
! 
! Reinstating the default IF.
! Inverting map 
dynam = dynam - 1.0
doamp = true
solint = solint * 1.5
flux_cutoff = imstat(rms) * dynam
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00956777 Jy
! Component: 100  -  total flux cleaned = 0.0158408 Jy
! Component: 150  -  total flux cleaned = 0.0198808 Jy
! Component: 200  -  total flux cleaned = 0.022735 Jy
! Total flux subtracted in 200 components = 0.022735 Jy
! Clean residual min=-0.023421 max=0.020057 Jy/beam
! Clean residual mean=-0.000009 rms=0.004999 Jy/beam
! Combined flux in latest and established models = 2.22651 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 42 model components to the UV plane model.
! The established model now contains 255 components and 2.22651 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.237042Jy  sigma=1.344024
! Fit after  self-cal, rms=0.228531Jy  sigma=1.285895
! 
! Reinstating the default IF.
! Inverting map and beam 
! Estimated beam: bmin=0.6911 mas, bmaj=1.211 mas, bpa=26.62 degrees
! Estimated noise=1.24199 mJy/beam.
! Added new window around map position (-10.8, 3.3).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0193866 Jy
! Component: 100  -  total flux cleaned = 0.0311017 Jy
! Component: 150  -  total flux cleaned = 0.0397654 Jy
! Component: 200  -  total flux cleaned = 0.0458794 Jy
! Total flux subtracted in 200 components = 0.0458794 Jy
! Clean residual min=-0.011240 max=0.014612 Jy/beam
! Clean residual mean=0.000000 rms=0.003232 Jy/beam
! Combined flux in latest and established models = 2.27239 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 37 model components to the UV plane model.
! The established model now contains 292 components and 2.27239 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.226670Jy  sigma=1.274078
! Fit after  self-cal, rms=0.225911Jy  sigma=1.270223
! 
! Reinstating the default IF.
! Inverting map and beam 
! Estimated beam: bmin=0.6913 mas, bmaj=1.211 mas, bpa=26.58 degrees
! Estimated noise=1.24132 mJy/beam.
! Added new window around map position (-4.2, 0.8).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0124672 Jy
! Component: 100  -  total flux cleaned = 0.0196916 Jy
! Component: 150  -  total flux cleaned = 0.0220574 Jy
! Component: 200  -  total flux cleaned = 0.0237124 Jy
! Total flux subtracted in 200 components = 0.0237124 Jy
! Clean residual min=-0.010682 max=0.012644 Jy/beam
! Clean residual mean=0.000005 rms=0.002922 Jy/beam
! Combined flux in latest and established models = 2.2961 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 24 model components to the UV plane model.
! The established model now contains 307 components and 2.2961 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.225035Jy  sigma=1.264578
! Fit after  self-cal, rms=0.224497Jy  sigma=1.262114
! 
! Reinstating the default IF.
! Inverting map and beam 
! Estimated beam: bmin=0.6913 mas, bmaj=1.211 mas, bpa=26.64 degrees
! Estimated noise=1.24092 mJy/beam.
! Added new window around map position (14.9, 13.4).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0117185 Jy
! Component: 100  -  total flux cleaned = 0.0188221 Jy
! Component: 150  -  total flux cleaned = 0.024726 Jy
! Component: 200  -  total flux cleaned = 0.0280667 Jy
! Total flux subtracted in 200 components = 0.0280667 Jy
! Clean residual min=-0.010210 max=0.009917 Jy/beam
! Clean residual mean=0.000009 rms=0.002629 Jy/beam
! Combined flux in latest and established models = 2.32416 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 45 model components to the UV plane model.
! The established model now contains 336 components and 2.32416 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.223755Jy  sigma=1.257322
! Fit after  self-cal, rms=0.223472Jy  sigma=1.255432
! 
! Reinstating the default IF.
! Inverting map and beam 
! Estimated beam: bmin=0.6914 mas, bmaj=1.211 mas, bpa=26.64 degrees
! Estimated noise=1.24118 mJy/beam.
doamp = false
solint = solint / 1.5
flux_cutoff = imstat(rms) * dynam
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal doamp, doflt, solint
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0051244 Jy
! Component: 100  -  total flux cleaned = 0.00785594 Jy
! Component: 150  -  total flux cleaned = 0.00981569 Jy
! Component: 200  -  total flux cleaned = 0.0108191 Jy
! Total flux subtracted in 200 components = 0.0108191 Jy
! Clean residual min=-0.009951 max=0.009834 Jy/beam
! Clean residual mean=0.000007 rms=0.002493 Jy/beam
! Combined flux in latest and established models = 2.33498 Jy
! Performing phase self-cal over 0.666667 minute time intervals
! Adding 45 model components to the UV plane model.
! The established model now contains 355 components and 2.33498 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.223179Jy  sigma=1.253576
! Fit after  self-cal, rms=0.223164Jy  sigma=1.253491
! 
! Reinstating the default IF.
! Inverting map 
print "***********************************"
! *********************************** 
print "*** Finished Amplitude Self-Cal ***"
! *** Finished Amplitude Self-Cal *** 
print "***********************************"
! *********************************** 
save 3C274_U
! Writing UV FITS file: 3C274_U.uvf
! Writing 355 model components to file: 3C274_U.mod
! wwins: Wrote 9 windows to 3C274_U.win
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6914 x 1.211 at 26.64 degrees (North through East)
! Clean map  min=-0.0097911  max=1.2034 Jy/beam
! Writing clean map to FITS file: 3C274_U.fits
! Writing difmap environment to: 3C274_U.par
print "****************************"
! **************************** 
print "*** Now, make some plots ***"
! *** Now, make some plots *** 
print "****************************"
! **************************** 
device 3C274_Uuv.ps/vps,1,2
! Attempting to open device: '3C274_Uuv.ps/vps'
radplot
! Using default options string "m1"
uvplot
device 3C274_Uim.ps/vps
! Attempting to open device: '3C274_Uim.ps/vps'
cmul = 1.0
uvweight 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.
loglevs (imstat(rms)*3),20000,2
! Inverting map and beam 
! Estimated beam: bmin=0.6914 mas, bmaj=1.211 mas, bpa=26.64 degrees
! Estimated noise=1.24118 mJy/beam.
! The new contour levels are:
!  -0.00748212 0.00748212 0.0149642 0.0299285 0.0598569 0.119714 0.239428 0.478855 0.957711 1.91542 3.83084 7.66169 15.3234 30.6467 61.2935 122.587 245.174 490.348 980.696 1961.39 3922.78 7845.57 15691.1
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.6914 x 1.211 at 26.64 degrees (North through East)
! Clean map  min=-0.0097911  max=1.2034 Jy/beam
quit
! Quitting program
! Log file 3C274_U.log closed on Mon Dec 16 17:03:31 1996