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@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
