diff --git a/src/darsia/measure/wasserstein.py b/src/darsia/measure/wasserstein.py
index 07d8a753..1c95e25f 100644
--- a/src/darsia/measure/wasserstein.py
+++ b/src/darsia/measure/wasserstein.py
@@ -896,7 +896,7 @@ def _compute_face_weight(self, flat_flux: np.ndarray) -> np.ndarray:
 
             # Combine weights**2 / |weight * flux| on faces
             face_weights = harm_avg_face_weights**2 / norm_weighted_face_flux
-            face_weights_inv = 1.0 / face_weights
+            face_weights_inv = norm_weighted_face_flux / harm_avg_face_weights**2
         else:
             raise ValueError(f"Mobility mode {self.mobility_mode} not supported.")