From 3e715ad7362805d6aac0f30435aab50e713d8ed6 Mon Sep 17 00:00:00 2001 From: jshipton Date: Mon, 19 Aug 2024 19:39:52 +0100 Subject: [PATCH] use correct beta_u in shallow water linear solvers --- gusto/solvers/linear_solvers.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) 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)