Date: 24th May 2001

The principal features of this version v5p1 of TELEMAC-2D are:

• New treatment of Shallow Water Equations as a wave equation
• Developments on porosity
• Development on the drag force of vertical structures.

Integer value. The default value is 1

1 : Classical option. Shallow Water equations are treated in their primitive form. Depth and velocity are coupled in the final linear system.

2 : Wave equation. The momentum equation is used to eliminate the velocity in the continuity equation. Automatically sets the following options:

> MASS-LUMPING ON THE VELOCITY (MASS-LUMPING SUR LA VITESSE) : 1.

> IMPLICITATION FOR DIFFUSION OF VELOCITY (IMPLICITATION POUR LA DIFFUSION DES VITESSES) : 0.

> IMPLICITATION FOR VELOCITY (IMPLICITATION POUR LA VITESSE) : 1.

> IMPLICITATION FOR DEPTH (IMPLICITATION FOR DEPTH) : 1.

> C-U PRECONDITIONING (PRECONDITIONNEMENT C-U) : NO

Moreover, the TYPE OF ADVECTION (FORME DE LA CONVECTION) for velocity must be : 1 or 6 (characteristics or PSI). 1 is recommended. These options do not change the results of a steady-state flow but may slightly smooth the results for a flow varying in time. With option 2, SOLVER : 1 (conjugate gradient) is recommended as the matrix of the final linear system is symmetric.

Logical value. The default value is Yes

If YES, output of initial conditions in the RESULTS FILE

Logical value. The default value is No

Drag forces from vertical structures are taken into account. SUBROUTINE DRAGFO must then be implemented.

Integer value. The default value is 1

Storage 2 has been suppressed (storage 3 is always better)

Storage 3 has been modified to work in parallel mode.

We now have :

3

Storage 1 : classical Element By Element storage Storage 3 : Edge-based storage

An example is given in the sources_util directory. The required parameters are:

The drag coefficient CD, the diameter of vertical structures DIAM, the number of vertical structures N, the definition of the area where the drag force will be applied. Only the sum of all CD*DIAM is used in the program so more complex situations with different drag coefficients and diameters may as well be treated.

Called only if TREATMENT OF TIDAL FLATS (TRAITEMENT DES BANCS DECOUVRANTS) : 3

An example is given in the sources_util directory. Can be used to set a porosity in an area of the computational domain. This porosity is given element by

element. Please note that the value of porosity is previously set by the treatment.