### Focus Manuals: Programing Guide: Building matrices and vectors: Available construction of Vectors

#### FORMULs of SUBROUTINE VECTOR

##### MASMAS

FORMUL = 'MASMAS '

Integrals of the bases, or product of a mass matrix by a vector with the value of 1 everywhere. ##### MASVEC

FORMUL = 'MASVEC '

Product of a mass matrix by a vector F. ##### SUPG

FORMUL = 'SUPG ' • is a vector with the components G and H. (this would have to be modified in dimension 3).
• is a vector with the components U, V and W. • is a vector with the components U and V.

VGRADP is the same formula, with corrections when the generalised sigma transformation is used.

##### FLUBOR

FORMUL = 'FLUBOR ' • is a vector with the components U, V and W.
• is the normal outer vector.
##### FLUBOR2

FORMUL = 'FLUBOR2 '

In 3D only, FLUBOR2 is like FLUBOR, but in the case of a generalised sigma transformation. • is a vector with the components U, V and W. • is a vector with the components U, V and W.
##### FLUBDF

FORMUL = 'FLUBDF ' • is a vector with the components U, V and W.
• is the normal outer vector. Beware the minus sign !!!!

If FORMUL(16:16) is equal to 'Y' or 'Z' instead of 'X' , the derivatives will be obtained according to y or z. If FORMUL(16:16) is equal to 'Y' or 'Z' instead of 'X' , the derivatives will be obtained according to y or z.

In 3 dimensions, variants are available:

• GRADF(X,Y) X and GRADF(X,Y) Y will consider only the gradient of a function

which does not depend on 'Z'.

• GRADF2 will take care of hydrostatic inconsistencies.
##### PRODF

FORMUL = 'PRODF' This vector is used in the calculation of the turbulent production with the model k-epsilon.

##### DIVQ

FORMUL = 'DIVQ ' • is a vector with the components U, V and W.
##### SUPGDIVU

FORMUL = 'SUPGDIVU ' • is a vector with the components U, V and W.
• is a vector with the components F, G and H.
##### FLUDIF

FORMUL = 'FLUDIF ' This is not currently used nor implemented. This is specifically for 3D computations with prisms, and unlike VGRADF, the test function is here a 2-dimensional test function (no dependency on z). This is used by Telemac-3D in subroutine WSTARW. This is used in 2D, mostly for computing fluxes. H in HUGRADP stands for the depth denoted h, which can be misleading as it does not refer to the function H which is an argument of subroutine VC19AA. A variant HUGRADP2 exists, in this case the velocity is not only of components (U,V), but . This is a way of treating the gradient of the free surface elevation as a piecewise constant function, which it is in reality when the depth is linear.

#### SUBROUTINE BY NAMES

The existing subroutines building vectors in version 6.+ are the following, their function can be deduced from the explanations above:

vc00aa.f vc00bb.f vc00cc.f vc00pp.f vc00pp2.f
vc00ft.f vc00ff.f vc00tt.f vc01aa.f vc01bb.f
vc01ff.f vc01ft.f vc01oo.f vc01pp.f vc01tt.f
vc01tt0.f vc03aa.f vc03bb.f vc04aa.f vc04pp.f
vc05oo.f vc05aa.f vc04tt.f vc05ff.f vc05ft.f
vc08aa.f vc08bb.f vc08cc.f vc08pp.f vc08tt.f
vc09aa.f vc10oo.f vc11aa.f vc11aa2.f vc11bb.f
vc11pp.f vc11tt.f vc11tt0.f vc13aa.f vc13bb.f
vc13cc.f vc13pp.f vc13pp2.f vc13tt.f vc14aa.f
vc15aa.f vc16aa.f vc18pp.f vc19aa.f