function Scut=coralTimeAlign(S);
%DESCRIPTION
%Scut=coralTimeAlign(S);
%
%This function takes a coral structure and all of the corresponding start
%times, finds their overlap, and returns the aligned structure.
%The only input argument is the structure S.  
%
%INPUT: coral structure
%OUTPUT: coral structure containing data temporally common to all
%seismograms.
%
%NOTE: This function uses coralCut.
%
% Josh Carmichael, 14.Aug.2007
%%
% GET MAX START TIME
[sM,sN]=size(S);
if(sM<sN); S=S'; end;
starts=[S.recStartTime]';
[N,M]=size(starts);
numDate=datenum(starts);
[maxStart,iMaxS]=max(numDate);
[minStart,iMinS]=min(numDate);
firstCut=starts(iMaxS,:); 
firstCut=firstCut';
lastCut=starts(iMinS,:);
lastCut=lastCut';
tbegin=firstCut; %to keep

%GET END TIME IN NUMBER OF SECONDS FROM START
% Assume a uniform sample rate!
numData=[S.recNumData]';
[minNum,iNum]=min(numData);
absNumData=numData(iNum);
sInts=[S.recSampInt]';
sInt=sInts(iNum);
numSeconds=absNumData*sInt;
DT=timediff(firstCut,lastCut);
timeRec=numSeconds-DT;
tend=timeadd(tbegin,timeRec); %to keep

%CUT THE OVERLAPPING TIME WINDOW OUT
opt.cutType='absTime';
opt.absStartTime=repmat(tbegin,1,N);
opt.absEndTime=repmat(tend,1,N);
Scut=coralCut(S,opt);
for(k=1:N)
    Scut(k).recEndTime=tend;
end;
if(sM<sN); Scut=Scut'; end;

%DISPLAY START AND END TIME
% disp(' ')
% sprintf('Start Time is:  %f', tbegin)
% disp(' ')
% sprintf('End Time is:  %f', tend)