Implementation of the convection-diffusion terms from the Navier-Stokes equations excluding the pressure related terms. It relies on tridiagonal solvers to obtain all the necessary derivative terms, including adding support for all possible options for the boundary conditions in the three spatial directions.