New features of Telemac-3D version 7.1
Estimated release date: end 2015
Latest article update: 05/01/2016
Main new features of this version 7.1
Boundary conditions on the bed
It is now possible to impose a flow rate on several liquid boundaries placed on the bed.
This is done throught the keywords:
"OPEN BOUNDARY CONDITIONS ON THE BED" and "PRESCRIBED FLOWRATES ON THE BED", "MAXIMUM NUMBER OF BOUNDARIES ON THE BED"
The liquid boundaries can be defined using subroutine "LIMI3D". An additional test case has been added, see steering file "./examples/telemac3d/bottom_bc/t3d_bottom_inlet.cas" and user fortran "./examples/telemac3d/bottom_bc/t3d_bottom_inlet.f".
Option for Dirac sources
It is now possible to define Dirac sources in Telemac-3D, as in Telemac-2D.
This is done throught the keyword TYPE OF SOURCES (TYPE DES SOURCES).
The option is set by default to the original formulation but the Dirac option is recommended with high numbers of sources.
Schematic tides for TELEMAC-2D and TELEMAC-3D
Modelling schematic tides (astronomical spring tide, mean spring tide, mean tide, mean neap tide, astronomical neap tide) to study schematic or average scenarii is available for the 4 tidal databases (JMJ already available since v6.1, and the other ones: TPXO-type from Oregon State University, NEA/FES, PREVIMER).
Keyword to activate OPTION FOR TIDAL BOUNDARY CONDITIONS with a value between 2 to 6.
Keywords to fill in: "GLOBAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER" with TPXO-type tidal database from Oregon State University (mandatory) and "LOCAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER" with JMJ, NEA/FES, PREVIMER tidal databases (if not filled in, a value is then automatically calculated).
New or modified key-words
"TURBULENCE REGIME FOR LATERAL SOLID BOUNDARIES" (REGIME DE TURBULENCE POUR LES PAROIS LATERALES), "TURBULENCE REGIME FOR THE BOTTOM" (REGIME DE TURBULENCE POUR LE FOND)
"MODEL" is replaced by "REGIME" in the English version of these 2 keywords.
"CONSTANT SEDIMENT SETTLING VELOCITY" (VITESSE DE CHUTE CONSTANTE)
with a final E in the French version (only), not CONSTANT.
"VELOCITIES OF THE SOURCES ALONG Z" (VITESSE DES SOURCES SELON Z)
Exactly like VELOCITIES OF THE SOURCES ALONG X, and Y. the Fortran variable WSCE already existed but was set to 0.
"COEFFICIENT FOR HORIZONTAL DIFFUSION OF TRACERS" (COEFFICIENT DE DIFFUSION HORIZONTAL DES TRACEURS), "COEFFICIENT FOR VERTICAL DIFFUSION OF TRACERS" (COEFFICIENT DE DIFFUSION VERTICAL DES TRACEURS)
One value per tracer can now be given.
"TYPE OF SOURCES" (TYPE DE SOURCES)
As in Telemac-2D, the user can set it to 1 (non Dirac sources) or 2 (Dirac sources). It is set by default to 1.
"MAXIMUM NUMBER OF BOUNDARIES" (NOMBRE MAXIMUM DE FRONTIERES)
Maximal number of boundaries in the mesh. Used for dimensioning arrays. Can be increased if needed.
"MAXIMUM NUMBER OF BOUNDARIES ON THE BED" (NOMBRE MAXIMUM DE FRONTIERES SUR LE FOND)
Maximal number of liquid boundaries on the bed. Used for dimensioning arrays. Can be increased if needed.
"MAXIMUM NUMBER OF SOURCES" (NOMBRE MAXIMUM DE SOURCES)
Maximal number of punctual sources in the mesh. Used for dimensioning arrays. Can be increased if needed.
"MAXIMUM NUMBER OF TRACERS" (NOMBRE MAXIMUM DE TRACEURS)
These 4 new keywords are used to give an a priori dimension to a number of allocatable files, they give the values of respectively MAXFRO, MAXBLB, MAXSCE and MAXTRA. A number of arrays with hardcoded dimension are now allocated with these dimensions. Changing these keyword will avoid changing the hardcoded values, which required recompiling the whole package. A number of user functions have been amended to cope with this change: Q3, TR3, SL3, VIT3, T3D_DEBSCE and T3D_TRSCE. To adapt a specific user function, either copy your personal bits of code into the new functions or do the following: remove the declarations and initialisation of DEJA and OK. Change the remaining occurences of OK into OKQ3 for function Q3, OKSL3 for function SL3, OKVIT3 for function VIT3, OKTR3 for function TR3, etc., OKTRSCE for function T3D_TRSCE and OKDEBSCE for function T3D_DEBSCE.
"OPEN BOUNDARY CONDITIONS ON THE BED" (CONDITIONS OUVERTES SUR LE FOND)
Defines if there are open boundary conditions on the bed. Default value = NO.
"PRESCRIBED FLOWRATES ON THE BED" (DEBITS IMPOSES SUR LE FOND)
Sets the value for flow rate on flow rate-imposed bed boundaries. Default value = 0.
"OPTION FOR TIDAL BOUNDARY CONDITIONS" (OPTION POUR LES CONDITIONS AUX LIMITES DE MAREE)
This keyword is now an array with a value given per liquid boundary file, separated by semicolons. Otherwise all liquid boundaries with prescribed velocities or depths would receive tidal conditions. This avoids a clash when using weirs in the domain. 0 codes for conditions other than tidal.
BEWARE !!! Old models must be changed if their tidal boundary is not number 1. In that case this keyword must be changed and more values given.
"LOCAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER" (NUMERO LOCAL DU POINT POUR CALER LA PLEINE MER)
Local number between 1 and the number of tidal boundary points (of the HARMONIC CONSTANTS FILE) where the tidal boundary conditions are computed with JMJ, NEA, FES, PREVIMER databases (except TPXO-type databases).
The tidal constituents have their phase shifted with respect to this point to start the simulation with a high water (for schematic tides only). Default value is 0, which means that the value is automatically calculated by TELEMAC.
"GLOBAL NUMBER OF THE POINT TO CALIBRATE HIGH WATER" (NUMERO GLOBAL DU POINT POUR CALER LA PLEINE MER)
Global number of the point with respect to which the tidal constituents have their phase shifted to start the calculation with a high water (for schematic tides only). Only harmonic constants databases like TPXO are concerned. The user has to give an existing value.
"OPTION FOR THE BOUNDARY CONDITIONS OF K-EPSILON" (OPTION POUR LES CONDITIONS AUX LIMITES DU K-EPSILON)
To compute the lateral boundary conditions of k and epsilon (KEPCL3 subroutine):
1 = no turbulence (= KMIN and EMIN defined in CSTKEP subroutine) which is the default value (same behaviour as before v7.0),
2 = Hans and Burchard's formula (introduced since v7.0).
Default value is now 1.E-10 rather than 1.E-20 previously.
Modifications in user subroutines or new subroutines
This subroutine has been changed to call noerod.f from Sisyphe library to get the non erodable bed. This is done to have the same value of ZR in both programmes in case of coupling. In further versions it will even be the same array in memory. It is important to understand that implementing noerod.f for specific cases is now necessary even when there is no coupling with Sisyphe. Note that this was already announced for version 7.0 but eventually not done.
As the number of arguments of CONDI_TPXO subroutine (in module TPXO) has changed, this subroutine has changed in version 7.1. The user is advised to update this subroutine with version 7.1.