should add background field and land mask in plot

[gaelforget]

https://github.com/ECCO-Hackweek/EH24-Drifters/blob/78bd420cf48e400c86be3c5c6598cb97071857f3/scripts/initial_simulations.jl#L41

For examples, see https://juliaclimate.github.io/MeshArrays.jl/dev/tutorials/geography.html

[kykinne]

import Pkg; Pkg.add("MeshArrays") using MeshArrays

fil=demo.download_polygons("countries.geojson") pol=MeshArrays.read_polygons(fil);

[lines!(a,l1,color = :white, linewidth = 0.5) for l1 in pol] #0-360 but countries are -180-180

Then also need to transform the particle lon positions -180 so that the grid are the same.

[YuanyuanSong99]

Perhaps the particle Lon positions should -360

[gaelforget]

import Pkg; Pkg.add("MeshArrays") using MeshArrays fil=demo.download_polygons("countries.geojson") pol=MeshArrays.read_polygons(fil); [lines!(a,l1,color = :white, linewidth = 0.5) for l1 in pol] #0-360 but countries are -180-180

Nice! Couple notes :

1. seems to need this before :

using MeshArrays, GeoJSON, CairoMakie, DataDeps
fi=Figure(); a=Axis(fi[1,1])

2. try lon180(x)=Float64(x>180.0 ? x-360.0 : x)

3. when you feel ready, please open a pull request to add this to the repository

[YuanyuanSong99]

I added examples_SouthAtlantic1.jl in the directory of scripts, with two ways to add land background. I recommend the second way. Here are figures by both ways.

[YuanyuanSong99]

Now we have a plot with depth colormap of particles

[YuanyuanSong99]

We also have a map using heatmap function