if (add==1) % add
  a=get(ha,'Currentpoint');
  a=a*scan_rate;
  a=round(a(1,1));
  findwin=10;
  [i1 i2]=max(data(currchannel,a+st-findwin:a+st+findwin));
  datind=a+st+i2-findwin-1;
  i=min(find(spks(ic,:)>datind));  
  if ~isempty(i)
    spks(ic,i+1:size(spks,2)+1)=spks(ic,i:size(spks,2));
    spks(ic,i)=datind;
  else
    spks(ic,max(find(spks(ic,:)))+1)=datind;
  end
  end
%keepaxes=1; ha=axes
  makind=spks(ic,find(spks(ic,:)));
plotspikes