perimeter/main_function.m

   1 wrfout='wrfout_d01_0001-01-01_00:00:00';
   2 interval=10;      % time step in wrfout in seconds
   3 count=1;          % update the wind every count time steps
   4 time=31;          % the number of the time step in wrfout closest to the perimeter time from above
   5 time_now=300.25;  % the perimeter time (s)
   6
   7 format long
   8
   9 % read burning/not burning map from wrfout, should not set A or tign_g
  10 % can be replaced by any other data acquisition routine
  11
  12 % read the fire map at perimeter time
  13
  14 [long,lat,fire_area]=read_file_perimeter(wrfout,time);
  15
  16 % read long, fire_area, lat from wrfout
  17 % interpolate/resample to the wrf fire grid
  18
  19 % now have data: long, lat, fire_area - burning or not burning (between 0 and 1, 0-1 OK)
  20 % wrfout - has ros
  21
  22 % JM should not take tign_g as input, only tign_g <= time_now,
  23
  24 tign=perimeter_in(long,lat,fire_area,wrfout,time_now,time,interval,count);
  25
  26 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  27 % The information below is not beig used for the current problem
  28 %data='data_perim.txt';
  29 %data='';            % contained the coordinates of the fire perimeter from the shapefile
  30 %wrfout='wrfout_d05_2012-09-12_21:15:01';
  31 %wrfout='wrfout_d05_2012-09-09_00:00:00'; - for the earlies Baker's fire
  32 %time=281;
  33 % time =49 for Witch fire;
  34 % time=100; for Baker's fire;
  35
  36
  37
  38
  39
  40
  41
  42