! Started logfile: difmap.log_2 on Mon May  9 13:01:38 2005
![@0552+398_K.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 0552+398_K
! Reading UV FITS file: 0552+398_K.uvf
! AN table 1: 16 integrations on 55 of 55 possible baselines.
! AN table 2: 16 integrations on 55 of 55 possible baselines.
! Apparent sampling: 0.211932 visibilities/baseline/integration-bin.
! *** This seems a bit low - see "help observe" on the binwid argument.
! Found source: 0552+398
! 
! There are 8 IFs, and a total of 8 channels:
! 
!  IF  Channel    Frequency  Freq offset  Number of   Overall IF
!       origin    at origin  per channel   channels    bandwidth
!  ------------------------------------------------------------- (Hz)
!  01        1  2.37245e+10        8e+06          1        8e+06
!  02        2  2.37375e+10        8e+06          1        8e+06
!  03        3  2.37765e+10        8e+06          1        8e+06
!  04        4  2.38545e+10        8e+06          1        8e+06
!  05        5  2.39975e+10        8e+06          1        8e+06
!  06        6  2.41015e+10        8e+06          1        8e+06
!  07        7  2.41665e+10        8e+06          1        8e+06
!  08        8  2.41925e+10        8e+06          1        8e+06
! 
! Polarization(s): RR
! 
! Read 788 lines of history.
! 
! Reading 2984 visibilities.
! A total of 114 model components were read from file 0552+398_K.mod
! rwins: Read 1 windows from 0552+398_K.win
! Multi-model file "0552+398_K.mtab" not available
mapsize 1024,0.03, 1024,0.03
! Map grid = 1024x1024 pixels with 0.0300x0.0300 milli-arcsec cellsize.
select RR, 1, 8
! Selecting polarization: RR,  channels: 1..8
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
! Reading IF 5 channels: 5..5
! Reading IF 6 channels: 6..6
! Reading IF 7 channels: 7..7
! Reading IF 8 channels: 8..8
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
![Exited script file: 0552+398_K.par]
select
! Selecting polarization: RR,  channels: 1..8
! Reading IF 1 channels: 1..1
! Reading IF 2 channels: 2..2
! Reading IF 3 channels: 3..3
! Reading IF 4 channels: 4..4
! Reading IF 5 channels: 5..5
! Reading IF 6 channels: 6..6
! Reading IF 7 channels: 7..7
! Reading IF 8 channels: 8..8
mappl
! Adding 114 model components to the UV plane model.
! The established model now contains 114 components and 1.93044 Jy
! Inverting map and beam 
! Estimated beam: bmin=0.5701 mas, bmaj=0.7372 mas, bpa=-51.78 degrees
! Estimated noise=2.44226 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.0469446 Jy
! Component: 100  -  total flux cleaned = 0.0626038 Jy
! Component: 150  -  total flux cleaned = 0.0641842 Jy
! Component: 200  -  total flux cleaned = 0.064819 Jy
! Total flux subtracted in 200 components = 0.064819 Jy
! Clean residual min=-0.026185 max=0.034674 Jy/beam
! Clean residual mean=-0.000057 rms=0.009778 Jy/beam
! Combined flux in latest and established models = 1.99526 Jy
selfcal
! Performing phase self-cal
! Adding 45 model components to the UV plane model.
! The established model now contains 159 components and 1.99526 Jy
! 
! Correcting IF 1.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 2.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 3.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 67 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 4.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 5.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 6.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 7.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 67 telescope corrections were flagged in sub-array 2.
! 
! Correcting IF 8.
!  A total of 176 telescope corrections were flagged in sub-array 1.
!  A total of 66 telescope corrections were flagged in sub-array 2.
! 
! Fit before self-cal, rms=0.280739Jy  sigma=1.349576
! Fit after  self-cal, rms=0.281169Jy  sigma=1.340294
clean
! Inverting map and beam 
! Estimated beam: bmin=0.5701 mas, bmaj=0.7372 mas, bpa=-51.78 degrees
! Estimated noise=2.44226 mJy/beam.
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00552109 Jy
! Component: 100  -  total flux cleaned = 0.00604967 Jy
! Component: 150  -  total flux cleaned = 0.00552587 Jy
! Component: 200  -  total flux cleaned = 0.00457158 Jy
! Total flux subtracted in 200 components = 0.00457158 Jy
! Clean residual min=-0.024357 max=0.029793 Jy/beam
! Clean residual mean=-0.000054 rms=0.008645 Jy/beam
! Combined flux in latest and established models = 1.99983 Jy
selfcal
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 179 components and 1.99983 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.278843Jy  sigma=1.330308
! Fit after  self-cal, rms=0.278681Jy  sigma=1.329684
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.5701 x 0.7372 at -51.78 degrees (North through East)
! Clean map  min=-0.023314  max=1.4943 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
clean
! Inverting map 
! clean: niter=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = -0.0015915 Jy
! Component: 100  -  total flux cleaned = -0.00202271 Jy
! Component: 150  -  total flux cleaned = -0.0024161 Jy
! Component: 200  -  total flux cleaned = -0.00316817 Jy
! Total flux subtracted in 200 components = -0.00316817 Jy
! Clean residual min=-0.023023 max=0.027956 Jy/beam
! Clean residual mean=-0.000046 rms=0.008030 Jy/beam
! Combined flux in latest and established models = 1.99666 Jy
selfcal
! Performing phase self-cal
! Adding 37 model components to the UV plane model.
! The established model now contains 190 components and 1.99666 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.277283Jy  sigma=1.323576
! Fit after  self-cal, rms=0.277204Jy  sigma=1.323250
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.5701 x 0.7372 at -51.78 degrees (North through East)
! Clean map  min=-0.022487  max=1.4928 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=200  gain=0.03  cutoff=0
! Component: 050  -  total flux cleaned = 0.00610237 Jy
! Component: 100  -  total flux cleaned = 0.000886566 Jy
! Component: 150  -  total flux cleaned = -0.00494153 Jy
! Component: 200  -  total flux cleaned = -0.00805092 Jy
! Total flux subtracted in 200 components = -0.00805092 Jy
! Clean residual min=-0.009482 max=0.009331 Jy/beam
! Clean residual mean=-0.000035 rms=0.004082 Jy/beam
! Combined flux in latest and established models = 1.98861 Jy
selfcal
! Performing phase self-cal
! Adding 60 model components to the UV plane model.
! The established model now contains 250 components and 1.98861 Jy
! 
! Correcting IF 1.
! 
! Correcting IF 2.
! 
! Correcting IF 3.
! 
! Correcting IF 4.
! 
! Correcting IF 5.
! 
! Correcting IF 6.
! 
! Correcting IF 7.
! 
! Correcting IF 8.
! 
! Fit before self-cal, rms=0.271760Jy  sigma=1.292977
! Fit after  self-cal, rms=0.271667Jy  sigma=1.292673
mappl
! Inverting map 
! 
! Move the cursor into the plot window and press 'H' for help
! Mean=-3.14416e-05  rms=0.00401677  min=-0.00994872  max=0.0101649 Jy/beam
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.5701 x 0.7372 at -51.78 degrees (North through East)
! Clean map  min=-0.015055  max=1.4906 Jy/beam
! 
! Move the cursor into the plot window and press 'H' for help
save 0552+398_K
! Writing UV FITS file: 0552+398_K.uvf
! Writing 250 model components to file: 0552+398_K.mod
! wwins: Wrote 2 windows to 0552+398_K.win
! Writing clean map to FITS file: 0552+398_K.fits
! Writing difmap environment to: 0552+398_K.par
device 0552+398_Kuv.ps/vps
! Attempting to open device: '0552+398_Kuv.ps/vps'
radpl
! Using default options string "m1"
uvplot
! Using default options string ""
device 0552+398_Kim.ps/vps
! Attempting to open device: '0552+398_Kim.ps/vps'
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
device /xw
! Attempting to open device: '/xw'
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
! 
! 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
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.5701 x 0.7372 at -51.78 degrees (North through East)
! Clean map  min=-0.015055  max=1.4906 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
device 0552+398_Kim.ps/vps
! Attempting to open device: '0552+398_Kim.ps/vps'
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.5701 x 0.7372 at -51.78 degrees (North through East)
! Clean map  min=-0.015055  max=1.4906 Jy/beam
exit
! Exiting program
! Log file difmap.log_2 closed on Mon May  9 13:08:31 2005