python-2.7 - 运行 flopy.modflow.HeadObservation 时出错:ValueError:不能从结构转换为非结构,除非结构只有一个字段
问题描述
我正在使用 Flopy 在 Python 2.7 中设置 MODFLOW 模型。我正在尝试通过 HOB 包添加头部观察。以下示例代码直接取自https://modflowpy.github.io/flopydoc/mfhob.html的函数文档:
import flopy
model = flopy.modflow.Modflow()
dis = flopy.modflow.ModflowDis(model, nlay=1, nrow=11, ncol=11,
nper=2, perlen=[1,1])
obs = flopy.modflow.mfhob.HeadObservation(model, layer=0, row=5,
column=5,
time_series_data=[[1.,54.4],
[2., 55.2]])
使用该函数的示例代码,我收到以下错误:
ValueError: Can't cast from structure to non-structure, except if the structure only has a single field.
当我尝试为我的模型创建头部观察时,我得到了同样的错误,这是稳态并且有一些不同的输入值。不幸的是,我无法找到一个可以与之比较的工作示例。有任何想法吗?
编辑:jdhughes 的代码就像一个魅力;但是我也忽略了将 Flopy 更新到最新版本 - 我尝试先更新 numpy,但直到我将 Flopy 从 3.2.8 更新到 3.2.9 才摆脱 ValueError。现在可以用了,谢谢!!!
解决方案
您需要创建 HeadObservation 类型的一个或多个实例并将其传递给ModflowHob
. 下面显示了具有两个观察位置的示例。
# create a new hob object
obs_data = []
# observation location 1
tsd = [[1., 1.], [87163., 2.], [348649., 3.],
[871621., 4.], [24439070., 5.], [24439072., 6.]]
names = ['o1.1', 'o1.2', 'o1.3', 'o1.4', 'o1.5', 'o1.6']
obs_data.append(flopy.modflow.HeadObservation(mf, layer=0, row=2, column=0,
time_series_data=tsd,
names=names, obsname='o1'))
# observation location 2
tsd = [[0., 126.938], [87163., 126.904], [871621., 126.382],
[871718.5943, 115.357], [871893.7713, 112.782]]
names = ['o2.1', 'o2.2', 'o2.3', 'o2.4', 'o2.5']
obs_data.append(flopy.modflow.HeadObservation(mf, layer=0, row=3, column=3,
time_series_data=tsd,
names=names, obsname='o2'))
hob = flopy.modflow.ModflowHob(mf, iuhobsv=51, obs_data=obs_data)
将提交问题以更新文档和文档字符串。
推荐阅读
- javascript - 网站在单个窗口中打开
- python - 复制数据框单元格并附加到所有行中数据框末尾的新列
- angular - Angular - 重新启动 CSS 动画
- magento - Magento 1.9 面包屑总是空的
- postgresql - PostgreSQL:硬 max_wal_size 或复制槽超时
- javascript - 谷歌图表时间线
- python - 使用带有请求和 lxml 的 xpath 抓取但有问题
- python - python多个和同时fifo
- vim - 如果
删除直到行尾。为什么不 拉到行尾? - windows - Capcom驱动漏洞:物理PC上的蓝屏,但虚拟机中没有