diff --git a/swiftest/sph_harmonics.py b/swiftest/sph_harmonics.py
index 25f74202a..0f1f0b8c9 100644
--- a/swiftest/sph_harmonics.py
+++ b/swiftest/sph_harmonics.py
@@ -82,7 +82,7 @@ def clm_from_ellipsoid(mass, density, a, b = None, c = None, lmax = 6, lref_radi
         if(lref_radius == True and ref_radius is None):
             ref_radius = shape_SH.expand(normalization = '4pi').coeffs[0, 0, 0]
             return clm, ref_radius
-        else if(lref_radius == True and ref_radius is not None):
+        elif(lref_radius == True and ref_radius is not None):
             clm = clm.change_ref(r0 = ref_radius)
             return clm, ref_radius
         else:
@@ -139,7 +139,7 @@ def clm_from_relief(mass, density, grid, lmax = 6, lref_radius = False, ref_radi
         if(lref_radius == True and ref_radius is None):
             ref_radius = shape_SH.expand(normalization = '4pi').coeffs[0, 0, 0]
             return clm, ref_radius
-        else if(lref_radius == True and ref_radius is not None):
+        elif(lref_radius == True and ref_radius is not None):
             clm = clm.change_ref(r0 = ref_radius)
             return clm, ref_radius
         else: