vis3d/verify_ros.m

   1 function err=check_ros(fuel)
   2 % err=check_ros(fuel)
   3 nfuels=length(fuel);
   4 for k=1:nfuels
   5     nwinds=length(fuel(k).wind);
   6     for j=1:nwinds
   7         speed=fuel(k).wind(j);
   8         ros=fire_ros(fuel(k),speed,0);
   9         err(k,j,1)=fuel(k).ros_wind(j)-ros;
  10     end
  11     nslopes=length(fuel(k).slope);
  12     for j=1:nslopes
  13         tanphi=fuel(k).slope(j);
  14         ros=fire_ros(fuel(k),0,tanphi);
  15         err(k,j,2)=fuel(k).ros_slope(j)-ros;
  16     end
  17 end