cycling/fuel_adjust.m

   1 function fuel_adjust(restart_file)
   2 %function adjust fuel moisture in a  wrf file used for input wrfinput or
   3 %wrfrst
   4
   5 fire_choice = input_num('which fire? Patch: [0], Camp: [1], Cougar: [3]',1)
   6 cycle = input_num('Which cycle? ',0);
   7 load_str = sprintf('p_%d.mat',cycle);
   8 load(load_str);
   9
  10 %nned to beef up what's in red, below
  11 w  = read_wrfout_tign(restart_file);
  12 %use analysis to et region for changing fuel moisture
  13 w.tign_g = p.analysis;
  14 w.tign_g = p.spinup;
  15 red = subset_domain(w);
  16
  17 %use the wrfout file to find subset of detections
  18
  19 time_bounds(2) = red.max_tign;
  20 time_bounds(1) = red.min_tign;
  21
  22 if fire_choice == 1
  23     fire_name = 'Camp fire';
  24     save_name = 'camp';
  25     prefix='../campTIFs/';
  26 elseif fire_choice == 0
  27     fire_name = 'Patch Springs fire';
  28     save_name = 'patch';
  29     prefix='../TIFs/';
  30 else
  31     fire_name = 'Cougar Creek fire';
  32     save_name = 'cougar';
  33     prefix = '../cougarTIFs/';
  34 end
  35
  36 fig.fig_map=0;
  37 fig.fig_3d=0;
  38 fig.fig_interp=0;
  39 if exist('g_match.mat','file')
  40     load('g_match.mat')
  41 else
  42     det_list=sort_rsac_files(prefix);
  43     g = load_subset_detections(prefix,det_list,red,time_bounds,fig);
  44     save('g_match.mat', 'g', '-v7.3')
  45 end
  46 %find list of detections
  47 min_con = 7;
  48 for i = 1:length(g)
  49     if i == 1
  50         fm = g(i).fxdata >=min_con;
  51     end
  52     fm = fm + (g(i).fxdata >=min_con);
  53 end %for i...
  54
  55 % flatten fm to height 1
  56 fm = fm > 0;
  57
  58 forecast_area = red.tign <= g(i).time;
  59
  60 mult = input_num('Enter fuel moisture multiplier: ',1)
  61 %mult=1.4;
  62 too_far = mult*(forecast_area-fm);
  63 new2far = ones(size(too_far));
  64 too_far_mask = too_far > 0;
  65 new2far(too_far_mask) = too_far(too_far_mask);
  66
  67 fast_fuel = red.nfuel_cat(too_far_mask);
  68 %figure,histogram(fast_fuel(:))
  69
  70 %go to w from red
  71 w_mask = ones(size(w.tign_g));
  72 w_mask(red.ispan,red.jspan) = new2far;
  73
  74
  75 fprintf('%s %s\n','Will write modified time into     ',restart_file)
  76 rewrite_bak=[restart_file,'.bak'];
  77 q=input_num(['1 to copy ',restart_file,' to ',rewrite_bak],1);
  78 if q,
  79     if system(['cp ',restart_file,' ',rewrite_bak]),
  80         warning('copy failed')
  81     end
  82 end
  83
  84 %%% change fuel moisture in restart file %%%
  85 new_moist = w.fmc_g;
  86 new_moist = new_moist.*w_mask;
  87 ncreplace(restart_file,'FMC_G',new_moist)
  88
  89
  90
  91 end % function