Dear all,
I'm new user of R package.
I'm trying to interpolate climate variables from ERA data.
Original data contains 6 column with 29,040 lines, and x, y dimension are
240x121.
I would like to get new data from column name: "lon", "lat" and "temp" by x, y
dimension are 3,586x1,800
So I used the code as below, but I got error message as "z[cbind(lx1 + 1, ly1)]
: out of range"
Help me please.
Thank you in advance.
==============================================================
library(fields)
d <- read.table("test.txt", dec=".", sep = "", nrows = 29040,
colClasses="numeric", col.names=c("lon", "lat", "geo_h", "temp", "rh", "hPa"))
vz=d$temp
z.temp <- matrix(vz, nrow=121, ncol=240, byrow=T)
obj <- list( x=seq( 0,358.5,,240), y=seq( -90,90,,121), z= z.temp[ 0:121,
0:240] )
grid.list<- list( x= seq( 0,358.5,,3586), y= seq( -90,90,,1800))
interp.surface.grid( obj, grid.list) -> look
==============================================================
Best Regards,
Parichat Wetchayont
----------------------------------------------------------------------------------
Radiation&Climate Physics Lab, Department of Geophysics,
Graduate School of Science, Tohoku University.
[email protected]
----------------------------------------------------------------------------------
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
