From 7fb0fd5fd2656e68db438cbb01b61c1857b5424c Mon Sep 17 00:00:00 2001
From: DankerMu <127004858+DankerMu@users.noreply.github.com>
Date: Fri, 13 Feb 2026 12:34:39 +0800
Subject: [PATCH] =?UTF-8?q?fix:=20SHUD=E6=9C=8D=E5=8A=A1=E5=99=A8=E7=8E=AF?=
 =?UTF-8?q?=E5=A2=83=20baseline=20Step1-3=20=E8=BF=90=E8=A1=8C=E9=AA=8C?=
 =?UTF-8?q?=E8=AF=81?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

- 更新 AutoSHUD 子模块 (Fun.Soil_Geol.R terra 兼容性修复)
- BUG_REPORT.md 追加 Appendix B: 出口边界条件标记 -4→-3 问题记录

验证环境: R 4.3.1 + terra 1.8.86, rSHUD 从子模块安装
验证结果: Step1-3 全部通过, qhh.sp.riv MD5 可复现 (78763f08)
---
 AutoSHUD                          |  2 +-
 testdata/mock_river/BUG_REPORT.md | 55 +++++++++++++++++++++++++++++++
 2 files changed, 56 insertions(+), 1 deletion(-)

diff --git a/AutoSHUD b/AutoSHUD
index dbf7b6c..09e8e2e 160000
--- a/AutoSHUD
+++ b/AutoSHUD
@@ -1 +1 @@
-Subproject commit dbf7b6c967f53bfb7224c83416ed9c689a4cc58d
+Subproject commit 09e8e2e4d3e146b8ca69da2b9faed43e17e5f142
diff --git a/testdata/mock_river/BUG_REPORT.md b/testdata/mock_river/BUG_REPORT.md
index 74d61ce..d30244d 100644
--- a/testdata/mock_river/BUG_REPORT.md
+++ b/testdata/mock_river/BUG_REPORT.md
@@ -696,3 +696,58 @@ Q_bug / Q_true = sqrt(S_bug / S_true) = sqrt(4e-4 / 0.01) = sqrt(0.04) = 0.2
 **@point 错误不影响计算结果**，因为 SHUD 求解器不读取该数据块。
 
 **建议**：修复 bug 后，对比修复前后的 QHH baseline 出口流量时序，量化 Slope 错误对水文模拟结果的实际影响。
+
+---
+
+## B. 出口河段边界条件标记被意外覆盖（`Down = -4` → `-3`）
+
+> 发现日期: 2026-02-13
+> 严重程度: **中等**（改变出口流量计算公式）
+> 状态: 待修复
+
+### B.1 问题概述
+
+Slope bug 修复时重新运行了 `shud.river()`，`sp.RiverDown()` 的默认出口标记为 `-3`，将 2022 年原始数据中有意设置的 `-4`（临界水深边界条件）覆盖为 `-3`（零水深梯度边界条件）。这是修复过程的**无意副作用**，不是 rSHUD 代码 bug。
+
+### B.2 证据链
+
+| 文件 | Down 值 | 来源 |
+|------|---------|------|
+| `runs/qhh/legacy/Calibration/input/qhh/qhh.sp.riv`（2022-06-08 原始） | `-4` | AutoSHUD Step3 产出后手动设置 |
+| `runs/qhh/baseline/input/qhh/qhh.sp.riv.bak.buggy`（buggy 备份） | `-4` | 与 legacy 文件 MD5 一致（`78f5afb3`） |
+| `/scratch/frd_muziyao/rSHUD/qhh/qhh.sp.riv`（旧版 rSHUD 目录） | `-3` | rSHUD 重新生成，默认值 |
+| `runs/qhh/baseline/input/qhh/qhh.sp.riv`（当前 baseline） | `-3` | Slope bug 修复时重跑 `shud.river()` |
+
+- rSHUD 全部 46 个历史版本中 `sp.RiverDown()` 的默认值均为 `rep(-3, nsp)`，从未出现 `-4`
+- AutoSHUD `Step3_BuidModel.R` 和 `terra_compat.R` 中也无 `-4` 赋值逻辑
+- 结论：2022 年原始文件中的 `-4` 是生成后**手动修改**的物理边界条件选择
+
+### B.3 SHUD solver 中 `-3` 与 `-4` 的语义差异
+
+代码位置：`SHUD/src/ModelData/MD_RiverFlux.cpp:36-58`
+
+```
+case -3:  /* 零水深梯度边界 */
+    s = BedSlope + 2·y/L
+    Q = Manning(A, n, R, s)
+
+case -4:  /* 临界水深边界 */
+    Q = A · √(g·y) · 60
+```
+
+此外 `MD_Lake.cpp:48` 中 `down <= -4` 会触发入湖路由（`toLake = (-3 - down) - 1`），但 QHH 案例中 `-4` 对应 `toLake = 0`，需确认是否有湖泊模块参与。
+
+### B.4 影响范围
+
+- 受影响河段：45 个出口河段（全部出口）
+- 影响方式：出口流量计算公式从临界水深变为零水深梯度
+  - `-4` 临界水深：`Q = A√(gy)·60`，与坡度无关，受水深控制
+  - `-3` 零水深梯度：`Q = Manning(S_bed + 2y/L)`，受坡度和水深共同控制
+- 潜在风险：若 QHH 案例启用了湖泊模块，`-4` 还会触发入湖路由逻辑，改为 `-3` 后湖泊路由断开
+
+### B.5 建议
+
+1. **确认原始意图**：核实 2022 年将出口标记设为 `-4` 的物理依据（临界水深 vs 入湖路由）
+2. **恢复 `-4` 标记**：若确认是有意选择，在当前 baseline 的 `qhh.sp.riv` 中将 45 个出口河段的 Down 恢复为 `-4`
+3. **工具链防护**：在 `shud.river()` 或 AutoSHUD Step3 中增加参数，允许指定出口边界条件类型，避免重新生成时丢失手动设置
+4. **回归对比**：分别用 `-3` 和 `-4` 运行 QHH baseline，对比出口流量时序差异