isort

Author: Joost Delsman

docs/api/changelog.rst

@@ -19,6 +19,8 @@ Fixed
 Changed
 ~~~~~~~
 
+- Removed ``meshzoo`` dependency.
+
 Added
 ~~~~~
 

examples/mf6/Henry.py

@@ -26,6 +26,9 @@
 with hydrostatic seawater with a density of 1025 kg m:sup:`-3`. This is
 represented by a General Head Boundary package.
 """
+
+# sphinx_gallery_thumbnail_number = -1
+
 # %%
 # We'll start with the usual imports. As this is a simple (synthetic)
 # structured model, we can make due with few packages.
@@ -264,11 +267,11 @@
 #
 # We'll open the head and concentration files.
 
-head = imod.mf6.open_hds(
+head = imod.mf6.out.open_hds(
     modeldir / "flow/flow.hds",
     modeldir / "flow/dis.dis.grb",
 )
-conc = imod.mf6.open_hds(
+conc = imod.mf6.out.open_conc(
     modeldir / "transport/transport.ucn",
     modeldir / "flow/dis.dis.grb",
 )
@@ -283,6 +286,6 @@
 # We can check the concentration to see that a fresh-saline interface has been
 # formed:
 
-conc.isel(y=0, time=-1).plot.contourf(yincrease=False)
+conc.isel(y=0, time=-1).plot.contourf(yincrease=False, cmap="RdYlBu_r")
 
 # %%

examples/mf6/circle.py

@@ -22,7 +22,6 @@
 # We'll start with the following imports:
 
 import matplotlib.pyplot as plt
-import meshzoo
 import numpy as np
 import xarray as xr
 import xugrid as xu
@@ -33,18 +32,19 @@
 # Create a mesh
 # -------------
 #
-# The first steps consists of generating a mesh. In this example, our geometry
-# is so simple that we don't need complex mesh generators. Instead we will use
-# ``meshzoo``, which can generate a number of meshes for simple domains such as
-# triangles, rectangles, regular polygons, or disks.
-#
-# We can use the nodes and the cell node connectivity directly to initialize a
-# xugrid Ugrid2d object, and we'll create a quick plot to take a look at the
-# grid.
+# The first steps consists of generating a mesh. In this example, we'll use data
+# included with iMOD Python for a circular mesh. Note that this is a `Ugrid2D
+# object. <https://deltares.github.io/xugrid/api/xugrid.Ugrid2d.html>`_
+# For more information on working with unstructured grids see the
+# `Xugrid documentation <https://deltares.github.io/xugrid/index.html>`_
+
+grid = imod.data.circle()
+
+grid
 
-nodes, triangles = meshzoo.disk(6, 6)
-nodes *= 1000.0
-grid = xu.Ugrid2d(*nodes.T, -1, triangles)
+
+# %%
+# We can plot this object as follows:
 
 fig, ax = plt.subplots()
 xu.plot.line(grid, ax=ax)
@@ -68,7 +68,7 @@
 # groundwater; with small flows in the center and larger flows as the recharge
 # accumulates while the groundwater flows towards the exterior boundary.
 
-nface = len(triangles)
+nface = grid.n_face
 nlayer = 2
 
 idomain = xu.UgridDataArray(