function wfSh = Shift(wf,dt,del);
%function wfSh = Shift(wf,dt,del);
%
%	Shifts series wf, with sample interval dt, by an amount del
%

numSam = length(wf);

ntf=2^ceil(log(numSam)/log(2));
jomega=sqrt(-1)*(0:ntf/2)*2*pi/(ntf*dt);

spwf=fft(wf,ntf);
spwfs=spwf(1:ntf/2+1).*exp(-(del)*jomega);
wfshift = flipud(real(ifft([spwfs,conj(fliplr(spwfs(2:ntf/2)))])));
wfSh = wfshift(1:numSam);

