diff --git a/gusto/solvers/linear_solvers.py b/gusto/solvers/linear_solvers.py
index 8214d2ba2..b482788ad 100644
--- a/gusto/solvers/linear_solvers.py
+++ b/gusto/solvers/linear_solvers.py
@@ -665,7 +665,7 @@ def _setup_solver(self):
 
         if 'coriolis' in equation.prescribed_fields._field_names:
             f = equation.prescribed_fields('coriolis')
-            eqn += beta_u_ * f * inner(w, equation.domain.perp(u)) * dx
+            eqn += beta_u * f * inner(w, equation.domain.perp(u)) * dx
 
         aeqn = lhs(eqn)
         Leqn = rhs(eqn)
@@ -881,7 +881,7 @@ def _setup_solver(self):
 
         if 'coriolis' in equation.prescribed_fields._field_names:
             f = equation.prescribed_fields('coriolis')
-            eqn += beta_u_ * f * inner(w, equation.domain.perp(u)) * dx
+            eqn += beta_u * f * inner(w, equation.domain.perp(u)) * dx
 
         aeqn = lhs(eqn)
         Leqn = rhs(eqn)