enable derivatives for 2d surfaces

[mkapur]

Ported from comments on website -- folks often do GAMs with a 2d smoother (i.e. lat an long). The Deriv function (outside of gratia) currently "handles" this by returning the same values for both smoother components. @gavinsimpson indicated developments to derivatives() function should include a tidy-based interface for this.