Displaying 1 result from an estimated 1 matches for "hcol".
Did you mean:
col
2010 Sep 03
1
calculate monthly mean
...open.ncdf("/srv/ccrc/data05/z3236814/mm5/co2/2000/q.200001.mon.nc")
v1 <- nc$var [[1]]
v2 <- nc$var [[2]]
v3 <- nc$var [[3]]
data1 <- get.var.ncdf(nc,v1)
data2 <- get.var.ncdf(nc,v2)
data3 <- get.var.ncdf(nc,v3)
coldat = data1[1:111,1:101,23,1:60]
## creat colour
hcol = cumsum(coldat)
coldat = hcol
hcol = hcol/max(hcol,na.rm=TRUE)
col <- hsv(h=hcol,s=1,v=1)
X <-scatterplot3d(data3[1:111,1:101],data2[1:111,1:101],data1[1:111,1:101,23,1:60],color=col)