@@ -89,10 +89,13 @@ KineticEnergy(model; kwargs...) = KineticEnergy(model, model.velocities...; kwar
8989 u_immersed_bc,
9090 v_immersed_bc,
9191 w_immersed_bc,
92+ buoyancy,
93+ background_fields,
94+ velocities,
9295 args... )
93- u∂ₜu = ℑxᶜᵃᵃ (i, j, k, grid, ψf, velocities. u, u_velocity_tendency, advection, coriolis, stokes_drift, closure, u_immersed_bc, args... )
94- v∂ₜv = ℑyᵃᶜᵃ (i, j, k, grid, ψf, velocities. v, v_velocity_tendency, advection, coriolis, stokes_drift, closure, v_immersed_bc, args... )
95- w∂ₜw = ℑzᵃᵃᶜ (i, j, k, grid, ψf, velocities. w, w_velocity_tendency, advection, coriolis, stokes_drift, closure, w_immersed_bc, args... )
96+ u∂ₜu = ℑxᶜᵃᵃ (i, j, k, grid, ψf, velocities. u, u_velocity_tendency, advection, coriolis, stokes_drift, closure, u_immersed_bc, buoyancy, background_fields, velocities, args... )
97+ v∂ₜv = ℑyᵃᶜᵃ (i, j, k, grid, ψf, velocities. v, v_velocity_tendency, advection, coriolis, stokes_drift, closure, v_immersed_bc, buoyancy, background_fields, velocities, args... )
98+ w∂ₜw = ℑzᵃᵃᶜ (i, j, k, grid, ψf, velocities. w, w_velocity_tendency, advection, coriolis, stokes_drift, closure, w_immersed_bc, buoyancy, background_fields, velocities, args... )
9699 return u∂ₜu + v∂ₜv + w∂ₜw
97100end
98101
0 commit comments