水库

该单元在reservoir.f子程序中处理。在前进模式下，hup的值可以为负或为正。如果为负，则没有计算洪水曲线，需要确定回水曲线（hdo设置为-1）。如果为正，则计算了回水曲线，hdo使用hns.f子程序计算为水跃后的高度。

进一步的计算取决于hdo的值：

如果hdo超过水库深度hr，则回水曲线导致水库中的跌落或水跃，下游深度设置为hr。
如果hdo小于或等于hr，则计算临界深度hk和正常深度he。如果hk < he，则从ndo处的hk开始计算回水曲线（nstack增加，istack(1,nstack)设置为nelem，istack(2,istack)设置为ndo）。如果hk ≥ he，则回水曲线不可行，未找到解。