max_hfx.m

   1 % function max_hfx(p,n)
   2 % p from drawfire, n=frame number
   3 [i,j,v]=find(p.fgrnhfx(:,:,n));
   4 [ni,nj]=size(p.fgrnhfx(:,:,n));
   5 [vm,km]=max(v);
   6 km=km(1);
   7 im=i(km);
   8 jm=j(km);
   9 fprintf('max fgrnhfx %f at (%i %i) = (%f %f) fxlong/fxlat\n',...
  10     p.fgrnhfx(im,jm,n),im,jm,p.fxlong(im,jm,n),p.fxlat(im,jm,n))
  11 r=5;
  12 is=max(1,im-r);
  13 ie=min(ni,im+r);
  14 js=max(1,jm-r);
  15 je=min(nj,jm+r);
  16 fprintf('displaying (%i:%i,%i:%i) of (1:%i,1:%i)\n',is,ie,js,ie,ni,nj)
  17 format short
  18 fgrnhfx=p.fgrnhfx(is:ie,js:je,n)
  19 fuel_frac=p.fuel_frac(is:ie,js:je,n)
  20 lfn=p.lfn(is:ie,js:je,n)
  21 tign_g=p.tign_g(is:ie,js:je,n)