Improved interpolation

Added a “method” argument to interpolation via grid expressions such as
gr1(g2,method = ‘linear’). Previously, only the default (linear) was
available.