diff --git a/src/parcels/interpolators/_xinterpolators.py b/src/parcels/interpolators/_xinterpolators.py
index d20313ffb..388925d1a 100644
--- a/src/parcels/interpolators/_xinterpolators.py
+++ b/src/parcels/interpolators/_xinterpolators.py
@@ -455,22 +455,22 @@ def is_land(ti: int, zi: int, yi: int, xi: int):
             f_v,
         )
         f_u = np.where(
-            is_land(0, 0, 0, 0) & is_land(0, 0, 0, 1) & is_land(0, 0, 1, 0 & is_land(0, 0, 1, 1) & (zeta > 0)),
+            is_land(0, 0, 0, 0) & is_land(0, 0, 0, 1) & is_land(0, 0, 1, 0) & is_land(0, 0, 1, 1) & (zeta > 0),
             f_u * (a + b * zeta) / zeta,
             f_u,
         )
         f_u = np.where(
-            is_land(0, 1, 0, 0) & is_land(0, 1, 0, 1) & is_land(0, 1, 1, 0 & is_land(0, 1, 1, 1) & (zeta < 1)),
+            is_land(0, 1, 0, 0) & is_land(0, 1, 0, 1) & is_land(0, 1, 1, 0) & is_land(0, 1, 1, 1) & (zeta < 1),
             f_u * (1 - b * zeta) / (1 - zeta),
             f_u,
         )
         f_v = np.where(
-            is_land(0, 0, 0, 0) & is_land(0, 0, 0, 1) & is_land(0, 0, 1, 0 & is_land(0, 0, 1, 1) & (zeta > 0)),
+            is_land(0, 0, 0, 0) & is_land(0, 0, 0, 1) & is_land(0, 0, 1, 0) & is_land(0, 0, 1, 1) & (zeta > 0),
             f_v * (a + b * zeta) / zeta,
             f_v,
         )
         f_v = np.where(
-            is_land(0, 1, 0, 0) & is_land(0, 1, 0, 1) & is_land(0, 1, 1, 0 & is_land(0, 1, 1, 1) & (zeta < 1)),
+            is_land(0, 1, 0, 0) & is_land(0, 1, 0, 1) & is_land(0, 1, 1, 0) & is_land(0, 1, 1, 1) & (zeta < 1),
             f_v * (1 - b * zeta) / (1 - zeta),
             f_v,
         )