function [out]=timeshft(in,d)
%
% Time shift a data matrix
%
% [out]=timeshft(in,d);
%
% out - timeshifted data
% in  - non-shifted data
% d   - matrix of time delays
%

[D L]=size(in);

%
% Find the maximum time delay.
% This defines where the new
% data matrix starts
%
max_delay=max(d);
start_pos=max_delay+1;
new_len=D-start_pos;

out=[];
for i = 1 : L
	s=start_pos-d(i);
	f=s+new_len;
	temp=in(s:f,i);
	out(:,i)=temp;
end