diff --git a/docs/examples/howto_conversions/plot_B_feflowlib_BC_mesh.py b/docs/examples/howto_conversions/plot_B_feflowlib_BC_mesh.py
index 3f7395f943f344a9ab4d7845df7b3deda4515a22..37178fc06a6708f6095c92945ff25670aa237cba 100644
--- a/docs/examples/howto_conversions/plot_B_feflowlib_BC_mesh.py
+++ b/docs/examples/howto_conversions/plot_B_feflowlib_BC_mesh.py
@@ -48,7 +48,7 @@ ogs_sim_res = ms.mesh(ms.timesteps[-1])
 ot.plot.contourf(ogs_sim_res, ot.variables.temperature)
 # %%
 # 6. The boundary meshes are manipulated to alter the model.
-# The original boundary conditions are shown in this example: :ref:`sphx_glr_auto_examples_howto_conversions_plot_F_HT_simulation.py`
+# The original boundary conditions are shown in this example: :ref:`sphx_glr_auto_examples_howto_conversions_plot_F_feflowlib_HT_simulation.py`
 # 6.1 The Dirichlet boundary conditions for the hydraulic head are set to 0. Therefore, no water flows from the left to the right edge.
 bc_flow = feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"]
 feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"][bc_flow == 10] = 0
diff --git a/docs/examples/howto_conversions/plot_F_feflowlib_HT_simulation.py b/docs/examples/howto_conversions/plot_F_feflowlib_HT_simulation.py
index 7928cf0ebcd869babdde81caeb77d651e0e47db1..836fb8ef53c567e80ebb25c99cb5d5676cb229b8 100644
--- a/docs/examples/howto_conversions/plot_F_feflowlib_HT_simulation.py
+++ b/docs/examples/howto_conversions/plot_F_feflowlib_HT_simulation.py
@@ -75,6 +75,7 @@ hydraulic_head_diff = ot.variables.Scalar(
     data_name="HEAD_difference", data_unit="m", output_unit="m"
 )
 ot.plot.contourf(diff_mesh, hydraulic_head_diff, vmin=-1.5e-9, vmax=1.5e-9)
+# %%
 # Plot differences in temperature.
 feflow_model.mesh["temperature"] = feflow_model.mesh["P_TEMP"]
 # Plot differences in temperature.