The discrete form of the total energy conservation, (C.3), is given by:
Substituting the discrete expression of the time derivative of the velocity either in vector invariant or in flux form, leads to the discrete equivalent of the
In flux from the vorticity term reduces to a Coriolis term in which the Coriolis parameter has been modified to account for the “metric” term. This altered Coriolis parameter is discretised at an f-point. It is given by:
Either the ENE or EEN scheme is then applied to obtain the vorticity term in flux form. It therefore conserves the total KE. The derivation is the same as for the vorticity term in the vector invariant form (§C.3.2).
The flux form operator of the momentum advection is evaluated using a centered second order finite difference scheme. Because of the flux form, the discrete operator does not contribute to the global budget of linear momentum. Because of the centered second order scheme, it conserves the horizontal kinetic energy, that is :
Let us first consider the first term of the scalar product ( just the the terms associated with the i-component of the advection) :
When the UBS scheme is used to evaluate the flux form momentum advection, the discrete operator does not contribute to the global budget of linear momentum (flux form). The horizontal kinetic energy is not conserved, but forced to decay ( the scheme is diffusive).
Gurvan Madec and the NEMO Team
NEMO European Consortium2017-02-17