not quite sure how this worked without the definition of the ke_form

firedrakeproject · Aug 26, 2023 · 2d9807d · 2d9807d
1 parent dda98e9
commit 2d9807d
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/gusto/equations.py b/gusto/equations.py
@@ -682,6 +682,7 @@ def __init__(self, domain, parameters, fexpr=None, bexpr=None,
         elif u_transport_option == "vector_advection_form":
             u_adv = prognostic(advection_form(w, u, u), 'u')
             if domain.move_mesh:
+                ke_form = prognostic(kinetic_energy_form(w, u, u), "u")
                 ke_form = transport.remove(ke_form)
                 ke_form = ke_form.label_map(
                     lambda t: t.has_label(transporting_velocity),