help

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20090202/ed653da5/attachment.pl>
I think you have to loop over time ( "i"  in your case, which
presumably is a variable in your environment ).

something like:

for (i in 1:length(dat[1,1,])) {
   vz<-as.vector(dat[,,i])
   oz<-interpp(vx,vy,vz,ox,oy,linear=TRUE)
   Kolkata[i]<-oz$z
}
hi;
i just started using R and need to do the following.
i want to plot weather data (which is in a .netCDF file).
i was trying to plot some climate data for westbengal region and then i
tried to interpolate the data for kolkata (which will be a subset of the
oroginal data for westbengal).

i tried the following commands
graphics.off()
library(ncdf)
library(akima)
ncdata<-open.ncdf("A2.TREFHT-gis.westbengal.nc")
lon<-get.var.ncdf(ncdata,"lon")
lat<-get.var.ncdf(ncdata,"lat")
time<-get.var.ncdf(ncdata,"time")
dat <- get.var.ncdf(ncdata,"TREFHT")
Kolkata <- array(dim=length(time))
 vx<-as.vector(matrix(lon,nrow=length(lon),ncol=length(lat)))
vy<-as.vector(matrix(lat,nrow=length(lon),ncol=length(lat),byrow=TRUE))
vz<-as.vector(dat[,,i])
 ox<-as.vector(88)
oy<-as.vector(22)
oz<-interpp(vx,vy,vz,ox,oy,linear=TRUE)
Kolkata[i]<-oz$z
 Kolkata<-Kolkata-273
plot(Kolkata,type="l")

now the command are working fine but producing the plot for a single day
only. whereas i want to plot the data for the entire time period and want a
line diagram.
so my question is how can i do that? please help.

       [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo