Fix the interpolation function
Fixed a bug in the diffusive flux
Added the vertical dimension to the 2nd order interpolation function
Keep the interpolations horizontal and let the statistical routines do
the vertical interpolation