Argo floats 3D

## =======================================================================================
## plotting argo float data
## implemented by Karline Soetaert
## =======================================================================================
# assumes data have been downloaded from e.g. "ftp://ftp.ifremer.fr/ifremer/argo/geo/atlantic_ocean/2012/01/"
## and put in a list and save in file argo.rda.

## ARGO
library(OceanView)

load(file = "argo.rda")

xlim <- c(-90, 50)
ylim <- c(-70, 70)
ii <- which (Hypsometry$x >= xlim[1] & Hypsometry$x <= xlim[2])
jj <- which (Hypsometry$y >= ylim[1] & Hypsometry$y <= ylim[2])

mx <- Hypsometry$x[ii]
my <- Hypsometry$y[jj]
mz <- Hypsometry$z[ii,jj]
mz [ mz > 0] <- 0

persp3D(x = mx, y = my, z = mz, col = "grey", plot = FALSE,
 main = "argo floats, januari 2012", cex.main = 2)

clim <- range(argo$temp, na.rm = TRUE)
with (argo, 
 for (i in 1:length(Time)) {
 pmin <- min(-press[[i]] , na.rm = TRUE) 
 scatter3D(longitude[[i]], latitude[[i]], -press[[i]] , colvar = temp[[i]] , 
 plot = FALSE, add = TRUE, pch = 16, cex = 0.6, clim = clim, 
 colkey = FALSE, NAcol = "transparent")
 }
)
# very slow:
plotdev(shade = 0.2)
colkey(clim = clim, clab = "dgC", add = TRUE, dist = -0.05, length = 0.5, width = 0.5)
#plotrgl(smooth = TRUE, lighting = TRUE)