软盘中压力期数据中的井 ID 和其他注释

Question

我想为我的井文件分配一个注释掉的井 ID 号（对 ghb 单元也这样做），但我找不到任何关于如何这样做的信息。我写了一些东西来创建我自己的 ghb 文件，但是如果我尝试将它加载回我的软盘 mf 类，然后用 mf.write_input() 将它写出我的其他包，注释不会保留，它会被覆盖。我知道在 mf.wrtie_input() 中我可以指定要写出的包，如果我拿走我之前制作的 ghb 文件（或井文件），那么原始文件不会被覆盖，这很好。但是我想知道是否有一种方法可以直接为每个包的 stress_period_data 添加注释，以便我可以将它们全部包含在 flopy 类中。

Answer 1

像在这个例子中一样，您可以扩展默认值dtype以包含MfList实例将携带到写入的额外属性：

well_dtype = [('k', '<i8'), ('i', '<i8'), ('j', '<i8'),('flux', '<f4'), ('wel_id', object)]

stress_period_data = np.zeros((3), dtype=well_dtype)

wel = flopy.modflow.ModflowWel(m, stress_period_data=stress_period_data, dtype=well_dtype)

我不确定加载wel具有额外属性的现有包的简单方法- 仅供参考