diff --git a/samples/python/kinetics/jet_stirred_reactor.py b/samples/python/kinetics/jet_stirred_reactor.py
index 1eac5505d0..c9692cd74f 100644
--- a/samples/python/kinetics/jet_stirred_reactor.py
+++ b/samples/python/kinetics/jet_stirred_reactor.py
@@ -41,7 +41,7 @@
 
 f, ax = plt.subplots(1, 3)
 plt.subplots_adjust(wspace=0.6)
-colours = ["xkcd:grey",'xkcd:purple']
+colors = ["xkcd:grey",'xkcd:purple']
 file = 'example_data/ammonia-CO-H2-Alzueta-2023.yaml'
 models = {'Original': 'baseline', 'LMR-R': 'linear-Burke'}
 
@@ -108,9 +108,9 @@ def getTemperatureDependence(gas, inputs):
     tempDependence = getTemperatureDependence(gas,inputs)
     ax[0].plot(tempDependence.index,
                np.subtract(tempDependence['temperature'],tempDependence.index),
-               color=colours[k],label=m)
-    ax[1].plot(tempDependence.index, tempDependence['O2']*100, color=colours[k])
-    ax[2].plot(tempDependence.index, tempDependence['H2']*100, color=colours[k])
+               color=colors[k],label=m)
+    ax[1].plot(tempDependence.index, tempDependence['O2']*100, color=colors[k])
+    ax[2].plot(tempDependence.index, tempDependence['H2']*100, color=colors[k])
 ax[0].plot(inputs['data']['T_range'], inputs['data']['deltaT'], 'o', fillstyle='none',
            color='k', label="Sabia et al.")
 ax[1].plot(inputs['data']['T_range'], inputs['data']['X_O2'], 'o', fillstyle='none',
diff --git a/samples/python/kinetics/shock_tube.py b/samples/python/kinetics/shock_tube.py
index dfc9c7a171..89ae072e6d 100644
--- a/samples/python/kinetics/shock_tube.py
+++ b/samples/python/kinetics/shock_tube.py
@@ -38,7 +38,7 @@
 
 file = 'example_data/ammonia-CO-H2-Alzueta-2023.yaml'
 models = {'Original': 'baseline', 'LMR-R': 'linear-Burke'}
-colours = {'Original': "xkcd:grey", 'LMR-R': 'xkcd:purple'}
+colors = {'Original': "xkcd:grey", 'LMR-R': 'xkcd:purple'}
 results = {}
 
 # %%
@@ -80,7 +80,7 @@
 fig, ax = plt.subplots()
 
 for m, timeHistory in results.items():
-    ax.plot(timeHistory.t*1e6, timeHistory('H2O').X*100, color=colours[m], label=m)
+    ax.plot(timeHistory.t*1e6, timeHistory('H2O').X*100, color=colors[m], label=m)
 
 ax.plot(expData['t'], np.array(expData['X_H2O'])*100, 'o', fillstyle='none', color='k',
         label='Shao et al.')
diff --git a/samples/python/onedim/flame_speed.py b/samples/python/onedim/flame_speed.py
index 4855995b6d..991b1a61c6 100644
--- a/samples/python/onedim/flame_speed.py
+++ b/samples/python/onedim/flame_speed.py
@@ -39,7 +39,7 @@
 
 file = 'example_data/ammonia-CO-H2-Alzueta-2023.yaml'
 models = {'Original': 'baseline', 'LMR-R': 'linear-Burke'}
-colours = {'Original': "xkcd:grey", 'LMR-R': 'xkcd:purple'}
+colors = {'Original': "xkcd:grey", 'LMR-R': 'xkcd:purple'}
 results = {}
 
 Tin = 296  # unburned gas temperature [K]
@@ -76,7 +76,7 @@
 
 fig, ax = plt.subplots()
 for m, vel_list in results.items():
-    ax.plot(phi_list, vel_list, color=colours[m], label=m)
+    ax.plot(phi_list, vel_list, color=colors[m], label=m)
 ax.plot(phi_data, expData['vel'], 'o', fillstyle='none', color='k', label='Ronney')
 ax.legend(frameon=False, loc='upper right')
 ax.set_ylabel(r'Burning velocity [cm $\rm s^{-1}$]')