function out=loadADCP(t,adcpdir,flg)
% ADCP=loadADCP(TIME,ADCPDIR) loads the data in the directory ADCPDIR
% between times TIME(1) and TIME(2) into memory, and returns that data to
% ADCP.  If there is no adcp data for that period, the empty matrix is
% returned.
% 
% if TIME='all', all adcp data in ADCPDIR is loaded.

if nargin<3
  flg=0;
end

% Load the data.
try
  mload([adcpdir '/adcpInfo.mat'])
catch
  mload([adcpdir '/adcp_info.mat'])
end
if isstr(t) && strcmp(t,'all')
  t=[-inf inf];
end
inds=find(t(1)<adcpinf.time,1,'first'):find(adcpinf.time<t(end),1,'last');
dr=dir(adcpdir);
for i0=1:length(dr)
  try
    [prfx]=getprefix(dr(i0).name);
    break
  end
end

[flnums i7 junk]=unique(adcpinf.trnum(inds));

out=[];
if isempty(flnums)
  return
end

fls=adcpinf.filename(inds(i7),:);
clear adcpinf
for i0=1:size(fls,1);
  tmp=load([adcpdir '/' fls(i0,:)]);
  if isfield(tmp,'adcp')
    out=mergefields(out,tmp.adcp);
  end
end
%keyboard
if isfield(out,'depth')
  out.depth=out.depth(:,1);
end
if isfield(out,'binpos')
  out.binpos=out.binpos(:,1);
end


fldnms=fieldnames(out);
if isfield(out,'time')
  inds=find(t(1)<out.time & out.time<t(end));
  nn=length(out.time);
else
  inds=find(within(out.navfirsttime,t));
  nn=length(out.navfirsttime);
end

for i0=1:length(fldnms)
  if size(out.(fldnms{i0}),2)==nn
    out.(fldnms{i0})=out.(fldnms{i0})(:,inds);
  else
    out.(fldnms{i0})=out.(fldnms{i0});
  end
end  

if nargout==0
  clear out
end

function dat = mergefields(dat1,dat2,dim2,direction);
% function dat = mergefields(dat1,dat2,dim2);
% merge structures dat1 and dat2 into one structure..
% Optional argument dim2 gives the dimension along which we should
% merge.  If dat1 is empty, then dat2 is returned.
%
% if there is a matrix size mismatch the data is made larger.
 
% $Revision: 1.4 $ $Date: 2003/04/25 20:39:01 $ $Author: aperlin $	

if nargin<3
  dim2=[];
end;

if isempty(dat1)
  dat=dat2;
  return;
end;

if nargin<4
  direction=2;
end;


names = fieldnames(dat1);
dat=dat1;
for i=1:length(names);
  data1 = getfield(dat1,names{i});
  data2 = getfield(dat2,names{i});
  [M1,N1]=size(data1);
  [M2,N2]=size(data2);
  if isempty(dim2)
    if direction==2
      if M1>M2
	data2 =[data2;NaN*ones(M1-M2,N2)];
      elseif M2>M1
	data1 =[data1;NaN*ones(M2-M1,N1)];
      end;
      dat = setfield(dat,names{i},[data1 data2]);
    else
      dat = setfield(dat,names{i},[data1; data2]);
    end;
    
  else
    if size(data2,2)==dim2
      if direction==2
	dat = setfield(dat,names{i},[data1 data2]);
      else
	dat = setfield(dat,names{i},[data1; data2]);
      end;
    end;
  end;    
end;