cycling/area_compare.m

   1 function [ac,gs] = area_compare(ps,tign_b)
   2 %compares areas of fires over time
   3 % outputs graph, areas and a mask
   4 % ps = cluster_paths(w) --- Forecast
   5 % tign_b = quish2(ps) --- "Analysis"
   6
   7 %tign_a = ps.tign_g;
   8 tign_a = ps.red.tign;
   9 %create a mask
  10 msk=tign_b-tign_a;
  11 msk(abs(msk)<0.2) = 0;
  12 msk(msk>0)=1; %add fuel moisture ==> slow fire
  13 msk(msk<0)=-1; %subtract fuel moist ==> speed up fire
  14
  15 %time +- 1/4 day at simulation start and end
  16 t_1 = max(min(tign_a(:)),min(tign_b(:)));%+0.25;
  17 t_2 = min(max(tign_a(:)),max(tign_b(:)))-0.1;
  18
  19 pts = 30;
  20 t = linspace(t_1,t_2,pts);
  21
  22 for i = 1:pts
  23     area_a(i) = sum(sum(tign_a < t(i)));
  24     area_b(i) = sum(sum(tign_b < t(i)));
  25 end
  26 %figure,plot(t(1:end-1),area_a(1:end-1),t(1:end-1),area_b(1:end-1))
  27 t_days = (t-ps.red.start_datenum);
  28 figure
  29 %plot(t_days,area_a,t_days,area_b)
  30 plot(t-t_1,area_a,t-t_1,area_b)
  31 %legend('wrf-a','wrf-b')
  32 legend('Forecast','Estimate')
  33 %legend('Wet fuel','Normal Fuel')
  34 xlabel('Simulation Time [days]')
  35 ylabel('Simulation Area [gride nodes]')
  36 title('Comparison of Areas')
  37 %xlim([0 6])
  38
  39 % figure,mesh(ps.red.fxlong,ps.red.fxlat,msk)
  40 % title('Mask for fuel adjustments')
  41 % xlabel('lon'),ylabel('lat')
  42 ac.area_a = area_a;
  43 ac.area_b = area_b;
  44 ac.msk = msk;
  45 gs = norm(area_a-area_b)/norm(area_a);
  46
  47 end
  48
  49
  50