Welford algorithm extension

This MR extends the Welford implementation to support the third-order central moments of a (vector) field.

It also adds a basic test case for different setups of the Welford (different spatial dimensions and different orders of validation).