getfile.m

function f=getfile(fin,f)
% 
% try to find data files given parts of directory file name 
% fin='file*'   part of a data file or directory
%

if nargin < 1, fin = pwd; end
if length(fin)<1, fin=pwd; end
disp([' GETFILE: search near file: ',fin])

fdo=0;
fup=0;

% directory up to this point
id=find(fin=='/' | fin=='\');
if length(id)>0
 fin1=fin(1:id(end));
else
 fin1='';
end

% look deeper 
a=dir([fin,'*']);
ia=1:size(a,1);

if length(ia)<1
 ii=id(end):length(fin);
 fin(ii)=lower(fin(ii));
 a=dir([fin,'*']);
 ia=1:size(a,1);
 disp([' try lower case ',fin])
end

if length(ia)<1
 ii=id(end):length(fin);
 fin(ii)=upper(fin(ii));
 disp([' try upper case ',fin])
 a=dir([fin,'*']);
 ia=1:size(a,1);
end


if length(ia)<1
 f.ladcpdo='';
 f.ladcpup='';
 disp(' no ADCP files found ')
 return
end

% check if directories are listed
if getfield(a,{1},'isdir')==1
 disp(' GETFILE: found directory/ies ')
 % parse out up/down looking directory
 ii=ia;
 for j=1:length(ii)
  i=ii(j);
  fin2=getfield(a,{i},'name');
  if length(fin2)>1
   if fin2(end)=='d' | fin2(end)=='D' | fin2(1:2)=='dn' | fin2(1:2)=='DN'
    f.ladcpdo=[fin1,fin2,filesep];
    fdo=i; ia(j)=NaN;
   elseif fin2(end)=='u' | fin2(end)=='U' | fin2(1:2)=='up' | fin2(1:2)=='UP'
    f.ladcpup=[fin1,fin2,filesep];
    fup=i; ia(j)=NaN;
   end
  end
 end

 ii=ia(isfinite(ia));
 for j=1:length(ii)
  i=ii(j);
  fin2=getfield(a,{i},'name');
  if fup==0
   if fin2(end)=='s' | fin2(end)=='S' 
    f.ladcpup=[fin1,fin2,filesep];
    fup=i; ia(j)=NaN;
   end
  end
 end

 ia=ia(isfinite(ia));
 if fdo==0 & ~isempty(ia)
  fin2=getfield(a,{ia(1)},'name');
  f.ladcpdo=[fin1,fin2,filesep];
  fdo=99;
 end

else
 % we are in directory
 % try to find down and up files
 f.ladcpdo=largefile(fin,'*d');
 if length(f.ladcpdo)<2
  f.ladcpdo=largefile(fin,'*D');
 end
 if length(f.ladcpdo)<2
  f.ladcpdo=largefile(fin,'*d0');
 end
 if length(f.ladcpdo)<2
  f.ladcpdo=largefile(fin,'*D0');
 end

 f.ladcpup=largefile(fin,'*u');
 if length(f.ladcpup)<2
  f.ladcpup=largefile(fin,'*U');
 end
 if length(f.ladcpup)<2
  f.ladcpup=largefile(fin,'*u0');
 end
 if length(f.ladcpup)<2
  f.ladcpup=largefile(fin,'*U0');
 end

 if length(f.ladcpdo)<2
  f.ladcpdo=largefile(fin);
 end


 return
end

if fdo>0
   f.ladcpdo=largefile(f.ladcpdo);
end

if fup>0
 f.ladcpup=largefile(f.ladcpup);
end

%------------------------------------------------------------------------------
function fout=largefile(fin,fin0)
% search for largest file in directory
if nargin<2, fin0=''; end
fout=' ';

% directory up to this point
ii=find(fin=='/' | fin=='\');
if length(ii)>0
 fin1=fin(1:ii(end));
else
 fin1='';
end

a=dir([fin,fin0,'*']);
ia=1:size(a,1);
% loop to get file size
if length(ia)>1
   ib = [];
   for i=ia
      fsize(i)=getfield(a,{i},'bytes');
      if fsize(i) > 10
	 if isbb(fopen([fin1, getfield(a,{i},'name')]));
	    ib = [ib i];
	 end
      end
   end
   if ~isempty(ib)
      [dum,ii]=max(fsize(ib));
      fout=[fin1,getfield(a,{ib(ii)},'name')];
   end
elseif length(ia)==1
  fout=[fin1,getfield(a,{ia},'name')];
end
%------------------------------------------------------------------------------
function i = isbb(fid)
%ISBB True if broad-band ADCP.

% check header and data source identification bytes
hid = 127;
sid = 127;
id = fread(fid,2,'uint8');
if length(id)<2
 err('ISBB: ****** can not read file id *****')
else
 i = id(1) == hid & id(2) == sid;
end

% rewind file
fseek(fid,0,'bof');