cycling/fmc_change.m

   1 function fmc_change(m,f)
   2 %perchentage change to add/subtract
   3 % msk - locations where to add subtract
   4 % f - string, path to a wrfinput or wrfrst file
   5 %f = 'wrfinput_d01';
   6 %msk = double(msk);
   7 %msk(msk>0) = 1;
   8
   9 %blur mask a little bit
  10 %msk= imgaussfilt(msk,1/2);
  11 %w = read_wrfout_tign(f);
  12 %load sm_mask.mat
  13 s = nc2struct(f,{'FMC_GC'},{})
  14 msk = ones(size(s.fmc_gc(:,:,1)));
  15 fprintf('Fuel levels 1--> 1hr, 2-->10hr 3-->100hr 4-->1000hr 5-->live \n')
  16 fprintf('Standard for now is [3,5]: 100hr, live fuels.\n')
  17 f_time = input_num('Which fuel levels? All = -1',[3,5]);
  18 if f_time < 0
  19     moist = m*msk + s.fmc_gc;
  20 else
  21     moist = s.fmc_gc;
  22     for i = 1:length(f_time)
  23         %mask area
  24         moist(:,:,f_time(i)) = m*msk + s.fmc_gc(:,:,f_time(i));
  25         %small adjustment globally
  26         moist(:,:,f_time(i)) = moist(:,:,f_time(i))+m/4;
  27     end
  28 end
  29 %moist2 = moist;
  30 rewrite_bak=[f,'.bak_before_fmc'];
  31 if system(['cp ',f,' ',rewrite_bak])
  32     fprintf('Error in copy \n')
  33 else
  34     fprintf('Copy ok. rewriting FMC_G \n')
  35     ncreplace(f,'FMC_GC',moist)
  36     %ncreplace(f,'FMC_G',moist)
  37 end
  38 end