diff --git a/py_wake/wind_turbines/_wind_turbines.py b/py_wake/wind_turbines/_wind_turbines.py
index be624ee0633b912729e897144f60210145a1a87a..5762170076512cd8bd5a0de05f93e148e9bc0aec 100644
--- a/py_wake/wind_turbines/_wind_turbines.py
+++ b/py_wake/wind_turbines/_wind_turbines.py
@@ -321,7 +321,7 @@ Use WindTurbines(names, diameters, hub_heights, power_ct_funcs) instead""", Depr
             'mode', [PowerCtTabular(wt['WindSpeed'], wt['PowerOutput'], power_unit, wt['ThrustCoEfficient'],
                                     ws_cutin=wt['LowSpeedCutIn'], ws_cutout=wt['HighSpeedCutOut'],
                                     ct_idle=wt['ct_idle'], additional_models=[]) for wt in wt_data],
-            default_value=0, additional_models=[SimpleYawModel()])
+            default_value=default_mode, additional_models=[SimpleYawModel()])
 
         char_data_tables = [np.array([pct.ws_tab, pct.power_ct_tab[0], pct.power_ct_tab[1]]).T
                             for pct in power_ct_funcs.windTurbineFunction_lst]
diff --git a/py_wake/wind_turbines/power_ct_functions.py b/py_wake/wind_turbines/power_ct_functions.py
index c6cd9aee40194593112db2460a381bba7357266f..49c09010abc13a611b029eea45d19513bed1190e 100644
--- a/py_wake/wind_turbines/power_ct_functions.py
+++ b/py_wake/wind_turbines/power_ct_functions.py
@@ -231,6 +231,7 @@ class PowerCtTabular(PowerCtFunction):
             ws = np.r_[ws, ws_cutout + eps, 100]
             power = np.r_[power, power_idle, power_idle]
             ct = np.r_[ct, ct_idle, ct_idle]
+        self.ws_cutin, self.ws_cutout, self.ct_idle = ws_cutin, ws_cutout, ct_idle
         idx = np.argsort(ws)
         ws, power, ct = np.asarray(ws)[idx], np.asarray(power)[idx], np.asarray(ct)[idx]