! Started logfile: J0403-3605_X_2003_05_07_pus_uvs.log on Tue Dec 12 14:00:06 2006
obs J0403-3605_X_2003_05_07_pus_uva.fits
! Reading UV FITS file: J0403-3605_X_2003_05_07_pus_uva.fits
! AN table 1: 79 integrations on 171 of 171 possible baselines.
! AN table 2: 40 integrations on 171 of 171 possible baselines.
! AN table 3: Unused.
! AN table 4: Unused.
! AN table 5: 60 integrations on 171 of 171 possible baselines.
! AN table 6: 28 integrations on 171 of 171 possible baselines.
! AN table 7: Unused.
! AN table 8: Unused.
! AN table 9: Unused.
! Apparent sampling: 0.0529988 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: J0403-3605
! 
! There are 4 IFs, and a total of 4 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  8.40997e+09        8e+06          1        8e+06
!  02        2  8.47997e+09        8e+06          1        8e+06
!  03        3  8.79497e+09        8e+06          1        8e+06
!  04        4  8.89997e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1735 lines of history.
! 
! Reading 7504 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
float field_size
field_size = 512
float field_cell
field_cell = 0.16
float taper_size
taper_size = 90.
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.160x0.160 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 1411 telescope corrections were flagged in sub-array 1.
!  A total of 494 telescope corrections were flagged in sub-array 2.
!  A total of 817 telescope corrections were flagged in sub-array 3.
!  A total of 496 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 2.
!  A total of 1411 telescope corrections were flagged in sub-array 1.
!  A total of 494 telescope corrections were flagged in sub-array 2.
!  A total of 817 telescope corrections were flagged in sub-array 3.
!  A total of 496 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 3.
!  A total of 1411 telescope corrections were flagged in sub-array 1.
!  A total of 494 telescope corrections were flagged in sub-array 2.
!  A total of 817 telescope corrections were flagged in sub-array 3.
!  A total of 496 telescope corrections were flagged in sub-array 4.
! 
! Correcting IF 4.
!  A total of 1411 telescope corrections were flagged in sub-array 1.
!  A total of 494 telescope corrections were flagged in sub-array 2.
!  A total of 817 telescope corrections were flagged in sub-array 3.
!  A total of 496 telescope corrections were flagged in sub-array 4.
! 
! Fit before self-cal, rms=1.609296Jy  sigma=0.012797
! Fit after  self-cal, rms=1.055413Jy  sigma=0.008561
! 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=0.61 mas, bmaj=0.9439 mas, bpa=13.04 degrees
! Estimated noise=2404.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 = 1.55974 Jy
! Component: 100  -  total flux cleaned = 1.9469 Jy
! Total flux subtracted in 100 components = 1.9469 Jy
! Clean residual min=-0.060596 max=0.141209 Jy/beam
! Clean residual mean=0.001366 rms=0.020059 Jy/beam
! Combined flux in latest and established models = 1.9469 Jy
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 9 components and 1.9469 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.419025Jy  sigma=0.002131
! Fit after  self-cal, rms=0.419130Jy  sigma=0.002131
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.140999 Jy
! Component: 100  -  total flux cleaned = 0.201503 Jy
! Total flux subtracted in 100 components = 0.201503 Jy
! Clean residual min=-0.041824 max=0.048116 Jy/beam
! Clean residual mean=-0.000027 rms=0.010611 Jy/beam
! Combined flux in latest and established models = 2.1484 Jy
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 14 components and 2.1484 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.378482Jy  sigma=0.001527
! Fit after  self-cal, rms=0.378603Jy  sigma=0.001527
! 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=0.7954 mas, bmaj=1.211 mas, bpa=9.12 degrees
! Estimated noise=1215.55 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.0144033 Jy
! Component: 100  -  total flux cleaned = 0.0144118 Jy
! Component: 150  -  total flux cleaned = 0.0139827 Jy
! Component: 200  -  total flux cleaned = 0.0143699 Jy
! Total flux subtracted in 200 components = 0.0143699 Jy
! Clean residual min=-0.031515 max=0.042125 Jy/beam
! Clean residual mean=-0.001888 rms=0.010264 Jy/beam
! Combined flux in latest and established models = 2.16277 Jy
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 18 components and 2.16277 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.377034Jy  sigma=0.001507
! Fit after  self-cal, rms=0.376974Jy  sigma=0.001507
! Inverting map 
print "************** FINISHED NATURAL WEIGHTING CLEAN **************"
! ************** FINISHED NATURAL WEIGHTING CLEAN ************** 
! 
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 90 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 3 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00*    FD        0.97     GG        1.00*    HN        1.00*    
!   KK        1.00*    KP        1.09     LA        0.98     MA        1.00*    
!   MC        1.00*    MK        1.00     NL        1.00*    ON        1.00*    
!   OV        1.00     PT        1.00     SC        0.98     TC        0.96     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00*    FD        1.00     GG        1.00*    HN        1.00*    
!   KK        0.96     KP        1.00*    LA        1.01     MA        1.00*    
!   MC        1.00*    MK        0.99     NL        1.00*    ON        1.00*    
!   OV        0.98     PT        1.04     SC        1.00*    TC        0.90     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
! 
! Correcting IF 2.
!  A total of 90 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 3 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00*    FD        0.97     GG        1.00*    HN        1.00*    
!   KK        1.00*    KP        1.02     LA        0.97     MA        1.00*    
!   MC        1.00*    MK        1.01     NL        1.00*    ON        1.00*    
!   OV        0.99     PT        1.04     SC        0.95     TC        1.07     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00*    FD        1.01     GG        1.00*    HN        1.00*    
!   KK        0.95     KP        1.00*    LA        1.02     MA        1.00*    
!   MC        1.00*    MK        0.99     NL        1.00*    ON        1.00*    
!   OV        0.97     PT        1.08     SC        1.00*    TC        0.94     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
! 
! Correcting IF 3.
!  A total of 90 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 3 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00*    FD        0.94     GG        1.00*    HN        1.00*    
!   KK        1.00*    KP        0.94     LA        0.97     MA        1.00*    
!   MC        1.00*    MK        0.98     NL        1.00*    ON        1.00*    
!   OV        1.01     PT        1.05     SC        0.96     TC        1.70     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00*    FD        0.98     GG        1.00*    HN        1.00*    
!   KK        1.00     KP        1.00*    LA        1.01     MA        1.00*    
!   MC        1.00*    MK        0.98     NL        1.00*    ON        1.00*    
!   OV        1.00     PT        1.09     SC        1.00*    TC        0.93     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
! 
! Correcting IF 4.
!  A total of 90 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
!  A total of 3 telescope corrections were flagged in sub-array 2.
!  Telescope amplitude corrections in sub-array 2:
!   BR        1.00*    FD        0.92     GG        1.00*    HN        1.00*    
!   KK        1.00*    KP        1.00     LA        0.95     MA        1.00*    
!   MC        1.00*    MK        1.01     NL        1.00*    ON        1.00*    
!   OV        1.00     PT        1.06     SC        1.00     TC        1.07     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 57 telescope corrections were flagged in sub-array 3.
!  Telescope amplitude corrections in sub-array 3:
!   BR        1.00*    FD        0.96     GG        1.00*    HN        1.00*    
!   KK        0.95     KP        1.00*    LA        0.99     MA        1.00*    
!   MC        1.00*    MK        1.02     NL        1.00*    ON        1.00*    
!   OV        0.98     PT        1.09     SC        1.00*    TC        0.97     
!   TS        1.00*    WF        1.00*    WZ        1.00*    
! 
!  A total of 36 telescope corrections were flagged in sub-array 4.
!  Telescope amplitude corrections in sub-array 4:
!   BR        0.00*    FD        0.00*    GG        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*    ON        0.00*    
!   OV        0.00*    PT        0.00*    SC        0.00*    TC        0.00*    
!   TS        0.00*    WF        0.00*    WZ        0.00*    
! 
! 
! Fit before self-cal, rms=0.376974Jy  sigma=0.001507
! Fit after  self-cal, rms=0.357911Jy  sigma=0.001227
dynam = dynam - 0.5
flux_cutoff = imstat(rms) * dynam
! Inverting map and beam 
! Estimated beam: bmin=0.7868 mas, bmaj=1.182 mas, bpa=5.648 degrees
! Estimated noise=1260.46 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.00713109 Jy
! Component: 100  -  total flux cleaned = 0.00771728 Jy
! Component: 150  -  total flux cleaned = 0.00824238 Jy
! Component: 200  -  total flux cleaned = 0.00870961 Jy
! Total flux subtracted in 200 components = 0.00870961 Jy
! Clean residual min=-0.015082 max=0.016779 Jy/beam
! Clean residual mean=0.000682 rms=0.003776 Jy/beam
! Combined flux in latest and established models = 2.17148 Jy
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 24 components and 2.17148 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.357248Jy  sigma=0.001221
! Fit after  self-cal, rms=0.357138Jy  sigma=0.001221
! 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.357138Jy  sigma=0.001221
! Fit after  self-cal, rms=0.356106Jy  sigma=0.001219
dynam = dynam - 0.75
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.7853 mas, bmaj=1.181 mas, bpa=5.671 degrees
! Estimated noise=1261.05 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00586585 Jy
! Component: 100  -  total flux cleaned = 0.00662864 Jy
! Component: 150  -  total flux cleaned = 0.00706775 Jy
! Component: 200  -  total flux cleaned = 0.00747462 Jy
! Total flux subtracted in 200 components = 0.00747462 Jy
! Clean residual min=-0.014080 max=0.015706 Jy/beam
! Clean residual mean=0.000658 rms=0.003527 Jy/beam
! Combined flux in latest and established models = 2.17896 Jy
selfcal
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 25 components and 2.17896 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.355677Jy  sigma=0.001215
! Fit after  self-cal, rms=0.355643Jy  sigma=0.001215
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 (1.92, 12.16).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0150913 Jy
! Component: 100  -  total flux cleaned = 0.0183115 Jy
! Component: 150  -  total flux cleaned = 0.0202204 Jy
! Component: 200  -  total flux cleaned = 0.0222847 Jy
! Total flux subtracted in 200 components = 0.0222847 Jy
! Clean residual min=-0.011605 max=0.009976 Jy/beam
! Clean residual mean=0.000209 rms=0.002612 Jy/beam
! Combined flux in latest and established models = 2.20124 Jy
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 43 components and 2.20124 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.354422Jy  sigma=0.001205
! Fit after  self-cal, rms=0.354398Jy  sigma=0.001205
! 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.354398Jy  sigma=0.001205
! Fit after  self-cal, rms=0.353834Jy  sigma=0.001203
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=0.7838 mas, bmaj=1.182 mas, bpa=5.602 degrees
! Estimated noise=1260.38 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00589656 Jy
! Component: 100  -  total flux cleaned = 0.00818235 Jy
! Component: 150  -  total flux cleaned = 0.00918652 Jy
! Component: 200  -  total flux cleaned = 0.0099133 Jy
! Total flux subtracted in 200 components = 0.0099133 Jy
! Clean residual min=-0.011445 max=0.008926 Jy/beam
! Clean residual mean=0.000110 rms=0.002289 Jy/beam
! Combined flux in latest and established models = 2.21116 Jy
selfcal
! Performing phase self-cal
! Adding 20 model components to the UV plane model.
! The established model now contains 49 components and 2.21116 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.353384Jy  sigma=0.001200
! Fit after  self-cal, rms=0.353372Jy  sigma=0.001200
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.000766324 Jy
! Component: 100  -  total flux cleaned = 0.00124202 Jy
! Component: 150  -  total flux cleaned = 0.00139185 Jy
! Component: 200  -  total flux cleaned = 0.00139809 Jy
! Total flux subtracted in 200 components = 0.00139809 Jy
! Clean residual min=-0.011084 max=0.008113 Jy/beam
! Clean residual mean=0.000082 rms=0.002210 Jy/beam
! Combined flux in latest and established models = 2.21256 Jy
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 52 components and 2.21256 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.353172Jy  sigma=0.001198
! Fit after  self-cal, rms=0.353167Jy  sigma=0.001198
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 90 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.353167Jy  sigma=0.001198
! Fit after  self-cal, rms=0.353231Jy  sigma=0.001197
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.977 mas, bmaj=5.173 mas, bpa=-2.529 degrees
! Estimated noise=1431.8 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00413467 Jy
! Component: 100  -  total flux cleaned = 0.00454518 Jy
! Component: 150  -  total flux cleaned = 0.00460145 Jy
! Component: 200  -  total flux cleaned = 0.00460221 Jy
! Total flux subtracted in 200 components = 0.00460221 Jy
! Clean residual min=-0.005954 max=0.006459 Jy/beam
! Clean residual mean=0.000124 rms=0.002395 Jy/beam
! Combined flux in latest and established models = 2.21716 Jy
selfcal
! Performing phase self-cal
! Adding 7 model components to the UV plane model.
! The established model now contains 54 components and 2.21716 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.353148Jy  sigma=0.001197
! Fit after  self-cal, rms=0.353133Jy  sigma=0.001197
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 = 4.47614e-06 Jy
! Component: 100  -  total flux cleaned = 0.000130034 Jy
! Component: 150  -  total flux cleaned = 0.000129407 Jy
! Component: 200  -  total flux cleaned = 0.000126276 Jy
! Total flux subtracted in 200 components = 0.000126276 Jy
! Clean residual min=-0.006014 max=0.006539 Jy/beam
! Clean residual mean=0.000120 rms=0.002385 Jy/beam
! Combined flux in latest and established models = 2.21728 Jy
! Performing phase self-cal
! Adding 9 model components to the UV plane model.
! The established model now contains 56 components and 2.21728 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.353170Jy  sigma=0.001196
! Fit after  self-cal, rms=0.353162Jy  sigma=0.001196
! 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=0.7827 mas, bmaj=1.182 mas, bpa=5.537 degrees
! Estimated noise=1260.76 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.7827 x 1.182 at 5.537 degrees (North through East)
! Clean map  min=-0.0085229  max=1.9877 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.353162Jy  sigma=0.001196
! Fit after  self-cal, rms=0.341678Jy  sigma=0.001047
! Inverting map and beam 
! Estimated beam: bmin=0.7822 mas, bmaj=1.181 mas, bpa=5.582 degrees
! Estimated noise=1261.22 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00311826 Jy
! Component: 100  -  total flux cleaned = 0.00328083 Jy
! Component: 150  -  total flux cleaned = 0.00342816 Jy
! Component: 200  -  total flux cleaned = 0.00356709 Jy
! Total flux subtracted in 200 components = 0.00356709 Jy
! Clean residual min=-0.009631 max=0.006817 Jy/beam
! Clean residual mean=0.000016 rms=0.001942 Jy/beam
! Combined flux in latest and established models = 2.22085 Jy
! Performing phase self-cal
! Adding 18 model components to the UV plane model.
! The established model now contains 65 components and 2.22085 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.341436Jy  sigma=0.001046
! Fit after  self-cal, rms=0.341432Jy  sigma=0.001046
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000116846 Jy
! Component: 100  -  total flux cleaned = 0.000235551 Jy
! Component: 150  -  total flux cleaned = 0.000237472 Jy
! Component: 200  -  total flux cleaned = 0.000344717 Jy
! Total flux subtracted in 200 components = 0.000344717 Jy
! Clean residual min=-0.009410 max=0.006441 Jy/beam
! Clean residual mean=0.000007 rms=0.001902 Jy/beam
! Combined flux in latest and established models = 2.22119 Jy
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 65 components and 2.2212 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.341294Jy  sigma=0.001045
! Fit after  self-cal, rms=0.341292Jy  sigma=0.001045
! 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.341292Jy  sigma=0.001045
! Fit after  self-cal, rms=0.341883Jy  sigma=0.001045
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=0.617 mas, bmaj=0.9316 mas, bpa=9.563 degrees
! Estimated noise=2443.43 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.5424 Jy
! Component: 100  -  total flux cleaned = 1.91836 Jy
! Component: 150  -  total flux cleaned = 2.05546 Jy
! Component: 200  -  total flux cleaned = 2.11818 Jy
! Total flux subtracted in 200 components = 2.11818 Jy
! Clean residual min=-0.023081 max=0.032870 Jy/beam
! Clean residual mean=0.000740 rms=0.006273 Jy/beam
! Combined flux in latest and established models = 2.11818 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 16 model components to the UV plane model.
! The established model now contains 16 components and 2.11818 Jy
! Inverting map and beam 
! Estimated beam: bmin=0.7815 mas, bmaj=1.182 mas, bpa=5.503 degrees
! Estimated noise=1261.51 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)
! Added new window around map position (0.96, 2.88).
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.045885 Jy
! Component: 100  -  total flux cleaned = 0.067163 Jy
! Component: 150  -  total flux cleaned = 0.0804671 Jy
! Component: 200  -  total flux cleaned = 0.0871113 Jy
! Total flux subtracted in 200 components = 0.0871113 Jy
! Clean residual min=-0.008410 max=0.008922 Jy/beam
! Clean residual mean=0.000391 rms=0.002268 Jy/beam
! Combined flux in latest and established models = 2.20529 Jy
! Adding 36 model components to the UV plane model.
! The established model now contains 46 components and 2.20529 Jy
! Inverting map 
uvtaper 0.3,taper_size
! Gaussian taper: value 0.3 at UV radius = 90 mega-wavelengths.
clean clean_niter,clean_gain
! Inverting map and beam 
! Estimated beam: bmin=1.974 mas, bmaj=5.166 mas, bpa=-2.532 degrees
! Estimated noise=1435.33 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00728344 Jy
! Component: 100  -  total flux cleaned = 0.0113459 Jy
! Component: 150  -  total flux cleaned = 0.0127732 Jy
! Component: 200  -  total flux cleaned = 0.0133991 Jy
! Total flux subtracted in 200 components = 0.0133991 Jy
! Clean residual min=-0.005192 max=0.005314 Jy/beam
! Clean residual mean=0.000232 rms=0.002203 Jy/beam
! Combined flux in latest and established models = 2.21869 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.000618153 Jy
! Component: 100  -  total flux cleaned = 0.00108904 Jy
! Component: 150  -  total flux cleaned = 0.00156032 Jy
! Component: 200  -  total flux cleaned = 0.00197555 Jy
! Total flux subtracted in 200 components = 0.00197555 Jy
! Clean residual min=-0.005082 max=0.005353 Jy/beam
! Clean residual mean=0.000155 rms=0.002106 Jy/beam
! Combined flux in latest and established models = 2.22066 Jy
! Adding 22 model components to the UV plane model.
! The established model now contains 64 components and 2.22066 Jy
! Inverting map 
uvtaper 0,0
! No UV-taper is currently set.
clean
! Inverting map and beam 
! Estimated beam: bmin=0.7815 mas, bmaj=1.182 mas, bpa=5.503 degrees
! Estimated noise=1261.51 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00113281 Jy
! Component: 100  -  total flux cleaned = -0.00115431 Jy
! Component: 150  -  total flux cleaned = -0.000575121 Jy
! Component: 200  -  total flux cleaned = -4.84928e-05 Jy
! Total flux subtracted in 200 components = -4.84928e-05 Jy
! Clean residual min=-0.006890 max=0.005952 Jy/beam
! Clean residual mean=0.000080 rms=0.001788 Jy/beam
! Combined flux in latest and established models = 2.22061 Jy
selfcal
! Performing phase self-cal
! Adding 17 model components to the UV plane model.
! The established model now contains 76 components and 2.22061 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.342439Jy  sigma=0.001042
! Fit after  self-cal, rms=0.342317Jy  sigma=0.001041
wmodel J0403-3605_X_2003_05_07_pus_map.mod
! Writing 76 model components to file: J0403-3605_X_2003_05_07_pus_map.mod
wobs J0403-3605_X_2003_05_07_pus_uvs.fits
! Writing UV FITS file: J0403-3605_X_2003_05_07_pus_uvs.fits
wwins J0403-3605_X_2003_05_07_pus_map.win
! wwins: Wrote 3 windows to J0403-3605_X_2003_05_07_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.00112377 Jy
! Component: 100  -  total flux cleaned = -0.00166187 Jy
! Component: 150  -  total flux cleaned = -0.00156407 Jy
! Component: 200  -  total flux cleaned = -0.00137056 Jy
! Component: 250  -  total flux cleaned = -0.00100441 Jy
! Component: 300  -  total flux cleaned = -0.000564978 Jy
! Component: 350  -  total flux cleaned = -5.4384e-05 Jy
! Component: 400  -  total flux cleaned = 0.000194404 Jy
! Component: 450  -  total flux cleaned = 0.000597307 Jy
! Component: 500  -  total flux cleaned = 0.000675263 Jy
! Component: 550  -  total flux cleaned = 0.000982299 Jy
! Component: 600  -  total flux cleaned = 0.00120851 Jy
! Component: 650  -  total flux cleaned = 0.00165097 Jy
! Component: 700  -  total flux cleaned = 0.00193906 Jy
! Component: 750  -  total flux cleaned = 0.00229402 Jy
! Component: 800  -  total flux cleaned = 0.00271134 Jy
! Component: 850  -  total flux cleaned = 0.00312228 Jy
! Component: 900  -  total flux cleaned = 0.00346055 Jy
! Component: 950  -  total flux cleaned = 0.00386014 Jy
! Component: 1000  -  total flux cleaned = 0.00431928 Jy
! Component: 1050  -  total flux cleaned = 0.00464128 Jy
! Component: 1100  -  total flux cleaned = 0.00489551 Jy
! Component: 1150  -  total flux cleaned = 0.00502168 Jy
! Component: 1200  -  total flux cleaned = 0.00520769 Jy
! Component: 1250  -  total flux cleaned = 0.00563586 Jy
! Component: 1300  -  total flux cleaned = 0.00575565 Jy
! Component: 1350  -  total flux cleaned = 0.00599482 Jy
! Component: 1400  -  total flux cleaned = 0.00623188 Jy
! Component: 1450  -  total flux cleaned = 0.00675849 Jy
! Component: 1500  -  total flux cleaned = 0.00693269 Jy
! Component: 1550  -  total flux cleaned = 0.00744845 Jy
! Component: 1600  -  total flux cleaned = 0.00761862 Jy
! Component: 1650  -  total flux cleaned = 0.00790047 Jy
! Component: 1700  -  total flux cleaned = 0.00829219 Jy
! Component: 1750  -  total flux cleaned = 0.00873552 Jy
! Component: 1800  -  total flux cleaned = 0.0090652 Jy
! Component: 1850  -  total flux cleaned = 0.00928328 Jy
! Component: 1900  -  total flux cleaned = 0.00960764 Jy
! Component: 1950  -  total flux cleaned = 0.00982221 Jy
! Component: 2000  -  total flux cleaned = 0.00992853 Jy
! Total flux subtracted in 2048 components = 0.0101938 Jy
! Clean residual min=-0.003153 max=0.003542 Jy/beam
! Clean residual mean=0.000183 rms=0.000935 Jy/beam
! Combined flux in latest and established models = 2.23081 Jy
keep
! Adding 118 model components to the UV plane model.
! The established model now contains 194 components and 2.23081 Jy
wmap "J0403-3605_X_2003_05_07_pus\_map.fits"
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.7815 x 1.182 at 5.503 degrees (North through East)
! Clean map  min=-0.0072789  max=1.979 Jy/beam
! Writing clean map to FITS file: J0403-3605_X_2003_05_07_pus_map.fits
float peak_flux
peak_flux = peak(flux,max)
print "The peak flux is",peak_flux
! The peak flux is 1.97905 
invert
! Inverting map 
float image_rms
image_rms = imstat(rms)
print "The final image rms is",image_rms
! The final image rms is 0.000940362 
device "J0403-3605_X_2003_05_07_pus\_map.eps/VPS"
! Attempting to open device: 'J0403-3605_X_2003_05_07_pus_map.eps/VPS'
loglevs 100*3*image_rms/peak_flux
! The new contour levels are:
!  -0.142548 0.142548 0.285096 0.570191 1.14038 2.28077 4.56153 9.12306 18.2461 36.4922 72.9845
print "clev=",int(3*image_rms*10000)/10000
! clev= 0.0028 
print "peak=",int(1000*peak_flux)/1000
! peak= 1.979 
mapcolor none
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
mapplot cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.7815 x 1.182 at 5.503 degrees (North through East)
! Clean map  min=-0.0072789  max=1.979 Jy/beam
mapcolor color, 0.9
! Mapplot colormap: rainbow, contrast: 0.9 brightness: 0.5.
device "J0403-3605_X_2003_05_07_pus\_map.gif/vgif"
! Attempting to open device: 'J0403-3605_X_2003_05_07_pus_map.gif/vgif'
mapplot cln
quit
! Quitting program
! Log file J0403-3605_X_2003_05_07_pus_uvs.log closed on Tue Dec 12 14:01:26 2006