pro testshift, shiftot

n1=640
n2=512
pwr=n2

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

   nwindow = 15

   nb = nwindow - 1

   name1 = '/home/serge/spicules/spi242'

   file = '/home/serge/spicules/spi24200.bin'

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

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

   xmaxshift = 30
   ymaxshift = 30
  

   for i = 1, 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)

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

; next file
endfor

return
end