使用软盘从MODFLOW-USG二进制文件中提取磁头错误

Question

我正试图从MODFLOW-USG二进制头部保存文件中提取模拟头部，因此我使用软盘如下:

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-USG模型加载磁头的正确方法是什么?

Answer 1

如果我将text='headu'传递给HeadFile构造函数，那么我可以读取头部保存文件，以获得一个非常简单的MODFLOW-USG模型:

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

Answer 2

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