0

我是 FloPy 的初学者，我正在使用 MODFLOW-USG 中的不饱和模型（使用理查兹方程），并且我有输出文件（.hds、.ddn、.cbb）。我需要导出不饱和单元格的饱和度（0 到 1 之间的参数），但我不知道该怎么做。我在一些论文中看到了参数“含水量”、“水饱和度”、“饱和度”，但我不知道如何得到它......

我一直在尝试从输出文件 .DDN 中使用代码在 FloPy 中对其进行后处理来获取饱和度：

flopy.utils.binaryfile.HeadFile(f"DRPA_TR_oc_.ddn",text="SATURATIONU", precision="auto")

和

flopy.utils.binaryfile.HeadFile(f"DRPA_TR_oc_.ddn",text="SATURATION", precision="auto")

但经过约 1 小时的计算，出现以下情况：

File "C:\Users\User1\Anaconda3\lib\site-packages\flopy\utils\binaryfile.py", line 369, in _get_header
    return header[0]
IndexError: index 0 is out of bounds for axis 0 with size 0

当我使用 .cbb 运行 CellBudgetFile 或使用 .hds 运行 HeadUFile 时，需要几分钟。所以，我不太明白错误在哪里以及如何修复它，以及为什么计算 .ddn 文件需要很多时间。

我希望您能帮助我修复错误或以其他方式获得饱和度（可能来自 .hds 或其他）:) 我会继续关注您的评论！提前致谢 ：）

4

0 回答 0