! Started logfile: 0108+388_X.log on Thu Jun 26 13:31:34 2008
observe 0108+388_X.edt
! Reading UV FITS file: 0108+388_X.edt
! AN table 1: 36 integrations on 190 of 190 possible baselines.
! AN table 2: 108 integrations on 190 of 190 possible baselines.
! Apparent sampling: 0.203947 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 0108+388
! 
! 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 1472 lines of history.
! 
! Reading 22320 visibilities.
![@gonzo6-0.1 0108+388_X]
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 = 50
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float win_mult
win_mult = 1.8
float solint
solint = 3.6/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.
select
! 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
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 430 telescope corrections were flagged in sub-array 1.
!  A total of 1544 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 430 telescope corrections were flagged in sub-array 1.
!  A total of 1544 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 430 telescope corrections were flagged in sub-array 1.
!  A total of 1549 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 470 telescope corrections were flagged in sub-array 1.
!  A total of 1549 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.829569Jy  sigma=13.282685
! Fit after  self-cal, rms=0.828739Jy  sigma=13.283783
! 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.4765 mas, bmaj=0.9186 mas, bpa=3.504 degrees
! Estimated noise=0.68978 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=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.141227 Jy
! Total flux subtracted in 50 components = 0.141227 Jy
! Clean residual min=-0.020724 max=0.077149 Jy/beam
! Clean residual mean=-0.000075 rms=0.006545 Jy/beam
! Combined flux in latest and established models = 0.141227 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 7 model components to the UV plane model.
! The established model now contains 7 components and 0.141227 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.339659Jy  sigma=4.484939
! Fit after  self-cal, rms=0.339574Jy  sigma=4.484423
! Inverting map 
! Added new window around map position (-5.1, -2.9).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0735739 Jy
! Total flux subtracted in 50 components = 0.0735739 Jy
! Clean residual min=-0.018256 max=0.044961 Jy/beam
! Clean residual mean=-0.000059 rms=0.005681 Jy/beam
! Combined flux in latest and established models = 0.214801 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 4 model components to the UV plane model.
! The established model now contains 9 components and 0.214801 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.321079Jy  sigma=4.062662
! Fit after  self-cal, rms=0.319369Jy  sigma=4.037379
! Inverting map 
! Added new window around map position (-0.9, -0.1).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0568408 Jy
! Total flux subtracted in 50 components = 0.0568408 Jy
! Clean residual min=-0.015758 max=0.043079 Jy/beam
! Clean residual mean=-0.000033 rms=0.005120 Jy/beam
! Combined flux in latest and established models = 0.271642 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 14 model components to the UV plane model.
! The established model now contains 20 components and 0.271642 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.305157Jy  sigma=3.671527
! Fit after  self-cal, rms=0.303796Jy  sigma=3.658428
! Inverting map 
! Added new window around map position (-1.6, -0.4).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0518823 Jy
! Total flux subtracted in 50 components = 0.0518823 Jy
! Clean residual min=-0.014462 max=0.033892 Jy/beam
! Clean residual mean=-0.000010 rms=0.004595 Jy/beam
! Combined flux in latest and established models = 0.323524 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 11 model components to the UV plane model.
! The established model now contains 28 components and 0.323524 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.291761Jy  sigma=3.334848
! Fit after  self-cal, rms=0.289777Jy  sigma=3.318198
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0432636 Jy
! Total flux subtracted in 50 components = 0.0432636 Jy
! Clean residual min=-0.013801 max=0.033487 Jy/beam
! Clean residual mean=0.000007 rms=0.004174 Jy/beam
! Combined flux in latest and established models = 0.366788 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 11 model components to the UV plane model.
! The established model now contains 31 components and 0.366788 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.280387Jy  sigma=3.072917
! Fit after  self-cal, rms=0.279099Jy  sigma=3.063881
! Inverting map 
! Added new window around map position (-4.1, -2.9).
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0390594 Jy
! Total flux subtracted in 50 components = 0.0390594 Jy
! Clean residual min=-0.012817 max=0.022867 Jy/beam
! Clean residual mean=0.000019 rms=0.003826 Jy/beam
! Combined flux in latest and established models = 0.405847 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 12 model components to the UV plane model.
! The established model now contains 35 components and 0.405847 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.270778Jy  sigma=2.859230
! Fit after  self-cal, rms=0.270101Jy  sigma=2.851713
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0327219 Jy
! Total flux subtracted in 50 components = 0.0327219 Jy
! Clean residual min=-0.012001 max=0.019543 Jy/beam
! Clean residual mean=0.000030 rms=0.003571 Jy/beam
! Combined flux in latest and established models = 0.438569 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 16 model components to the UV plane model.
! The established model now contains 38 components and 0.438569 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.263756Jy  sigma=2.692532
! Fit after  self-cal, rms=0.263396Jy  sigma=2.688142
! Inverting map 
! clean: niter=50  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0279523 Jy
! Total flux subtracted in 50 components = 0.0279523 Jy
! Clean residual min=-0.011953 max=0.017221 Jy/beam
! Clean residual mean=0.000039 rms=0.003376 Jy/beam
! Combined flux in latest and established models = 0.466522 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 17 model components to the UV plane model.
! The established model now contains 44 components and 0.466522 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.258514Jy  sigma=2.560592
! Fit after  self-cal, rms=0.258392Jy  sigma=2.558321
! 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.6
clean_niter = clean_niter * 2
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.581 mas, bmaj=1.107 mas, bpa=3.773 degrees
! Estimated noise=0.397034 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)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0517423 Jy
! Component: 100  -  total flux cleaned = 0.0891947 Jy
! Total flux subtracted in 100 components = 0.0891947 Jy
! Clean residual min=-0.015460 max=0.028909 Jy/beam
! Clean residual mean=0.000610 rms=0.006094 Jy/beam
! Combined flux in latest and established models = 0.555716 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 40 model components to the UV plane model.
! The established model now contains 66 components and 0.555716 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.246028Jy  sigma=2.241784
! Fit after  self-cal, rms=0.245907Jy  sigma=2.232460
! 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 21 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.88     FD        0.85     GC        0.88     HH        1.00*    
!   HN        1.10     KK        0.99     KP        0.94     LA        1.00     
!   MA        1.00*    MC        1.00*    MK        1.01     NL        0.86     
!   NY        0.86     ON        1.00*    OV        0.86     PT        0.88     
!   SC        1.07     TS        1.00     WF        0.95     WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.51     FD        0.93     GC        1.01     HH        1.00*    
!   HN        1.00*    KK        1.00*    KP        0.97     LA        0.86     
!   MA        0.68     MC        0.82     MK        1.00*    NL        1.02     
!   NY        0.51     ON        1.00*    OV        1.83     PT        0.89     
!   SC        1.01     TS        1.00*    WF        1.47     WZ        1.08     
! 
! 
! Correcting IF 2.
!  A total of 21 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.85     FD        0.95     GC        0.95     HH        1.00*    
!   HN        1.03     KK        1.01     KP        1.00     LA        1.01     
!   MA        1.00*    MC        1.00*    MK        1.06     NL        0.93     
!   NY        1.03     ON        1.00*    OV        1.01     PT        0.89     
!   SC        1.15     TS        0.99     WF        0.88     WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.43     FD        0.92     GC        0.78     HH        1.00*    
!   HN        1.00*    KK        1.00*    KP        0.92     LA        0.85     
!   MA        0.67     MC        0.86     MK        1.00*    NL        0.99     
!   NY        0.42     ON        1.00*    OV        1.90     PT        0.87     
!   SC        1.00     TS        1.00*    WF        1.73     WZ        1.02     
! 
! 
! Correcting IF 3.
!  A total of 21 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.85     FD        1.05     GC        0.97     HH        1.00*    
!   HN        1.15     KK        1.07     KP        0.96     LA        1.01     
!   MA        1.00*    MC        1.00*    MK        1.06     NL        0.79     
!   NY        1.12     ON        1.00*    OV        0.99     PT        0.90     
!   SC        1.48     TS        0.94     WF        0.66     WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.50     FD        0.98     GC        0.95     HH        1.00*    
!   HN        1.00*    KK        1.00*    KP        0.92     LA        0.88     
!   MA        0.82     MC        0.89     MK        1.00*    NL        1.02     
!   NY        0.43     ON        1.00*    OV        1.83     PT        0.76     
!   SC        1.06     TS        1.00*    WF        1.26     WZ        1.00     
! 
! 
! Correcting IF 4.
!  A total of 30 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.88     FD        1.01     GC        1.09     HH        1.00*    
!   HN        0.98     KK        0.96     KP        0.97     LA        0.97     
!   MA        1.00*    MC        1.00*    MK        1.03     NL        0.92     
!   NY        1.00*    ON        1.00*    OV        0.99     PT        0.92     
!   SC        1.22     TS        0.97     WF        0.80     WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.39     FD        1.03     GC        0.94     HH        1.00*    
!   HN        1.00*    KK        1.00*    KP        0.91     LA        0.89     
!   MA        0.83     MC        1.15     MK        1.00*    NL        1.06     
!   NY        0.57     ON        1.00*    OV        1.91     PT        0.86     
!   SC        1.18     TS        1.00*    WF        0.88     WZ        1.01     
! 
! 
! Fit before self-cal, rms=0.245907Jy  sigma=2.232460
! Fit after  self-cal, rms=0.167261Jy  sigma=1.656728
dynam = dynam - 0.5
doamp = false
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.5735 mas, bmaj=1.117 mas, bpa=4.084 degrees
! Estimated noise=0.347361 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)
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00351604 Jy
! Component: 100  -  total flux cleaned = 0.00711602 Jy
! Total flux subtracted in 100 components = 0.00711602 Jy
! Clean residual min=-0.005645 max=0.007637 Jy/beam
! Clean residual mean=0.000182 rms=0.001188 Jy/beam
! Combined flux in latest and established models = 0.562832 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 13 model components to the UV plane model.
! The established model now contains 78 components and 0.562832 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.165984Jy  sigma=1.633855
! Fit after  self-cal, rms=0.166147Jy  sigma=1.629243
! Inverting map 
dynam = dynam - 1.0
doamp = true
solint = solint * (6.0/3.6)
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=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00154753 Jy
! Component: 100  -  total flux cleaned = 0.00316888 Jy
! Total flux subtracted in 100 components = 0.00316888 Jy
! Clean residual min=-0.004945 max=0.006081 Jy/beam
! Clean residual mean=0.000177 rms=0.001092 Jy/beam
! Combined flux in latest and established models = 0.566001 Jy
! Performing amp+phase self-cal over 1 minute time intervals
! Adding 14 model components to the UV plane model.
! The established model now contains 87 components and 0.566001 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.165506Jy  sigma=1.619469
! Fit after  self-cal, rms=0.246591Jy  sigma=1.575781
! Inverting map and beam 
! Estimated beam: bmin=0.5695 mas, bmaj=1.119 mas, bpa=4.121 degrees
! Estimated noise=0.345107 mJy/beam.
doamp = false
solint = solint / (6.0/3.6)
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=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000672119 Jy
! Component: 100  -  total flux cleaned = 0.00183643 Jy
! Total flux subtracted in 100 components = 0.00183643 Jy
! Clean residual min=-0.004411 max=0.005469 Jy/beam
! Clean residual mean=0.000142 rms=0.000908 Jy/beam
! Combined flux in latest and established models = 0.567838 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 14 model components to the UV plane model.
! The established model now contains 93 components and 0.567838 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.246293Jy  sigma=1.568227
! Fit after  self-cal, rms=0.247266Jy  sigma=1.566244
! Inverting map 
! Added new window around map position (0.6, 0.1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00200231 Jy
! Component: 100  -  total flux cleaned = 0.00283542 Jy
! Total flux subtracted in 100 components = 0.00283542 Jy
! Clean residual min=-0.003965 max=0.004124 Jy/beam
! Clean residual mean=0.000135 rms=0.000858 Jy/beam
! Combined flux in latest and established models = 0.570673 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 13 model components to the UV plane model.
! The established model now contains 98 components and 0.570673 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.247025Jy  sigma=1.560387
! Fit after  self-cal, rms=0.247089Jy  sigma=1.560826
! Inverting map 
! Added new window around map position (-6.2, -2.4).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000863361 Jy
! Component: 100  -  total flux cleaned = 0.00140885 Jy
! Total flux subtracted in 100 components = 0.00140885 Jy
! Clean residual min=-0.003491 max=0.004029 Jy/beam
! Clean residual mean=0.000132 rms=0.000817 Jy/beam
! Combined flux in latest and established models = 0.572082 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 17 model components to the UV plane model.
! The established model now contains 106 components and 0.572082 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.246913Jy  sigma=1.556451
! Fit after  self-cal, rms=0.246920Jy  sigma=1.556671
! Inverting map 
! Added new window around map position (7.1, 18.7).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00175296 Jy
! Component: 100  -  total flux cleaned = 0.00207847 Jy
! Total flux subtracted in 100 components = 0.00207847 Jy
! Clean residual min=-0.003286 max=0.003321 Jy/beam
! Clean residual mean=0.000128 rms=0.000777 Jy/beam
! Combined flux in latest and established models = 0.57416 Jy
! Performing phase self-cal over 0.6 minute time intervals
! Adding 21 model components to the UV plane model.
! The established model now contains 112 components and 0.57416 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.246785Jy  sigma=1.553271
! Fit after  self-cal, rms=0.246790Jy  sigma=1.553308
! Inverting map 
print "***********************************"
! *********************************** 
print "*** Finished Amplitude Self-Cal ***"
! *** Finished Amplitude Self-Cal *** 
print "***********************************"
! *********************************** 
save 0108+388_X
! Writing UV FITS file: 0108+388_X.uvf
! Writing 112 model components to file: 0108+388_X.mod
! wwins: Wrote 8 windows to 0108+388_X.win
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.5695 x 1.119 at 4.121 degrees (North through East)
! Clean map  min=-0.0035397  max=0.14973 Jy/beam
! Writing clean map to FITS file: 0108+388_X.fits
! Writing difmap environment to: 0108+388_X.par
print "****************************"
! **************************** 
print "*** Now, make some plots ***"
! *** Now, make some plots *** 
print "****************************"
! **************************** 
![@gonzo_plot 0108+388_X]
device 0108+388_Xuv.ps/vps,1,2
! Attempting to open device: '0108+388_Xuv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device 0108+388_Xim.ps/vps
! Attempting to open device: '0108+388_Xim.ps/vps'
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 -1,2048,2
! The new contour levels are:
!  -1 1 2 4 8 16 32 64 128 256 512 1024 2048
cmul = imstat(rms)*3
! Inverting map and beam 
! Estimated beam: bmin=0.5695 mas, bmaj=1.119 mas, bpa=4.121 degrees
! Estimated noise=0.345107 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.5695 x 1.119 at 4.121 degrees (North through East)
! Clean map  min=-0.0035397  max=0.14973 Jy/beam
![Exited script file: gonzo_plot]
![Exited script file: gonzo6-0.1]
quit
! Quitting program
! Log file 0108+388_X.log closed on Thu Jun 26 13:31:50 2008