! Started logfile: difmap.log on Mon Feb 20 09:47:52 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
![@1342+663_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 1342+663_S
! Reading UV FITS file: 1342+663_S.uvf
! AN table 1: 31 integrations on 171 of 171 possible baselines.
! AN table 2: 30 integrations on 171 of 171 possible baselines.
! AN table 3: 81 integrations on 171 of 171 possible baselines.
! AN table 4: 38 integrations on 171 of 171 possible baselines.
! Apparent sampling: 0.172027 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 1342+663
! 
! 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 21180 visibilities.
! A total of 54 model components were read from file 1342+663_S.mod
! rwins: Read 2 windows from 1342+663_S.win
! Multi-model file "1342+663_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: 1342+663_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 54 model components to the UV plane model.
! The established model now contains 54 components and 0.822176 Jy
! Inverting map and beam 
! Estimated beam: bmin=4.028 mas, bmaj=4.763 mas, bpa=7.844 degrees
! Estimated noise=361.035 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: 4.028 x 4.763 at 7.844 degrees (North through East)
! Clean map  min=-0.0020326  max=0.79178 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.00226188 Jy
! Component: 100  -  total flux cleaned = 0.00356762 Jy
! Total flux subtracted in 100 components = 0.00356762 Jy
! Clean residual min=-0.002009 max=0.002029 Jy/beam
! Clean residual mean=-0.000013 rms=0.000460 Jy/beam
! Combined flux in latest and established models = 0.825743 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.000845094 Jy
! Component: 100  -  total flux cleaned = 0.00139066 Jy
! Total flux subtracted in 100 components = 0.00139066 Jy
! Clean residual min=-0.001964 max=0.001969 Jy/beam
! Clean residual mean=-0.000011 rms=0.000453 Jy/beam
! Combined flux in latest and established models = 0.827134 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000315031 Jy
! Component: 100  -  total flux cleaned = 0.000525818 Jy
! Total flux subtracted in 100 components = 0.000525818 Jy
! Clean residual min=-0.001943 max=0.001935 Jy/beam
! Clean residual mean=-0.000010 rms=0.000449 Jy/beam
! Combined flux in latest and established models = 0.82766 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000196315 Jy
! Component: 100  -  total flux cleaned = 0.000344256 Jy
! Total flux subtracted in 100 components = 0.000344256 Jy
! Clean residual min=-0.001924 max=0.001911 Jy/beam
! Clean residual mean=-0.000010 rms=0.000446 Jy/beam
! Combined flux in latest and established models = 0.828004 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.000175496 Jy
! Component: 100  -  total flux cleaned = 0.000306815 Jy
! Total flux subtracted in 100 components = 0.000306815 Jy
! Clean residual min=-0.001910 max=0.001887 Jy/beam
! Clean residual mean=-0.000010 rms=0.000444 Jy/beam
! Combined flux in latest and established models = 0.828311 Jy
clean
! clean: niter=100  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 9.43957e-05 Jy
! Component: 100  -  total flux cleaned = 0.000156467 Jy
! Total flux subtracted in 100 components = 0.000156467 Jy
! Clean residual min=-0.001900 max=0.001877 Jy/beam
! Clean residual mean=-0.000010 rms=0.000442 Jy/beam
! Combined flux in latest and established models = 0.828467 Jy
mappl
! 
! 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: 4.028 x 4.763 at 7.844 degrees (North through East)
! Clean map  min=-0.0018204  max=0.79167 Jy/beam
! Clean map  min=-0.0018204  max=0.79209 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
mappl
! Adding 55 model components to the UV plane model.
! The established model now contains 109 components and 0.828467 Jy
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
cmul=3*imstata(rms)
! Syntax error in line:
! cmul=3*imstata(rms)
!              ^
! Unable to identify keyword "imstata"
! 
cmul=3*imstat(rms)
mappl cln
! restore: Substituting estimate of restoring beam from last 'invert'.
! Restoring with beam: 4.028 x 4.763 at 7.844 degrees (North through East)
! Clean map  min=-0.0018173  max=0.79209 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
xyran -30,30,-30,30
! The mapplot display area will be limited to:
!  Relative  RA: -30 -> 30 (milli-arcsec)
!  Relative DEC: -30 -> 30 (milli-arcsec)
mappl cln
! 
! Move the cursor into the plot window and press 'H' for help
mapcol non
! Mapplot colormap: none, contrast: 1 brightness: 0.5.
mappl cln
! 
! Move the cursor into the plot window and press 'H' for help
exit
! Exiting program
! Log file difmap.log closed on Mon Feb 20 09:57:30 2006