scipy.interpolate.griddata is slow

`export_gridded_model()` currently uses `scipy.interpolate.griddata`. However, this is slow for high-resolution gridding, for models with lots of elements, or for models with lots of timesteps.

Is there a better interpolation scheme for moving from and unstructured mesh to a structured grid?