! Started logfile: difmap.log on Wed Apr 30 13:28:57 2008
observe 1333-337_S.edt
! Reading UV FITS file: 1333-337_S.edt
! AN table 1: 69 integrations on 136 of 136 possible baselines.
! Apparent sampling: 0.125639 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 1333-337
! 
! 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.23697e+09        8e+06          1        8e+06
!  02        2  2.26697e+09        8e+06          1        8e+06
!  03        3  2.35697e+09        8e+06          1        8e+06
!  04        4  2.37697e+09        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1141 lines of history.
! 
! Reading 4716 visibilities.
select rr
! 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
uvweight 2,-1
! Uniform weighting binwidth: 2 (pixels).
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
mapsize 2048,0.3
! Map grid = 2048x2048 pixels with 0.300x0.300 milli-arcsec cellsize.
peakwin 0.3
! Inverting map and beam 
! Estimated beam: bmin=2.996 mas, bmaj=15.28 mas, bpa=-17.08 degrees
! Estimated noise=2720.64 mJy/beam.
! Added new window around map position (0, 0).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.140574 Jy
! Component: 100  -  total flux cleaned = 0.154354 Jy
! Total flux subtracted in 100 components = 0.154354 Jy
! Clean residual min=-0.012179 max=0.014909 Jy/beam
! Clean residual mean=-0.000050 rms=0.003101 Jy/beam
! Combined flux in latest and established models = 0.154354 Jy
selfcal
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 12 components and 0.154354 Jy
! 
! Correcting IF 1.
!  A total of 802 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 805 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 753 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 863 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.107620Jy  sigma=0.001093
! Fit after  self-cal, rms=0.107283Jy  sigma=0.001092
clean
! Inverting map and beam 
! Estimated beam: bmin=3.053 mas, bmaj=15.32 mas, bpa=-17.3 degrees
! Estimated noise=2662.33 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.00214952 Jy
! Component: 100  -  total flux cleaned = 0.00214183 Jy
! Total flux subtracted in 100 components = 0.00214183 Jy
! Clean residual min=-0.012029 max=0.015308 Jy/beam
! Clean residual mean=-0.000046 rms=0.003108 Jy/beam
! Combined flux in latest and established models = 0.156496 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 17 components and 0.156496 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.107036Jy  sigma=0.001090
! Fit after  self-cal, rms=0.107033Jy  sigma=0.001090
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 8.53517e-05 Jy
! Component: 100  -  total flux cleaned = 7.72786e-05 Jy
! Total flux subtracted in 100 components = 7.72786e-05 Jy
! Clean residual min=-0.012019 max=0.015302 Jy/beam
! Clean residual mean=-0.000045 rms=0.003109 Jy/beam
! Combined flux in latest and established models = 0.156573 Jy
selfcal
! Performing phase self-cal
! Adding 6 model components to the UV plane model.
! The established model now contains 18 components and 0.156573 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.107027Jy  sigma=0.001090
! Fit after  self-cal, rms=0.107023Jy  sigma=0.001090
mapl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
uvweight 0,-2
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -2.
! Radial weighting is not currently selected.
clean
! Inverting map and beam 
! Estimated beam: bmin=5.586 mas, bmaj=22.71 mas, bpa=-15.86 degrees
! Estimated noise=1064.87 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0079489 Jy
! Component: 100  -  total flux cleaned = 0.00817018 Jy
! Total flux subtracted in 100 components = 0.00817018 Jy
! Clean residual min=-0.010529 max=0.014645 Jy/beam
! Clean residual mean=-0.000145 rms=0.003069 Jy/beam
! Combined flux in latest and established models = 0.164743 Jy
selfcal
! Performing phase self-cal
! Adding 3 model components to the UV plane model.
! The established model now contains 18 components and 0.164744 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.106406Jy  sigma=0.001086
! Fit after  self-cal, rms=0.106407Jy  sigma=0.001086
clean
! Inverting map 
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 7.49261e-05 Jy
! Component: 100  -  total flux cleaned = 8.57514e-05 Jy
! Total flux subtracted in 100 components = 8.57514e-05 Jy
! Clean residual min=-0.010510 max=0.014611 Jy/beam
! Clean residual mean=-0.000145 rms=0.003068 Jy/beam
! Combined flux in latest and established models = 0.164829 Jy
selfcal
! Performing phase self-cal
! Adding 3 model components to the UV plane model.
! The established model now contains 18 components and 0.164829 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.106390Jy  sigma=0.001086
! Fit after  self-cal, rms=0.106390Jy  sigma=0.001086
mapl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00     HN        1.00*    KB        1.00*    
!   KK        1.00*    KP        0.82     LA        0.95     MK        1.00*    
!   NL        1.02     NY        1.00*    OV        1.18     PT        0.99     
!   SC        1.00*    TS        1.00*    WF        1.00*    WZ        1.00*    
!   ZC        1.00*    
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        0.99     HN        1.00*    KB        1.00*    
!   KK        1.00*    KP        0.92     LA        0.87     MK        1.00*    
!   NL        1.08     NY        1.00*    OV        1.23     PT        0.96     
!   SC        1.00*    TS        1.00*    WF        1.00*    WZ        1.00*    
!   ZC        1.00*    
! 
! 
! Correcting IF 3.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.01     HN        1.49     KB        1.00*    
!   KK        1.00*    KP        0.83     LA        0.93     MK        1.17     
!   NL        1.05     NY        1.00*    OV        1.21     PT        0.98     
!   SC        0.73     TS        1.00*    WF        1.00*    WZ        1.00*    
!   ZC        1.00*    
! 
! 
! Correcting IF 4.
!  A total of 3 telescope corrections were flagged in sub-array 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00*    FD        1.00*    HN        0.75     KB        1.00*    
!   KK        1.00*    KP        0.79     LA        0.92     MK        0.77     
!   NL        1.07     NY        1.00*    OV        1.19     PT        0.98     
!   SC        0.81     TS        1.00*    WF        1.00*    WZ        1.00*    
!   ZC        1.00*    
! 
! 
! Fit before self-cal, rms=0.106390Jy  sigma=0.001086
! Fit after  self-cal, rms=0.102253Jy  sigma=0.001052
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 6 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=6.367 mas, bmaj=25.76 mas, bpa=-15.49 degrees
! Estimated noise=1060.49 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -0.000234448 Jy
! Component: 100  -  total flux cleaned = -0.000233705 Jy
! Total flux subtracted in 100 components = -0.000233705 Jy
! Clean residual min=-0.005228 max=0.004601 Jy/beam
! Clean residual mean=-0.000060 rms=0.001316 Jy/beam
! Combined flux in latest and established models = 0.164596 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 23 components and 0.164596 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101417Jy  sigma=0.001049
! Fit after  self-cal, rms=0.101239Jy  sigma=0.001046
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 6.367 x 25.76 at -15.49 degrees (North through East)
! Clean map  min=-0.0051948  max=0.16266 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 0.3
! Inverting map 
! Added new window around map position (7.8, -18.3).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.000149094 Jy
! Component: 100  -  total flux cleaned = 0.000143278 Jy
! Total flux subtracted in 100 components = 0.000143278 Jy
! Clean residual min=-0.005226 max=0.004837 Jy/beam
! Clean residual mean=-0.000064 rms=0.001330 Jy/beam
! Combined flux in latest and established models = 0.164739 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 24 components and 0.164739 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101237Jy  sigma=0.001046
! Fit after  self-cal, rms=0.101236Jy  sigma=0.001046
selfcal true,true,120
! 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.101236Jy  sigma=0.001046
! Fit after  self-cal, rms=0.101941Jy  sigma=0.001045
clean
! Inverting map and beam 
! Estimated beam: bmin=6.448 mas, bmaj=25.95 mas, bpa=-15.65 degrees
! Estimated noise=1064.79 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.000673202 Jy
! Component: 100  -  total flux cleaned = 0.000669999 Jy
! Total flux subtracted in 100 components = 0.000669999 Jy
! Clean residual min=-0.005089 max=0.004406 Jy/beam
! Clean residual mean=-0.000057 rms=0.001283 Jy/beam
! Combined flux in latest and established models = 0.165409 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 25 components and 0.165409 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101935Jy  sigma=0.001045
! Fit after  self-cal, rms=0.101952Jy  sigma=0.001044
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
peakwin 0.3
! Inverting map 
! Added new window around map position (-71.1, 127.8).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = -4.90733e-06 Jy
! Component: 100  -  total flux cleaned = -1.15456e-05 Jy
! Total flux subtracted in 100 components = -1.15456e-05 Jy
! Clean residual min=-0.005097 max=0.004403 Jy/beam
! Clean residual mean=-0.000057 rms=0.001283 Jy/beam
! Combined flux in latest and established models = 0.165397 Jy
selfcal
! Performing phase self-cal
! Adding 5 model components to the UV plane model.
! The established model now contains 25 components and 0.165397 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101950Jy  sigma=0.001044
! Fit after  self-cal, rms=0.101949Jy  sigma=0.001044
mapl cln
! Inverting map 
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 6.448 x 25.95 at -15.65 degrees (North through East)
! Clean map  min=-0.0050708  max=0.16346 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
loglevs 2
! The new contour levels are:
!  -2 2 4 8 16 32 64
mapl cln
! 
! Move the cursor into the plot window and press 'H' for help
loglevs 3
! The new contour levels are:
!  -3 3 6 12 24 48 96
mapl cln
! 
! Move the cursor into the plot window and press 'H' for help
selfcal true,true,1
! Performing amp+phase self-cal over 1 minute time intervals
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.101949Jy  sigma=0.001044
! Fit after  self-cal, rms=0.102501Jy  sigma=0.001020
loglevs 2
! The new contour levels are:
!  -2 2 4 8 16 32 64
mapl cln
! Inverting map and beam 
! Estimated beam: bmin=6.374 mas, bmaj=25.69 mas, bpa=-16.68 degrees
! Estimated noise=1065.4 mJy/beam.
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 6.374 x 25.69 at -16.68 degrees (North through East)
! Clean map  min=-0.0042666  max=0.16469 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
radpl
! Using default options string "m1"
! Move the cursor into the plot window and press 'H' for help
! Applying 14 buffered edits.
clean
! Inverting map and beam 
! Estimated beam: bmin=6.395 mas, bmaj=25.69 mas, bpa=-16.7 degrees
! Estimated noise=1066.17 mJy/beam.
! clean: niter=100  gain=0.05  cutoff=0
! Component: 050  -  total flux cleaned = 0.0012161 Jy
! Component: 100  -  total flux cleaned = 0.0012874 Jy
! Total flux subtracted in 100 components = 0.0012874 Jy
! Clean residual min=-0.004099 max=0.003763 Jy/beam
! Clean residual mean=-0.000047 rms=0.001100 Jy/beam
! Combined flux in latest and established models = 0.166685 Jy
selfcal
! Performing phase self-cal
! Adding 4 model components to the UV plane model.
! The established model now contains 27 components and 0.166685 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.100284Jy  sigma=0.001014
! Fit after  self-cal, rms=0.100222Jy  sigma=0.001013
peakwin 0.4
! Inverting map 
! Added new window around map position (10.5, -24.3).
mapl
! 
! Move the cursor into the plot window and press 'H' for help
mapl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 6.395 x 25.69 at -16.7 degrees (North through East)
! Clean map  min=-0.0040278  max=0.16452 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
save 1333-337_S
! Writing UV FITS file: 1333-337_S.uvf
! Writing 27 model components to file: 1333-337_S.mod
! wwins: Wrote 1 windows to 1333-337_S.win
! Writing clean map to FITS file: 1333-337_S.fits
! Writing difmap environment to: 1333-337_S.par
![@gonzo_plot 1333-337_S]
device 1333-337_Suv.ps/vps,1,2
! Attempting to open device: '1333-337_Suv.ps/vps'
radplot
! Using default options string "m1"
uvplot
! Using default options string ""
device 1333-337_Sim.ps/vps
! Attempting to open device: '1333-337_Sim.ps/vps'
uvweight 0,-1
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
loglevs -1,2048,2
! The new contour levels are:
!  -1 1 2 4 8 16 32 64 128 256 512 1024 2048
cmul = imstat(rms)*3
! Inverting map and beam 
! Estimated beam: bmin=5.94 mas, bmaj=24.89 mas, bpa=-16.62 degrees
! Estimated noise=1092.83 mJy/beam.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 5.94 x 24.89 at -16.62 degrees (North through East)
! Clean map  min=-0.004303  max=0.16366 Jy/beam
![Exited script file: gonzo_plot]
quit
! Quitting program
! Log file difmap.log closed on Wed Apr 30 13:41:12 2008