standalone/make_rain_input.m

   1 function make_rain_input(n,mmph,hours)
   2 % create moisture_input.txt for moisture_test.exe
   3 % usage: make_rain_input(n,mmph)
   4 % input:
   5 %    n        number of time steps
   6 %    mmph     rain intensity mm/h
   7 %    hours    length of simulation in hours
   8 % example
   9 % make_rain_input(n,100,14) should give 1.723 for 14-hour wetting time lag
  10 %
  11 % note: the last line in the output from moisture_test.exe
  12 %       should not depend on n
  13
  14 mm=[hours*[0:n]'/n,ones(n+1,1)*[300,1e5,0.01],mmph*hours*[0:n]'/n];
  15 save('moisture_input.txt','mm','-ascii')
  16
  17 end