detection/readmod14files.m

   1 function r=readmod14files(file_search,ax)
   2 % r=readmod14files(file_search,ax)
   3 %   file_search   search string
   4 %   ax            bounds [min_lon, max_lon, min_lat, max_lat]
   5
   6 % find the files
   7 d=dir(file_search);d={d.name};
   8 if(isempty(d)), error(['No files found for ',file_search]),end
   9 prefix=regexp(file_search,'^.*/','match');  % all up to the last /
  10 if isempty(prefix),  % comes as cell, convert to string
  11     prefix='';
  12 else
  13     prefix=prefix{1};
  14 end
  15
  16 % order the files in time
  17 nfiles=length(d);
  18 t=zeros(1,nfiles);
  19 for i=1:nfiles
  20     t(i)=rsac2time(d{i});
  21 end
  22 [t,i]=sort(t);
  23 d={d{i}};
  24
  25 k=0;
  26 for i=1:nfiles,
  27     file=d{i};
  28     v=readmod14([prefix,file]);
  29      % select fire detection within the domain
  30     xj=find(v.lon >= ax(1) & v.lon <= ax(2));
  31     xi=find(v.lat >= ax(3) & v.lat <= ax(4));
  32     if ~isempty(xi) & ~isempty(xj)
  33         k=k+1;
  34         v.data=v.data(xi,xj);    % subset data
  35         v.lon=v.lon(xj);
  36         v.lat=v.lat(xi);
  37         r.time(k)=t(i);
  38         v.axis=ax;
  39         r.x{k}=v;
  40     end
  41 end
  42 end