perimeter_new/main_function_ak3.m

   1 %clear
   2 clc
   3 addpath ../netcdf
   4 addpath ../util1_jan
   5 addpath ..
   6 %addpath Output
   7 if ~exist('saved_data','var'),
   8      global saved_data
   9      saved_data=0;
  10 end
  11
  12
  13 %wrfout{1}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-11_00:00:00';
  14 % wrfout{2}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-12_00:00:00';
  15 % wrfout{1}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-13_00:00:00';
  16 % wrfout{4}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-14_00:00:00';
  17 % wrfout{5}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-15_00:00:00';
  18 % wrfout{6}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-16_00:00:00';
  19 % wrfout{7}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-17_00:00:00';
  20 % wrfout{8}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-18_00:00:00';
  21 % wrfout{9}='/share_home/akochans/WRF341F/wrf-fire/WRFV3/test/em_utfire_1d_med_4km_nofire/wrfout_d01_2013-08-19_00:00:00';
  22
  23 wrfout{1}='wrfout_d01_2013-08-11_00:00:00';
  24 wrfout{2}='wrfout_d01_2013-08-12_00:00:00';
  25 wrfout{3}='wrfout_d01_2013-08-13_00:00:00';
  26 wrfout{4}='wrfout_d01_2013-08-14_00:00:00';
  27 wrfout{5}='wrfout_d01_2013-08-15_00:00:00';
  28 wrfout{6}='wrfout_d01_2013-08-16_00:00:00';
  29 wrfout{7}='wrfout_d01_2013-08-17_00:00:00';
  30 wrfout{8}='wrfout_d01_2013-08-18_00:00:00';
  31 wrfout{9}='wrfout_d01_2013-08-19_00:00:00';
  32
  33
  34 input_type=2;     % Type of the input, 1- data file; 0 - wrfout file; 2- fire_area - is given in input_file;
  35 input_file='perimeter.txt';
  36 %input_file='Input_data_Patch_Springs_8-12-2013_2123.dat';
  37 % File that contains perimeter data
  38                   % It is being used only when input_type=1;;
  39 num_wrf=9;        % The total number of wrfouts that are being used
  40 time_step=48;      % The amount of time_steps in each wrfout
  41 interval=1800;      % time step in wrfout in seconds (every 15 min=900 sec)
  42 time=48;           % the number of the time step in the latest wrfout
  43                   %  closest to the perimeter time from below
  44 format long
  45
  46 % read the fire map at perimeter time
  47 [long,lat,fire_area]=read_file_perimeter(wrfout{num_wrf},wrfout{num_wrf}, time,input_type,input_file);
  48 % now have data: long, lat,
  49 % fire_area - (input_type=0)- burning or not burning (between 0 and 1, 0-1 OK)
  50 %             (input_type=1)- set of ordered points of the boundary 1st=last;
  51 %                  bound(i,1)-horisontal; bound(i,1)-vertical coordinate
  52 tign=perimeter_in(long,lat,fire_area,wrfout,time,interval,time_step,num_wrf, input_type);
  53
  54 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  55 % The information below is not being used for the current problem
  56 %data='data_perim.txt';
  57 %data='';            % contained the coordinates of the fire perimeter from the shapefile
  58 %wrfout='wrfout_d05_2012-09-12_21:15:01';
  59 %wrfout='wrfout_d05_2012-09-09_00:00:00'; - for the earlies Baker's fire
  60 %time=281;
  61 % time =49 for Witch fire;
  62 % time=100; for Baker's fire;
  63
  64
  65
  66
  67
  68
  69
  70