MODFLOW 6数据集的VTK输出？

Question

我知道Flopy可以从MF2005数据集(例如export.ipynb)导出VTK文件，但是我试图让这个示例处理MODFLOW 6数据集和一些代码编辑，但没有成功。

Flopy实际上是否支持从MODFLOW 6数据集导出VTK，还是这一功能还没有实现？

回应Sarah的回复：

我尝试使用vtk.export_heads函数，发现模型对象似乎期待着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)

...and我得到一个‘属性错误:详细’，因为在mf6.mf仿真中不存在详细信息(但是在从MF2005导出时使用的等效modflow.Modflow中)。我强迫详细的错误只是为了让代码进展，然后就出现了一个‘属性错误: modelGrid'，它看起来是一个更实质性和更复杂的缺失属性。

Answer 1

我还想导出vtk格式的modflow 6模型的结果。例程存在(vtk.export_heads、vtk.export_model等)。您可以检查这些页面：https://flopy.readthedocs.io/en/3.3.2/source/flopy.export.vtk.html https://flopy.readthedocs.io/en/3.3.2/source/flopy.export.utils.html

不幸的是，我面临着非常长的计算时间。我有一个两百万网格模型，它需要几个小时。现在，它是残废的

如果你想办法减少出口时间，请告诉我！

Answer 2

# 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)

看起来您正在导出模拟，尝试导出模拟，您可以使用如下一行：

vtk.export_model(loaded_sim.get_model(), head_output_folder)