VTK使用Flopy导出MODFLOW 6数据集?

VTK使用Flopy导出MODFLOW 6数据集?,vtk,flopy,Vtk,Flopy,我知道Flopy可以从MF2005数据集导出VTK文件(例如),但是我尝试让该示例使用MODFLOW 6数据集和一些代码编辑,但没有成功 Flopy是否真的支持MODFLOW 6数据集的VTK导出,或者该功能还不可用 对于莎拉·勒雷的回答: 我尝试使用vtk.export_heads函数,发现model对象似乎需要一些MF6公式中不存在的属性。你有没有面对这个问题?如果有,你是如何解决的?修改源代码 import os import sys try: import flopy excep

我知道Flopy可以从MF2005数据集导出VTK文件(例如),但是我尝试让该示例使用MODFLOW 6数据集和一些代码编辑,但没有成功

Flopy是否真的支持MODFLOW 6数据集的VTK导出,或者该功能还不可用

对于莎拉·勒雷的回答:
我尝试使用vtk.export_heads函数,发现model对象似乎需要一些MF6公式中不存在的属性。你有没有面对这个问题?如果有,你是如何解决的?修改源代码

import os
import sys
try:
    import flopy
except:
    fpth = os.path.abspath(os.path.join('..', '..'))
    sys.path.append(fpth)
    import flopy

from flopy.export import vtk

sim_name = "example_model"
sim_path = os.path.join("data", "flopy_mf6_tutorial")
heads_file_name = sim_name + ".hds"
heads_file = os.path.join(sim_path, heads_file_name)
head_output_folder = os.path.join(sim_path, "heads")

# load the simulation
loaded_sim = flopy.mf6.MFSimulation.load(sim_name, 'mf6', 'mf6', sim_path)

#export VTK heads
vtk.export_heads(loaded_sim, heads_file, head_output_folder, smooth=False, kstpkper=[(0,0)], point_scalars=False, nanval=-999.99)

…我得到一个“attribute error:verbose”,因为verbose在mf6.MFSimulation中不存在(但在从MF2005等导出时使用的等效modflow.modflow中存在)。我强制verbose为False只是为了让代码继续进行,然后出现了一个“属性错误:modelGrid”,它看起来像是一个更重要、更复杂的缺失属性。

我还想以vtk格式导出modflow 6模型的结果。存在例程(vtk.export\u头、vtk.export\u模型等)。您可以查看以下页面:

不幸的是,我面临着非常长的计算时间。我有一个两百万网格的模型,需要几个小时。就目前而言,这是一个严重的问题


如果您找到减少导出时间的方法,请告诉我

谢谢莎拉,非常感谢你的回复!我尝试使用vtk.export_heads函数,发现model对象似乎需要一些MF6公式中不存在的属性。你有没有面对这个问题?如果有,你是如何解决的?修改源代码?我已经编辑了我的问题,包括与此相关的详细信息!