! Started logfile: difmap.log on Mon Aug 12 10:08:06 2002
mapunits mas
! Distances in the map plane now have units: milli-arcsec.
! Distances in the UV plane now have units: mega-wavelengths.
multi_model false
! Use one model for all channel/polarization selections.
get 0048-097_Q
! Reading UV FITS file: 0048-097_Q.uvf
! AN table 1: 87 integrations on 45 of 45 possible baselines.
! Apparent sampling: 0.820434 visibilities/baseline/integration-bin.
! Found source: 0048-097
! 
! 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  4.29545e+10        8e+06          1        8e+06
!  02        2  4.30345e+10        8e+06          1        8e+06
!  03        3  4.32145e+10        8e+06          1        8e+06
!  04        4  4.33545e+10        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 1889 lines of history.
! 
! Reading 12848 visibilities.
! A total of 119 model components were read from file 0048-097_Q.mod
! rwins: Read 4 windows from 0048-097_Q.win
! Multi-model file "0048-097_Q.mtab" not available
mapsize 1024,0.02, 1024,0.02
! Map grid = 1024x1024 pixels with 0.0200x0.0200 milli-arcsec cellsize.
select RR, 1, 4
! 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
uvtaper 0, 0
! No UV-taper is currently set.
uvrange 0, 0
! The full UV range of the data is currently selected for gridding.
uvweight 0, -1, false
! Uniform weighting is not currently selected.
! Gridding weights will be scaled by errors raised to the power -1.
! Radial weighting is not currently selected.
selftaper 0, 0
! No selfcal UV-taper is currently set.
selflims 0, 0
! Selfcal amplitude corrections will not be limited.
! Selfcal phase corrections will not be limited.
xyrange 0, 0, 0, 0
! The mapplot display area will be limited to:
!  Relative  RA: (The whole available range)
!  Relative DEC: (The whole available range)
beamloc 0, 0, 0.01, 0.3
! Mapplot clean beam ellipse center: 0,0. Size range: 0.01 -> 0.3.
polvec 0, 0, 0, 1, 1
! Give polarization vectors lengths of 0 milli-arcsec/Jy.
! Draw vectors where unpolarized flux > 0 Jy and polarized flux > 0 Jy.
! Draw polarization vectors in every X pixel and in every Y pixel.
integer niter
niter=200
float gain
gain=0.03
float cutoff
cutoff=0
float cmul
cmul=0
logical docont
docont=true
mapcolor grey, 1.000000, 0.500000
! Mapplot colormap: grey, contrast: 1 brightness: 0.5.
mapfunc linear, 0.000000, 0.000000
! Mapplot transfer-function = linear, Data range = data min -> data max.
string vflags
vflags=""
string rflags
rflags=""
string pflags
pflags=""
string tflags
tflags=""
selfflag true, 3, 4
! - In phase-only self-cal, good data on baselines that are not in closed
!   arrays of at least 3 telescopes will be flagged.
! - In amplitude self-cal, good data on baselines that are not in closed
!   arrays of at least 4 telescopes will be flagged.
scangap 3600
! Delimiting interscan gap changed to 3600 seconds in all sub-arrays.
float levs(1)
levs(1:1) = 0
mappl
! Adding 119 model components to the UV plane model.
! The established model now contains 119 components and 0.876581 Jy
! Inverting map and beam 
! Estimated beam: bmin=0.1522 mas, bmaj=0.3878 mas, bpa=-1.836 degrees
! Estimated noise=1.14114 mJy/beam.
! 
! Move the cursor into the plot window and press 'H' for help
clrmod true,true,true
! clrmod: Cleared the established, tentative and continuum models.
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
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.
uncal true,true,true
! uncal: All telescope phase corrections have been un-done.
! uncal: All telescope amplitude corrections have been un-done.
! uncal: All telescope correction flags have been un-done.
mappl
! Inverting map and beam 
! Estimated beam: bmin=0.1353 mas, bmaj=0.3547 mas, bpa=-1.033 degrees
! Estimated noise=2.13053 mJy/beam.
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.625159 Jy
! Component: 100  -  total flux cleaned = 0.761731 Jy
! Component: 150  -  total flux cleaned = 0.800878 Jy
! Component: 200  -  total flux cleaned = 0.81983 Jy
! Total flux subtracted in 200 components = 0.81983 Jy
! Clean residual min=-0.017273 max=0.024567 Jy/beam
! Clean residual mean=0.000004 rms=0.004660 Jy/beam
! Combined flux in latest and established models = 0.81983 Jy
clrmod
! clrmod: Cleared the tentative model.
clean 100,0.03
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.625159 Jy
! Component: 100  -  total flux cleaned = 0.761731 Jy
! Total flux subtracted in 100 components = 0.761731 Jy
! Clean residual min=-0.021572 max=0.041172 Jy/beam
! Clean residual mean=-0.000007 rms=0.005768 Jy/beam
! Combined flux in latest and established models = 0.761731 Jy
selfcal false,false,0.5
! Performing phase self-cal over 0.5 minute time intervals
! Adding 4 model components to the UV plane model.
! The established model now contains 4 components and 0.761731 Jy
! 
! Correcting IF 1.
!  A total of 163 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 2.
!  A total of 99 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 3.
!  A total of 100 telescope corrections were flagged in sub-array 1.
! 
! Correcting IF 4.
!  A total of 108 telescope corrections were flagged in sub-array 1.
! 
! Fit before self-cal, rms=0.282423Jy  sigma=1.568113
! Fit after  self-cal, rms=0.282255Jy  sigma=1.566866
clean
! Inverting map and beam 
! Estimated beam: bmin=0.1353 mas, bmaj=0.3547 mas, bpa=-1.033 degrees
! Estimated noise=2.13053 mJy/beam.
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0378943 Jy
! Component: 100  -  total flux cleaned = 0.0567863 Jy
! Total flux subtracted in 100 components = 0.0567863 Jy
! Clean residual min=-0.017174 max=0.025377 Jy/beam
! Clean residual mean=0.000004 rms=0.004742 Jy/beam
! Combined flux in latest and established models = 0.818517 Jy
selfcal
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 17 components and 0.818517 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.277433Jy  sigma=1.537519
! Fit after  self-cal, rms=0.277654Jy  sigma=1.536057
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 1.5
! Added new window around map position (-0.28, 0.54).
mappl
! 
! Move the cursor into the plot window and press 'H' for help
! Mean=0.00663957  rms=0.00710841  min=0.00169308  max=0.0134651 Jy/beam
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0223539 Jy
! Component: 100  -  total flux cleaned = 0.0335782 Jy
! Total flux subtracted in 100 components = 0.0335782 Jy
! Clean residual min=-0.016532 max=0.017690 Jy/beam
! Clean residual mean=-0.000001 rms=0.004268 Jy/beam
! Combined flux in latest and established models = 0.852095 Jy
selfcal
! Performing phase self-cal
! Adding 8 model components to the UV plane model.
! The established model now contains 22 components and 0.852095 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276478Jy  sigma=1.529226
! Fit after  self-cal, rms=0.276397Jy  sigma=1.528606
clean
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00864672 Jy
! Component: 100  -  total flux cleaned = 0.0118097 Jy
! Total flux subtracted in 100 components = 0.0118097 Jy
! Clean residual min=-0.016494 max=0.014719 Jy/beam
! Clean residual mean=-0.000001 rms=0.004100 Jy/beam
! Combined flux in latest and established models = 0.863905 Jy
selfcal
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 26 components and 0.863905 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276266Jy  sigma=1.527966
! Fit after  self-cal, rms=0.276242Jy  sigma=1.527764
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
peakwin 1.5
! Added new window around map position (-0.12, -3.02).
mappl
! 
! Move the cursor into the plot window and press 'H' for help
mappl
! 
! Move the cursor into the plot window and press 'H' for help
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.0133454 Jy
! Component: 100  -  total flux cleaned = 0.0170238 Jy
! Total flux subtracted in 100 components = 0.0170238 Jy
! Clean residual min=-0.016009 max=0.014166 Jy/beam
! Clean residual mean=0.000001 rms=0.003868 Jy/beam
! Combined flux in latest and established models = 0.880929 Jy
selfcal
! Performing phase self-cal
! Adding 15 model components to the UV plane model.
! The established model now contains 39 components and 0.880929 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276038Jy  sigma=1.526835
! Fit after  self-cal, rms=0.276003Jy  sigma=1.526573
clean
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00391713 Jy
! Component: 100  -  total flux cleaned = 0.00368127 Jy
! Total flux subtracted in 100 components = 0.00368127 Jy
! Clean residual min=-0.015927 max=0.014163 Jy/beam
! Clean residual mean=0.000001 rms=0.003781 Jy/beam
! Combined flux in latest and established models = 0.88461 Jy
selfcal
! Performing phase self-cal
! Adding 16 model components to the UV plane model.
! The established model now contains 48 components and 0.88461 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276024Jy  sigma=1.526755
! Fit after  self-cal, rms=0.276020Jy  sigma=1.526714
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
gscale true
! Performing overall amplitude self-cal
! 
! Correcting IF 1.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        1.00     HN        1.00     KP        1.01     
!   LA        1.00     MK        1.00     NL        1.00     OV        1.00     
!   PT        1.00     SC        1.00     
! 
! 
! Correcting IF 2.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        1.00     HN        0.99     KP        1.01     
!   LA        1.00     MK        1.00     NL        1.00     OV        1.00     
!   PT        1.01     SC        1.00     
! 
! 
! Correcting IF 3.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.01     FD        1.00     HN        0.99     KP        1.01     
!   LA        1.00     MK        1.00     NL        1.00     OV        1.00     
!   PT        1.01     SC        1.00     
! 
! 
! Correcting IF 4.
!  Telescope amplitude corrections in sub-array 1:
!   BR        1.00     FD        1.00     HN        0.99     KP        1.01     
!   LA        1.00     MK        1.00     NL        1.00     OV        1.00     
!   PT        1.01     SC        1.00     
! 
! 
! Fit before self-cal, rms=0.276020Jy  sigma=1.526714
! Fit after  self-cal, rms=0.276913Jy  sigma=1.526399
doamp false
! Syntax error in line:
! doamp false
!     ^
! Unable to identify keyword "doamp"
! 
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.
clean
! Inverting map and beam 
! Estimated beam: bmin=0.1521 mas, bmaj=0.3877 mas, bpa=-1.88 degrees
! Estimated noise=1.14566 mJy/beam.
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.00293691 Jy
! Component: 100  -  total flux cleaned = -0.00245424 Jy
! Total flux subtracted in 100 components = -0.00245424 Jy
! Clean residual min=-0.007705 max=0.007509 Jy/beam
! Clean residual mean=-0.000004 rms=0.001827 Jy/beam
! Combined flux in latest and established models = 0.882156 Jy
selfcal
! Performing phase self-cal
! Adding 11 model components to the UV plane model.
! The established model now contains 59 components and 0.882156 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276781Jy  sigma=1.525528
! Fit after  self-cal, rms=0.276783Jy  sigma=1.525492
clean
! Inverting map 
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000192934 Jy
! Component: 100  -  total flux cleaned = -0.000164323 Jy
! Total flux subtracted in 100 components = -0.000164323 Jy
! Clean residual min=-0.007703 max=0.007495 Jy/beam
! Clean residual mean=-0.000004 rms=0.001800 Jy/beam
! Combined flux in latest and established models = 0.881992 Jy
selfcal
! Performing phase self-cal
! Adding 12 model components to the UV plane model.
! The established model now contains 63 components and 0.881991 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Fit before self-cal, rms=0.276729Jy  sigma=1.525185
! Fit after  self-cal, rms=0.276729Jy  sigma=1.525177
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
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
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 0.1521 x 0.3877 at -1.88 degrees (North through East)
! Clean map  min=-0.0074127  max=0.79735 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
save 0048-097_Q
! Writing UV FITS file: 0048-097_Q.uvf
! Writing 63 model components to file: 0048-097_Q.mod
! wwins: Wrote 3 windows to 0048-097_Q.win
! Writing clean map to FITS file: 0048-097_Q.fits
! Writing difmap environment to: 0048-097_Q.par
wlog 0048-097_Q.log
! Syntax error in line:
! wlog 0048-097_Q.log
!    ^
! Unable to identify keyword "wlog"
! 
help log
help logfile
help write
help wlog
! help: no such symbol 'wlog'
help difmap
help
help file_io
A
! Syntax error in line:
! A
! ^
! "a" is ambiguous and could match any of:
! Function: abs()
! Function: acos()
! Function: addcmp()
! Function: addhist()
! Function: addwin()
! Function: advance()
! Variable: amplitude
! Help topic: antenna_names
! Function: apropos()
! Function: asin()
! Function: atan()
! Function: atan2()
! 
help
help general
apropos log
device 0048-097_Qim.ps
! Attempting to open device: '0048-097_Qim.ps'
! Error occured in command: device
! 
device 0048-097_Qim.ps/vps
! Attempting to open device: '0048-097_Qim.ps/vps'
mappl cln
quit
! Quitting program
! Log file difmap.log closed on Mon Aug 12 10:21:00 2002