使用Flopy创建vtk输出

Question

我正在尝试将模型网格从Modflow模型写入vtk输出，以在Paraview中打开。

我使用了地下水远景，并创建了MODFLOW-USG（非结构化网格）

当我尝试使用导出功能时，出现以下错误：

追踪（最近通话）：

文件“ C：\ Program Files \ JetBrains \ PyCharm Community Edition 2019.2.5 \ helpers \ pydev \ pydevd.py”，行1415，在_exec中 pydev_imports.execfile（file，globals，locals）＃执行脚本

文件“ C：\ Program Files \ JetBrains \ PyCharm Community Edition 2019.2.5 \ helpers \ pydev_pydev_imps_pydev_execfile.py”，第18行，在execfile中 exec（compile（contents +“ \ n”，file，'exec'），glob，loc）

文件“ C：/ Work / Models / Work / Projects / Job 302 / Export grid.py”，第16行，在 vtkobj = fp.export.vtk.Vtk（mf）

文件“ C：\ ProgramData \ Anaconda3 \ lib \ site-packages \ flopy \ export \ vtk.py”，第206行，在init中 self.shape =（self.modelgrid.nlay，self.modelgrid.nrow，

AttributeError：“网格”对象没有属性“ nlay”

下面是我使用的代码：

import os
import flopy as fp
from flopy.export import vtk as fv

model_path = r'C:\Work\Models\302'
os.chdir(model_path)

# MODFLOW name file:
namefile = 'westend_Alluv3_11yr_A.nam'

mf = fp.modflow.Modflow.load(namefile,verbose=False,check=False)

# create the vtk object
vtkobj = fp.export.vtk.Vtk(mf)

谢谢