! Started logfile: J1230+1223_S_2001_04_09_pus_uvs.log on Fri Oct  5 17:16:00 2007
obs J1230+1223_S_2001_04_09_pus_uva.fits
! Reading UV FITS file: J1230+1223_S_2001_04_09_pus_uva.fits
! AN table 1: 20 integrations on 190 of 190 possible baselines.
! AN table 2: 156 integrations on 190 of 190 possible baselines.
! AN table 3: 156 integrations on 190 of 190 possible baselines.
! Apparent sampling: 0.13142 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J1230+1223
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  2.22997e+09        8e+06          1        8e+06
!  02        2  2.25997e+09        8e+06          1        8e+06
!  03        3  2.34997e+09        8e+06          1        8e+06
!  04        4  2.36997e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1949 lines of history.
! 
! Reading 33160 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 = 1024 
print "field_cell =",field_cell
! field_cell = 0.5 
print "taper_size =",taper_size
! taper_size = 20 
![@muppet J1230+1223_S_2001_04_09_pus]
integer clean_niter
clean_niter = 100
float clean_gain
clean_gain = 0.03
float dynam
dynam = 6.0
float soltime1
soltime1 = 120
float soltime2
soltime2 = 20
float thresh
thresh = 0.5
float win_mult
win_mult = 1.8
float time_av
time_av = 16
float old_peak
float new_peak
float flux_cutoff
mapsize field_size, field_cell
! Map grid = 1024x1024 pixels with 0.500x0.500 milli-arcsec cellsize.
startmod
! Applying default point source starting model.
! Performing phase self-cal
! Adding 1 model components to the UV plane model.
! The established model now contains 1 components and 1 Jy
! 
! Correcting IF 1.
!  A total of 370 telescope corrections were flagged in sub-array 1.
!  A total of 2493 telescope corrections were flagged in sub-array 2.
!  A total of 1926 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 2.
!  A total of 370 telescope corrections were flagged in sub-array 1.
!  A total of 2493 telescope corrections were flagged in sub-array 2.
!  A total of 1926 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 3.
!  A total of 370 telescope corrections were flagged in sub-array 1.
!  A total of 2493 telescope corrections were flagged in sub-array 2.
!  A total of 1926 telescope corrections were flagged in sub-array 3.
! 
! Correcting IF 4.
!  A total of 370 telescope corrections were flagged in sub-array 1.
!  A total of 2493 telescope corrections were flagged in sub-array 2.
!  A total of 1926 telescope corrections were flagged in sub-array 3.
! 
! Fit before self-cal, rms=1.020209Jy  sigma=0.006678
! Fit after  self-cal, rms=0.714491Jy  sigma=0.004315
! clrmod: Cleared the established, tentative and continuum models.
! Redundant starting model cleared.
uvw 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=1.815 mas, bmaj=3.391 mas, bpa=8.184 degrees
! Estimated noise=1086.31 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.617397 Jy
! Component: 100  -  total flux cleaned = 0.891952 Jy
! Total flux subtracted in 100 components = 0.891952 Jy
! Clean residual min=-0.040902 max=0.152508 Jy/beam
! Clean residual mean=-0.000153 rms=0.010989 Jy/beam
! Combined flux in latest and established models = 0.891952 Jy
! Performing phase self-cal
! Adding 10 model components to the UV plane model.
! The established model now contains 10 components and 0.891952 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.609358Jy  sigma=0.003935
! Fit after  self-cal, rms=0.606516Jy  sigma=0.003928
! Inverting map 
! Added new window around map position (-2.5, 0.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.165398 Jy
! Component: 100  -  total flux cleaned = 0.269221 Jy
! Total flux subtracted in 100 components = 0.269221 Jy
! Clean residual min=-0.029093 max=0.091790 Jy/beam
! Clean residual mean=-0.000071 rms=0.007871 Jy/beam
! Combined flux in latest and established models = 1.16117 Jy
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 24 components and 1.16117 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.495871Jy  sigma=0.003118
! Fit after  self-cal, rms=0.494031Jy  sigma=0.003109
! Inverting map 
! Added new window around map position (-6, 1).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0938757 Jy
! Component: 100  -  total flux cleaned = 0.15813 Jy
! Total flux subtracted in 100 components = 0.15813 Jy
! Clean residual min=-0.025024 max=0.068768 Jy/beam
! Clean residual mean=-0.000021 rms=0.006331 Jy/beam
! Combined flux in latest and established models = 1.3193 Jy
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 36 components and 1.3193 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.440308Jy  sigma=0.002710
! Fit after  self-cal, rms=0.438924Jy  sigma=0.002702
! Inverting map 
! Added new window around map position (3, -0.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0635573 Jy
! Component: 100  -  total flux cleaned = 0.110362 Jy
! Total flux subtracted in 100 components = 0.110362 Jy
! Clean residual min=-0.021854 max=0.046689 Jy/beam
! Clean residual mean=0.000015 rms=0.005418 Jy/beam
! Combined flux in latest and established models = 1.42967 Jy
! Performing phase self-cal
! Adding 29 model components to the UV plane model.
! The established model now contains 58 components and 1.42967 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.407322Jy  sigma=0.002464
! Fit after  self-cal, rms=0.407195Jy  sigma=0.002463
! Inverting map 
! Added new window around map position (-9, 2.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0488958 Jy
! Component: 100  -  total flux cleaned = 0.0847598 Jy
! Total flux subtracted in 100 components = 0.0847598 Jy
! Clean residual min=-0.021139 max=0.034034 Jy/beam
! Clean residual mean=0.000042 rms=0.004881 Jy/beam
! Combined flux in latest and established models = 1.51443 Jy
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 81 components and 1.51443 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.386033Jy  sigma=0.002300
! Fit after  self-cal, rms=0.385610Jy  sigma=0.002297
! Inverting map 
! Added new window around map position (-12, 4.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0375679 Jy
! Component: 100  -  total flux cleaned = 0.0667784 Jy
! Total flux subtracted in 100 components = 0.0667784 Jy
! Clean residual min=-0.019760 max=0.030459 Jy/beam
! Clean residual mean=0.000063 rms=0.004503 Jy/beam
! Combined flux in latest and established models = 1.5812 Jy
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 100 components and 1.5812 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.370769Jy  sigma=0.002180
! Fit after  self-cal, rms=0.370544Jy  sigma=0.002179
! Inverting map 
! Added new window around map position (-22.5, 10.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0330228 Jy
! Component: 100  -  total flux cleaned = 0.0585291 Jy
! Total flux subtracted in 100 components = 0.0585291 Jy
! Clean residual min=-0.018094 max=0.027638 Jy/beam
! Clean residual mean=0.000076 rms=0.004188 Jy/beam
! Combined flux in latest and established models = 1.63973 Jy
! Performing phase self-cal
! Adding 30 model components to the UV plane model.
! The established model now contains 115 components and 1.63973 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.358998Jy  sigma=0.002087
! Fit after  self-cal, rms=0.358622Jy  sigma=0.002085
! Inverting map 
! Added new window around map position (-19, 6.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0309539 Jy
! Component: 100  -  total flux cleaned = 0.0544116 Jy
! Total flux subtracted in 100 components = 0.0544116 Jy
! Clean residual min=-0.017690 max=0.020113 Jy/beam
! Clean residual mean=0.000086 rms=0.003884 Jy/beam
! Combined flux in latest and established models = 1.69414 Jy
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 139 components and 1.69414 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.347561Jy  sigma=0.001997
! Fit after  self-cal, rms=0.347240Jy  sigma=0.001996
! Inverting map 
print "************** FINISHED UNIFORM WEIGHTING CLEAN **************"
! ************** FINISHED UNIFORM WEIGHTING CLEAN ************** 
! 
uvw 0,-1
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
win_mult = win_mult * 1.6
clean_niter = clean_niter * 2
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.001 mas, bmaj=3.951 mas, bpa=16.58 degrees
! Estimated noise=646.746 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.0966525 Jy
! Component: 100  -  total flux cleaned = 0.168037 Jy
! Component: 150  -  total flux cleaned = 0.222485 Jy
! Component: 200  -  total flux cleaned = 0.264635 Jy
! Total flux subtracted in 200 components = 0.264635 Jy
! Clean residual min=-0.024699 max=0.041150 Jy/beam
! Clean residual mean=0.002222 rms=0.010704 Jy/beam
! Combined flux in latest and established models = 1.95878 Jy
! Performing phase self-cal
! Adding 109 model components to the UV plane model.
! The established model now contains 212 components and 1.95878 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.311909Jy  sigma=0.001711
! Fit after  self-cal, rms=0.310042Jy  sigma=0.001698
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 30 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GC        0.00*    GG        0.00*    
!   HH        0.00*    HN        0.00*    KK        0.00*    KP        0.00*    
!   LA        0.00*    MA        0.00*    MC        0.00*    MK        0.00*    
!   NL        0.00*    NY        0.00*    OV        0.00*    PT        0.00*    
!   SC        0.00*    TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 111 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        1.03     GC        1.06     GG        1.00*    
!   HH        1.06     HN        0.93     KK        1.04     KP        1.01     
!   LA        0.96     MA        0.99     MC        0.86     MK        1.02     
!   NL        1.01     NY        1.07     OV        1.03     PT        0.96     
!   SC        1.06     TS        1.14     WF        0.95     WZ        0.97     
! 
!  A total of 42 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.04     FD        0.91     GC        1.00*    GG        1.00*    
!   HH        1.04     HN        1.01     KK        1.00*    KP        1.02     
!   LA        0.90     MA        1.07     MC        1.04     MK        1.00*    
!   NL        1.04     NY        0.98     OV        1.03     PT        0.92     
!   SC        1.01     TS        1.13     WF        0.84     WZ        0.98     
! 
! 
! Correcting IF 2.
!  A total of 30 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GC        0.00*    GG        0.00*    
!   HH        0.00*    HN        0.00*    KK        0.00*    KP        0.00*    
!   LA        0.00*    MA        0.00*    MC        0.00*    MK        0.00*    
!   NL        0.00*    NY        0.00*    OV        0.00*    PT        0.00*    
!   SC        0.00*    TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 111 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.98     FD        1.01     GC        1.04     GG        1.00*    
!   HH        1.16     HN        0.96     KK        1.05     KP        0.90     
!   LA        0.98     MA        0.82     MC        0.94     MK        1.05     
!   NL        1.02     NY        1.11     OV        0.98     PT        0.99     
!   SC        1.06     TS        1.18     WF        0.93     WZ        0.95     
! 
!  A total of 42 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.02     FD        0.91     GC        1.00*    GG        1.00*    
!   HH        0.99     HN        0.98     KK        1.00*    KP        1.00     
!   LA        0.93     MA        0.91     MC        1.17     MK        1.00*    
!   NL        1.03     NY        1.01     OV        1.00     PT        0.96     
!   SC        1.03     TS        0.99     WF        0.91     WZ        0.96     
! 
! 
! Correcting IF 3.
!  A total of 30 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GC        0.00*    GG        0.00*    
!   HH        0.00*    HN        0.00*    KK        0.00*    KP        0.00*    
!   LA        0.00*    MA        0.00*    MC        0.00*    MK        0.00*    
!   NL        0.00*    NY        0.00*    OV        0.00*    PT        0.00*    
!   SC        0.00*    TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 111 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.04     FD        1.03     GC        1.00     GG        1.00*    
!   HH        0.99     HN        0.98     KK        1.06     KP        1.01     
!   LA        0.95     MA        0.90     MC        0.88     MK        0.99     
!   NL        1.01     NY        1.05     OV        1.00     PT        0.98     
!   SC        1.02     TS        1.19     WF        0.92     WZ        0.96     
! 
!  A total of 42 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.07     FD        0.91     GC        1.00*    GG        1.00*    
!   HH        1.00     HN        1.02     KK        1.00*    KP        1.03     
!   LA        0.89     MA        1.00     MC        0.96     MK        1.00*    
!   NL        1.03     NY        0.92     OV        1.00     PT        0.97     
!   SC        0.98     TS        1.16     WF        0.91     WZ        0.96     
! 
! 
! Correcting IF 4.
!  A total of 30 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GC        0.00*    GG        0.00*    
!   HH        0.00*    HN        0.00*    KK        0.00*    KP        0.00*    
!   LA        0.00*    MA        0.00*    MC        0.00*    MK        0.00*    
!   NL        0.00*    NY        0.00*    OV        0.00*    PT        0.00*    
!   SC        0.00*    TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 111 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.99     FD        1.02     GC        1.02     GG        1.00*    
!   HH        0.87     HN        1.01     KK        1.02     KP        1.00     
!   LA        0.98     MA        1.01     MC        1.81     MK        0.97     
!   NL        1.01     NY        1.12     OV        0.99     PT        1.09     
!   SC        1.00     TS        1.10     WF        0.95     WZ        1.23     
! 
!  A total of 42 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.02     FD        0.92     GC        1.00*    GG        1.00*    
!   HH        1.09     HN        1.01     KK        1.00*    KP        1.02     
!   LA        0.92     MA        0.87     MC        0.77     MK        1.00*    
!   NL        1.03     NY        1.14     OV        1.00     PT        0.95     
!   SC        0.99     TS        0.53     WF        0.91     WZ        1.42     
! 
! 
! Fit before self-cal, rms=0.310042Jy  sigma=0.001698
! Fit after  self-cal, rms=0.281918Jy  sigma=0.001580
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.004 mas, bmaj=4.141 mas, bpa=15.87 degrees
! Estimated noise=638.198 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.0212413 Jy
! Component: 100  -  total flux cleaned = 0.0307277 Jy
! Component: 150  -  total flux cleaned = 0.0333649 Jy
! Component: 200  -  total flux cleaned = 0.0355521 Jy
! Total flux subtracted in 200 components = 0.0355521 Jy
! Clean residual min=-0.015304 max=0.036941 Jy/beam
! Clean residual mean=0.001185 rms=0.007677 Jy/beam
! Combined flux in latest and established models = 1.99433 Jy
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 234 components and 1.99433 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.280192Jy  sigma=0.001566
! Fit after  self-cal, rms=0.280517Jy  sigma=0.001563
! Inverting map 
selfcal true, true, soltime1
! Performing amp+phase self-cal over 120 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.280517Jy  sigma=0.001563
! Fit after  self-cal, rms=0.284118Jy  sigma=0.001550
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.004 mas, bmaj=4.141 mas, bpa=15.35 degrees
! Estimated noise=636.354 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0130412 Jy
! Component: 100  -  total flux cleaned = 0.0177523 Jy
! Component: 150  -  total flux cleaned = 0.0203546 Jy
! Component: 200  -  total flux cleaned = 0.0222927 Jy
! Total flux subtracted in 200 components = 0.0222927 Jy
! Clean residual min=-0.016598 max=0.034902 Jy/beam
! Clean residual mean=0.000927 rms=0.007460 Jy/beam
! Combined flux in latest and established models = 2.01662 Jy
selfcal
! Performing phase self-cal
! Adding 27 model components to the UV plane model.
! The established model now contains 242 components and 2.01662 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.283334Jy  sigma=0.001544
! Fit after  self-cal, rms=0.284972Jy  sigma=0.001543
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! Added new window around map position (-50.5, 20).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0380465 Jy
! Component: 100  -  total flux cleaned = 0.0606965 Jy
! Component: 150  -  total flux cleaned = 0.0751628 Jy
! Component: 200  -  total flux cleaned = 0.085176 Jy
! Total flux subtracted in 200 components = 0.085176 Jy
! Clean residual min=-0.013947 max=0.026534 Jy/beam
! Clean residual mean=0.000631 rms=0.005528 Jy/beam
! Combined flux in latest and established models = 2.1018 Jy
! Performing phase self-cal
! Adding 60 model components to the UV plane model.
! The established model now contains 300 components and 2.1018 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.280029Jy  sigma=0.001499
! Fit after  self-cal, rms=0.279750Jy  sigma=0.001497
! Inverting map 
! Added new window around map position (-15.5, -40.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0265618 Jy
! Component: 100  -  total flux cleaned = 0.0418581 Jy
! Component: 150  -  total flux cleaned = 0.0526096 Jy
! Component: 200  -  total flux cleaned = 0.0613736 Jy
! Total flux subtracted in 200 components = 0.0613736 Jy
! Clean residual min=-0.011289 max=0.017429 Jy/beam
! Clean residual mean=0.000424 rms=0.004417 Jy/beam
! Combined flux in latest and established models = 2.16317 Jy
! Performing phase self-cal
! Adding 48 model components to the UV plane model.
! The established model now contains 343 components and 2.16317 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.277322Jy  sigma=0.001473
! Fit after  self-cal, rms=0.277136Jy  sigma=0.001472
! Inverting map 
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.277136Jy  sigma=0.001472
! Fit after  self-cal, rms=0.307144Jy  sigma=0.001467
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.006 mas, bmaj=4.145 mas, bpa=15.26 degrees
! Estimated noise=634.143 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0105736 Jy
! Component: 100  -  total flux cleaned = 0.0188896 Jy
! Component: 150  -  total flux cleaned = 0.0258582 Jy
! Component: 200  -  total flux cleaned = 0.0315512 Jy
! Total flux subtracted in 200 components = 0.0315512 Jy
! Clean residual min=-0.010311 max=0.015974 Jy/beam
! Clean residual mean=0.000484 rms=0.003899 Jy/beam
! Combined flux in latest and established models = 2.19472 Jy
selfcal
! Performing phase self-cal
! Adding 42 model components to the UV plane model.
! The established model now contains 364 components and 2.19472 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.306593Jy  sigma=0.001462
! Fit after  self-cal, rms=0.307506Jy  sigma=0.001461
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0050988 Jy
! Component: 100  -  total flux cleaned = 0.00838676 Jy
! Component: 150  -  total flux cleaned = 0.0115165 Jy
! Component: 200  -  total flux cleaned = 0.0138819 Jy
! Total flux subtracted in 200 components = 0.0138819 Jy
! Clean residual min=-0.009926 max=0.015218 Jy/beam
! Clean residual mean=0.000438 rms=0.003772 Jy/beam
! Combined flux in latest and established models = 2.20861 Jy
! Performing phase self-cal
! Adding 53 model components to the UV plane model.
! The established model now contains 378 components and 2.20861 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.307309Jy  sigma=0.001460
! Fit after  self-cal, rms=0.307300Jy  sigma=0.001460
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
win_mult = win_mult * 1.2
selfcal true, true, soltime2
! Performing amp+phase self-cal over 20 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.307300Jy  sigma=0.001460
! Fit after  self-cal, rms=0.348799Jy  sigma=0.001458
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.494 mas, bmaj=10.09 mas, bpa=-7.273 degrees
! Estimated noise=826.11 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0158458 Jy
! Component: 100  -  total flux cleaned = 0.0261106 Jy
! Component: 150  -  total flux cleaned = 0.0335075 Jy
! Component: 200  -  total flux cleaned = 0.036893 Jy
! Total flux subtracted in 200 components = 0.036893 Jy
! Clean residual min=-0.017025 max=0.029207 Jy/beam
! Clean residual mean=0.000865 rms=0.007732 Jy/beam
! Combined flux in latest and established models = 2.2455 Jy
selfcal
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 390 components and 2.2455 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.348767Jy  sigma=0.001458
! Fit after  self-cal, rms=0.349037Jy  sigma=0.001458
flux_cutoff = imstat(rms) * dynam
! Inverting map 
repeat
  if (peak(flux) > flux_cutoff) peakwin win_mult
  clean clean_niter,clean_gain
  flux_cutoff = imstat(rms) * dynam
  selfcal
  new_peak = peak(flux)
  until(new_peak<=flux_cutoff)
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00300303 Jy
! Component: 100  -  total flux cleaned = 0.00561293 Jy
! Component: 150  -  total flux cleaned = 0.00702558 Jy
! Component: 200  -  total flux cleaned = 0.00788838 Jy
! Total flux subtracted in 200 components = 0.00788838 Jy
! Clean residual min=-0.016694 max=0.028813 Jy/beam
! Clean residual mean=0.000796 rms=0.007644 Jy/beam
! Combined flux in latest and established models = 2.25339 Jy
! Performing phase self-cal
! Adding 23 model components to the UV plane model.
! The established model now contains 395 components and 2.25339 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.349052Jy  sigma=0.001458
! Fit after  self-cal, rms=0.349055Jy  sigma=0.001458
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
win_mult = win_mult / 1.2
print "************** FINISHED AMPLITUDE SELF-CAL **************"
! ************** FINISHED AMPLITUDE SELF-CAL ************** 
restore
! Inverting map and beam 
! Estimated beam: bmin=2.007 mas, bmaj=4.146 mas, bpa=15.21 degrees
! Estimated noise=634.139 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.007 x 4.146 at 15.21 degrees (North through East)
! Clean map  min=-0.0092532  max=0.73378 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.349055Jy  sigma=0.001458
! Fit after  self-cal, rms=0.882508Jy  sigma=0.001351
! Inverting map and beam 
! Estimated beam: bmin=1.999 mas, bmaj=4.141 mas, bpa=15.29 degrees
! Estimated noise=631.899 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000427726 Jy
! Component: 100  -  total flux cleaned = -0.000612275 Jy
! Component: 150  -  total flux cleaned = -0.000903136 Jy
! Component: 200  -  total flux cleaned = -0.00117051 Jy
! Total flux subtracted in 200 components = -0.00117051 Jy
! Clean residual min=-0.009583 max=0.014110 Jy/beam
! Clean residual mean=0.000301 rms=0.003434 Jy/beam
! Combined flux in latest and established models = 2.25222 Jy
! Performing phase self-cal
! Adding 26 model components to the UV plane model.
! The established model now contains 408 components and 2.25222 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.882435Jy  sigma=0.001350
! Fit after  self-cal, rms=0.882387Jy  sigma=0.001350
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000115126 Jy
! Component: 100  -  total flux cleaned = -0.000350481 Jy
! Component: 150  -  total flux cleaned = -0.000126694 Jy
! Component: 200  -  total flux cleaned = 8.07313e-05 Jy
! Total flux subtracted in 200 components = 8.07313e-05 Jy
! Clean residual min=-0.009505 max=0.014074 Jy/beam
! Clean residual mean=0.000302 rms=0.003435 Jy/beam
! Combined flux in latest and established models = 2.2523 Jy
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 416 components and 2.2523 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.882361Jy  sigma=0.001349
! Fit after  self-cal, rms=0.882362Jy  sigma=0.001349
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.882362Jy  sigma=0.001349
! Fit after  self-cal, rms=0.968361Jy  sigma=0.001348
print "************** CLEARING MODEL AND STARTING OVER **********"
! ************** CLEARING MODEL AND STARTING OVER ********** 
clrmod true
! clrmod: Cleared the established and tentative models.
uvw 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.798 mas, bmaj=3.467 mas, bpa=7.545 degrees
! Estimated noise=1100.71 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.619855 Jy
! Component: 100  -  total flux cleaned = 0.900739 Jy
! Component: 150  -  total flux cleaned = 1.06946 Jy
! Component: 200  -  total flux cleaned = 1.19061 Jy
! Total flux subtracted in 200 components = 1.19061 Jy
! Clean residual min=-0.024914 max=0.069194 Jy/beam
! Clean residual mean=-0.000106 rms=0.007317 Jy/beam
! Combined flux in latest and established models = 1.19061 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 25 model components to the UV plane model.
! The established model now contains 25 components and 1.19061 Jy
! Inverting map and beam 
! Estimated beam: bmin=1.999 mas, bmaj=4.143 mas, bpa=15.24 degrees
! Estimated noise=631.529 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.22059 Jy
! Component: 100  -  total flux cleaned = 0.375076 Jy
! Component: 150  -  total flux cleaned = 0.488868 Jy
! Component: 200  -  total flux cleaned = 0.573657 Jy
! Total flux subtracted in 200 components = 0.573657 Jy
! Clean residual min=-0.023573 max=0.053744 Jy/beam
! Clean residual mean=0.002074 rms=0.011800 Jy/beam
! Combined flux in latest and established models = 1.76427 Jy
! Adding 95 model components to the UV plane model.
! The established model now contains 105 components and 1.76427 Jy
! Inverting map 
! Added new window around map position (-15.5, 5.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0647679 Jy
! Component: 100  -  total flux cleaned = 0.114923 Jy
! Component: 150  -  total flux cleaned = 0.15736 Jy
! Component: 200  -  total flux cleaned = 0.194098 Jy
! Total flux subtracted in 200 components = 0.194098 Jy
! Clean residual min=-0.016329 max=0.032108 Jy/beam
! Clean residual mean=0.001332 rms=0.007617 Jy/beam
! Combined flux in latest and established models = 1.95837 Jy
! Adding 90 model components to the UV plane model.
! The established model now contains 179 components and 1.95837 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 20 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=8.482 mas, bmaj=10.1 mas, bpa=-7.202 degrees
! Estimated noise=825.401 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0715367 Jy
! Component: 100  -  total flux cleaned = 0.123707 Jy
! Component: 150  -  total flux cleaned = 0.16323 Jy
! Component: 200  -  total flux cleaned = 0.193905 Jy
! Total flux subtracted in 200 components = 0.193905 Jy
! Clean residual min=-0.022184 max=0.036795 Jy/beam
! Clean residual mean=0.001629 rms=0.009823 Jy/beam
! Combined flux in latest and established models = 2.15227 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.024205 Jy
! Component: 100  -  total flux cleaned = 0.0434883 Jy
! Component: 150  -  total flux cleaned = 0.0589433 Jy
! Component: 200  -  total flux cleaned = 0.0713702 Jy
! Total flux subtracted in 200 components = 0.0713702 Jy
! Clean residual min=-0.018430 max=0.028651 Jy/beam
! Clean residual mean=0.000970 rms=0.007881 Jy/beam
! Combined flux in latest and established models = 2.22364 Jy
! Adding 160 model components to the UV plane model.
! The established model now contains 317 components and 2.22364 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=1.999 mas, bmaj=4.143 mas, bpa=15.24 degrees
! Estimated noise=631.529 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000967279 Jy
! Component: 100  -  total flux cleaned = 0.00339777 Jy
! Component: 150  -  total flux cleaned = 0.00649467 Jy
! Component: 200  -  total flux cleaned = 0.009311 Jy
! Total flux subtracted in 200 components = 0.009311 Jy
! Clean residual min=-0.010805 max=0.014321 Jy/beam
! Clean residual mean=0.000341 rms=0.003548 Jy/beam
! Combined flux in latest and established models = 2.23295 Jy
selfcal
! Performing phase self-cal
! Adding 27 model components to the UV plane model.
! The established model now contains 336 components and 2.23295 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.968035Jy  sigma=0.001355
! Fit after  self-cal, rms=0.967752Jy  sigma=0.001353
wmodel J1230+1223_S_2001_04_09_pus_map.mod
! Writing 336 model components to file: J1230+1223_S_2001_04_09_pus_map.mod
wobs J1230+1223_S_2001_04_09_pus_uvs.fits
! Writing UV FITS file: J1230+1223_S_2001_04_09_pus_uvs.fits
wwins J1230+1223_S_2001_04_09_pus_map.win
! wwins: Wrote 11 windows to J1230+1223_S_2001_04_09_pus_map.win
x = (field_size-8) * field_cell / 4
addwin -x,x,-x,x
clean (field_size*4),0.01
! Inverting map 
! clean: niter=4096  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.00662054 Jy
! Component: 100  -  total flux cleaned = 0.012523 Jy
! Component: 150  -  total flux cleaned = 0.0180297 Jy
! Component: 200  -  total flux cleaned = 0.0232439 Jy
! Component: 250  -  total flux cleaned = 0.0280391 Jy
! Component: 300  -  total flux cleaned = 0.0322644 Jy
! Component: 350  -  total flux cleaned = 0.0365282 Jy
! Component: 400  -  total flux cleaned = 0.0406603 Jy
! Component: 450  -  total flux cleaned = 0.0444952 Jy
! Component: 500  -  total flux cleaned = 0.048391 Jy
! Component: 550  -  total flux cleaned = 0.0521777 Jy
! Component: 600  -  total flux cleaned = 0.0553839 Jy
! Component: 650  -  total flux cleaned = 0.0588218 Jy
! Component: 700  -  total flux cleaned = 0.0618699 Jy
! Component: 750  -  total flux cleaned = 0.0646981 Jy
! Component: 800  -  total flux cleaned = 0.066884 Jy
! Component: 850  -  total flux cleaned = 0.0684546 Jy
! Component: 900  -  total flux cleaned = 0.0699972 Jy
! Component: 950  -  total flux cleaned = 0.0715102 Jy
! Component: 1000  -  total flux cleaned = 0.0727261 Jy
! Component: 1050  -  total flux cleaned = 0.073923 Jy
! Component: 1100  -  total flux cleaned = 0.0752307 Jy
! Component: 1150  -  total flux cleaned = 0.0765172 Jy
! Component: 1200  -  total flux cleaned = 0.0777818 Jy
! Component: 1250  -  total flux cleaned = 0.0789018 Jy
! Component: 1300  -  total flux cleaned = 0.0802501 Jy
! Component: 1350  -  total flux cleaned = 0.0812153 Jy
! Component: 1400  -  total flux cleaned = 0.0821696 Jy
! Component: 1450  -  total flux cleaned = 0.0829937 Jy
! Component: 1500  -  total flux cleaned = 0.0838046 Jy
! Component: 1550  -  total flux cleaned = 0.0842632 Jy
! Component: 1600  -  total flux cleaned = 0.0848319 Jy
! Component: 1650  -  total flux cleaned = 0.0855039 Jy
! Component: 1700  -  total flux cleaned = 0.0860578 Jy
! Component: 1750  -  total flux cleaned = 0.0867153 Jy
! Component: 1800  -  total flux cleaned = 0.0873637 Jy
! Component: 1850  -  total flux cleaned = 0.0877917 Jy
! Component: 1900  -  total flux cleaned = 0.0885315 Jy
! Component: 1950  -  total flux cleaned = 0.0888454 Jy
! Component: 2000  -  total flux cleaned = 0.0891565 Jy
! Component: 2050  -  total flux cleaned = 0.0896705 Jy
! Component: 2100  -  total flux cleaned = 0.0901784 Jy
! Component: 2150  -  total flux cleaned = 0.0906806 Jy
! Component: 2200  -  total flux cleaned = 0.0909787 Jy
! Component: 2250  -  total flux cleaned = 0.0913723 Jy
! Component: 2300  -  total flux cleaned = 0.0917619 Jy
! Component: 2350  -  total flux cleaned = 0.0920513 Jy
! Component: 2400  -  total flux cleaned = 0.0922431 Jy
! Component: 2450  -  total flux cleaned = 0.0928121 Jy
! Component: 2500  -  total flux cleaned = 0.0929048 Jy
! Component: 2550  -  total flux cleaned = 0.092999 Jy
! Component: 2600  -  total flux cleaned = 0.0932763 Jy
! Component: 2650  -  total flux cleaned = 0.0935515 Jy
! Component: 2700  -  total flux cleaned = 0.093642 Jy
! Component: 2750  -  total flux cleaned = 0.0938214 Jy
! Component: 2800  -  total flux cleaned = 0.0940885 Jy
! Component: 2850  -  total flux cleaned = 0.0940004 Jy
! Component: 2900  -  total flux cleaned = 0.0944384 Jy
! Component: 2950  -  total flux cleaned = 0.0944382 Jy
! Component: 3000  -  total flux cleaned = 0.0946106 Jy
! Component: 3050  -  total flux cleaned = 0.0946959 Jy
! Component: 3100  -  total flux cleaned = 0.0949499 Jy
! Component: 3150  -  total flux cleaned = 0.0951174 Jy
! Component: 3200  -  total flux cleaned = 0.0950341 Jy
! Component: 3250  -  total flux cleaned = 0.0950342 Jy
! Component: 3300  -  total flux cleaned = 0.0950338 Jy
! Component: 3350  -  total flux cleaned = 0.0950335 Jy
! Component: 3400  -  total flux cleaned = 0.0950336 Jy
! Component: 3450  -  total flux cleaned = 0.0951142 Jy
! Component: 3500  -  total flux cleaned = 0.0951141 Jy
! Component: 3550  -  total flux cleaned = 0.0948762 Jy
! Component: 3600  -  total flux cleaned = 0.0947984 Jy
! Component: 3650  -  total flux cleaned = 0.0947989 Jy
! Component: 3700  -  total flux cleaned = 0.0947991 Jy
! Component: 3750  -  total flux cleaned = 0.0947223 Jy
! Component: 3800  -  total flux cleaned = 0.0947984 Jy
! Component: 3850  -  total flux cleaned = 0.0949494 Jy
! Component: 3900  -  total flux cleaned = 0.0945722 Jy
! Component: 3950  -  total flux cleaned = 0.0947972 Jy
! Component: 4000  -  total flux cleaned = 0.0947966 Jy
! Component: 4050  -  total flux cleaned = 0.0947227 Jy
! Total flux subtracted in 4096 components = 0.0946493 Jy
! Clean residual min=-0.003707 max=0.004806 Jy/beam
! Clean residual mean=0.000017 rms=0.001435 Jy/beam
! Combined flux in latest and established models = 2.3276 Jy
keep
! Adding 423 model components to the UV plane model.
! The established model now contains 754 components and 2.3276 Jy
![Exited script file: muppet]
wmap "J1230+1223_S_2001_04_09_pus\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 1.999 x 4.143 at 15.24 degrees (North through East)
! Clean map  min=-0.0058008  max=0.72783 Jy/beam
! Writing clean map to FITS file: J1230+1223_S_2001_04_09_pus_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 0.727826 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.00144785 
device "J1230+1223_S_2001_04_09_pus\_map.eps/VPS"
! Attempting to open device: 'J1230+1223_S_2001_04_09_pus_map.eps/VPS'
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.596783 0.596783 1.19357 2.38713 4.77426 9.54852 19.097 38.1941 76.3882
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0043 
print "peak=",int(1000*peak_flux)/1000
! peak= 0.727 
mapcolor none
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
mapplot cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 1.999 x 4.143 at 15.24 degrees (North through East)
! Clean map  min=-0.0058008  max=0.72783 Jy/beam
mapcolor color, 0.9
! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5.
device "J1230+1223_S_2001_04_09_pus\_map.gif/vgif"
! Attempting to open device: 'J1230+1223_S_2001_04_09_pus_map.gif/vgif'
mapplot cln
quit
! Quitting program
! Log file J1230+1223_S_2001_04_09_pus_uvs.log closed on Fri Oct  5 17:16:47 2007