diff --git a/dumux/porenetwork/util/extract_pore_network_with_porespy.py b/dumux/porenetwork/util/extract_pore_network_with_porespy.py
index 872f3c1d83a1ac783e3feaf2f01e401e7150cd94..2800855233321148a88dd4d966dd1de2997201e3 100755
--- a/dumux/porenetwork/util/extract_pore_network_with_porespy.py
+++ b/dumux/porenetwork/util/extract_pore_network_with_porespy.py
@@ -202,11 +202,13 @@ def extractNetwork(
     # use PoreSpy to sanitize some parameters (might become obsolete as some point)
     porenetwork = op.io.network_from_porespy(snowOutput.network)  # porenetwork in OpenPNM format
 
-    # trimming pore network to avoid singularity (remove isolated and disconnected pores)
+    # trimming pore network to avoid singularity (remove pores in disconnected clusters,
+    # this includes isolated pores).
+    # isolated pores could be removed separately by using
+    # op.topotools.trim(network=porenetwork, pores=health["isolated_pores"])
     print("Number of pores before trimming: ", porenetwork.Np)
     health = op.utils.check_network_health(porenetwork)
     op.topotools.trim(network=porenetwork, pores=health["disconnected_pores"])
-    op.topotools.trim(network=porenetwork, pores=health["isolated_pores"])
     print("Number of pores after trimming: ", porenetwork.Np)
 
     filename = outputName + ("_dual" if dualSnow else "")