sortSites.m

function [SitesOrdered,UseSite] = sortSites(Sites,DataFiles,Tlims)
%
%  Usage: [SitesOrdered,UseSite] = sortSites(Sites,DataFiles,Tlims);
%
%  Sort sites by start time of first data;
%  ALSO: if optional third argument Tlims is present, return
%     an indicator array, UseSite = 1 if any segment for this
%     site overlaps time limits given in Tlims.  If Tlims
%     is not input, UseSite is identically 1 for all sites

if nargin < 3
   Tlims = [-10^20,10^20];
end
nSites = length(Sites);
UseSite = ones(nSites,1);
SitesOrdered = cell(nSites,1);
t1 = zeros(nSites,1);
t2 = t1;
for k = 1:nSites
   temp = DataFiles(Sites{k}.inds);
   nFiles = length(temp);
   x = zeros(nFiles,1);
   for l = 1:nFiles
      x(l) = datenum(temp{l}.start);
   end
   t1(k) = min(x);
   t2(k) = max(x); 
end
[dum,ind] = sort(t1);
for k = 1:nSites
   SitesOrdered{k} = Sites{ind(k)};
   UseSite(k) = t1(ind(k)) >= Tlims(1) & t2(ind(k)) <= Tlims(2);
end