aml 做的　scs曲线计算径流模型

/*accept the arguments
&s oddem05 = [response 'Enter your dem' oddem05]
&s odlu05 = [response 'enter your lu' odlu05]
&s FC = [response '田间持水量状态下的土壤含水量' 0.25]
&S SAT = [response '饱和状态下的土壤含水量' 0.25]
&s SW = [response '降雨前土壤含水量' 0.2]
&s Rday = [response '日降雨量(mm)' 200]
/*define the internal var
&s slp = slp
&s cn2 = cn2
&s cn3 = cn3
&s cn2s = cn2s
&s cn3s = cn3s
&s cn1s = cn1s
&s smax = smax
&s smin = smin
&s w2 = w2
&s w1 = w1
&s s = s
/*define the result val
&s Qsurf = Qsurf
/* the process of the computing
%slp% = slope(oddem05, percentrise)
items %odlu05%.vat
list %odlu05%.vat
if(%odlu05% == 1)
%cn2% = 83
else if(%odlu05% == 2)
%cn2% = 61
else if(%odlu05% == 3)
%cn2% = 48
else if(%odlu05% == 4)
%cn2% = 55
else
%cn2% = -9999
endif
%cn3% = %cn2% * exp(0.00673 * (100 - %cn2%))
%cn2s% = (%cn3% - %cn2%) / 3 * (1 - 2 * exp(- 13.86 * %slp%))
%cn3s% = %cn2s% * exp(0.00673 * (100 - %cn2s%))
%cn1s% = %cn2s% - 20 * (100 - %cn2s%) / (100 - %cn2s% + exp(2.533 - 0.0636 * (100 - %

cn2s%)))
%smax% = 25.4 * (1000 / %cn1s% - 10)
%smin% = 25.4 * (1000 / %cn3s% - 10)
%w2% = (ln(%FC% / (1 - %smax% / %smin%)) - ln(%SAT% / (1 - 2.54 / %smax%) - %SAT%)) / (%SAT%

- %FC%)
%W1% = ln(%FC% / (1 - %smax% / %smin%) - %FC%) + %w2% * %FC%
%s% = %smax% * (1 - %SW%) / (%SW% + exp(%w1% - %w2% * %SW%))
%Qsurf% = (%Rday% - 0.2 * %s%) / (%Rday% + 0.8 * %s%)
/* delete all internal layers
kill %slp% all
kill %cn2% all
kill %cn3% all
kill %cn2s% all
kill %cn3s% all
kill %cn1s% all
kill %smax% all
kill %smin% all
kill %w2% all
kill %w1% all
kill %s% all