Add test case for global models

Summary

This issue is part of a larger effort of supporting global models in TELEMAC (see #985).

In this context, we need a test case that would enable bench-marking of openTELEMAC for the following processes and modules:

  • global surges 2D / 3D
  • global tides + surges 2D / 3D (barotropic)
  • global 3D baroclinic
  • global waves (and ideally coupling with T2D or T3D)
  • global tsunamis (#1027)

Why is this feature useful?

Ability to test continuously the features important global models, like:

  • Nearest Neighbor Interpolation for 1D TS (already implemented !203 (merged))
  • Self Attraction and Loading #1110
  • global Multi fault tsunamis (#1027)

How to implement it?

Mesh

A very simple implementation is FESOM's pimesh:

pimesh

The PIMESH is 3140 nodes. We could also use a larger mesh.

FIY Here are the meshes size for the Malpasset test case:

  • the small version is 13541 nodes.
  • the large version is 53081 nodes.

Meteo

We are thinking about putting a subset of ERA5 for our own CI pipelines on Zenodo. They guarantee the same DOI for 25 years We have 50GB allowed per repo. Here are the sizes for one months of data:

  • 3 variables: u, v,p -> 4GB
  • 4 variables (In the future we'll need surface temp for better baroclinic simulations) -> 4.3GB ERA5 forcing (u10, v10, pmsl - and ideally tmsl for baroclinic implementation)

The vnv would need to:

  1. fetch the data on Zenodo
  2. convert it to a Selafin file (using xarray-selafin)
  3. interpolate using converter.py the meteo onto the mesh.

If (2) is too complicated to do, we can directly put a Selafin grid on Zenodo too.

Edited by Thomas SAILLOUR