New vegetation friction law in TELEMAC-2D based on Cui et al. (2023)

See merge request !243

Improvements to TELEMAC-2D friction laws

See merge request !240

(bottom and free surface planes) with mesh transformation = 3 (fixed z)
Check introduced in MR !233

Definition of the vertical mesh by keywords in the TELEMAC-3D steering file

See merge request !233

Alexander BREUGEM authored 1 month ago and Chi-Tuân Pham committed 1 month ago

New dimension for keyword MESH TRANSFORMATION (array with 1 value per plane,
not a single value)
2 new keywords:
- SIGMA DISTRIBUTION OF PLANES
- PRESCRIBED ELEVATION OF PLANES
They allow to avoid to implement any Fortran file (USER_MESH_TRANSF)

Fix the interpolation of wave directions in time series files

Closes #1071

See merge request !236

[t2d][src] better time step managment in the finite volume kernel.

See merge request !225

New method Nearest Neighbor Search to interpolate time series files

Closes #1033

See merge request !203

Adding new Saint-Venant-Exner coupled finite volume solver for sediment transport in 1D: Arrakis

See merge request !234

[mascaret][courlis] Bug fix in mascaret.dico (following merge !238)

See merge request !239

This fixes a bug on Windows that caused the contour of some meshes to be
extracted incorrectly.

[VnV] Changes for the IMDC system

See merge request !237

Alexander BREUGEM authored 2 months ago and Chi-Tuân Pham committed 2 months ago

+ add missing med tag for the T3D lock-exchange sensitivity study

(1 argument forgotten for GAIA + arguments for T2D and T3D called with optional
syntax)

[fix] Time Series extraction for longitude/latitude meshes across IDL
(International Date Line)

Closes #1034

See merge request !202

Thomas SAILLOUR authored 2 months ago and Chi-Tuân Pham committed 2 months ago

(International Date Line)

MASCARET new keyword VARIABLES_PRECISION

Closes #1078

See merge request !238

- NERD for the advection of velocities (as for tracers) rather than
characteristics to decrease differences between sequential and parallel runs
- use of the non-hydrostatic version (but more CPU time consuming)
and conjugate gradient + preconditiong = 34 = 2*17 to solve PPE (more efficient
choice than GMRES here)
- use of preconditioning 34 = 2*17 for the diffusion of tracers + velocities
so that accuracy for the diffusion of tracers can be 1.E-15 without significant
increase of CPU time
- increase of MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES to 100

Update of epsilons (1 order of magnitude less than before for tracers)
Smaller GRAPHIC PRINTOUT PERIOD to draw time series for figure

(Intel Turbo Boost on EDF clusters)

- TREATMENT OF NEGATIVE DEPTHS = 2 (new default value since release 9.0)
to improve mass conservation (water depth + temperature, water depth up to
machine accuracy)
- solver and accuracy for propagation set to default (GMRES + 1.E-8)
- accuracy for diffusion of tracers can be set to 1.E-15 without significant
increase of CPU time
- preconditioning for diffusion = 34 = 2*17 for tracers + velocities for
more efficiency (less iterations)
- solver and preconditioning for PPE = conjugate gradient (1) and 34
- non-hydrostatic version (new default since 8.0)

- default options to solve propagation (Krylov space dimensions, accuracy)
without significant increase of CPU time
- no diffusion along x and y with keyword
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 0. rather than
SCHEME FOR DIFFUSION OF VELOCITIES = 0

This allows CPUs to operate at a higher clock rate, to speed up CPU-bound calculations

- SOLVER FOR PROPAGATION = 7 (GMRES = default) as quite more efficient than
old choice 2 (conjugate residual) or 1 (conjugate gradient)
This choice enables to set ACCURACY FOR PROPAGATION to its default value 1.E-8
rather than 1.E-6 with quite the same CPU times

- TIDAL FLATS = YES (default value) to improve mass conservation (water depth
+ tracer, water depth up to machine accuracy)
- SOLVER FOR PROPAGATION = 7 (GMRES = default) as the most efficient
- with theses choices of solvers, accuracies for propagation and PPE can be set
to default = 1.E-8
- ==> less differences between sequential and parallel runs

- TREATMENT OF NEGATIVE DEPTHS = 2 (new default value since release 9.0)
to improve mass conservation (water depth + temperature, water depth up to
machine accuracy)
- SCHEME FOR ADVECTION OF VELOCITIES = 5 (PSI) rather than 1 (characteristics)
for better mass conservation + reduced differences between seq/par runs
- use of gradient conjugate + preconditioning 34 = 2*17 to solve PPE, more
efficient than GMRES for this example
- accuracies for diffusion of velocities/tracers/k-epsilon + PPE set to
1.E-15 as easily and efficiently reached with this set-up