Compute mesh data attributes on demand

Separate the notion of Nlr (Eucclhyd normals) of the calculation of
Cjr in 3D. It remain to implement its calculation in 2D to finish this
part.

This should allow a trivial implementation of the Eucclhyd scheme and
simplify coding off pressure boundary conditions.