Hirlam · roelstappers · Oct 18, 2023 · Oct 16, 2023 · Oct 16, 2023 · Oct 16, 2023
diff --git a/Project.toml b/Project.toml
@@ -5,6 +5,7 @@ version = "0.1.1"
 
 [deps]
 CommonDataModel = "1fbeeb36-5f17-413c-809b-666fb144f157"
+Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
 
 [compat]

diff --git a/docs/make.jl b/docs/make.jl
@@ -11,11 +11,22 @@ makedocs(;
         prettyurls=get(ENV, "CI", "false") == "true",
         canonical="https://Hirlam.github.io/NCPlots.jl",
         edit_link="main",
+        collapselevel = 1,
         assets=String[],
     ),
     pages=[
         "Home" => "index.md",
-        "Examples" => "example.md",
+        "Examples" => [
+            "Basic" => "examples/era5_pv.md",
+            "Light" => "examples/light.md",
+            "Animation" => "examples/animations.md",
+            "Multiple plot in axis" => "examples/single_axis.md",
+            "Geopotential height" => "examples/geopotential.md",
+            "Limited Area plots" => "examples/lam.md",
+            "Metcoop Increments" => "examples/metcoop.md"
+            "Sliders and Menus" => "examples/sliders_menus.md",
+
+        ], 
         "Getting data" => "cds.md",
         "Keyboard control" => "keyboard.md",
         "API reference" => "references.md"

diff --git a/docs/src/animations.md b/docs/src/animations.md
diff --git a/docs/src/assets/metcoop_menu.png b/docs/src/assets/metcoop_menu.png
diff --git a/docs/src/example.md b/docs/src/example.md
diff --git a/docs/src/examples/animations.md b/docs/src/examples/animations.md
@@ -0,0 +1,31 @@
+# Animations
+
+
+
+## Animated plots
+
+To make these animations use Observables and `@lift` to lift the dataset view e.g.
+
+```julia 
+ds = Dataset(...) 
+t = Observable(1)
+pv = @lift(view(ds,time=$t)["pv"]) 
+plot(pv) 
+```
+
+Updating `t` wil update the plot e.g.
+
+```julia
+for i=1:100
+   t[] = i
+   sleep(0.01)
+end 
+```
+
+![](../assets/pv_era5_2.gif)
+
+You can access eyeposition in `ax.scene` and update in the for loop  
+
+![](../assets/test.mp4)
+
+## Recording animations
diff --git a/docs/src/examples/era5_pv.md b/docs/src/examples/era5_pv.md
@@ -0,0 +1,19 @@
+
+# [Example](@id example) 
+
+
+
+These examples use NetCDF files from the `docs/src/assets` directory. For the animations download larger datasets as exlained in the Getting data section. 
+
+
+## ERA5 potential vorticity 
+
+```julia   
+using NCPlots, GLMakie, NCDatasets
+ds = Dataset("assets/era5_pv_z_500hPa.nc")
+pv = view(ds, time=1)["pv"]
+fig, ax, plt = plot(pv, colormap=:RdBu, colorrange=(-3e-6,3e-6))
+save("assets/era5_pv_docs.png", fig) # hide
+```
+
+![](../assets/era5_pv_docs.png) 
diff --git a/docs/src/examples/geopotential.md b/docs/src/examples/geopotential.md
@@ -0,0 +1,18 @@
+
+# Plot on geopotential height surface 
+
+Multiply `x`, `y`, `z` by geopotential height 
+
+![](../assets/logo2.png)
+
+
+## Plot on orography
+
+See plot geopotential height. Color show temperature increments from data assimilation  
+
+![](../assets/envar_2019081803.png) 
+
+
+
+
+
diff --git a/docs/src/examples/lam.md b/docs/src/examples/lam.md
@@ -0,0 +1,21 @@
+# Limited area plots
+
+## CARRA 
+
+
+![](../assets/east_domain.png)
+![](../assets/east_west_domain.png)
+
+## Metcoop
+
+```julia
+archive="/lustre/storeB/immutable/archive/projects/metproduction/MEPS/"
+ds = Dataset("$archive/2023/10/01/meps_det_2_5km_20231001T00Z.nc")
+field = view(ds,hybrid=65,time=1)["air_temperature_ml"]
+plot(field)
+```
+
+
+
+![](../assets/metcoop.png) 
+
diff --git a/docs/src/examples/light.md b/docs/src/examples/light.md
@@ -0,0 +1,7 @@
+# Control light
+
+Access properties `specular` `diffuse` etc. in the returned `plt`.
+
+![](../assets/marble2.png)
+
+Explain how to update camera, change lightposition  
diff --git a/docs/src/examples/metcoop.md b/docs/src/examples/metcoop.md
@@ -0,0 +1,42 @@
+# Metcoop 
+
+Example reading from MEPS archive at MET Norway  
+
+```julia
+using Glob, NCDatasets, NCPlots, GLMakie
+
+archive="/lustre/storeB/immutable/archive/projects/metproduction/MEPS/2023/01/30"
+
+files = glob("meps_det_2_5km*.nc", archive)   # use e.g.  "....*00.nc" to only match files with 00:00 UTC
+
+ds = NCDatasets.Dataset(files, aggdim="cycle", isnewdim=true, constvars=["x","y","hybrid","time"],deferopen=false) 
+
+fig = Figure() 
+fcint = 3; llmax=66; nlev=65
+
+
+tsl, lsl,csl = SliderGrid(fig[2,1], 
+  (label="LL",range=1:(llmax-fcint),startvalue=1), 
+  (label="Level",range=nlev:-1:1,startvalue=nlev),
+   (label="Cycle",range=2:length(files),startvalue=2),
+   tellheight=false,
+   width=350
+)
+ax= Axis(fig[1,2])
+
+varname = "air_temperature_ml"
+var = ds[varname]
+
+dsv = @lift(var[:,:,$(lsl.value),$(tsl.value),$(csl.value)] -
+            var[:,:,$(lsl.value),$(tsl.value)+3,$(csl.value)-1]
+            )
+
+title = @lift(string(Analysis increment, " " ,varname,  basename(files[$(cycle_sl.value)]))
+subtitle = @lift(string(" Level ", $(lsl.value), "forecast", $(tsl.value)-1, " Hour") )
+
+heatmap!(ax,dsv,colormap=:RdBu,colorrange=(-2,2))
+Label(fig[0,:],label,tellwidth=false)
+Label(fig[1,:],label,tellwidth=false)
+```
+
+![](../assets/metcoop_menu.png)
diff --git a/docs/src/examples/single_axis.md b/docs/src/examples/single_axis.md
@@ -0,0 +1,6 @@
+
+# Mutliple plots in single axis 
+
+See `docs/plotlogo.jl` for an example 
+
+![](../assets/logo.png)
diff --git a/docs/src/examples/sliders_menus.md b/docs/src/examples/sliders_menus.md
@@ -0,0 +1,33 @@
+# Sliders Menus
+
+
+
+```julia
+
+crangedict = Dict(
+    "pv" => (-2e-6,2e-6),
+    "z" => (20000.0,50000.0)
+)
+
+fig = Figure()
+
+menu = Menu(fig[1,1],options=setdiff(keys(ds),dimnames(ds)))
+menu2 = Menu(fig[1,2], options=colorschemes)
+sl = SliderGrid(fig[2,1],
+    (label="time", range=1:length(ds["time"]),startvalue=1,format = x-> Dates.format(ds["time"][x],"yyyymmdd-HH"))
+    )
+ax = LScene(fig[3,1], show_axis=false)
+colorrange=@lift(crangedict($(menu.selection)))
+colormap=@lift(colorschem)
+
+ds = Dataset("...")
+dsv = @lift(view(ds,time=$(sl.sliders[1].value))[$(menu.selection)])
+plotvar!(ax,dsv,colorrange=colorrange,colormap=:RdBu)
+
+```
+
+
+## Level slider
+
+
+
diff --git a/docs/src/index.md b/docs/src/index.md
@@ -3,10 +3,9 @@
 NCPlots is a Julia package for plotting (meteorological) data on the sphere. Datasets should conform to the CF convention and implement the [CommonDataModel.jl](https://github.com/JuliaGeo/CommonDataModel.jl) interface like [NCDatasets](https://github.com/Alexander-Barth/NCDatasets.jl) and [GRIBDatasets](https://github.com/JuliaGeo/GRIBDatasets.jl)  
 
 
-See [here](@ref Examples) for some examples 
-See [Getting data](@ref getting_data) how to obtain ERA5 and CARRA data. 
+See [here](@ref example) for a basic example and [Getting data](@ref getting_data) how to obtain ERA5 and CARRA data. 
 
-## Installation
+### Installation
 
 NCPlots is in the Harmonie registry. 
 Add the Harmonie registry (hit `]` in the Julia REPL to enter package mode).  
@@ -22,7 +21,7 @@ Then install NCPlots with
 pkg> add NCPlots 
 ```
 
-## Data requirements 
+### Data requirements 
 
 Variables should only have dimensions `longitude` and `latitude`. If there are additional dimensions create a view
 e.g. for a dataset `ds` that contains a variable `pv` where `pv` in addition to `longitude` and `latitude` also has dimension  `time` do