diff --git a/docs/photonics/examples/calculate_GVD.py b/docs/photonics/examples/calculate_GVD.py
index 813e5284..d4abacac 100644
--- a/docs/photonics/examples/calculate_GVD.py
+++ b/docs/photonics/examples/calculate_GVD.py
@@ -98,7 +98,7 @@ def n_SiO2(wavelength):
 
 # %% [markdown]
 # Create the mesh, and sweep wavelength using the same mesh. The target mode is te mode, so the mode is selected by highest te fraction
-# %%
+# %% tags=["remove-stderr"]
 neff_list = []
 aeff_list = []
 basis0 = Basis(mesh, ElementTriP0())