Fix cornercases in Lagrange-interpolation subgrids

[cschwan]

We should not allow incorrectly set parameters, which will lead to problems similar to the one reported here: #75 (comment)

[cschwan]

This has been fixed in #299 by the introduction of the struct Interp, whose constructor carefully checks the values.