! Started logfile: J0136+4751_S_1999_12_20_pus_uvs.log on Fri Dec 21 12:47:24 2007
obs J0136+4751_S_1999_12_20_pus_uva.fits
! Reading UV FITS file: J0136+4751_S_1999_12_20_pus_uva.fits
! AN table 1: 125 integrations on 190 of 190 possible baselines.
! AN table 2: 160 integrations on 190 of 190 possible baselines.
! AN table 3: 11 integrations on 190 of 190 possible baselines.
! AN table 4: 148 integrations on 190 of 190 possible baselines.
! AN table 5: 69 integrations on 190 of 190 possible baselines.
! Apparent sampling: 0.258767 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J0136+4751
! 
! 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.22497e+09        8e+06          1        8e+06
!  02        2  2.24497e+09        8e+06          1        8e+06
!  03        3  2.33497e+09        8e+06          1        8e+06
!  04        4  2.36497e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 2302 lines of history.
! 
! Reading 100888 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 = 512 
print "field_cell =",field_cell
! field_cell = 0.5 
print "taper_size =",taper_size
! taper_size = 20 
![@muppet J0136+4751_S_1999_12_20_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 = 512x512 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 1596 telescope corrections were flagged in sub-array 1.
!  A total of 2259 telescope corrections were flagged in sub-array 2.
!  A total of 220 telescope corrections were flagged in sub-array 3.
!  A total of 1603 telescope corrections were flagged in sub-array 4.
!  A total of 580 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 2.
!  A total of 1597 telescope corrections were flagged in sub-array 1.
!  A total of 2262 telescope corrections were flagged in sub-array 2.
!  A total of 220 telescope corrections were flagged in sub-array 3.
!  A total of 1603 telescope corrections were flagged in sub-array 4.
!  A total of 582 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 3.
!  A total of 1596 telescope corrections were flagged in sub-array 1.
!  A total of 2284 telescope corrections were flagged in sub-array 2.
!  A total of 220 telescope corrections were flagged in sub-array 3.
!  A total of 1682 telescope corrections were flagged in sub-array 4.
!  A total of 604 telescope corrections were flagged in sub-array 5.
! 
! Correcting IF 4.
!  A total of 1597 telescope corrections were flagged in sub-array 1.
!  A total of 2258 telescope corrections were flagged in sub-array 2.
!  A total of 220 telescope corrections were flagged in sub-array 3.
!  A total of 1602 telescope corrections were flagged in sub-array 4.
!  A total of 585 telescope corrections were flagged in sub-array 5.
! 
! Fit before self-cal, rms=1.080695Jy  sigma=0.008052
! Fit after  self-cal, rms=0.381661Jy  sigma=0.002941
! 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.944 mas, bmaj=2.479 mas, bpa=-2.991 degrees
! Estimated noise=573.608 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.839177 Jy
! Component: 100  -  total flux cleaned = 1.0742 Jy
! Total flux subtracted in 100 components = 1.0742 Jy
! Clean residual min=-0.013721 max=0.104607 Jy/beam
! Clean residual mean=0.000088 rms=0.004175 Jy/beam
! Combined flux in latest and established models = 1.0742 Jy
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 5 components and 1.0742 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.340661Jy  sigma=0.002571
! Fit after  self-cal, rms=0.340331Jy  sigma=0.002570
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.1151 Jy
! Component: 100  -  total flux cleaned = 0.176251 Jy
! Total flux subtracted in 100 components = 0.176251 Jy
! Clean residual min=-0.012339 max=0.032290 Jy/beam
! Clean residual mean=0.000048 rms=0.002619 Jy/beam
! Combined flux in latest and established models = 1.25045 Jy
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 16 components and 1.25045 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.282479Jy  sigma=0.001901
! Fit after  self-cal, rms=0.281656Jy  sigma=0.001899
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0421343 Jy
! Component: 100  -  total flux cleaned = 0.0719973 Jy
! Total flux subtracted in 100 components = 0.0719973 Jy
! Clean residual min=-0.011181 max=0.020080 Jy/beam
! Clean residual mean=0.000031 rms=0.002269 Jy/beam
! Combined flux in latest and established models = 1.32244 Jy
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 23 components and 1.32244 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.270324Jy  sigma=0.001725
! Fit after  self-cal, rms=0.270057Jy  sigma=0.001725
! Inverting map 
! Added new window around map position (-2.5, 2.5).
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.024589 Jy
! Component: 100  -  total flux cleaned = 0.0433986 Jy
! Total flux subtracted in 100 components = 0.0433986 Jy
! Clean residual min=-0.010703 max=0.011701 Jy/beam
! Clean residual mean=0.000020 rms=0.002136 Jy/beam
! Combined flux in latest and established models = 1.36584 Jy
! Performing phase self-cal
! Adding 14 model components to the UV plane model.
! The established model now contains 30 components and 1.36584 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.265695Jy  sigma=0.001650
! Fit after  self-cal, rms=0.265443Jy  sigma=0.001649
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0159495 Jy
! Component: 100  -  total flux cleaned = 0.0284529 Jy
! Total flux subtracted in 100 components = 0.0284529 Jy
! Clean residual min=-0.010303 max=0.011128 Jy/beam
! Clean residual mean=0.000013 rms=0.002066 Jy/beam
! Combined flux in latest and established models = 1.3943 Jy
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 38 components and 1.3943 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.263548Jy  sigma=0.001611
! Fit after  self-cal, rms=0.263408Jy  sigma=0.001611
! 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.465 mas, bmaj=3.091 mas, bpa=-5.229 degrees
! Estimated noise=307.039 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.0298693 Jy
! Component: 100  -  total flux cleaned = 0.0485016 Jy
! Component: 150  -  total flux cleaned = 0.0608016 Jy
! Component: 200  -  total flux cleaned = 0.0691013 Jy
! Total flux subtracted in 200 components = 0.0691013 Jy
! Clean residual min=-0.007853 max=0.016189 Jy/beam
! Clean residual mean=-0.000090 rms=0.002267 Jy/beam
! Combined flux in latest and established models = 1.4634 Jy
! Performing phase self-cal
! Adding 28 model components to the UV plane model.
! The established model now contains 53 components and 1.4634 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.261583Jy  sigma=0.001556
! Fit after  self-cal, rms=0.261147Jy  sigma=0.001553
! Inverting map 
! Added new window around map position (-1.5, 7.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0173034 Jy
! Component: 100  -  total flux cleaned = 0.0275908 Jy
! Component: 150  -  total flux cleaned = 0.0348315 Jy
! Component: 200  -  total flux cleaned = 0.0400701 Jy
! Total flux subtracted in 200 components = 0.0400701 Jy
! Clean residual min=-0.007874 max=0.008975 Jy/beam
! Clean residual mean=-0.000277 rms=0.001968 Jy/beam
! Combined flux in latest and established models = 1.50347 Jy
! Performing phase self-cal
! Adding 28 model components to the UV plane model.
! The established model now contains 80 components and 1.50347 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.260221Jy  sigma=0.001537
! Fit after  self-cal, rms=0.260083Jy  sigma=0.001536
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
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        0.99     FD        0.99     GC        1.03     HH        1.00*    
!   HN        1.02     KK        0.98     KP        1.08     LA        1.00     
!   MC        1.00*    MK        0.98     NY        0.99     ON        1.03     
!   OV        0.94     PT        1.01     SC        1.12     TS        1.01     
!   WF        1.01     WZ        1.16     GN        0.99     NL        1.00     
! 
!  A total of 27 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.02     FD        1.04     GC        1.04     HH        1.00*    
!   HN        1.00     KK        1.01     KP        1.14     LA        1.03     
!   MC        1.00     MK        0.99     NY        0.97     ON        1.24     
!   OV        0.99     PT        1.07     SC        1.07     TS        1.08     
!   WF        1.00     WZ        0.98     GN        1.07     NL        1.05     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.00*    FD        0.00*    GC        0.00*    HH        0.00*    
!   HN        0.00*    KK        0.00*    KP        0.00*    LA        0.00*    
!   MC        0.00*    MK        0.00*    NY        0.00*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TS        0.00*    
!   WF        0.00*    WZ        0.00*    GN        0.00*    NL        0.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.96     FD        0.98     GC        1.04     HH        1.00*    
!   HN        1.08     KK        1.06     KP        0.99     LA        0.98     
!   MC        0.98     MK        0.98     NY        1.00     ON        1.06     
!   OV        0.96     PT        1.00     SC        1.02     TS        1.05     
!   WF        1.02     WZ        0.95     GN        0.99     NL        0.98     
! 
!  A total of 9 telescope corrections were flagged in sub-array 5.
!  Telescope amplitude corrections in sub-array 5:
!   BR        0.95     FD        0.99     GC        1.01     HH        1.00*    
!   HN        1.07     KK        1.01     KP        0.95     LA        0.95     
!   MC        0.90     MK        1.03     NY        0.99     ON        1.02     
!   OV        1.00     PT        1.01     SC        1.02     TS        1.00*    
!   WF        1.15     WZ        0.96     GN        0.98     NL        0.95     
! 
! 
! Correcting IF 2.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.99     FD        1.04     GC        1.02     HH        1.00*    
!   HN        1.08     KK        0.95     KP        1.01     LA        0.99     
!   MC        1.00*    MK        0.96     NY        0.97     ON        1.07     
!   OV        0.94     PT        0.99     SC        1.12     TS        1.04     
!   WF        1.04     WZ        1.22     GN        1.00     NL        1.00     
! 
!  A total of 36 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.03     FD        1.08     GC        1.04     HH        1.00*    
!   HN        1.06     KK        0.95     KP        1.08     LA        1.02     
!   MC        1.27     MK        0.95     NY        0.96     ON        1.25     
!   OV        0.97     PT        1.05     SC        1.06     TS        1.08     
!   WF        1.01     WZ        1.00     GN        1.07     NL        1.03     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.00*    FD        0.00*    GC        0.00*    HH        0.00*    
!   HN        0.00*    KK        0.00*    KP        0.00*    LA        0.00*    
!   MC        0.00*    MK        0.00*    NY        0.00*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TS        0.00*    
!   WF        0.00*    WZ        0.00*    GN        0.00*    NL        0.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.96     FD        1.03     GC        1.02     HH        1.00*    
!   HN        1.18     KK        0.95     KP        0.98     LA        0.96     
!   MC        1.07     MK        0.95     NY        0.99     ON        1.07     
!   OV        0.94     PT        0.98     SC        1.03     TS        1.05     
!   WF        1.01     WZ        0.98     GN        1.00     NL        0.94     
! 
!  A total of 9 telescope corrections were flagged in sub-array 5.
!  Telescope amplitude corrections in sub-array 5:
!   BR        0.96     FD        1.04     GC        0.99     HH        1.00*    
!   HN        1.03     KK        0.99     KP        0.98     LA        0.93     
!   MC        1.08     MK        1.01     NY        1.00     ON        1.06     
!   OV        0.98     PT        0.97     SC        0.97     TS        1.00*    
!   WF        1.09     WZ        1.00     GN        0.98     NL        0.93     
! 
! 
! Correcting IF 3.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.98     FD        1.01     GC        1.02     HH        1.00*    
!   HN        1.04     KK        0.97     KP        0.98     LA        1.00     
!   MC        1.00*    MK        0.94     NY        0.96     ON        0.97     
!   OV        0.95     PT        0.97     SC        1.15     TS        0.98     
!   WF        1.00     WZ        1.17     GN        1.00     NL        0.99     
! 
!  A total of 42 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.01     FD        1.05     GC        1.05     HH        1.00*    
!   HN        1.01     KK        0.98     KP        1.05     LA        1.04     
!   MC        1.00*    MK        0.96     NY        0.96     ON        1.28     
!   OV        0.99     PT        1.02     SC        1.12     TS        1.02     
!   WF        0.99     WZ        1.02     GN        1.12     NL        1.01     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.00*    FD        0.00*    GC        0.00*    HH        0.00*    
!   HN        0.00*    KK        0.00*    KP        0.00*    LA        0.00*    
!   MC        0.00*    MK        0.00*    NY        0.00*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TS        0.00*    
!   WF        0.00*    WZ        0.00*    GN        0.00*    NL        0.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.95     FD        1.00     GC        1.01     HH        1.00*    
!   HN        1.08     KK        0.98     KP        0.97     LA        0.99     
!   MC        1.00*    MK        0.95     NY        1.00     ON        1.04     
!   OV        0.98     PT        0.95     SC        1.04     TS        1.00     
!   WF        1.01     WZ        0.99     GN        1.02     NL        0.96     
! 
!  A total of 9 telescope corrections were flagged in sub-array 5.
!  Telescope amplitude corrections in sub-array 5:
!   BR        0.94     FD        1.00     GC        1.01     HH        1.00*    
!   HN        0.99     KK        1.03     KP        0.98     LA        0.95     
!   MC        1.00*    MK        1.00     NY        1.00     ON        1.08     
!   OV        1.00     PT        0.95     SC        1.02     TS        1.00*    
!   WF        1.06     WZ        1.03     GN        0.97     NL        0.93     
! 
! 
! Correcting IF 4.
!  A total of 18 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.96     FD        1.00     GC        1.02     HH        1.00*    
!   HN        1.04     KK        0.97     KP        1.00     LA        1.01     
!   MC        1.00*    MK        0.97     NY        0.96     ON        0.98     
!   OV        0.92     PT        1.03     SC        1.14     TS        1.00     
!   WF        1.00     WZ        1.16     GN        1.03     NL        0.97     
! 
!  A total of 24 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        0.99     FD        1.04     GC        1.05     HH        1.00*    
!   HN        1.02     KK        0.96     KP        1.07     LA        1.04     
!   MC        0.99     MK        0.98     NY        0.95     ON        1.29     
!   OV        0.95     PT        1.07     SC        1.10     TS        1.05     
!   WF        1.01     WZ        1.02     GN        1.10     NL        1.01     
! 
!  Telescope amplitude corrections in sub-array 3:
!   BR        0.00*    FD        0.00*    GC        0.00*    HH        0.00*    
!   HN        0.00*    KK        0.00*    KP        0.00*    LA        0.00*    
!   MC        0.00*    MK        0.00*    NY        0.00*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TS        0.00*    
!   WF        0.00*    WZ        0.00*    GN        0.00*    NL        0.00*    
! 
!  A total of 15 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.93     FD        0.99     GC        1.02     HH        1.00*    
!   HN        1.16     KK        1.01     KP        0.98     LA        0.99     
!   MC        0.91     MK        0.99     NY        0.99     ON        1.04     
!   OV        0.94     PT        1.00     SC        1.04     TS        1.00     
!   WF        1.01     WZ        0.98     GN        1.01     NL        0.95     
! 
!  A total of 9 telescope corrections were flagged in sub-array 5.
!  Telescope amplitude corrections in sub-array 5:
!   BR        0.93     FD        1.00     GC        1.00     HH        1.00*    
!   HN        1.05     KK        1.00     KP        1.00     LA        0.96     
!   MC        1.19     MK        1.04     NY        0.98     ON        1.14     
!   OV        0.96     PT        0.96     SC        1.02     TS        1.00*    
!   WF        1.18     WZ        1.02     GN        0.98     NL        0.93     
! 
! 
! Fit before self-cal, rms=0.260083Jy  sigma=0.001536
! Fit after  self-cal, rms=0.261927Jy  sigma=0.001406
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=2.477 mas, bmaj=3.118 mas, bpa=-6.026 degrees
! Estimated noise=308.379 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.00412771 Jy
! Component: 100  -  total flux cleaned = -0.00365935 Jy
! Component: 150  -  total flux cleaned = -0.00311767 Jy
! Component: 200  -  total flux cleaned = -0.0026466 Jy
! Total flux subtracted in 200 components = -0.0026466 Jy
! Clean residual min=-0.005431 max=0.005446 Jy/beam
! Clean residual mean=0.000014 rms=0.001225 Jy/beam
! Combined flux in latest and established models = 1.50082 Jy
! Performing phase self-cal
! Adding 21 model components to the UV plane model.
! The established model now contains 94 components and 1.50082 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.261545Jy  sigma=0.001402
! Fit after  self-cal, rms=0.261659Jy  sigma=0.001401
! 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.261659Jy  sigma=0.001401
! Fit after  self-cal, rms=0.267540Jy  sigma=0.001340
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.478 mas, bmaj=3.125 mas, bpa=-6.182 degrees
! Estimated noise=308.407 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00112311 Jy
! Component: 100  -  total flux cleaned = -0.000226773 Jy
! Component: 150  -  total flux cleaned = 0.000169089 Jy
! Component: 200  -  total flux cleaned = 0.000408983 Jy
! Total flux subtracted in 200 components = 0.000408983 Jy
! Clean residual min=-0.003125 max=0.005792 Jy/beam
! Clean residual mean=0.000296 rms=0.000768 Jy/beam
! Combined flux in latest and established models = 1.50123 Jy
selfcal
! Performing phase self-cal
! Adding 19 model components to the UV plane model.
! The established model now contains 104 components and 1.50123 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267415Jy  sigma=0.001338
! Fit after  self-cal, rms=0.267634Jy  sigma=0.001338
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 (-6, 2.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0050404 Jy
! Component: 100  -  total flux cleaned = 0.00563925 Jy
! Component: 150  -  total flux cleaned = 0.00575822 Jy
! Component: 200  -  total flux cleaned = 0.00607052 Jy
! Total flux subtracted in 200 components = 0.00607052 Jy
! Clean residual min=-0.002843 max=0.004169 Jy/beam
! Clean residual mean=0.000268 rms=0.000715 Jy/beam
! Combined flux in latest and established models = 1.5073 Jy
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 126 components and 1.5073 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267456Jy  sigma=0.001336
! Fit after  self-cal, rms=0.267431Jy  sigma=0.001336
! Inverting map 
! Added new window around map position (4, 12).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00342381 Jy
! Component: 100  -  total flux cleaned = 0.00423476 Jy
! Component: 150  -  total flux cleaned = 0.00473847 Jy
! Component: 200  -  total flux cleaned = 0.00511049 Jy
! Total flux subtracted in 200 components = 0.00511049 Jy
! Clean residual min=-0.002566 max=0.003488 Jy/beam
! Clean residual mean=0.000246 rms=0.000688 Jy/beam
! Combined flux in latest and established models = 1.51241 Jy
! Performing phase self-cal
! Adding 33 model components to the UV plane model.
! The established model now contains 142 components and 1.51241 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267294Jy  sigma=0.001335
! Fit after  self-cal, rms=0.267285Jy  sigma=0.001335
! Inverting map 
! Added new window around map position (1.5, 23).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00323114 Jy
! Component: 100  -  total flux cleaned = 0.0045274 Jy
! Component: 150  -  total flux cleaned = 0.00523601 Jy
! Component: 200  -  total flux cleaned = 0.00581431 Jy
! Total flux subtracted in 200 components = 0.00581431 Jy
! Clean residual min=-0.002470 max=0.003121 Jy/beam
! Clean residual mean=0.000225 rms=0.000657 Jy/beam
! Combined flux in latest and established models = 1.51822 Jy
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 168 components and 1.51822 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267177Jy  sigma=0.001334
! Fit after  self-cal, rms=0.267173Jy  sigma=0.001334
! Inverting map 
! Added new window around map position (-5, 16).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00319177 Jy
! Component: 100  -  total flux cleaned = 0.0035412 Jy
! Component: 150  -  total flux cleaned = 0.00400664 Jy
! Component: 200  -  total flux cleaned = 0.00444197 Jy
! Total flux subtracted in 200 components = 0.00444197 Jy
! Clean residual min=-0.002427 max=0.002853 Jy/beam
! Clean residual mean=0.000208 rms=0.000634 Jy/beam
! Combined flux in latest and established models = 1.52267 Jy
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 183 components and 1.52267 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267098Jy  sigma=0.001334
! Fit after  self-cal, rms=0.267095Jy  sigma=0.001334
! Inverting map 
! Added new window around map position (21, -6.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00192059 Jy
! Component: 100  -  total flux cleaned = 0.0028133 Jy
! Component: 150  -  total flux cleaned = 0.00336847 Jy
! Component: 200  -  total flux cleaned = 0.00389672 Jy
! Total flux subtracted in 200 components = 0.00389672 Jy
! Clean residual min=-0.002494 max=0.002681 Jy/beam
! Clean residual mean=0.000190 rms=0.000630 Jy/beam
! Combined flux in latest and established models = 1.52656 Jy
! Performing phase self-cal
! Adding 51 model components to the UV plane model.
! The established model now contains 200 components and 1.52656 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267029Jy  sigma=0.001333
! Fit after  self-cal, rms=0.267029Jy  sigma=0.001333
! Inverting map 
! Added new window around map position (8, -51).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00209375 Jy
! Component: 100  -  total flux cleaned = 0.00255302 Jy
! Component: 150  -  total flux cleaned = 0.00286605 Jy
! Component: 200  -  total flux cleaned = 0.00322737 Jy
! Total flux subtracted in 200 components = 0.00322737 Jy
! Clean residual min=-0.002429 max=0.002478 Jy/beam
! Clean residual mean=0.000177 rms=0.000626 Jy/beam
! Combined flux in latest and established models = 1.52979 Jy
! Performing phase self-cal
! Adding 54 model components to the UV plane model.
! The established model now contains 213 components and 1.52979 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.266970Jy  sigma=0.001333
! Fit after  self-cal, rms=0.266968Jy  sigma=0.001333
! 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.266968Jy  sigma=0.001333
! Fit after  self-cal, rms=0.267483Jy  sigma=0.001323
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=2.477 mas, bmaj=3.135 mas, bpa=-6.336 degrees
! Estimated noise=308.679 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000185517 Jy
! Component: 100  -  total flux cleaned = 0.000499468 Jy
! Component: 150  -  total flux cleaned = 0.00120302 Jy
! Component: 200  -  total flux cleaned = 0.0016851 Jy
! Total flux subtracted in 200 components = 0.0016851 Jy
! Clean residual min=-0.002111 max=0.002571 Jy/beam
! Clean residual mean=0.000171 rms=0.000547 Jy/beam
! Combined flux in latest and established models = 1.53148 Jy
selfcal
! Performing phase self-cal
! Adding 31 model components to the UV plane model.
! The established model now contains 232 components and 1.53148 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267423Jy  sigma=0.001322
! Fit after  self-cal, rms=0.267454Jy  sigma=0.001322
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 (-21.5, 5.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00147053 Jy
! Component: 100  -  total flux cleaned = 0.00244798 Jy
! Component: 150  -  total flux cleaned = 0.00292629 Jy
! Component: 200  -  total flux cleaned = 0.00337729 Jy
! Total flux subtracted in 200 components = 0.00337729 Jy
! Clean residual min=-0.002117 max=0.002403 Jy/beam
! Clean residual mean=0.000155 rms=0.000534 Jy/beam
! Combined flux in latest and established models = 1.53485 Jy
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 248 components and 1.53485 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267399Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267403Jy  sigma=0.001321
! Inverting map 
! Added new window around map position (-26, 29.5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00160634 Jy
! Component: 100  -  total flux cleaned = 0.00243823 Jy
! Component: 150  -  total flux cleaned = 0.00297374 Jy
! Component: 200  -  total flux cleaned = 0.00336839 Jy
! Total flux subtracted in 200 components = 0.00336839 Jy
! Clean residual min=-0.002162 max=0.002316 Jy/beam
! Clean residual mean=0.000141 rms=0.000524 Jy/beam
! Combined flux in latest and established models = 1.53822 Jy
! Performing phase self-cal
! Adding 44 model components to the UV plane model.
! The established model now contains 253 components and 1.53822 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267357Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267362Jy  sigma=0.001321
! Inverting map 
! Added new window around map position (9, 35).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00188719 Jy
! Component: 100  -  total flux cleaned = 0.00282469 Jy
! Component: 150  -  total flux cleaned = 0.00326102 Jy
! Component: 200  -  total flux cleaned = 0.00362718 Jy
! Total flux subtracted in 200 components = 0.00362718 Jy
! Clean residual min=-0.002141 max=0.002318 Jy/beam
! Clean residual mean=0.000127 rms=0.000517 Jy/beam
! Combined flux in latest and established models = 1.54185 Jy
! Performing phase self-cal
! Adding 52 model components to the UV plane model.
! The established model now contains 270 components and 1.54185 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267304Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267304Jy  sigma=0.001321
! Inverting map 
! Added new window around map position (-13, -36).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00180924 Jy
! Component: 100  -  total flux cleaned = 0.00218079 Jy
! Component: 150  -  total flux cleaned = 0.00252898 Jy
! Component: 200  -  total flux cleaned = 0.00276798 Jy
! Total flux subtracted in 200 components = 0.00276798 Jy
! Clean residual min=-0.002137 max=0.002101 Jy/beam
! Clean residual mean=0.000116 rms=0.000508 Jy/beam
! Combined flux in latest and established models = 1.54462 Jy
! Performing phase self-cal
! Adding 56 model components to the UV plane model.
! The established model now contains 287 components and 1.54462 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267260Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267257Jy  sigma=0.001321
! 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.267257Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267831Jy  sigma=0.001320
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=7.967 mas, bmaj=8.477 mas, bpa=48.81 degrees
! Estimated noise=446.907 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00297079 Jy
! Component: 100  -  total flux cleaned = 0.00490594 Jy
! Component: 150  -  total flux cleaned = 0.00609204 Jy
! Component: 200  -  total flux cleaned = 0.00689968 Jy
! Total flux subtracted in 200 components = 0.00689968 Jy
! Clean residual min=-0.002924 max=0.002621 Jy/beam
! Clean residual mean=0.000306 rms=0.000803 Jy/beam
! Combined flux in latest and established models = 1.55152 Jy
selfcal
! Performing phase self-cal
! Adding 26 model components to the UV plane model.
! The established model now contains 304 components and 1.55152 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267797Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267793Jy  sigma=0.001320
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.000839756 Jy
! Component: 100  -  total flux cleaned = 0.00120754 Jy
! Component: 150  -  total flux cleaned = 0.00133584 Jy
! Component: 200  -  total flux cleaned = 0.00149914 Jy
! Total flux subtracted in 200 components = 0.00149914 Jy
! Clean residual min=-0.002906 max=0.002635 Jy/beam
! Clean residual mean=0.000286 rms=0.000788 Jy/beam
! Combined flux in latest and established models = 1.55301 Jy
! Performing phase self-cal
! Adding 38 model components to the UV plane model.
! The established model now contains 321 components and 1.55301 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.267786Jy  sigma=0.001321
! Fit after  self-cal, rms=0.267786Jy  sigma=0.001321
! 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.478 mas, bmaj=3.138 mas, bpa=-6.401 degrees
! Estimated noise=308.838 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.478 x 3.138 at -6.401 degrees (North through East)
! Clean map  min=-0.0018691  max=1.1573 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.267786Jy  sigma=0.001321
! Fit after  self-cal, rms=1.130541Jy  sigma=0.001227
! Inverting map and beam 
! Estimated beam: bmin=2.479 mas, bmaj=3.139 mas, bpa=-6.607 degrees
! Estimated noise=308.81 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00266545 Jy
! Component: 100  -  total flux cleaned = -0.00366459 Jy
! Component: 150  -  total flux cleaned = -0.00445451 Jy
! Component: 200  -  total flux cleaned = -0.00498216 Jy
! Total flux subtracted in 200 components = -0.00498216 Jy
! Clean residual min=-0.001876 max=0.001941 Jy/beam
! Clean residual mean=0.000097 rms=0.000456 Jy/beam
! Combined flux in latest and established models = 1.54803 Jy
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 325 components and 1.54803 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.130531Jy  sigma=0.001227
! Fit after  self-cal, rms=1.130507Jy  sigma=0.001227
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000573321 Jy
! Component: 100  -  total flux cleaned = -0.000724948 Jy
! Component: 150  -  total flux cleaned = -0.000770512 Jy
! Component: 200  -  total flux cleaned = -0.00076966 Jy
! Total flux subtracted in 200 components = -0.00076966 Jy
! Clean residual min=-0.001837 max=0.001936 Jy/beam
! Clean residual mean=0.000099 rms=0.000451 Jy/beam
! Combined flux in latest and established models = 1.54726 Jy
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 338 components and 1.54726 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.130503Jy  sigma=0.001227
! Fit after  self-cal, rms=1.130510Jy  sigma=0.001227
! Inverting map 
! Added new window around map position (-21.5, 39).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00143789 Jy
! Component: 100  -  total flux cleaned = 0.00157459 Jy
! Component: 150  -  total flux cleaned = 0.00170207 Jy
! Component: 200  -  total flux cleaned = 0.00174235 Jy
! Total flux subtracted in 200 components = 0.00174235 Jy
! Clean residual min=-0.001799 max=0.001826 Jy/beam
! Clean residual mean=0.000091 rms=0.000443 Jy/beam
! Combined flux in latest and established models = 1.549 Jy
! Performing phase self-cal
! Adding 46 model components to the UV plane model.
! The established model now contains 350 components and 1.549 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.130500Jy  sigma=0.001226
! Fit after  self-cal, rms=1.130484Jy  sigma=0.001226
! Inverting map 
! Performing amp+phase self-cal
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=1.130484Jy  sigma=0.001226
! Fit after  self-cal, rms=0.948281Jy  sigma=0.001226
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.956 mas, bmaj=2.502 mas, bpa=-4.895 degrees
! Estimated noise=582.442 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.85101 Jy
! Component: 100  -  total flux cleaned = 1.07958 Jy
! Component: 150  -  total flux cleaned = 1.191 Jy
! Component: 200  -  total flux cleaned = 1.25068 Jy
! Total flux subtracted in 200 components = 1.25068 Jy
! Clean residual min=-0.005126 max=0.030967 Jy/beam
! Clean residual mean=0.000076 rms=0.001762 Jy/beam
! Combined flux in latest and established models = 1.25068 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 14 model components to the UV plane model.
! The established model now contains 14 components and 1.25068 Jy
! Inverting map and beam 
! Estimated beam: bmin=2.481 mas, bmaj=3.142 mas, bpa=-6.695 degrees
! Estimated noise=308.866 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.0787319 Jy
! Component: 100  -  total flux cleaned = 0.124787 Jy
! Component: 150  -  total flux cleaned = 0.155353 Jy
! Component: 200  -  total flux cleaned = 0.177269 Jy
! Total flux subtracted in 200 components = 0.177269 Jy
! Clean residual min=-0.002423 max=0.012774 Jy/beam
! Clean residual mean=0.000605 rms=0.001695 Jy/beam
! Combined flux in latest and established models = 1.42795 Jy
! Adding 81 model components to the UV plane model.
! The established model now contains 87 components and 1.42795 Jy
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0169209 Jy
! Component: 100  -  total flux cleaned = 0.0303239 Jy
! Component: 150  -  total flux cleaned = 0.0411479 Jy
! Component: 200  -  total flux cleaned = 0.0500004 Jy
! Total flux subtracted in 200 components = 0.0500004 Jy
! Clean residual min=-0.001938 max=0.005597 Jy/beam
! Clean residual mean=0.000384 rms=0.000918 Jy/beam
! Combined flux in latest and established models = 1.47795 Jy
! Adding 72 model components to the UV plane model.
! The established model now contains 134 components and 1.47795 Jy
! Inverting map 
! Added new window around map position (3, 5).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00738177 Jy
! Component: 100  -  total flux cleaned = 0.0135679 Jy
! Component: 150  -  total flux cleaned = 0.0188214 Jy
! Component: 200  -  total flux cleaned = 0.0233506 Jy
! Total flux subtracted in 200 components = 0.0233506 Jy
! Clean residual min=-0.001683 max=0.003402 Jy/beam
! Clean residual mean=0.000284 rms=0.000645 Jy/beam
! Combined flux in latest and established models = 1.5013 Jy
! Adding 62 model components to the UV plane model.
! The established model now contains 166 components and 1.5013 Jy
! Inverting map 
! Added new window around map position (0, 18).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00413807 Jy
! Component: 100  -  total flux cleaned = 0.0076954 Jy
! Component: 150  -  total flux cleaned = 0.010906 Jy
! Component: 200  -  total flux cleaned = 0.01383 Jy
! Total flux subtracted in 200 components = 0.01383 Jy
! Clean residual min=-0.001555 max=0.002605 Jy/beam
! Clean residual mean=0.000226 rms=0.000542 Jy/beam
! Combined flux in latest and established models = 1.51513 Jy
! Adding 71 model components to the UV plane model.
! The established model now contains 200 components and 1.51513 Jy
! Inverting map 
! Added new window around map position (-8, 9).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00289913 Jy
! Component: 100  -  total flux cleaned = 0.00541941 Jy
! Component: 150  -  total flux cleaned = 0.00775511 Jy
! Component: 200  -  total flux cleaned = 0.00993272 Jy
! Total flux subtracted in 200 components = 0.00993272 Jy
! Clean residual min=-0.001566 max=0.002059 Jy/beam
! Clean residual mean=0.000184 rms=0.000494 Jy/beam
! Combined flux in latest and established models = 1.52507 Jy
! Adding 69 model components to the UV plane model.
! The established model now contains 234 components and 1.52506 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=7.973 mas, bmaj=8.484 mas, bpa=48.61 degrees
! Estimated noise=445.701 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00468875 Jy
! Component: 100  -  total flux cleaned = 0.00805403 Jy
! Component: 150  -  total flux cleaned = 0.0107478 Jy
! Component: 200  -  total flux cleaned = 0.0130309 Jy
! Total flux subtracted in 200 components = 0.0130309 Jy
! Clean residual min=-0.001799 max=0.002949 Jy/beam
! Clean residual mean=0.000440 rms=0.000717 Jy/beam
! Combined flux in latest and established models = 1.5381 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.00199141 Jy
! Component: 100  -  total flux cleaned = 0.0037414 Jy
! Component: 150  -  total flux cleaned = 0.00528476 Jy
! Component: 200  -  total flux cleaned = 0.0066475 Jy
! Total flux subtracted in 200 components = 0.0066475 Jy
! Clean residual min=-0.001854 max=0.002865 Jy/beam
! Clean residual mean=0.000344 rms=0.000678 Jy/beam
! Combined flux in latest and established models = 1.54474 Jy
! Adding 131 model components to the UV plane model.
! The established model now contains 352 components and 1.54474 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=2.481 mas, bmaj=3.142 mas, bpa=-6.695 degrees
! Estimated noise=308.866 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.000478998 Jy
! Component: 100  -  total flux cleaned = 0.000135157 Jy
! Component: 150  -  total flux cleaned = 0.000753745 Jy
! Component: 200  -  total flux cleaned = 0.00117514 Jy
! Total flux subtracted in 200 components = 0.00117514 Jy
! Clean residual min=-0.001957 max=0.001873 Jy/beam
! Clean residual mean=0.000095 rms=0.000443 Jy/beam
! Combined flux in latest and established models = 1.54592 Jy
selfcal
! Performing phase self-cal
! Adding 32 model components to the UV plane model.
! The established model now contains 369 components and 1.54592 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.948325Jy  sigma=0.001227
! Fit after  self-cal, rms=0.948886Jy  sigma=0.001226
wmodel J0136+4751_S_1999_12_20_pus_map.mod
! Writing 369 model components to file: J0136+4751_S_1999_12_20_pus_map.mod
wobs J0136+4751_S_1999_12_20_pus_uvs.fits
! Writing UV FITS file: J0136+4751_S_1999_12_20_pus_uvs.fits
wwins J0136+4751_S_1999_12_20_pus_map.win
! wwins: Wrote 17 windows to J0136+4751_S_1999_12_20_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=2048  gain=0.01  cutoff=0
! Component: 050  -  total flux cleaned = 0.000690112 Jy
! Component: 100  -  total flux cleaned = 0.00129289 Jy
! Component: 150  -  total flux cleaned = 0.00180938 Jy
! Component: 200  -  total flux cleaned = 0.00230871 Jy
! Component: 250  -  total flux cleaned = 0.00270995 Jy
! Component: 300  -  total flux cleaned = 0.00313016 Jy
! Component: 350  -  total flux cleaned = 0.00356907 Jy
! Component: 400  -  total flux cleaned = 0.00399944 Jy
! Component: 450  -  total flux cleaned = 0.00426393 Jy
! Component: 500  -  total flux cleaned = 0.00452494 Jy
! Component: 550  -  total flux cleaned = 0.00480853 Jy
! Component: 600  -  total flux cleaned = 0.00508909 Jy
! Component: 650  -  total flux cleaned = 0.00539168 Jy
! Component: 700  -  total flux cleaned = 0.00569098 Jy
! Component: 750  -  total flux cleaned = 0.00598707 Jy
! Component: 800  -  total flux cleaned = 0.0062314 Jy
! Component: 850  -  total flux cleaned = 0.00647355 Jy
! Component: 900  -  total flux cleaned = 0.0066659 Jy
! Component: 950  -  total flux cleaned = 0.00690396 Jy
! Component: 1000  -  total flux cleaned = 0.00702187 Jy
! Component: 1050  -  total flux cleaned = 0.00711587 Jy
! Component: 1100  -  total flux cleaned = 0.0073026 Jy
! Component: 1150  -  total flux cleaned = 0.00744143 Jy
! Component: 1200  -  total flux cleaned = 0.00757942 Jy
! Component: 1250  -  total flux cleaned = 0.00767073 Jy
! Component: 1300  -  total flux cleaned = 0.00782952 Jy
! Component: 1350  -  total flux cleaned = 0.00782974 Jy
! Component: 1400  -  total flux cleaned = 0.00803158 Jy
! Component: 1450  -  total flux cleaned = 0.00809845 Jy
! Component: 1500  -  total flux cleaned = 0.00820944 Jy
! Component: 1550  -  total flux cleaned = 0.00829775 Jy
! Component: 1600  -  total flux cleaned = 0.00847313 Jy
! Component: 1650  -  total flux cleaned = 0.0085602 Jy
! Component: 1700  -  total flux cleaned = 0.00871205 Jy
! Component: 1750  -  total flux cleaned = 0.0087336 Jy
! Component: 1800  -  total flux cleaned = 0.00892692 Jy
! Component: 1850  -  total flux cleaned = 0.00896973 Jy
! Component: 1900  -  total flux cleaned = 0.00916111 Jy
! Component: 1950  -  total flux cleaned = 0.0091611 Jy
! Component: 2000  -  total flux cleaned = 0.00928771 Jy
! Total flux subtracted in 2048 components = 0.00939262 Jy
! Clean residual min=-0.001225 max=0.001121 Jy/beam
! Clean residual mean=0.000054 rms=0.000394 Jy/beam
! Combined flux in latest and established models = 1.55531 Jy
keep
! Adding 157 model components to the UV plane model.
! The established model now contains 526 components and 1.55531 Jy
![Exited script file: muppet]
wmap "J0136+4751_S_1999_12_20_pus\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.481 x 3.142 at -6.695 degrees (North through East)
! Clean map  min=-0.0016151  max=1.161 Jy/beam
! Writing clean map to FITS file: J0136+4751_S_1999_12_20_pus_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 1.16104 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000393814 
device "J0136+4751_S_1999_12_20_pus\_map.eps/VPS"
! Attempting to open device: 'J0136+4751_S_1999_12_20_pus_map.eps/VPS'
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.101757 0.101757 0.203515 0.40703 0.81406 1.62812 3.25624 6.51248 13.025 26.0499 52.0998
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0011 
print "peak=",int(1000*peak_flux)/1000
! peak= 1.161 
mapcolor none
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
mapplot cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 2.481 x 3.142 at -6.695 degrees (North through East)
! Clean map  min=-0.0016151  max=1.161 Jy/beam
mapcolor color, 0.9
! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5.
device "J0136+4751_S_1999_12_20_pus\_map.gif/vgif"
! Attempting to open device: 'J0136+4751_S_1999_12_20_pus_map.gif/vgif'
mapplot cln
quit
! Quitting program
! Log file J0136+4751_S_1999_12_20_pus_uvs.log closed on Fri Dec 21 12:48:30 2007