drawfire.m

   1 function p=drawfire(f)
   2 % p=drawfire(f)
   3 % show heatflux stored in wrfout
   4 vars={'FIRE_AREA','FGRNHFX','LFN','TIGN_G','FUEL_FRAC','FXLONG','FXLAT','NFUEL_CAT','Times'};
   5 p=nc2struct(f,vars,{'DX','DY'});
   6 nframes=size(p.fgrnhfx,3);
   7 [i,j,a]=find(p.fgrnhfx(:,:,end));
   8 is=min(i);
   9 ie=max(i);
  10 js=min(j);
  11 je=max(j);
  12 r=max(ie-is,je-js)/2;
  13 ci = (is+ie)/2;
  14 cj = (js+je)/2;
  15 is=round(ci-r);
  16 ie=round(ci+r);
  17 js=round(cj-r);
  18 je=round(cj+r);
  19 for n=2:nframes
  20     figure(1)
  21     subplot(1,2,1)
  22     mesh(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),p.fgrnhfx(is:ie,js:je,n)),
  23     title(['FGRNHFX frame ',num2str(n)])
  24     subplot(1,2,2)
  25     %{
  26     mesh(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),p.fire_area(is:ie,js:je,n))
  27     hold on
  28     contour(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),p.lfn(is:ie,js:je,n),[0 0],'k')
  29     hold off
  30     view(0.9,90)
  31     title(['FIRE\_AREA frame ',num2str(n)])
  32     subplot(2,2,3)
  33     mesh(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),p.fire_area(is:ie,js:je,n-1))
  34     view(0.9,90)
  35     title(['FIRE\_AREA frame ',num2str(n-1)])
  36     subplot(2,2,4)
  37     %}
  38     mesh(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),...
  39         p.fuel_frac(is:ie,js:je,n-1)-p.fuel_frac(is:ie,js:je,n))
  40     title(['FUEL\_FRAC frame ',num2str(n),' decrement'])
  41     drawnow
  42     pause(0.3)
  43     % figure(2)
  44     % mesh(p.fxlong(is:ie,js:je,n),p.fxlat(is:ie,js:je,n),p.fire_area(is:ie,js:je,n))
  45     % view(0.9,90)
  46 end