vis3d/find_ignition.m

   1 function r=find_ignition(f)
   2 % r=find_ignition(f)
   3 % finds the first frame with nodes where fire_area > 0
   4 % arguments:
   5 %     f wrfout file name
   6 % returns:
   7 %     structure with ignited nodes and their coordinates
   8
   9 p = nc2struct(f,{'Times'},{});
  10 times=char(p.times)';
  11 steps=size(times,1);
  12 r = [];
  13 for step=1:steps
  14     p=nc2struct(f,{'FIRE_AREA'},{},step)
  15     [i,j,a]=find(p.fire_area);
  16     n = length(i);
  17     disp([num2str(step),' ',times(step,:),' ignited ',num2str(n)])
  18     if n
  19         q=nc2struct(f,{'FXLONG','FXLAT'},{},step);
  20         r.file=f;
  21         r.frame=step;
  22         r.i=i';
  23         r.j=j';
  24         r.times=times(step,:);
  25         for k=n:-1:1
  26             r.fxlong(k)=q.fxlong(i(k),j(k));
  27             r.fxlat(k)=q.fxlat(i(k),j(k));
  28             r.fire_area(k)=p.fire_area(i(k),j(k));
  29         end
  30         break
  31     end
  32 end
  33 end
  34
  35
  36
  37
  38