问题标签 [flopy]

我正在尝试学习 flopy，并尝试加载 MODFLOW-USG 测试用例之一。测试用例是一个带有嵌套网格和 GNC 包的 1 层模型（示例 01B）。我没有修改与 MODFLOW-USG 一起分发的测试用例文件，并且测试用例在 run.bat 文件中运行良好。

首先，我尝试按照示例使用flopy.modflow.mf 模块加载名称文件：

但我得到了以下信息ValueError：

然后我尝试按照flopy.modflow.mfdisu 模块的示例加载 DISU 包：

但这给了我以下信息TypeError：

用软盘加载现有 MODFLOW-USG 模型的正确方法是什么？

我在形状文件中有河流的线要素类，并且在线的开始和结束处有河流包的属性。例如：线要素类起点和终点的水位。

现在，我想在网格上导入河流包属性，以便对河流线重叠的所有节点的值进行插值。谁能帮我在 gis、flopy 或 python 中的任何其他模块中做到这一点？或者帮我解决这个问题？

问题如灰色箭头所示：

我正在设置我的模型空间参考如下：

我相信属性设置正确，因为我导出的 geotiff 栅格旋转正确。但是，我的 get_rc() 结果并不正确。我通过的 x,y 对位于单元格中心，未旋转网格的结果非常完美。

我正在使用的方法调用（伪代码）是：

其中 X 和 Y 是 numpy 数组。

我可以使用 flopy.mf6.MFSimulation.load 读取现有的 MODFLOW 6 模拟。现在我想找出它有多少个压力周期，作为一个整数，由 tdis 包中的 nper 定义。最简单的方法是什么？

几周以来，我一直在用 Python 中的 Flopy 编写 MODFLOW 模型。由于 Python 的透明性，我选择用 Flopy 编写模型。但是，有时我的模型没有运行，但它没有告诉我哪里出了问题，这使得错误处理变得困难。

目前我的模型在运行模型时出错。它使用我手动添加的消息引发错误（并且很常用）：

我得到的错误是：

此外，所有 MODFLOW 文件都已创建，但 .hds 和 .cbc 文件包含零字节。

我的问题：是否有人有技巧可以有效而智能地追踪这些错误。

我建立了一个耦合的 MODFLOW-MT3DMS-SEAWAT 模型。根据软件包和可执行文件 swt_v4 与此示例类似（当然我使用了更多的输入变量）。

但是，当我尝试运行模型时，出现以下错误：

如果我运行该命令print(swt.namfile)，通常会提示为swt_model_0605.nam

我使用以下软件包：DIS、BAS、LPF、PCG、OC、GHB、WEL、ADV、BTN、DSP、GCG、SSM 和 VDF。

有没有办法追踪这个错误？

我正在尝试在“FloPy”（v'3.2.12'）中设置模型，但我在处理 modflow 输出文件时遇到了一些问题。当关注开发者网站 ( https://modflowpy.github.io/flopydoc/tutorial2.html ) 上的“turorial2”时，似乎一切正常。但是，cbc 文件不包含有关一般头部边界条件的信息。我在排水边界条件方面也遇到了同样的问题。在运行模型之前，我需要更改或设置一些其他设置吗？还是有另一种方法可以从 GHB 边界条件中读取速率信息？

仅供参考-到目前为止，我已经尝试更改 ModflowOc 设置，并且还使用不同的求解器运行模型。

谢谢您的帮助，

有关设置的一些附加信息（应与教程相同）：

我正在使用 FloPy 加载现有的 MODFLOW-NWT 模型。UPW 软件包依赖于外部文件来获取水力特性。我注意到无法使用常规加载命令加载 UPW。

我也试过单独加载UPW，还是不行。我在存储库中找不到此示例。

我有一个带输出的工作模型。我试图利用 flopy 来询问结果。具体来说，我希望能够（即时）为特定单元生成区域预算信息。我当然可以不厌其烦地构建我感兴趣的区域文件，但我希望能够根据需要提取结果，而无需经过中间步骤。

我正在寻找一个可以在 Python 中使用的工作流，以仅在 Python 控制台中为特定单元格从 CBB（已生成）中提取净流量信息（宁愿避免生成区域文件、导入然后提取的传统样式） .

编辑

试图让该实用程序完全正常工作时遇到了问题。这是最新的尝试。无法判断 flopy 是否在 CBB 本身上窒息，或者我提供的区域是否给它带来了问题。

该模型为 7 层，196 行，241 列。我正在尝试提取 Layer=3、Row=58、Col=30。我为区域创建的列表对象在这里：

然后我使用以下内容将其写入文件：

然后我为区域预算类/方法构建我的 modflow 模型：

所有这些都运行良好，直到最后一个命令，我得到以下错误：

请记住，我仍然对跳过文件写入部分感兴趣。我把它包括在内是为了展示我到目前为止的工作

我尝试在软盘中加载、更新和写入新的输入文件。我尝试了很多事情，但我做不到。这是我的代码：

modflowModel是和flopy.modflow.Modflow对象。代码末尾的注释是我尝试编写更新的新输入但不起作用的行。