modflowpy · mjreno · Nov 9, 2025 · Nov 6, 2025 · Nov 7, 2025 · Nov 7, 2025
diff --git a/.gitignore b/.gitignore
@@ -1,5 +1,7 @@
 build/
 docs/examples/*.ipynb
+docs/examples/quickstart/*
+docs/examples/twri/*
 venv/
 .DS_Store
 .ruff_cache

diff --git a/docs/examples/advanced_packages.py b/docs/examples/advanced_packages.py
diff --git a/docs/examples/image/qs.png b/docs/examples/image/qs.png
diff --git a/docs/examples/image/qsprj.png b/docs/examples/image/qsprj.png
diff --git a/docs/examples/image/quickstart.png b/docs/examples/image/quickstart.png
diff --git a/docs/examples/qsplot.py b/docs/examples/qsplot.py
diff --git a/docs/examples/quickstart.py b/docs/examples/quickstart.py
@@ -11,6 +11,7 @@
 
 name = "quickstart"
 workspace = Path(__file__).parent / name
+workspace.mkdir(exist_ok=True)
 time = Time(perlen=[1.0], nstp=[1])
 grid = StructuredGrid(
     nlay=1,
@@ -25,7 +26,7 @@
 gwf_name = "mymodel"
 ims = Ims(parent=sim, models=[gwf_name])  # temporary hack
 gwf = Gwf(parent=sim, name=gwf_name, save_flows=True, dis=grid)
-npf = Npf(parent=gwf, save_specific_discharge=True)
+npf = Npf(parent=gwf, print_flows=True, save_flows=True, save_specific_discharge=True)
 chd = Chd(
     parent=gwf,
     head={0: {(0, 0, 0): 1.0, (0, 9, 9): 0.0}},

diff --git a/docs/examples/quickstart/mfsim.lst b/docs/examples/quickstart/mfsim.lst