diff --git a/py_wake/rotor_avg_models/gaussian_overlap_model.py b/py_wake/rotor_avg_models/gaussian_overlap_model.py
index 436424fc7251922a3b77cb138d4c68496e8c5e74..1c0d574ae960ff5e481a2609101192e65a84a5ff 100644
--- a/py_wake/rotor_avg_models/gaussian_overlap_model.py
+++ b/py_wake/rotor_avg_models/gaussian_overlap_model.py
@@ -9,7 +9,7 @@ from py_wake.utils.grid_interpolator import GridInterpolator
 
 class GaussianOverlapAvgModel(RotorAvgModel):
     def __init__(self, filename=os.path.dirname(__file__) + f'/gaussian_overlap_.02_.02_128_512.nc'):
-        table = xr.load_dataarray(filename)
+        table = xr.load_dataarray(filename, engine='h5netcdf')
         R_sigma = np.arange(0, 20.001, 0.01)
         CW_sigma = np.arange(0, 10.01, 0.01)
         dat = table.interp(R_sigma=R_sigma, CW_sigma=CW_sigma, method='cubic')