Replacing matplotlib.mlab.griddata by scipy.interpolate.griddata

As discussed in issue #146, I'd some problems when I tried to run `PrismMesh.carvetopo` because of `matplotlib.mlab.griddata`. So, @leouieda suggested to replace the `matplotlib.mlab.griddata` functions  by `scipy.interpolate.griddata`.

@leouieda I don't know what's the best interpolation method to set as default. With `scipy.interpolate.griddata` our options are `linear`, `nearest` and `cubic`. Here I'm using `cubic`.

## Checklist:

- [x] Make tests for new code
- [x] Create/update docstrings
- [x] Code follows PEP8 style conventions
- [x] Code and docs have been spellchecked
- [x] Changelog entry
- [x] Include new dependencies in docs, requirements.txt, README
- [x] Documentation builds properly
- [x] All tests pass
- [x] Can be merged