Open mdsumner opened 4 years ago
still doesn't work very well for mesh_plot
f <- "tos_Omon_IPSL-CM5A-MR_rcp85_r1i1p1_200601-210012.nc"
x <- tidync(f)
lon <- hyper_array(x)[[2]]
lat <- hyper_array(x)[[1]]
time_idx <- 1
qm <- rgl::qmesh3d(vertices = rbind(x = as.vector(lon),
y = as.vector(lat),
z = rep(c(hyper_array(x %>% activate("D2,D1,D0"))[["tos"]][,,time_idx]), each = 4L), h = 1),
indices = matrix(seq_along(lon), nrow = 4L))
regrid to normal is best I think