Interpolation routines for loading fields #11
Description
A priori, the isostasy and ice domain can be distinct. To handle this, we need efficient interpolation routines. @alex-robinson has something up his sleeve for this but we still need to include it nicely to the code. This might also be helpful for loading fields and reduce the size of IsostasyData.