Flopy Inactive Cell in DISU grid

[Javier González]

Hi,

Is there a quick command to plot the inactive cells of a MODFLOW-USG unstructured grid in a plan view figure?

I was looking at this notebook but seems that the plot_inactive()  method used for MF-6 does not work with USG unstructured grids.

https://github.com/modflowpy/flopy/blob/develop/examples/Notebooks/flopy3.3_PlotMapView.ipynb

I am doing this, after I create my UnstructuredGrid Object (ug)

mview = flopy.plot.PlotMapView(modelgrid = ug)

mview.plot_grid(color='k', lw=0.5)

mview.plot_inactive()

It tells me :

AssertionError: Ibound/Idomain array must be provided

I've tried to get the basic package from the model with bas = model.get_package('bas') but it gives me a NoneType object.

I am a bit lost.

Javier.

[Javier González]

Got it. it is a BAS6 package.

This worked:

bas6  = model.get_package('bas6')
ibound_array = bas6.ibound.array
mview.plot_inactive(ibound=ibound_array)