diff --git a/py_wake/deficit_models/fuga.py b/py_wake/deficit_models/fuga.py
index 0cbab4919bb665a3b136166734b991290870aa37..17ad3abb1c8cd35ee829a5c89312f34795376324 100644
--- a/py_wake/deficit_models/fuga.py
+++ b/py_wake/deficit_models/fuga.py
@@ -261,7 +261,9 @@ class FugaMultiLUTDeficit(XRLUTDeficitModel, FugaDeficit):
         if z_lst is None:
             z_lst = np.sort(np.unique([da.z for da in da_lst]))
         x_lst = np.sort(np.unique([da.x for da in da_lst]))
+        x_lst = x_lst[(x_lst >= np.max([da.x[0] for da in da_lst])) & (x_lst <= np.min([da.x[-1] for da in da_lst]))]
         y_lst = np.sort(np.unique([da.y for da in da_lst]))
+        y_lst = y_lst[(y_lst >= np.max([da.y[0] for da in da_lst])) & (y_lst <= np.min([da.y[-1] for da in da_lst]))]
         da_lst = [da.interp(z=z_lst, x=x_lst, y=y_lst) for da in da_lst]
 
         # combine_by_coords does not always merge attributes correctly