Multidimensional linear interpolation with modern Fortran
