涉及计算公式：
模型一:Shuttleworth-Wallace (S-W)：
（1）raa计算：
hc=10.5, K:0.41， za = hc + 2， d0 = hc * 2 / 3， Zw=15，n = 4.25，U已知输入值。
U* = K * Ua / Log((za - d0) / z0c)
Ua = U * Log((za - d0) / z0c) / Log((Zw - d0) / z0c)
Kh = K * U*
* (hc - d0)
raa = Log((za - d0) / (hc - d0)) / (K * Ustar) + hc / (n * Kh) * (Exp(n * (1 - (z0c + d0) / hc)) - 1)
（2）rac计算：
W= 0.05，
rb = 100 / n * (w / Uh) ^ 0.5 / (1 - Exp(-n / 2))
Uh = U * Log((hc - d0) / z0c) / Log((Zw - d0) / z0c)
rac = rb * σ / LAI （LAI已知输入值）
（3）rsc计算：
1、rstmin最小气孔阻力，rstmax最大气孔阻力/
2、LAI <= 2 时： LAIe = LAI
3、LAI＞2时：LAIe = 2 + (LAI - 2) / 3
4、F1 = (rstmin / rstmax + ff) / (1 + ff)
5、ff = 0.55 * Rs / Rscri * 2 / LAI（Rs < 0 则 Rs = 0，Rs为已知输入值）
6、如果：Thet >= Thetf 则 F2 = 1
如果：Thet <= Thetr 则F2 = 0
其他： F2 = (Thet - Thetr) / (Thetf - Thetr)
7、F3 = 1 - 0.061 * VPD
VPD = est – ea
est = 0.6108 * Exp(17.27 * Ta / (Ta + 237.3))（Ta已知输入值）
ea = est * RH / 100（RH已知输入值）
8、F4 = 1 - 0.0016 * (295-Tak)2
Tak = Ta + 273
9、如果 F1 <= 0 Or F2 <= 0 Or F3 <= 0 Or F4 <= 0 则 rsc = 50000:
其他则： rsc = rstmin /(LAIe * F1 * F2 * F3 * F4)
（4）ras计算：
1、z0g = 0.02，
2、ras = hc * Exp(n) / (n * Kh) * (Exp(-n * z0g / hc) - Exp(-n * (z0c + d0) / hc))
（5）rss计算：
1、b1 = 35: b2 = 2.3: b3 = 100
2、θ_s=0.4296，θ=0.4229
3、rss = b1 * (Thetf / Thets) b2 + b3
（5）蒸散（ET）计算：
1、Δ= Delta = 4098 * 0.6108 * Exp(17.27 * Ta / (Ta + 237.3)) / (Ta + 237.3) ^ 2
2、γ= Gamma = Cp * P / (epsilon * Lambda)
3、λ=Lambda = 2.501 - 0.002361 * Ta
4、ρ=rho = P / (Tvk * R)
5、ε=epsilon= 0.622
6、Gsc = 0.082
7、Cp = 0.001013
8、rra = (Delta + Gamma) * raa
9、rrc = (Delta + Gamma) * rac + Gamma * rsc
10、rrs = (Delta + Gamma) * ras + Gamma * rss
11、Cc = 1 / (1 + (rrc * rra) / (rrs * (rrc + rra)))
12、Cs = 1 / (1 + (rrs * rra) / (rrc * (rrs + rra)))
13、LEc = (Delta * (Rn - G) + (rho * Cp * (est - ea) * 1000000 - Delta * rac * (Rnss - G)) / (raa + rac)) / (Delta + Gamma * (1 + rsc / (raa + rac)))
Rn、G已知输入值
Rnss = Rn * Exp(-Cr * LAI)
Cr = 0.5
14、LEs = (Delta * (Rn - G) + (rho * Cp * (est - ea) * 1000000 - Delta * ras * (Rn - Rnss)) / (raa + ras)) / (Delta + Gamma * (1 + rss / (raa + rac)))
15、LE = Cc * LEc + Cs * Les
16、ET = 3.6 * (Cc * LEc + Cs * LEs) / Lambda
17、T = 3.6 * (Cc * LEc) / Lambda
18、E = 3.6 * (Cs * LEs) / Lambda
模型二、Penman-Monteith 模型
LE=(Δ（R_n-G）+ρC_p （e_s
- e_a
）/r_a)/(Δ+γ(1+r_sc)/r_a )
（1）
（1）ra计算：
如果：Ua > 2则：
ra = (Log((za - d0) / z0c)) ^ 2 / (K ^ 2 * Ua)
其他：
ra = 4.72 * (Log((za - d0) / z0c)) ^ 2 / (1 + 0.54 * Ua)
（2）rsc计算：
rsc = rstmin / (LAIe * F1 * F2 * F3 * F4)（具体算法同上个模型中）
（3）蒸散计算
1、LE = (Delta * (Rn - G) + rho * Cp * 1000000 * VPD / ra) / (Delta + Gamma * (1 + rsc / ra))
2、ET = 3.6 * LE / Lambda（具体算法同上个模型中）
