pro testshift1, shiftot

n1=640
n2=512
pwr=n2

   im_init   = bytarr(n1,n2)
   im        = bytarr(pwr,pwr)
   refimage  = bytarr(pwr,pwr)
   shiftot   = fltarr(2,15)

   kmax=0
   kkmax=0
   lmax=0
   llmax=0
  
   testfor0,kmax,kkmax,lmax,llmax

   nwindow = 15

   nb = nwindow - 1

   name1 = 'spi_c242'

   file = 'spi_c24200.bin'

   openr, 1, file
   readu, 1, im_init
   close,1

   refimage = im_init(0:511,0:511)

   imcut = refimage(lmax:511-llmax,kmax:511-kkmax)
   refimage = congrid( imcut, 512,512,/interp)

   xmaxshift = 30
   ymaxshift = 30
  

   for i = 0, nb do begin
   ;=====================

   number = string( format='(2I)', i)
   no     = strcompress( number,/remove_all)
   ss =  no
   if (i lt 10) then ss = '0' + no

   file = name1 + ss + '.bin'

   openr, 1, file
   readu, 1, im_init
   close,1

   im = im_init(0:511,0:511)

   imcut = im( lmax:511-llmax,kmax:511-kkmax)
   im    = congrid( imcut,512,512,/interp)

   shift = corr_2d( refimage, im, xmaxshift, ymaxshift )
   shiftot(*,i) = shift
   print,shift

; next file
endfor

return
end