问题标签 [flopy]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
9 浏览

flopy - 如何在不饱和模型 MODFLOW-USG 中获得饱和?

我是 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 小时的计算,出现以下情况:

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

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