! Started logfile: difmap.log on Thu Feb 16 11:04:55 2006
![@/home/corey/survey/vsop_difmap_macros]
! Unable to open command file: /home/corey/survey/vsop_difmap_macros
! The file name was taken from environment variable DIFMAP_LOGIN
![@0234+285_S.par]
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 0234+285_S
! Reading UV FITS file: 0234+285_S.uvf
! AN table 1: 91 integrations on 171 of 171 possible baselines.
! AN table 2: 18 integrations on 171 of 171 possible baselines.
! AN table 3: 17 integrations on 171 of 171 possible baselines.
! AN table 4: 37 integrations on 171 of 171 possible baselines.
! Apparent sampling: 0.116851 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 0234+285
! 
! 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 1569 lines of history.
! 
! Reading 13028 visibilities.
! A total of 136 model components were read from file 0234+285_S.mod
! rwins: Read 4 windows from 0234+285_S.win
! Multi-model file "0234+285_S.mtab" not available
mapsize 1024,0.4, 1024,0.4
! Map grid = 1024x1024 pixels with 0.400x0.400 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=100
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
![Exited script file: 0234+285_S.par]
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
cmul=3*imstat(rms)
! Adding 136 model components to the UV plane model.
! The established model now contains 136 components and 1.41652 Jy
! Inverting map and beam 
! Estimated beam: bmin=3.67 mas, bmaj=5.327 mas, bpa=-9.61 degrees
! Estimated noise=593.709 mJy/beam.
loglevs 1,4096
! The new contour levels are:
!  -1 1 2 4 8 16 32 64 128 256 512 1024 2048
mapcol col
! Mapplot colormap: rainbow, contrast: 1 brightness: 0.5.
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.67 x 5.327 at -9.61 degrees (North through East)
! Clean map  min=-0.0060712  max=0.82963 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mappl
! Inverting map 
! 
! 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.00361547 Jy
! Component: 100  -  total flux cleaned = 0.00491424 Jy
! Total flux subtracted in 100 components = 0.00491424 Jy
! Clean residual min=-0.006397 max=0.005566 Jy/beam
! Clean residual mean=0.000178 rms=0.001283 Jy/beam
! Combined flux in latest and established models = 1.42144 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000642638 Jy
! Component: 100  -  total flux cleaned = 0.000784191 Jy
! Total flux subtracted in 100 components = 0.000784191 Jy
! Clean residual min=-0.006317 max=0.005352 Jy/beam
! Clean residual mean=0.000173 rms=0.001253 Jy/beam
! Combined flux in latest and established models = 1.42222 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000128165 Jy
! Component: 100  -  total flux cleaned = 0.000123788 Jy
! Total flux subtracted in 100 components = 0.000123788 Jy
! Clean residual min=-0.006201 max=0.005210 Jy/beam
! Clean residual mean=0.000172 rms=0.001235 Jy/beam
! Combined flux in latest and established models = 1.42235 Jy
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.000112373 Jy
! Component: 100  -  total flux cleaned = -4.5135e-06 Jy
! Total flux subtracted in 100 components = -4.5135e-06 Jy
! Clean residual min=-0.006106 max=0.005277 Jy/beam
! Clean residual mean=0.000172 rms=0.001222 Jy/beam
! Combined flux in latest and established models = 1.42234 Jy
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.00444955 Jy
! Component: 100  -  total flux cleaned = 0.00445305 Jy
! Total flux subtracted in 100 components = 0.00445305 Jy
! Clean residual min=-0.005889 max=0.004837 Jy/beam
! Clean residual mean=0.000143 rms=0.001183 Jy/beam
! Combined flux in latest and established models = 1.42679 Jy
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.0018674 Jy
! Component: 100  -  total flux cleaned = 0.00163235 Jy
! Total flux subtracted in 100 components = 0.00163235 Jy
! Clean residual min=-0.005875 max=0.004718 Jy/beam
! Clean residual mean=0.000132 rms=0.001162 Jy/beam
! Combined flux in latest and established models = 1.42843 Jy
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.000109165 Jy
! Component: 100  -  total flux cleaned = -0.000108806 Jy
! Total flux subtracted in 100 components = -0.000108806 Jy
! Clean residual min=-0.005829 max=0.004637 Jy/beam
! Clean residual mean=0.000133 rms=0.001150 Jy/beam
! Combined flux in latest and established models = 1.42832 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -1.96946e-06 Jy
! Component: 100  -  total flux cleaned = -1.81218e-06 Jy
! Total flux subtracted in 100 components = -1.81218e-06 Jy
! Clean residual min=-0.005758 max=0.004565 Jy/beam
! Clean residual mean=0.000133 rms=0.001140 Jy/beam
! Combined flux in latest and established models = 1.42832 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 1.09205e-06 Jy
! Component: 100  -  total flux cleaned = 1.58061e-06 Jy
! Total flux subtracted in 100 components = 1.58061e-06 Jy
! Clean residual min=-0.005698 max=0.004523 Jy/beam
! Clean residual mean=0.000133 rms=0.001130 Jy/beam
! Combined flux in latest and established models = 1.42832 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 9.00036e-08 Jy
! Component: 100  -  total flux cleaned = 4.00505e-08 Jy
! Total flux subtracted in 100 components = 4.00505e-08 Jy
! Clean residual min=-0.005649 max=0.004476 Jy/beam
! Clean residual mean=0.000133 rms=0.001122 Jy/beam
! Combined flux in latest and established models = 1.42832 Jy
mappl
! 
! Move the cursor into the plot window and press 'H' for help
cmu
! Syntax error in line:
! cmu
!   ^
! Missing '=' in variable assignment
! 
cmul=3*imstat(rms)
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.67 x 5.327 at -9.61 degrees (North through East)
! Clean map  min=-0.0054231  max=0.83233 Jy/beam
! Clean map  min=-0.0054234  max=0.82944 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mappl
! Adding 75 model components to the UV plane model.
! The established model now contains 211 components and 1.42832 Jy
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 3.67 x 5.327 at -9.61 degrees (North through East)
! Clean map  min=-0.0054233  max=0.82944 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mapcol non
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
xyran -20,40,-30,30
! The mapplot display area will be limited to:
!  Relative  RA: -20 -> 40 (milli-arcsec)
!  Relative DEC: -30 -> 30 (milli-arcsec)
mappl cln
! 
! Move the cursor into the plot window and press 'H' for help
xyran -30,30,-20,40
! The mapplot display area will be limited to:
!  Relative  RA: -30 -> 30 (milli-arcsec)
!  Relative DEC: -20 -> 40 (milli-arcsec)
mappl cln
! 
! Move the cursor into the plot window and press 'H' for help
dev 0234+285_S.ps/vps
! Attempting to open device: '0234+285_S.ps/vps'
mappl cln
xyran 0,0,0,0
! The mapplot display area will be limited to:
!  Relative  RA: (The whole available range)
!  Relative DEC: (The whole available range)
dev 0234+285_Sim.ps/vps
! Attempting to open device: '0234+285_Sim.ps/vps'
mappl cln
dev 0234+285_Suv.ps
! Attempting to open device: '0234+285_Suv.ps'
! Error occured in command: device
! 
dev 0234+285_Suv.ps/vps
! Attempting to open device: '0234+285_Suv.ps/vps'
rad
! Using default options string "m1"
save 0234+285_S
! Writing UV FITS file: 0234+285_S.uvf
! Writing 211 model components to file: 0234+285_S.mod
! wwins: Wrote 10 windows to 0234+285_S.win
! Writing clean map to FITS file: 0234+285_S.fits
! Writing difmap environment to: 0234+285_S.par
exit
! Exiting program
! Log file difmap.log closed on Thu Feb 16 11:09:40 2006