fields package question

Thanks, Eric, for looking into that.
The values are below and since I subset the new abcissa  is smaller range
grid_new.l <- list(abcissa=c(-15.0,-14.),ordinate=y)
I am emailing form gmail - don't know why is using html to format when all
is in ascii

x
 [1] -15.20180 -15.01948 -14.86533 -14.73180 -14.61402 -14.50866 -14.41335
 [8] -14.32634 -14.24629 -14.17219
y
 [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
z
[5,] 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642
 [6,] 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143
 [7,] 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133
 [8,] 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326
 [9,] 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149
[10,] 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031
           [,8]      [,9]     [,10]
 [1,] 1.1900951 1.1900951 1.1900951
 [2,] 1.0636935 1.0636935 1.0636935
 [3,] 0.8927228 0.8927228 0.8927228
 [4,] 0.7554456 0.7554456 0.7554456
 [5,] 0.6467642 0.6467642 0.6467642
 [6,] 0.5597143 0.5597143 0.5597143
 [7,] 0.4854133 0.4854133 0.4854133
 [8,] 0.4278326 0.4278326 0.4278326
 [9,] 0.3834149 0.3834149 0.3834149
[10,] 0.3433031 0.3433031 0.3433031



On Tue, Dec 25, 2018 at 12:45 AM Eric Berger <ericjberger at gmail.com>
wrote:

Since you don't provide lambda, rh or qext it is impossible to reproduce
what you are seeing.
Also note that in this mailing list HTML formatted emails are not passed
along.



On Tue, Dec 25, 2018 at 4:13 AM M P <mzp3769 at gmail.com> wrote:

Hello,
I used commands below to obtain a surface, can plot it and all looks as
expected.
How do I evaluate values at new point. I tried as below but that produces
errors.
Thanks for suggestions/help.

x <- log(lambda)
y <- rh
z <- qext[,,2]

grid.l <- list(abcissa=x,ordinate=y)
xg <- make.surface.grid(grid.l)
out.p <- as.surface(xg,z)
plot.surface(out.p,type="p")

tried:
grid_new.l <- list(abcissa=c(-15.0,-10.),ordinate=y)
xg_new <- make.surface.grid(grid_new.l)

out_new.p <- predict.surface(out.p,xg_new)

results in this prompt:
predict.surface is now the function predictSurface>

        [[alternative HTML version deleted]]