1 % This code is done for the real case for the witch_fire
2 Please see directions about how to run the code:
4 A) 1) You just want to run witch example from the test file:
5 Put the file test_ignition.m (it is located at wrf-fire/other/Matlab/ignition/) in the folder
9 4) Output will be in data_out1.txt, created especially for wrfinput.
12 B) You want to run the code manually by parts and play with coefficients:
14 1a) Run matlab and go to the folder wrf-fire/WRFV3/run and then do
15 addpath ../../other/Matlab/ignition
18 1b) Put the files: ignition.m; line_dist.m; line_sign.m; read_file_ignition.m in the folder
21 2) Make sure your file, where your perimeter is specified and your wrfout file lie in the same folder.
25 4) Run [unit_long,unit_lat,long,lat,time_now,mesh_size,ign_pnt,bound]=read_file_ignition('data.txt','wrfout_d01_2007-10-21_12:00:00_real_case');
27 5) Run B=ignition(unit_long,unit_lat,long,lat,time_now,mesh_size,ign_pnt,bound);
29 6) Matrix of time of ignition made for wrf_output will be printed to
33 6) Please use read_array_2d to interpret or graph results.
38 7) If you use different from my example perimeter data and wrfout, then
39 specify their names when you run ignition7.m
41 First three rows of your perimeter file should be:
42 1rt row - time_now (second number is not needed, is set to 0);
43 2nd row - size of the mesh;
44 3rd row - coordinates of ignition point;
45 All next rows - coordinates of all the
46 points on the boundary (lon,lat).