从MODFLOW-USG二进制文件中用flopy提取头错误

Question

我正在尝试从MODFLOW二进制头保存文件中提取模拟头，因此我使用flopy如下所示：

modelname = 'Model'

import flopy.utils.binaryfile as bf

# Create the headfile object
headobj = bf.HeadFile(modelname+'.hds')
head = headobj.get_data()

但是，对于最后一行，我得到了以下错误：

ValueError: cannot reshape array of size 17349 into shape (17338,8670)

记录列表是

(1, 1, 1., 1., b'           HEADU', 1, 8669, 1)
(1, 1, 1., 1., b'           HEADU', 8670, 17338, 2)

什么是正确的方式，从一个现有的MODFLOW模型装载头部与flopy？

Answer 1

如果我将text='headu'传递给HeadFile构造函数，那么我可以读取相当简单的MODFLOW模型的头保存文件：

headobj = bf.HeadFile(modelname+'.hds',text='headu')
head = headobj.get_data()

Answer 2

hdobj = bf.HeadUFile(modelname+'.hds')
head = hdobj.get_data()