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



我知道Flopy可以从MF2005数据集导出VTK文件(例如flopy3_VTK_export.ipynb(,但我试图让该示例用于MODFLOW 6数据集和一些代码编辑,但没有成功。

Flopy真的支持从MODFLOW 6数据集导出VTK吗?还是该功能还不可用?

在回应Sarah Leray的回复时:
我尝试使用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)

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

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

不幸的是,我面临着非常长的计算时间。我有一个200万网格的模型,它需要几个小时。目前,它正在削弱

如果你能找到减少出口时间的方法,请告诉我!

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

最新更新