Article: @@@
Summary:

数值实现植物水分势

该模型采用电路模拟方法来计算四个植物节点的水分势（ψ）：阳光叶（ψsunleaf）、阴暗叶（ψshadeleaf）、茎（ψstem）和根（ψroot）。这些势能是根据土壤基质势和非胁迫性蒸腾作用计算得出的，进而用于计算气孔导度和水分胁迫因子（βt）。

关键步骤：
1. 根据水分势定义蒸腾流（Esun, Eshade）和水供应流（q1a, q1b, q2, qsoil）。
2. 使用线性化显式方法结合牛顿法迭代求解满足水流连续性方程的水分势向量（ψ）。
3. 通过解矩阵方程 AΔψ = b 计算水分势的变化（Δψ），其中 A 和 b 由流方程导出。
4. 使用 ψm+1 = ψm + Δψ 更新水分势，迭代直至水流平衡（b → 0）。
5. 利用最终的蒸腾解计算气孔导度和水分胁迫因子（βt,sun, βt,shade），这些因子随后用于光合作用模块。
6. 将 ψ 的解保存为植物水分势，土壤到根部的流（q3,i）作为理查兹方程地下流方程中的土壤蒸腾汇。
@@@