import numpy as np
import flopy
# open the binary file
sobj = flopy.utils.HeadFile('model.zdisplacement.bin',
text='Z DISPLACEMENT')
# get all of the available times in the file
times = sobj.get_times()
# extract the data for the last time in the file
zd = sobj.get_data(totim=times[-1])
# save the z-displacement for the first layer (layer 0) to an ascii file
# zd is a 3D numpy array with a shape of (nlay, nrow, ncol)
np.savetxt('layer0.zdisplacement.txt', zd[0])
如果有多个图层,则需要保存每个图层的数据。你知道吗
您可以使用以下命令输出文件中的所有数据:
for t in sobj.get_times():
zd = sobj.get_data(totim=t)
for k in range(nlay):
fpth = 'layer{}_{}.zdisplacement.txt'.format(k, t)
np.savetxt(fpth, zd[k])
子包的二进制输出与MODFLOW二进制头文件的格式相同。您需要知道写入二进制文件的输出文本字符串的名称。有关确定给定子包二进制文件的文本字符串的子包,请参见MODFLOW-2005 online documentation中的表1。你知道吗
下面演示如何使用
flopy
和numpy
将二进制文件中的Z DISPLACEMENT
数据转换为ascii文件:如果有多个图层,则需要保存每个图层的数据。你知道吗
您可以使用以下命令输出文件中的所有数据:
相关问题 更多 >
编程相关推荐