Fix floating point error in 2d/3d coordinate interpolation

[w-k-jones]

Resolves #433. Clipping of hdim1/2 values to size of dimensions (or size + 1 for periodic coordinates) in feature_location, and extrapolation enabled for 2d/3d coordinates as with 1d coords.

• [x] Have you followed our guidelines in CONTRIBUTING.md?
• [x] Have you self-reviewed your code and corrected any misspellings?
• [x] Have you written documentation that is easy to understand?
• [x] Have you written descriptive commit messages?
• [x] Have you added NumPy docstrings for newly added functions?
• [x] Have you formatted your code using black?
• [ ] If you have introduced a new functionality, have you added adequate unit tests?
• [x] Have all tests passed in your local clone?
• [ ] If you have introduced a new functionality, have you added an example notebook?
• [x] Have you kept your pull request small and limited so that it is easy to review?
• [x] Have the newest changes from this branch been merged?

[github-actions[bot]]

Linting results by Pylint:

Your code has been rated at 8.70/10 (previous run: 8.70/10, +0.00) _{The linting score is an indicator that reflects how well your code version follows Pylint’s coding standards and quality metrics with respect to the RC_v1.5.x branch. A decrease usually indicates your new code does not fully meet style guidelines or has potential errors.}

[freemansw1]

@w-k-jones, are you happy to merge? Only issue seems to be codecov.