declarations_telemac3d.f

Go to the documentation of this file.

!                   *****************************
                    MODULE declarations_telemac3d
!                   *****************************
!
!
!***********************************************************************
! TELEMAC3D   V8P2
!***********************************************************************
!
!brief    DECLARATION OF PRINCIPAL TELEMAC3D VARIABLES
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  G. ANTOINE (EDF LAB, LNHE)
!+        19/09/2014
!+        V7P0
!+   Adding variables for mixed sediment
!
!history  A. JOLY (EDF LAB, LNHE)
!+        27/08/2015
!+        V7P1
!+   Imposed flowrates on the bed.
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        24/03/2016
!+        V7P2
!+   New keywords for distributive schemes. Adding Z at time t(n), as ZN
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief_def
!
!       NOTE: THIS MODULE IS ORGANISED IN 10 PARTS
!
!       (1) VECTORS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!       (2) MATRICES (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!       (3) BLOCKS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!       (4) INTEGERS
!       (5) LOGICAL VALUES
!       (6) REALS
!       (7) STRINGS
!       (8) SLVCFG STRUCTURES
!       (9) MESH STRUCTURE
!      (10) ALIASES
!
!-----------------------------------------------------------------------
! (1) VECTORS (REAL AND INTEGER)
!-----------------------------------------------------------------------
!
!
!     3D VELOCITY COMPONENTS AT PREVIOUS TIMESTEP (TIME N)
!
      TYPE(bief_obj), TARGET :: un,vn,wn
!
!     3D VELOCITY COMPONENTS AFTER ADVECTION
!
      TYPE(bief_obj), TARGET :: uc,vc,wc
!
!     3D VELOCITY COMPONENTS AFTER DIFFUSION
!
      TYPE(bief_obj), TARGET :: ud,vd,wd
!
!     3D VELOCITY COMPONENTS AT CURRENT TIMESTEP (TIME N+1)
!
      TYPE(bief_obj), TARGET :: u,v,w
!
!     EXPLICIT SOURCE TERMS ON VELOCITIES
!
      TYPE(bief_obj), TARGET :: s0u,s0v,s0w
!
!     IMPLICIT SOURCE TERMS ON VELOCITIES U
!
      TYPE(bief_obj), TARGET :: s1u,s1v,s1w
!
!     RIGHT-HAND SIDE ON VELOCITIES EQUATIONS
!
      TYPE(bief_obj), TARGET :: smu,smv
!
!     PRESCRIBED VELOCITY COMPONENTS ON THE BOTTOM
!
      TYPE(bief_obj), TARGET :: uborf,vborf,wborf
!
!     PRESCRIBED VELOCITY COMPONENTS ON THE LATERAL BOUNDARY
!
      TYPE(bief_obj), TARGET :: uborl,vborl,wborl
!
!     SAVED PRESCRIBED VELOCITY COMPONENTS ON THE LATERAL BOUNDARY FROM BC FILE
!
      TYPE(bief_obj), TARGET :: uborcli,vborcli
!
!     PRESCRIBED VELOCITY COMPONENTS AT THE FREE SURFACE
!
      TYPE(bief_obj), TARGET :: ubors,vbors,wbors
!
!     ORIGINAL PRESCRIBED VALUES AT THE FREE SURFACE
!
      TYPE(bief_obj), TARGET :: trborsave,uborsave,vborsave
      TYPE(bief_obj), TARGET :: wborsave,kborsave,eborsave
      TYPE(bief_obj), TARGET :: nuborsave
!
!     (FRICTION VELOCITY)**2 FOR BOTTOM AND LATERAL BOUNDARIES
!
      TYPE(bief_obj), TARGET :: uetcar,uetcal
!
!     LOGARITHMIC LAW FOR U-VELOCITY COMPONENT, ON THE BOTTOM: NU*DU/DN = AUBORF*U + BUBORF
!
      TYPE(bief_obj), TARGET :: auborf,buborf
!
!     LOGARITHMIC LAW FOR U-VELOCITY COMPONENT, ON THE LATERAL BOUNDARIES: NU*DU/DN = AUBORL*U + BUBORL
!
      TYPE(bief_obj), TARGET :: auborl,buborl
!
!     LOGARITHMIC LAW FOR U-VELOCITY COMPONENT, AT THE FREE SURFACE: NU*DU/DN = AUBORS*U + BUBORS
!
      TYPE(bief_obj), TARGET :: aubors,bubors
!
!     LOGARITHMIC LAW FOR V-VELOCITY COMPONENT, ON THE BOTTOM: NU*DV/DN = AVBORF*V + BVBORF
!
      TYPE(bief_obj), TARGET :: avborf,bvborf
!
!     LOGARITHMIC LAW FOR V-VELOCITY COMPONENT, ON THE LATERAL BOUNDARIES: NU*DV/DN = AVBORL*V + BVBORL
!
      TYPE(bief_obj), TARGET :: avborl,bvborl
!
!     LOGARITHMIC LAW FOR V-VELOCITY COMPONENT, AT THE FREE SURFACE: NU*DV/DN = AVBORS*V + BVBORS
!
      TYPE(bief_obj), TARGET :: avbors,bvbors
!
!     LOGARITHMIC LAW FOR W-VELOCITY COMPONENT, ON THE BOTTOM: NU*DW/DN = AWBORF*W + BWBORF
!
      TYPE(bief_obj), TARGET :: awborf,bwborf
!
!     LOGARITHMIC LAW FOR W-VELOCITY COMPONENT, ON THE LATERAL BOUNDARIES: NU*DW/DN = AWBORL*W + BWBORL
!
      TYPE(bief_obj), TARGET :: awborl,bwborl
!
!     LOGARITHMIC LAW FOR W-VELOCITY COMPONENT, AT THE FREE SURFACE: NU*DW/DN = AWBORS*W + BWBORS
!
      TYPE(bief_obj), TARGET :: awbors,bwbors
!
!     TYPES OF BOUNDARY CONDITIONS FOR VELOCITY COMPONENTS ON THE BOTTOM
!
      TYPE(bief_obj), TARGET :: liubof,livbof,liwbof
!
!     TYPES OF BOUNDARY CONDITIONS FOR U ON THE LATERAL BOUNDARIES
!
      TYPE(bief_obj), TARGET :: liubol,livbol,liwbol
      TYPE(bief_obj), TARGET :: liubol_user,livbol_user
!
!     TYPES OF BOUNDARY CONDITIONS FOR VELOCITY COMPONENTS AT THE FREE SURFACE
!
      TYPE(bief_obj), TARGET :: liubos,livbos,liwbos
!
!     'COLOUR' OF BOUNDARY NODES (TAKEN IN BOUNDARY CONDITIONS FILE)
!
      TYPE(bief_obj), TARGET :: boundary_colour
!
!     PLANE NUMBER OF LAST PLANE WHICH HAS A CRUSHED LAYER ABOVE IT
!
      TYPE(bief_obj), TARGET :: ipbot
!
!     SIGMA-TRANSFORMED VERTICAL VELOCITY COMPONENT
!
      TYPE(bief_obj), TARGET :: ws
!
!     DYNAMIC PRESSURE AT TIME N+1 AND TIME N (NON-HYDROSTATIC)
!
      TYPE(bief_obj), TARGET :: dp,dpn
!
!     HYDROSTATIC PRESSURE (NON-HYDROSTATIC)
!
      TYPE(bief_obj), TARGET :: ph
!
!     PRESCRIBED DYNAMIC PRESSURE ON BOTTOM, LATERAL BOUNDARIES AND FREE SURFACE
!
      TYPE(bief_obj), TARGET :: pborf, pborl, pbors
!
!     TYPE OF BOUNDARY CONDITIONS FOR DYNAMIC PRESSURE
!
      TYPE(bief_obj), TARGET :: lipbof, lipbol, lipbos
!
!     K OF K-EPSILON MODEL AT PREVIOUS TIMESTEP (TIME N)
!
      TYPE(bief_obj), TARGET :: akn
!
!     K OF K-EPSILON MODEL AFTER ADVECTION
!
      TYPE(bief_obj), TARGET :: akc
!
!     K OF K-EPSILON MODEL AT TIME N+1
!
      TYPE(bief_obj), TARGET :: ak
!
!     EXPLICIT SOURCE TERMA FOR K IN K-EPSILON MODEL
!
      TYPE(bief_obj), TARGET :: s0ak
!
!     IMPLICIT SOURCE TERMA FOR K IN K-EPSILON MODEL
!
      TYPE(bief_obj), TARGET :: s1ak
!
!     PRESCRIBED K OF K-EPSILON MODEL ON: LATERAL BOUNDARY
!                                         FREE SURFACE
!                                         BOTTOM
!
      TYPE(bief_obj), TARGET :: kborl,kbors,kborf
!
!     COEFFICIENTS OF LOGARITHMIC LAW FOR K OF K-EPSILON MODEL, ON THE BOTTOM
!                                                               FREE SURFACE
!
      TYPE(bief_obj), TARGET :: akborf,bkborf,akborl
      TYPE(bief_obj), TARGET :: bkborl,akbors,bkbors
!
!     TYPES OF BOUNDARY CONDITIONS FOR K OF K-EPSILON MODEL ON THE BOTTOM
!                                                        LATERAL BOUNDARY
!                                                            FREE SURFACE
      TYPE(bief_obj), TARGET :: likbof,likbol,likbos
!
!     EPSILON OF K-EPSILON MODEL AT PREVIOUS TIMESTEP (TIME N)
!
      TYPE(bief_obj), TARGET :: epn
!
!     EPSILON OF K-EPSILON MODEL AFTER ADVECTION
!
      TYPE(bief_obj), TARGET :: epc
!
!     EPSILON OF K-EPSILON MODEL AT TIME N+1
!
      TYPE(bief_obj), TARGET :: ep
!
!     EXPLICIT/IMPLICIT SOURCE TERM OF EPSILON IN K-EPSILON MODEL
!
      TYPE(bief_obj), TARGET :: s0ep,s1ep
!
!     PRESCRIBED EPSILON OF K-EPSILON MODEL ON THE BOTTOM
!                                        LATERAL BOUNDARY
!                                            FREE SURFACE
!
      TYPE(bief_obj), TARGET :: eborf,eborl,ebors
!
!     COEFFICIENTS IN LOGARITHMIC LAW FOR EPSILON OF K-EPSILON MODEL, ON THE BOTTOM
!                                                                  LATERAL BOUNDARY
!                                                                      FREE SURFACE
      TYPE(bief_obj), TARGET :: aeborf,aeborl,aebors
      TYPE(bief_obj), TARGET :: beborf,beborl,bebors
!
!     TYPES OF BOUNDARY CONDITIONS FOR EPSILON OF K-EPSILON MODEL ON THE BOTTOM
!                                                              LATERAL BOUNDARY
!                                                                  FREE SURFACE
      TYPE(bief_obj), TARGET :: liebof,liebol,liebos
!
!     NUTILD OF S-A MODEL AT PREVIOUS TIMESTEP (TIME N)
!
      TYPE(bief_obj), TARGET :: nun
!
!     NUTILD OF S-A MODEL AFTER ADVECTION
!
      TYPE(bief_obj), TARGET :: nuc
!
!     NUTILD OS S-A MODEL AT TIME N+1
!
      TYPE(bief_obj), TARGET :: nu
!
!     EXPLICIT/IMPLICIT SOURCE TERM OF NUTILD IN S-A MODEL
!
      TYPE(bief_obj), TARGET :: s0nu,s1nu
!
!     PRESCRIBED NUTILD OF S-A MODEL ON THE BOTTOM
!                                 LATERAL BOUNDARY
!                                 FREE SURFACE
      TYPE(bief_obj), TARGET :: nuborf,nuborl,nubors
!
!     COEFFICIENTS IN LOGARITHMIC LAW FOR NUTIL OF S-A MODEL, ON THE BOTTOM
!                                                          LATERAL BOUNDARY
!                                                              FREE SURFACE
      TYPE(bief_obj), TARGET :: anuborf,anuborl,anubors
      TYPE(bief_obj), TARGET :: bnuborf,bnuborl,bnubors
!
!     TYPES OF BOUNDARY CONDITIONS FOR NUTILD OF S-A MODEL ON THE BOTTOM
!                                                        LATERAL BOUNDARY
!                                                            FREE SURFACE
      TYPE(bief_obj), TARGET :: linubof,linubol,linubos
!
!
!     VERTICALLY INTEGRATED VELOCITY COMPONENTS
!
      TYPE(bief_obj), TARGET :: u2d,v2d
!
!     VERTICALLY INTEGRATED VELOCITY COMPONENTS ON BOUNFDARIES
!
      TYPE(bief_obj), TARGET :: ubor2d,vbor2d
!
!     VERTICALLY INTEGRATED VELOCITY COMPONENTS AT PREVIOUS TIME STEP
!
      TYPE(bief_obj), TARGET :: un2d,vn2d
!
!     SOURCE TERMS ON VELOCITIES
!
      TYPE(bief_obj), TARGET :: fu,fv
!
!     WAVE STRESSES FROM ARTEMIS OR TOMAWAC
!
      TYPE(bief_obj), TARGET :: fxh,fyh
!
!     FLUXES ON BOUNDARIES
!
      TYPE(bief_obj), TARGET :: flbor,flblim
!
!     ADVECTION FIELD
!
      TYPE(bief_obj), TARGET :: uconv,vconv,wconv
!
!     VERTICAL VELOCITY IN TRANSFORMED MESH
!
      TYPE(bief_obj), TARGET :: wsconv
!
!     ADVECTION FIELD FOR CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: uconvc,vconvc
!
!     PSEUDO-VISCOSITY IN WAVE EQUATION
!
      TYPE(bief_obj), TARGET :: nuwave
!
!     PIECE-WISE LINEAR FUNCTION, AND ITS GRADIENT
!     USED WITH "COMPATIBILITY OF FREE SURFACE GRADIENT"
!
      TYPE(bief_obj), TARGET :: zconv,grazco
!
!     D**-1 IN WAVE EQUATION (SEE BOOK)
!
      TYPE(bief_obj), TARGET :: dm1
!
!     EDGE BY EDGE FLUXES (AND ASSEMBLED IN PARALLEL VERSION)
!
      TYPE(bief_obj), TARGET :: flodel,flopar
!
!     LIMITATION OF FLUXES
!
      TYPE(bief_obj), TARGET :: flulim
!
!     VOLUME AROUND POINTS AT TIME N+1, NOT ASSEMBLED (VOLU)
!                                       ASSEMBLED IN PARALLEL (VOLUPAR)
!
      TYPE(bief_obj), TARGET :: volu,volupar
!
!     VOLUME AROUND POINTS AT TIME N
!
      TYPE(bief_obj), TARGET :: volun,volunpar
!
!     VARIOUS VOLUMES (AT AN INTERMEDIATE TIME STEP)
!
      TYPE(bief_obj), TARGET :: volut,volu3d,volu3dpar
!
!     INTEGRAL OF 2D BASES (SCALAR AND ASSEMBLED IN PARALLEL)
!
      TYPE(bief_obj), TARGET :: volu2d,v2dpar
!
!     INVERSE OF INTEGRAL OF BASES IN 2D, IN 3D
!
      TYPE(bief_obj), TARGET :: unsv2d,unsv3d
!
!     INTERNAL FLUXES, PER NODE
!
      TYPE(bief_obj), TARGET :: fluint
!
!     FLUX AT BOUNDARIES, AND ASSEMBLED FORM IN PARALLEL
!
      TYPE(bief_obj), TARGET :: fluext,fluextpar
!
!     FLUXES INTEGRATED ALONG THE VERTICAL
!
      TYPE(bief_obj), TARGET :: flint2
!
!     ABSOLUTE VALUE OF THE VORTICITY
!
      TYPE(bief_obj), TARGET :: rotat
!
!     ABSOLUTE DENSITY
!
      TYPE(bief_obj), TARGET :: rho
!
!     POTENTIAL DENSITY
!
      TYPE(bief_obj), TARGET :: rhopot
!
!     RELATIVE DENSITY = (RHO-RHO0)/RHO0
!
      TYPE(bief_obj), TARGET :: deltar
!
!     RICHARDSON NUMBER
!
      TYPE(bief_obj), TARGET :: ri
!
!     FRICTION COEFFICIENT ON THE BOTTOM
!
      TYPE(bief_obj), TARGET :: rugof
!
!     FRICTION COEFFICIENT ON THE LATERAL BOUNDARY
!
      TYPE(bief_obj), TARGET :: rugol
!
!     ADIMENSIONAL FRICTION COEFFICIENT
!
      TYPE(bief_obj), TARGET :: cf
!
!     WIND VELOCITY
!
!     TYPE(BIEF_OBJ), TARGET :: WIND
!
!     MAGNITUDE OF WIND VELOCITY
!
      TYPE(bief_obj), TARGET :: windspd
!
!     DENSITY, TEMPERATURE, SALINITY FOR COUPLING WITH AED2
!
      TYPE(bief_obj), TARGET :: rho4aed2,temp4aed2,sali4aed2
!
!     THICKNESS OF CELLS FOR COUPLING WITH AED2
!
      TYPE(bief_obj), TARGET :: thick4aed2
!
!     NUMBER OF THE BENTHIC LAYER FOR COUPLING WITH AED2
!
      TYPE(bief_obj), TARGET :: benth4aed2
!
!     ATMOSPHERIC PRESSURE
!
!     TYPE(BIEF_OBJ), TARGET :: PATMOS
!
!     PARAMETERS FOR GLOBAL MASS AND FLUX BALANCES
!
      TYPE(bief_obj), TARGET :: masini,masse,massen
!
!     FLUX (ACROSS THE BOUNDARIES)
!
      TYPE(bief_obj), TARGET :: flux
!
!     FLUX OF TRACERS (ACROSS THE BOUNDARIES, SOURCES, BED)
!
      TYPE(bief_obj), TARGET :: fluxb
!
!     CUMULATED FLUXES OF TRACERS SINCE BEGINNING OF COMPUTATION
!
      TYPE(bief_obj), TARGET :: flucum
!
!     DEPTH AT TIME N+1, AT TIME TN
!
      TYPE(bief_obj), TARGET :: h,hn
!
!     PROPAGATION DEPTH, I.E. DEPTH INTO THE TERM DIV(HU)
!
      TYPE(bief_obj), TARGET :: hprop
!
!     HIGH WATER MARK, HIGH WATER TIME
!
      TYPE(bief_obj), TARGET :: maxz,tmaxz
!
!     VARIATION OF DEPTH : H(N+1)-H(N)
!
      TYPE(bief_obj), TARGET :: dh
!
!     RIGHT-HAND SIDE FOR WATER DEPTH
!
      TYPE(bief_obj), TARGET :: smh
!
!     PRESCRIBED DEPTH ON LATERAL BOUNDARIES
!
      TYPE(bief_obj), TARGET :: hbor
!
!     SAVED PRESCRIBED DEPTH ON LATERAL BOUNDARIES FROM BC FILE
!
      TYPE(bief_obj), TARGET :: hborcli
!
!     RAIN IN M/S MULTIPLIED BY VOLU2D, AND VERSION ASSEMBLED IN PARALLEL
!
      TYPE(bief_obj), TARGET :: pluie,parapluie
!
!     FOR CULVERTS/TUBES/BRIDGES MANAGEMENT
      INTEGER, POINTER, DIMENSION(:) :: entbus,sorbus
      DOUBLE PRECISION, POINTER, DIMENSION(:) :: ubus,vbus
      TYPE(bief_obj),TARGET :: lrgbus, haubus, altbus
      TYPE(bief_obj),TARGET :: csbus , cebus , delbus, angbus, lbus
      TYPE(bief_obj),TARGET :: dbus, tbus, clpbus
      TYPE(bief_obj),TARGET :: secbus, longbus,aabus
      TYPE(bief_obj),TARGET :: cv, c56, cv5, c5, ctrash, fricbus, circ
      DOUBLE PRECISION relaxb
!
!     BLOCK OF VARIABLES TO BE ADVECTED BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: fn3d,fc3d
!
!     BLOCK OF VARIOUS BIEF_OBJ STRUCTURES FOR VARIABLES TO BE ADVECTED
!
      TYPE(bief_obj), TARGET :: bl_fn,bl_fc,bl_s0f,bl_fsc
      TYPE(bief_obj), TARGET :: bl_bol,bl_borl
!
!     VERTICAL COORDINATES FOR PROPAGATION STEP
!
      TYPE(bief_obj), TARGET :: zprop
!
!     BOTTOM ELEVATION
!
      TYPE(bief_obj), TARGET :: zf
!
!     BOTTOM COORDINATES GIVEN PER ELEMENT
!
      TYPE(bief_obj), TARGET :: zfe
!
!     PERCENTUAL MESH PLANES
!
      TYPE(bief_obj), TARGET :: zstar
!
!     Z ELEVATIONS USED ONLY WITH SUPG ADVECTION
!
      TYPE(bief_obj), TARGET :: zt
!
!     DH/DT
!
      TYPE(bief_obj), TARGET :: dssudt
!
!     PRESCRIBED ELEVATIONS OF PLANES
!
      TYPE(bief_obj), TARGET :: zplane
!
!     TYPE OF SIGMA TRANSFORMATION
!
      TYPE(bief_obj), TARGET :: transf_plane
!
!     COORDINATES IN THE TRANSFORMED MESH FOR THE METHOD OF CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: zchar
!
!     MASK OF ELEMENTS
!
      TYPE(bief_obj), TARGET :: maskel
!
!     MASK OF NODES
!
      TYPE(bief_obj), TARGET :: maskpt
!
!     3D MASK OF BOUNDARY ELEMENTS (CASE WITH MASKING)
!
      TYPE(bief_obj), TARGET :: maskbr
!
!     3D MASK OF BOUNDARY ELEMENTS (WORKING ARRAY FOR DEBIMP_3D)
!
      TYPE(bief_obj), TARGET :: masktr
!
!     POSITIONS OF FLOATING BODIES (X,Y, Z AND Z IN TRANSFORMED MESH)
!
      TYPE(bief_obj), TARGET :: xflot,yflot,zflot
!
!     BARYCENTRIC COORDINATES OF FLOATING BODIES, AT THE FOOT OF CHARACTERISTICS
!     HORIZONTAL AND VERTICAL
!
      TYPE(bief_obj), TARGET :: shpflo,shzflo
!
!     TAGS FOR FLOATING BODIES
!
      TYPE(bief_obj), TARGET :: tagflo
!
!     CLASS FOR FLOATING BODIES
!
      TYPE(bief_obj), TARGET :: clsflo
!
!     STARTING ELEMENTS OF FLOATING BODIES FOR CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: eltflo
!
!     STARTING LEVEL OF FLOATING BODIES FOR CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: etaflo
!
!     TYPE OF BOUNDARY CONDITIONS ON DEPTH (AND VALUE OF USER SAVED IF THOMPSON)
!
      TYPE(bief_obj), TARGET :: lihbor,lihbor_user
!
!     LIQUID BOUNDARY NUMNER OF BOUNDARY POINTS, OF BOUNDARY ELEMENTS
!
      TYPE(bief_obj), TARGET :: numliq,numliq_elm
!
!     LIQUID BED BOUNDARY NUMBER OF BOUNDARY POINTS
!
      TYPE(bief_obj), TARGET :: nliqbed
!
!     PROPAGATION BC TYPES (TELEMAC2D'S PROPAG)
!
      TYPE(bief_obj), TARGET :: limpro
!
!     RIGHT HAND SIDE FOR THE LINEAR EQUATIONS IN 3D
!
      TYPE(bief_obj), TARGET :: sem3d
!
!     RIGHT HAND SIDE FOR THE LINEAR EQUATIONS IN 2D
!
      TYPE(bief_obj), TARGET :: sem2d
!
!     ELEMENT-ORIENTED WORKING ARRAYS
!
      TYPE(bief_obj), TARGET :: te1,te2,te3
!
!     PIECE-WISE LINEAR FREE SURFACE
!
      TYPE(bief_obj), TARGET :: zflats
!
!     VOID VECTOR STRUCTURE
!
      TYPE(bief_obj), TARGET :: svide
!
!     RIGHT HAND SIDE OF CONTINUITY EQUATION WHEN SOURCES
!
      TYPE(bief_obj), TARGET :: sources
!
!     RIGHT HAND SIDE OF CONTINUITY EQUATION WHEN DEALING WITH OPEN BOUNDARIES ON THE BED
!
      TYPE(bief_obj), TARGET :: bedflu
!
!     SETTLING VELOCITY
!
      TYPE(bief_obj), TARGET :: wchu
!
!     SETTLING VELOCITY FOR SEDI3D
!
      TYPE(bief_obj), TARGET :: s3d_wchu
!
!     THICKNESS OF SOLID FRACTION OF THE BED LAYER ( S3D_EPAI=DZ/(1+S3D_IVIDE), DZ BED LAYER THICKNESS )
!
!      TYPE(BIEF_OBJ), TARGET :: S3D_EPAI
!
!     VOID RATIO (GIBSON MODEL ONLY)
!
      TYPE(bief_obj), TARGET :: s3d_ivide
!
!     TIME COUNTER FOR CONSOLIDATION MODEL (MULTILAYER MODEL)
!
      TYPE(bief_obj), TARGET :: s3d_temp
!
!     ELEVATIONS OF TIME T(N) SAVED
!
      TYPE(bief_obj), TARGET :: zn
!
!     CONCENTRATION OF MUD BED LAYER (MULTILAYER MODEL)
!
!     TYPE(BIEF_OBJ), TARGET :: S3D_CONC
!
      DOUBLE PRECISION,DIMENSION(:,:), TARGET, ALLOCATABLE :: s3d_conc
!
! TBE FLOCULATION AND HINDERED SETTLING MODEL
!     HINDERED SETTLING SWITCH
      LOGICAL s3d_hinder
!     HINDERED SETTLING METHOD
      INTEGER s3d_hind_type
!
!     SEDIMENT CLASSES TO CONSIDER FOR SUSPENSION
      INTEGER isusp
!
      DOUBLE PRECISION s3d_cgel, s3d_cini
!
!     FLOCULATION SWITCH
      LOGICAL s3d_floc
!     FLOCULATION METHOD
      INTEGER s3d_floc_type
!
! BLOCK OF DATA FOR SETTLING AND DIFFUSION 1D PROFILE MODEL
      TYPE(bief_obj), TARGET :: trav1
! CV...
! uniformisation avec SISYPHE
!
!     MAXIMUM NUMBER OF LAYERS ON THE MESH
!
      INTEGER, PARAMETER :: s3d_nlaymax= 20
!
!     BED SHEAR STRENGTH OF THE MUD BED LAYER (MULTILAYER MODEL)
!
!      TYPE(BIEF_OBJ), TARGET :: S3D_TOCE

      DOUBLE PRECISION,DIMENSION(:,:), TARGET, ALLOCATABLE :: s3d_toce
!
      TYPE(bief_obj), TARGET :: s3d_laytoce
!
      DOUBLE PRECISION,DIMENSION(:,:), TARGET, ALLOCATABLE :: s3d_epai
!
      TYPE(bief_obj), TARGET :: s3d_laythi
!
!     MUD CONCENTRATION FOR EACH LAYER (Constante for initialisation)
!
      DOUBLE PRECISION :: s3d_conc_layer(s3d_nlaymax)
!
!     CRITICAL SHEAR STRESS FOR EROSION (Initialisation)
!
      DOUBLE PRECISION :: s3d_toce_layer(s3d_nlaymax)
!
!     SEDIMENT BED LAYER THICKNESS (Initialisation)
!
      DOUBLE PRECISION :: s3d_es_layer(s3d_nlaymax)
!
!  Mass bilan
! MASSINI IN SISYPHE
!
      DOUBLE PRECISION :: s3d_massusp0, s3d_massusp
!
!      DOUBLE PRECISION :: MASTEN, MASTOU
!
!      DOUBLE PRECISION :: MASTCP, MASFIN
!
      DOUBLE PRECISION :: s3d_masdep
!
      DOUBLE PRECISION :: s3d_masbed0,s3d_masbed

! ...CV
!
!     THICKNESS OF FRESH DEPOSIT (FLUID MUD LAYER)
!
      TYPE(bief_obj), TARGET :: s3d_hdep
!
!     EROSION FLUX
!
      TYPE(bief_obj), TARGET :: s3d_fluer
!
!     PROBABILITY OF DEPOSIT
!
      TYPE(bief_obj), TARGET :: s3d_pdepo
!
! CV...
!       Harmonisation avec Sisyphe
!
!     DEPOSITION FLUX
!
      TYPE(bief_obj), TARGET :: s3d_fludp
!
!     DEPOSITION FLUX (Implict part : can replace S3D_PDEPO)
!
      TYPE(bief_obj), TARGET :: s3d_fludpt
!
!     BED EVOLUTION
!
      TYPE(bief_obj), TARGET :: s3d_zf_s
!
!     CUMULATED BED EVOLUTION
!
      TYPE(bief_obj), TARGET :: s3d_esomt
!
! ...CV
!
!     RIDIG BED ELEVATION
!
      TYPE(bief_obj), TARGET :: zr
!
!     EQUILIBRIUM CONCENTRATION
!
      TYPE(bief_obj), TARGET :: s3d_cref
!
!     AVERAGE DIAMETER OF GRAINS
!
      TYPE(bief_obj), TARGET :: s3d_dmoy
!
!     NUMBER OF POINTS WITHIN THE BED ALONG THE VERTICAL
!
      TYPE(bief_obj), TARGET :: npf
!
!     VECTOR OF SIZE NDP*NELEM (NUMBER OF POINTS IN THE ELEMENT * NELEM)
!
      TYPE(bief_obj), TARGET :: wel
!
!     3D VECTOR TO STORE DH AND HN IN 3D MESH (RESTART PROCEDURE)
!
      TYPE(bief_obj), TARGET :: dhhn
!
!     FOR TIDAL BOUNDARY CONDITIONS
!
      TYPE(bief_obj),TARGET :: hbtide,ubtide,vbtide,numtide
!
!     VARIABLES TRANSMITTED FROM TOMAWAC TO SISYPHE
!
      TYPE(bief_obj),TARGET :: dirmoy,hm0,tpr5,orbvel
!
!     VARIABLES USED IN CASE OF MIXED SEDIMENT
!
      TYPE(bief_obj),TARGET :: s3d_epaico,s3d_epainco,s3d_pvsco
      TYPE(bief_obj),TARGET :: s3d_pvsnco,s3d_fluerc
      TYPE(bief_obj),TARGET :: s3d_fludptc,s3d_fludptnc,s3d_fludpc
      TYPE(bief_obj),TARGET :: s3d_fludpnc,s3d_fluernc
!
!      DISTANCE FROM WALLS
!
      TYPE(bief_obj), TARGET :: wdist
!
!     STRAIN NORM
!
      TYPE(bief_obj), TARGET :: strain
!
!     ROTATION NORM
!
      TYPE(bief_obj), TARGET :: rotan
!
!-----------------------------------------------------------------------
! (2) MATRICES
!-----------------------------------------------------------------------
!
!
!     SUPG NON SYMMETRIC MATRIX
!
      TYPE(bief_obj), TARGET :: msupg
!
!     MURD NON SYMMETRIC MATRIX
!
      TYPE(bief_obj), TARGET :: mmurd
!
!     DIFFUSION MATRIX
!
      TYPE(bief_obj), TARGET :: mdiff
!
!     MURD MATRIX FOR TIDAL FLATS VERSION
!
      TYPE(bief_obj), TARGET :: murd_tf
!
!     WORKING MATRICES
!
      TYPE(bief_obj), TARGET :: mtra1,mtra2
!
!     2D BOUNDARY MATRIX
!
      TYPE(bief_obj), TARGET :: mbor2d
!
!     WORKING 2D MATRIX
!
      TYPE(bief_obj), TARGET :: matr2h
!
!-----------------------------------------------------------------------
! (3) BLOCKS
!-----------------------------------------------------------------------
!
!
!     TRACERS AT PREVIOUS TIME STEP (TAN IN VERSIONS PRIOR TO 6.3
!                                    WHICH WAS A CONFLICT WITH TANgent)
!
      TYPE(bief_obj), TARGET :: trn
!
!     TRACERS AFTER ADVECTION
!
      TYPE(bief_obj), TARGET :: tac
!
!     TRACERS AT CURRENT TIME STEP
!
      TYPE(bief_obj), TARGET :: ta
!
!     EXPLICIT SOURCE TERM FOR TRACERS
!
      TYPE(bief_obj), TARGET :: s0ta
!
!     IMPLICIT SOURCE TERM FOR TRACERS
!
      TYPE(bief_obj), TARGET :: s1ta
!
!     PRESCRIBED TRACERS ON THE BOTTOM
!
      TYPE(bief_obj), TARGET :: taborf
!
!     PRESCRIBED TRACERS ON THE LATERAL BOUNDARY
!
      TYPE(bief_obj), TARGET :: taborl
!
!     PRESCRIBED TRACERS AT THE FREE SURFACE
!
      TYPE(bief_obj), TARGET :: tabors
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: atabof
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: btabof
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: atabol
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: btabol
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: atabos
!
!     LOGARITHMIC LAW FOR TRACERS: NU*DTA/DN = ATABO*TA + BTABO
!
      TYPE(bief_obj), TARGET :: btabos
!
!     TYPES OF BOUNDARY CONDITIONS FOR TRACERS ON THE BOTTOM
!
      TYPE(bief_obj), TARGET :: litabf
!
!     TYPES OF BOUNDARY CONDITIONS FOR TRACERS ON THE LATERAL BOUNDARY
!
      TYPE(bief_obj), TARGET :: litabl
!
!     TYPES OF BOUNDARY CONDITIONS FOR TRACERS AT THE FREE SURFACE
!
      TYPE(bief_obj), TARGET :: litabs
!
!     VISCOSITY
!
      TYPE(bief_obj), TARGET :: viscvi
!
!     DIFFUSIVITY FOR TRACERS
!
      TYPE(bief_obj), TARGET :: viscta
!
!     VISCOSITY IN SPALART_ALLMARAS SCHEME
!
      TYPE(bief_obj), TARGET :: viscnu
!
!     BOTTOM GRADIENTS
!
      TYPE(bief_obj), TARGET :: gradzf
!
!     FREE SURFACE GRADIENT (BLOCK OF 2 COMPONENTS)
!
      TYPE(bief_obj), TARGET :: gradzs
!
!     OLD FREE SURFACE GRADIENT (BLOCK OF 2 COMPONENTS)
!
      TYPE(bief_obj), TARGET :: gradzn
!
!     2D MASK, FOR TELEMAC2D COMPATIBILITY
!
      TYPE(bief_obj), TARGET :: mask,mask_3d
!
!     BLOCKS OF ARRAYS FOR THE USER IN 3D AND 2D
!
      TYPE(bief_obj), TARGET :: prive,prive2d
!
!     BLOCK OF DIFFERENTIATED VECTORS
!
      TYPE(bief_obj), TARGET :: advar,advar2d
!
!     BLOCKS OF ARRAYS FOR THE USER
!
      TYPE(bief_obj), TARGET :: mat2d, tm1
!
!     BLOCKS THAT WILL CONTAIN WORK STRUCTURES
!
      TYPE(bief_obj), TARGET :: trav3,trav2,itrav3
!
      TYPE(bief_obj), TARGET :: varsor, varcl
!
      TYPE(bief_obj), TARGET :: varso3
!
!     VALUES OF ADVECTED VARIABLES AT SOURCES
!
      TYPE(bief_obj), TARGET :: u_sce,v_sce,w_sce,ak_sce,ep_sce,ta_sce
      TYPE(bief_obj), TARGET :: nu_sce
!
!     BLOCK OF WORk ARRAYS FOR LIPS ADVECTION SCHEME
!
      TYPE(bief_obj), TARGET :: tb2
!
!
!     CALCULATED VARIABLES FOR THE COUPLING TOMAWAC/TELEMAC3D
!
!
!     WAVE INDUCED PRESSURE AND ITS DERIVATIVES
!
      TYPE(bief_obj), TARGET ::  wip, wipdx, wipdy
!
!     STOKES VELOCITY
!
      TYPE(bief_obj), TARGET :: ustokes, vstokes, wstokes
!
!     2D PART OF WSTOKES
!
      TYPE(bief_obj), TARGET :: wst1, wst2
!
!     STOKES VELOCITY INTEGRATED OVER DEPTH AND ITS DERIVATIVES
!
      TYPE(bief_obj), TARGET :: us2d, vs2d
      TYPE(bief_obj), TARGET :: dus2d, dvs2d
      TYPE(bief_obj), TARGET :: dus2dx, dvs2dy
!
!     DERIVATIVES OF STOKES VELOCITY
!
      TYPE(bief_obj), TARGET :: ustx, vsty
!
!     SURFACES STRESS DUE TO DEPTH INDUCED WAVE BREAKING
!
      TYPE(bief_obj), TARGET :: fdx, fdy
!
!     DISSIPATION DUE TO BOTTOM FRICTION
!
      TYPE(bief_obj), TARGET :: fbx, fby
!
!     MODIFIED FRICTION COEFFICIENT DUE TO WAVES + CURRENTS
!
      TYPE(bief_obj), TARGET :: cfwc
!
!     MODIFIED TURBULENT WAVES VISCOSITY DUE TO BREAKING
!
      TYPE(bief_obj), TARGET :: fdk
!
!     FORCES DUE TO EFFECT OF WIND ON WAVES
!
      TYPE(bief_obj), TARGET :: fwx, fwy
!
!     BOTTOM MOMENTUM SWITCH
!
      LOGICAL bot_moment
!
!-----------------------------------------------------------------------
! (4) INTEGERS
!-----------------------------------------------------------------------
!
!     MAXIMUM OF SOURCE POINTS
!
      INTEGER, TARGET :: maxsce
!
!     MAXIMUM OF LIQUID BOUNDARIES
!
      INTEGER maxfro
!
!     MAXIMUM OF LIQUID BOUNDARIES ON THE BED
!
      INTEGER maxblb
!
!     MAXIMUM OF TRACERS
!
      INTEGER, TARGET :: maxtra
!
!     MAXIMUM OF OUTPUT VARIABLES
!
      INTEGER, PARAMETER :: maxvar = 140
      INTEGER :: nvar_hyd
!
!     MAXIMUM OF OUTPUT VARIABLES FOR 3D
!
      INTEGER, PARAMETER :: maxva3 = 140
      INTEGER :: nvar_t3d
!
!     ADDRESS OF TRACERS IN VARSO3D
!     IT IS COMPUTED IN NOMVAR_TELEMAC3D
!
      INTEGER :: adr_trac
!
!     ADDRESS OF TRACERS IN VARSO2D
!     IT IS COMPUTED IN NOMVAR_2D_IN_3D
!
      INTEGER :: adr_trac_2d
!
!     ADDRESS OF LAST 2D WORK ARRAY BEFORE TRACERS
!
      INTEGER, PARAMETER :: adr_trav2 = 17
!
!     NUMBER OF TIME STEPS
!
      INTEGER, TARGET :: nit
!
!     TIME STEP
!
      INTEGER, TARGET :: lt
!
!     NUMBER OF HORIZONTAL PLANES
!
      INTEGER, TARGET :: nplan
!
!     NUMBER OF TRACERS
!
      INTEGER, TARGET :: ntrac
!
!     NUMBER OF PRESCRIBED TRACERS GIVEN
!
      INTEGER ntracer
!
!     PRINTOUT PERIOD FOR FLOATING BODIES
!
      INTEGER floprd
!
!     NUMBER OF DROGUES (CONSIDERED AS THE MAXIMUM)
!
      INTEGER nflot_max,nflot
!
!     PRINTOUT PERIOD FOR RESULT FILE
!
      INTEGER, TARGET :: graprd
!
!     LISTING PRINTOUT PERIOD
!
      INTEGER, TARGET :: lisprd
!
!     NUMBER OF FIRST TIME STEP FOR GRAPHIC PRINTOUTS
!
      INTEGER gradeb
!
!     NUMBER OF FIRST TIME STEP FOR LISTING PRINTOUTS
!
      INTEGER lisdeb
!
!     NUMBER OF BOTTOM SMOOTHINGS
!
      INTEGER lisfon
!
!     NUMBER OF SUB ITERATIONS FOR NON LINEARITIES
!
      INTEGER nsousi
!
!     HORIZONTAL TURBULENCE MODEL
!
      INTEGER iturbh
!
!     VERTICAL TURBULENCE MODEL
!
      INTEGER iturbv
!
!     TURBULENCE MODEL FOR THE BOTTOM
!
      INTEGER lisruf
!
!     TURBULENCE MODEL FOR LATERAL SOLID BOUNDARIES
!
      INTEGER lisrul
!
!     INITIAL GUESS FOR DEPTH
!
      INTEGER iordrh
!
!     SPATIAL PROJECTION TYPE AND SAVE OF ITS ORIGINAM VALUE
!
      INTEGER protyp,protyp0
!
!     TOTAL NUMBER OF SOURCES
!     (including source points and culverts)
!
      INTEGER nsce
!
!     NUMBER OF SOURCE POINTS
!
      INTEGER nptsce
!
!     NUMBER OF SINK/SOURCES GIVEN BY THEIR NODE NUMBERS IN THE 2D MESH
!
      INTEGER nrejeti
!
!     ADDRESSES OF SOURCE POINTS IN THE 2D MESH
!
      INTEGER, ALLOCATABLE, TARGET :: isce(:)
!
!     PLANE NUMBERS OF SOURCE POINTS IN THE 2D MESH
!
      INTEGER, ALLOCATABLE, TARGET :: ksce(:)
!
!     ENTRY CULVERT AND OUTPUT CULVERT ARRAYS
!
      INTEGER, POINTER, DIMENSION(:) :: kentbus,ksorbus
!
!     NUMBER OF VELOCITIES OF SOURCES GIVEN
!
      INTEGER nrejeu
!
!     SCHEME FOR ADVECTION OF VELOCITIES
!
      INTEGER schcvi
!
!     SCHEME FOR ADVECTION OF TRACERS
!
      INTEGER, ALLOCATABLE :: schcta(:)
!
!     SCHEME FOR ADVECTION OF K-EPSILON
!
      INTEGER schcke
!
!     SCHEME FOR ADVECTION OF DEPTH
!
      INTEGER schch
!
!     SCHEME FOR DIFFUSION OF VELOCITIES
!
      INTEGER schdvi
!
!     SCHEME FOR DIFFUSION OF TRACERS
!
      INTEGER schdta
!
!     SCHEME FOR DIFFUSION OF K-EPSILON
!
      INTEGER schdke
!
!     TREATMENT ON TIDAL FLATS FOR VELOCITIES
!
      INTEGER trbavi
!
!     TREATMENT ON TIDAL FLATS FOR TRACERS
!
      INTEGER trbata
!
!     TREATMENT ON TIDAL FLATS FOR K-EPSILON
!
      INTEGER trbake
!
!     TREATMENT ON TIDAL FLATS FOR S-A
!
      INTEGER trbasa
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED DISCHARGE
!
      INTEGER ndebit
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED ELEVATION
!
      INTEGER ncote
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED VELOCITY
!
      INTEGER nvit
!
!     NUMBER OF BED BOUNDARIES WITH PRESCRIBED DISCHARGE
!
      INTEGER nbedflo
!
!     ORIGINAL DATE OF TIME
!
      INTEGER mardat(3)
!
!     ORIGINAL HOUR OF TIME
!
      INTEGER martim(3)
!
!     VECTOR LENGTH
!
      INTEGER lvmac
!
!     NUMBER OF PRIVATE ARRAYS, NUMBER OF PRIVATE ARRAYS WITH GIVEN NAME
!
      INTEGER npriv,n_names_priv2d,npriv2d
!
!     NUMBER OF DIFFERENTIATING ARRAYS, AND THOSE WITH A GIVEN NAME
!
      INTEGER nadvar,n_names_advar2d,nadvar2d
!
!     RANK OF TEMPERATURE IN TRACERS
!
      INTEGER ind_t
!
!     RANK OF SALINITY IN TRACERS
!
      INTEGER ind_s
!
!     RANK OF SEDIMENTS IN TRACERS
!
      INTEGER ind_sed
!
!     NUMBER OF POINTS IN THE ELEMENT
!
      INTEGER ndp
!
!     LAW OF BOTTOM FRICTION
!
      INTEGER kfrot
!
!     LAW OF FRICTION ON LATERAL BOUNDARIES
!
      INTEGER kfrotl
!
!     MATRIX STORAGE
!
      INTEGER optass
!
!     MATRIX-VECTOR PRODUCT
!
      INTEGER produc
!
!     OPTION FOR THE TREATMENT OF TIDAL FLATS
!
      INTEGER optban
!
!     TREATMENT OF NEGATIVE DEPTHS
!
      INTEGER opt_hneg
!
!     OPTION FOR THE DIFFUSION OF VELOCITIES
!
      INTEGER opdvit
!
!     OPTION FOR THE TREATMENT OF SOURCES
!
      INTEGER optsou
!
!     OPTION FOR THE TREATMENT OF CULVERTS
!
      INTEGER optbuse
!
!     SUPG OPTION
!
      INTEGER optsup(4)
!
!     MATRIX STORAGE FOR 2D MATRICES
!
      INTEGER optass2d
!
!     3D DISCRETISATION TYPE
!
      INTEGER ielm3
!
!     2D DISCRETISATION TYPES, HORIZONTAL, VERTICAL
!
      INTEGER ielm2h,ielm2v
!
!     DISCRETISATION TYPES
!
      INTEGER ielm0, ielmh, ielmu, ielm1, ielmx
!
!     NUMBER OF LAYERS OF 3D ELEMENTS (NPLAN - 1)
!
      INTEGER netage
!
!     NUMBER OF VARIABLES TREATED IN MASS BALANCE
!
      INTEGER nvbil
!
!     MAXIMUM NUMBER OF HORIZONTAL PLANES WITHIN THE BED (GIBSON MODEL)
! OBSOLETE, replaced by S3D_NCOUCHou S3D_NLAYMAX
      INTEGER s3d_npfmax

! NEW
!     NUMBER OF LAYERS WITHIN THE BED (MULTILAYER MODEL)
!
      INTEGER s3d_ncouch
!
!     MIXING LENGTH MODEL
!
      INTEGER mixing
!
!     DAMPING FUNCTION
!
      INTEGER damping
!
!     VELOCITY PROFILES
!
      INTEGER, ALLOCATABLE :: profvel(:)
!
!     TREATMENT OF FLUXES AT THE BOUNDARIES
!
      INTEGER, ALLOCATABLE :: dirflu(:)
!
!     VELOCITY VERTICAL PROFILES
!
      INTEGER, ALLOCATABLE :: verprovel(:)
!
!     TRACERS VERTICAL PROFILES
!
      INTEGER, ALLOCATABLE :: verprotra(:)
!
!     NUMBER OF LIQUID BOUNDARIES
!
      INTEGER nfrliq
!
!     CHOICE OF MESH TRANSFORMATION
!
      INTEGER transf
!
!     CHOICE OF DENSITY LAW
!
      INTEGER denlaw
!
!     DELWAQ PRINTOUT PERIOD
!
      INTEGER waqprd
!
!     ORIGIN COORDINATES
!
      INTEGER i_orig,j_orig
!
!     NUMBER OF POINTS GIVEN FOR EACH DISCHARGE-ELEVATIONS CURVES
!
      INTEGER, ALLOCATABLE :: pts_curves(:)
!
!     STAGE-DISCHARGE CURVES
!
      INTEGER, ALLOCATABLE :: sta_dis_curves(:)
!
!     KEYWORD DEBUGGER
!
      INTEGER, TARGET :: debug
!
!     SKIN FRICTION
!
      INTEGER s3d_icr
!
!     Equilibrium transport formula
!
      INTEGER s3d_icq
!
!     RECORD NUMBER IN THE WAVE DRIVEN CURRENTS FILE
!
      INTEGER npth
!
!     NUMBER OF CULVERTS/TUBES/BRIDGES
!
      INTEGER nbuse
!
!     GEOMETRY FILE NUMBER
!
      INTEGER, TARGET :: t3dgeo
!
!     BOUNDARY CONDITIONS FILE NUMBER
!
      INTEGER, TARGET :: t3dcli
!
!     PREVIOUS COMPUTATION FILE NUMBER
!
      INTEGER t3dpre
!
!     RESULT FILE NUMBER
!
      INTEGER, TARGET :: t3dres
!
!     BOTTOM TOPOGRAPHY FILE NUMBER
!
      INTEGER t3dfon
!
!     FILE NUMBER FOR SCOPE
!
!     INTEGER T3DSCO
!
!     2D RESULT FILE NUMBER
!
      INTEGER, TARGET :: t3dhyd
!
!     FORMATTED DATA FILE 1
!
      INTEGER t3dfo1
!
!     FORMATTED DATA FILE 2
!
      INTEGER t3dfo2
!
!     BINARY DATA FILE 1
!
      INTEGER t3dbi1
!
!     BINARY DATA FILE 2
!
      INTEGER t3dbi2
!
!     ASCII ATMOSPHERIC DATA FILE
!
      INTEGER t3atma
!
!     BINARY ATMOSPHERIC DATA FILE
!
      INTEGER t3atmb
!
!     BINARY BOUNDARY DATA FILE
!
      INTEGER t3dbnd
!
!     SEDIMENTOLOGICAL RESULT FILE
!
      INTEGER s3d_t3dsed
!
!     PREVIOUS COMPUTATION SEDIMENTOLOGICAL FILE
!
      INTEGER t3dsus
!
!     REFERENCE FILE NUMBER
!
      INTEGER t3dref
!
!     CULVERTS/BRIDGES DATA FILE
!
      INTEGER t3dbus
!
!     LIQUID BOUNDARIES FILE NUMBER
!
      INTEGER t3dimp
!
!     FILES FOR DELWAQ
!
      INTEGER t3ddl1,t3ddl2,t3ddl3,t3ddl4,t3ddl5,t3ddl6,t3ddl7,t3ddl8
      INTEGER t3ddl9,t3dl10,t3dl11
!
!     STAGE-DISCHARGE CURVES FILE NUMBER
!
      INTEGER t3dpar
!
!     SOURCES FILE NUMBER
!
      INTEGER t3dvef
!
!     BINARY RESULTS FILE NUMBER
!
      INTEGER t3drbi
!
!     FORMATTED RESULTS FILE NUMBER
!
      INTEGER t3drfo
!
!     MIGRHYCAR STEERING FILE
!
      INTEGER t3dmig
!
!     RESTART FILE
!
      INTEGER t3drst
!
!     HARMONIC CONSTANTS FILE NUMBER
!
      INTEGER t3dhar
!
!     TIDAL MODEL FILE NUMBER
!
      INTEGER t3dtid
!
!     ASCII TIDAL MODEL DATABASE FILE NUMBER
!
      INTEGER t3dbdd
!
!     BINARY TIDAL MODEL DATABASES 1 AND 2 FILE NUMBER
!
      INTEGER t3dbb1,t3dbb2
!
!     ASCII AND BINARY DROGUES FILE
!
      INTEGER t3dflo,t3dblo
!
!     FILE FOR 2D CONTINUATION
!
      INTEGER t3ds2d
!
!     WAQ STEERING FILE (INTERNAL WAQ NOT THOSE OF DELWAQ)
!
      INTEGER t3dwaq
!
!     WAQ DICO (INTERN WAQ NOT THOSE OF DELWAQ)
!
      INTEGER t3dwqd
!
!     LAMBERT 93 CONVERSION FILE
!
      INTEGER t3dl93
!
!     NUMBER OF VARIABLES TO BE ADVECTED BY A GIVEN SCHEME
!     AND THEIR LIST IN THE LIST OF ADVECTED VARIABLES
!     S_ADV: ASSOCIATED ADVECTION SCHEME
!     NOM_ADV: RANK TO FIND THE NAME IN TEXT3
!
      INTEGER n_adv(0:15),list_adv(100,0:15)
      INTEGER, ALLOCATABLE :: nom_adv(:),s_adv(:)
!
!     COUPLING PERIODS FOR SISYPHE AND TOMAWAC
!
      INTEGER percou_sis,percou_wac
!
!     BOUNDARY CONDITION ON THE BOTTOM
!
      INTEGER bc_bottom
!
!     THOMPSON BOUNDARY CONDITIONS
!
      INTEGER, ALLOCATABLE :: frtype(:)
!
!     TIDAL TYPE OF BOUNDARY CONDITIONS
!
      INTEGER, ALLOCATABLE, TARGET :: bnd_tide(:)
!
!     VERTICAL VELOCITY DERIVATIVES
!
      INTEGER linlog
!
!     RECORD NUMBER FOR RESTART
!
      INTEGER start_record
!
!     OPTION FOR TIDAL BOUNDARY CONDITIONS
!
      INTEGER tidaltype
!
!     TIDAL DATA BASE
!
      INTEGER tidaldb
!
!     GEOGRAPHIC SYSTEM
!
      INTEGER geosyst
!
!     ZONE NUMBER IN GEOGRAPHIC SYSTEM
!
      INTEGER numzone
!
!     NUMBERS OF THE POINT TO CALIBRATE HIGH WATER (SCHEMATIC TIDES)
!
      INTEGER icalhwb,icalhwg
!
!     FOR MONITORING K-EPSILON MODEL (SEE SUBROUTINE CSTKEP)
!
      INTEGER optprod,limkf,limef,limks,limes
!
!     FOR MONITORING S-A MODEL (SEE SUBROUTINE CSTSA)
!
      INTEGER limnuf,limnus
!
!     MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES
!
      INTEGER maxadv
!
!     OPTION FOR CHARACTERISTICS,
!     NUMBER OF GAUSS POINTS FOR WEAK CHARACTERISTICS
!
      INTEGER optcha,ngauss
!
!     OPTION FOR THE COMPUTATION OF THE LATERAL BOUNDARY CONDITIONS
!     OF K AND EPSILON
!
      INTEGER optbcke
!
!     SCHEME OPTION FOR ADVECTION OF VELOCITIES, K-EPSILON
!
      INTEGER optadv_vi,optadv_ke
!
!     OPTION FOR WIND
!
      INTEGER optwind
!
!     SCHEME OPTION FOR ADVECTION OF TRACERS
!
      INTEGER, ALLOCATABLE :: optadv_tr(:)
!
!     NUMBER OF CORRECTIONS OF DISTRIBUTIVE SCHEMES
!
      INTEGER nco_dist
!
!     NUMBER OF SUB-STEPS OF DISTRIBUTIVE SCHEMES
!
      INTEGER nsp_dist
!
!     NUMBER OF DIRECTIONS FOR DIFFERENTIATING in VECTOR MODES
!
      INTEGER ad_numofdir
!
!     GOTM STEERING FILE
!
      INTEGER t3dgtm
!
!     NUMBER OF TRACER FOR AMR MESH TRANSFORMATION
!
      INTEGER itrac_amr
!
!-----------------------------------------------------------------------
! (5) LOGICAL VALUES
!-----------------------------------------------------------------------
!
!     WILL BE USED BY FUNCTIONS Q3, SL3, TR3 AND VIT3
!
      LOGICAL, ALLOCATABLE :: okq3(:),oksl3(:),oktr3(:,:),okvit3(:)
      LOGICAL, ALLOCATABLE :: okdebsce(:),oktrsce(:,:)
!
!     LOGICAL STEERING PARAMETERS
!
      LOGICAL debu
!
!     IF YES, CORIOLIS
!
      LOGICAL coriol
!
!     IF YES, WIND TAKEN INTO ACCOUNT
!
      LOGICAL vent
!
!     AIR PRESSURE
!
      LOGICAL atmos
!
!     IF YES, SEDIMENT
!
      LOGICAL s3d_sedi
!
!     IF YES, TIDAL FLATS
!
      LOGICAL bandec
!
!     IF YES, PROPAGATION LINEARISED
!
      LOGICAL prolin
!
!     IF YES, MASS BALANCE
!
      LOGICAL bilmas
!
!     INFORMATION ABOUT MASS-BALANCE FOR EACH LISTING PRINTOUT
!
      LOGICAL infmas
!
!     IF YES, MESH WITH SIGMA TRANSFORMATION
!
      LOGICAL sigmag
!
!     IF YES, HYDROSTATIC INCONSISTENCY FILTER
!
      LOGICAL inchyd
!
!     ARRAYS SAYING IF A VARIABLE IS TO BE PRINTED IN RESULTS FILES
!
      LOGICAL sorg2d(maxvar),sorg3d(maxva3),sorest(maxva3)
!
!     ARRAYS SAYING IF A VARIABLE IS TO BE PRINTED IN THE LISTING
!
      LOGICAL sorimp(maxvar),sorim3(maxva3),soris3(maxva3)
!
!     IF YES, CONTINUATION FROM A 2D FILE
!
      LOGICAL suit2
!
!     IF YES, RAIN OR EVAPORATION
!
      LOGICAL rain
!
!     IF YES, STABILISED INITIAL CONDITION (OBSOLETE...)
!
!     LOGICAL CONSOL
!
!     IF YES, ELEMENTS MASKED BY USER
!
      LOGICAL mskuse
!
!     DIF(I) WILL SAY IF DIFFUSION SOLVER I IS USED FOR AT LEAST ONE OF THE VARIABLES
!
      LOGICAL dif(0:2)
!
!     IF YES, SPHERICAL COORDINATES (NOT IMPLEMENTED)
!
      LOGICAL spheri
!
!     IF YES, THERE ARE MASKED ELEMENTS
!
      LOGICAL msk
!
!     IF YES, LISTING PRINTOUT
!
      LOGICAL listin
!
!     IF YES, INFORMATION PRINTED ON LISTING
!
      LOGICAL infogr
!
!     IF YES, VALIDATION
!
      LOGICAL valid
!
!     IF YES, CONSOLIDATION MODEL
!
      LOGICAL s3d_tasse
!
!     CHOICE OF MODEL
!
      INTEGER s3d_itass
!
!     IF YES, S3D_GIBSONCONSOLIDATION MODEL
!
      LOGICAL s3d_gibson
!
!     IF YES, INFLUENCE OF TURBULENCE ON SETTLING VELOCITY
!
      LOGICAL s3d_turbwc
!
!     IF YES, COHESIVE SEDIMENT
!
      LOGICAL s3d_sedco
!
!     IF YES, NON-HYDROSTATIC VERSION
!
      LOGICAL nonhyd
!
!     CONSISTENT PROJECTION (OBSOLETE...)
!
!     LOGICAL CONPRO
!
!     FOR INITIALISATION OF K-EPSILON (SET IN CONDIM)
!
      LOGICAL akep
!
!     FOR INITIALISATION OF K-OMEGA (SET IN CONDIM)
!
      LOGICAL akom
!
!     FOR INITIALISATION OF S-A (SET IN CONDIM)
!
      LOGICAL aksa
!
!     IF YES, INITIAL TIME RESET TO ZERO IN A COMPUTATION CONTINUED
!
      LOGICAL raztim
!
!     IF YES, DYNAMIC PRESSURE IN WAVE EQUATION
!
      LOGICAL dpwaveq
!
!     IF YES, DYNAMIC BOUNDARY CONDITION
!
      LOGICAL cldyn
!
!     IF YES, SALINITY FOR DELWAQ
!
      LOGICAL sali_del
!
!     IF YES, TEMPERATURE FOR DELWAQ
!
      LOGICAL temp_del
!
!     IF YES, VELOCITY FOR DELWAQ
!
      LOGICAL velo_del
!
!     IF YES, DIFFUSION FOR DELWAQ
!
      LOGICAL diff_del
!
!     IF YES, WAVE DRIVEN CURRENTS
!
      LOGICAL courou
!
!     IF YES, BYPASS VOID VOLUMES
!
      LOGICAL bypass
!
!     PROJECTION OF VELOCITY ON LATERAL SOLID BOUNDARIES
!
      LOGICAL velprolat
!
!     PROJECTION OF VELOCITY ON BOTTOM
!
      LOGICAL velprobot
!
!     FOR CALCULATING FLUXES OF ADVECTED VARIABLES
!
      LOGICAL,ALLOCATABLE :: calcflu(:)
!
!     FOR TAKING INTO ACCOUNT RAIN IN ADVECTION OF VARIABLES
!
      LOGICAL,ALLOCATABLE :: calcrain(:)
!
!     OIL SPILL MODEL
!
      LOGICAL spill_model
!
!     THOMFR (TRUE IF ONE OF THE VALUES IN FRTYPE IS 2)
!
      LOGICAL thomfr
!
!     RESTART MODE (TO GENERATE A RESTART FILE)
!
      LOGICAL restart_mode
!
!     IF YES, INFERENCE OF MINOR CONSTITUENTS
!
      LOGICAL intmicon
!
!     FOR MONITORING K-EPSILON MODEL (SEE SUBROUTINE CSTKEP)
!
      LOGICAL clipk,clipe,yap,wsik,clipnu
!
!     TIDE GENERATING FORCE
!
      LOGICAL maree
!
!     OPEN BOUNDARY CONDITIONS ON THE BED
!
      LOGICAL bedbou
!
!     EXPLICIT SCHEME FOR SETTLING
      LOGICAL  s3d_read_toce
      INTEGER s3d_setdep
!
!     MIXED SEDIMENT, NON COHESIVE SEDIMENT
!
      LOGICAL s3d_mixte,s3d_sednco
!
!     CONTINUITY CORRECTION ON OPEN BOUNDARIES
!
      LOGICAL concor
!
!     IF YES KEEP LONG LAT IN RESULT FILE
!
      LOGICAL keep_lonlat
!
!     TO CHANGE LATER
!
      LOGICAL charr_tel, susp_tel
!
!     IF YES, BOTTOM SMOOTHINGS ARE DONE AFTER MODIFICATIONS BY THE USER
!
      LOGICAL lisfon_after
!
!     IF YES, INITIAL VELOCITIES ARE COMPUTED FROM OSU TIDAL SOLUTIONS
!     E.G. TPXO. OTHERWISE = 0.
!
      LOGICAL vitini_tpxo
!
!-----------------------------------------------------------------------
! (6) REALS
!-----------------------------------------------------------------------
!
!     TIME
!
      DOUBLE PRECISION, TARGET :: at
!
!     TIMESTEP
!
      DOUBLE PRECISION, TARGET :: dt
!
!     DURATION
!
      DOUBLE PRECISION, TARGET :: duree
!
!     GRAVITY ACCELERATION
!
      DOUBLE PRECISION grav
!
!     CORIOLIS COEFFICIENT
!
      DOUBLE PRECISION fcor
!
!     DRAG COEFFICIENT OF WIND
!
      DOUBLE PRECISION fair
!
!     IF YES, THE DRAG COEFFICIENT OF WIND IS COMPUTED
!     WITH THE WIND VELOCITY
!
      LOGICAL fairaccu
!
!     WIND VELOCITY ALONG X, AND Y REPLACED BY CST_WINDX,CST_WINDY
!     IN METEO_TELEMAC MODULE
!
!     DOUBLE PRECISION FUAIR,FVAIR
!
!     WATER DENSITY AT REFERENCE CONCENTRATION
!
      DOUBLE PRECISION, TARGET :: rho0
!
!     FRICTION COEFFICIENT FOR THE BOTTOM
!
      DOUBLE PRECISION rugof0
!
!     FRICTION COEFFICIENT FOR LATERAL SOLID BOUNDARIES
!
      DOUBLE PRECISION rugol0
!
!     ZERO
!
      DOUBLE PRECISION zero
!
!     MINIMAL VALUE FOR DEPTH
!
      DOUBLE PRECISION hmin
!
!     MEAN DEPTH FOR LINEARIZATION
!
      DOUBLE PRECISION haulin
!
!     COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES
!
      DOUBLE PRECISION dnuvih
!
!     COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES
!
      DOUBLE PRECISION dnuviv
!
!     COEFFICIENT FOR HORIZONTAL DIFFUSION OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE :: dnutah(:)
!
!     COEFFICIENT FOR VERTICAL DIFFUSION OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE :: dnutav(:)
!
!     INITIAL DEPTH
!
      DOUBLE PRECISION hautin
!
!     INITIAL ELEVATION
!
      DOUBLE PRECISION cotini
!
!     RAIN OR EVAPORATION IN MM PER DAY
!     REPLACED BY CST_RAINFALL IN METEO_TELEMAC MODULE)
!
!     DOUBLE PRECISION RAIN_MMPD
!
!     IMPLICITATION FOR DEPTH
!
      DOUBLE PRECISION tetah
!
!     IMPLICITATION FOR VELOCITIES
!
      DOUBLE PRECISION tetau
!
!     IMPLICITATIONS FOR DIFFUSION ?
!
      DOUBLE PRECISION tetad,tetadi
!
!     MASS-LUMPING FOR DEPTH, DIFFUSION, VELOCITIES
!
      DOUBLE PRECISION aggloh,agglod,agglou
!
!     ARRAY OF PRESCRIBED FLOWRATES, ELEVATIONS, VELOCITIES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: debimp(:), cotimp(:)
!
!     ARRAY OF PRESCRIBED  VELOCITIES
!
      DOUBLE PRECISION, POINTER :: vitimp(:)
!
!     ARRAY OF PRESCRIBED FLOWRATES, SURFACE, VELOCITIES FOR LIQUID BOUNDARIES IMPOSED ON THE BED
!
      DOUBLE PRECISION, ALLOCATABLE :: bedflo(:),bedqarea(:)
!
!     BETA EXPANSION COEFFICIENT FOR TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: betac(:)
!
!     REFERENCE CONCENTRATION OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: t0ac(:)
!
!     INITIAL VALUES OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: trac0(:)
!
!     VALUES OF TRACERS IN THE RAIN
!
      DOUBLE PRECISION, ALLOCATABLE :: train(:)
!
!     DENSITY OF THE SEDIMENT
!
      DOUBLE PRECISION s3d_rhos
!
!     CRITICAL SHEAR STRESS FOR DEPOSITION
!
      DOUBLE PRECISION s3d_tocd
!
!     CONCENTRATION (G/L) OF FRESH DEPOSITS
!
      DOUBLE PRECISION s3d_cfdep
!
! CV: non cohesive bed porosity
!
      DOUBLE PRECISION s3d_xkv
!
!     REFERENCE BED LAYER THICKNESS FOR NEW LAYER CREATION
!
      DOUBLE PRECISION s3d_epai0
!
!     TIMESTEP FOR CONSOLIDATION
!
      DOUBLE PRECISION s3d_dtc
!
!     CONCENTRATION (G/L) OF THE CONSOLIDATED MUD
!
      DOUBLE PRECISION s3d_cfmax
!
!     EROSION COEFFICIENT (EMPIRICAL PARTHENIADES COEFFICIENT)
!
      DOUBLE PRECISION s3d_mpart
!
!     CRITICAL SHEAR STRESS FOR EROSION (FRESH DEPOSIT)
! CV : now an array
!      DOUBLE PRECISION S3D_TOCE
!
!     FLOCULATION COEFFICIENT
!
      DOUBLE PRECISION s3d_turba
!
!     COEFFICIENT RELATIVE TO FLOC DESTRUCTION
!
      DOUBLE PRECISION s3d_turbb
!
!     CONSOLIDATION TIME SCALE (ONLY FOR MULTILAYER MODEL)
!
      DOUBLE PRECISION s3d_trest(30)
!
!     CONSTANT SEDIMENT SETTLING VELOCITY (M/S)
!
      DOUBLE PRECISION s3d_wchu0
!
!     MEAN DIAMETER OF THE SEDIMENT
!
      DOUBLE PRECISION s3d_d50
!
!     GEOGRAPHICAL LATITUDE IN GRAD, POSITIVE FOR NORTHERN AND NEGATIVE ON SOUTHERN HEMISPHERE
!     NOT USED
!
!     DOUBLE PRECISION PHILAT
!
!     NOT USED
!
!     DOUBLE PRECISION DELTA
!
!     LATITUDE AND LONGITUDE OF THE ORIGIN POINT
!
      DOUBLE PRECISION latit,longit
!
!     NORTH
!
      DOUBLE PRECISION nord
!
!     FREE SURFACE GRADIENT COMPATIBILITY IN WAVE EQUATION
!
      DOUBLE PRECISION tetazcomp
!
!     ABSCISSAE, ORDINATES AND ELEVATIONS OF SOURCES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: xsce(:)
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: ysce(:)
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: zsce(:)
!
!     VELOCITIES OF THE SOURCES ALONG X, Y, Z
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: usce(:)
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: vsce(:)
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: wsce(:)
!
!     WATER DISCHARGE OF SOURCES TAKEN FROM STEERING FILE
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: qsce(:)
!
!     WATER DISCHARGE OF SOURCES COMPUTED WITH T3D_DEBSCE (VARIATIONS IN TIME)
!
      DOUBLE PRECISION, ALLOCATABLE :: qsce2(:)
!
!     VALUE OF THE TRACERS AT THE SOURCES TAKEN FROM STEERING FILE
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: tasce(:,:)
!
!     PRESCRIBED VALUES OF TRACERS AT LIQUID BOUNDARIES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: tracer(:)
!
!     RATIO BETWEEN SKIN FRICTION AND MEAN DIAMETER
!
      DOUBLE PRECISION s3d_kspratio
!
!     SHIELDS PARAMETER
!
      DOUBLE PRECISION s3d_ac
!
!     FLUX AT BOUNDARIES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: flux_boundaries(:)
!
!     FOR STORING DISCHARGE-ELEVATION CURVES
!
      DOUBLE PRECISION, ALLOCATABLE :: qz(:,:,:)
!
!     THRESHOLD DEPTH FOR WIND
!
      DOUBLE PRECISION hwind
!
!     CUMULATED FLUXES OF WATER SINCE BEGINNING OF COMPUTATION
!
      DOUBLE PRECISION fluxtotcum
!
!     INITIAL, PAST AND PRESENT MASS OF WATER
!
      DOUBLE PRECISION masini_water, massen_water,masse_water
!
!     THRESHOLD FOR VISCOSITY CORRECTION ON TIDAL FLATS
!
      DOUBLE PRECISION hlim
!
!     THRESHOLD FOR VISCOSITY CORRECTION ON TIDAL FLATS
!
      DOUBLE PRECISION s3d_hsed
!
!     COEFFICIENT TO CALIBRATE TIDAL RANGE
!
      DOUBLE PRECISION, TARGET :: ctide
!
!     COEFFICIENT TO CALIBRATE TIDAL VELOCITIES
!
      DOUBLE PRECISION, TARGET :: ctidev
!
!     COEFFICIENT TO CALIBRATE SEA LEVEL
!
      DOUBLE PRECISION, TARGET :: msl
!
!     MINIMUM DEPTH TO COMPUTE TIDAL VELOCITIES BOUNDARY CONDITIONS
!
      DOUBLE PRECISION hmin_vit_bc
!
!     MINIMUM DEPTH TO COMPUTE TIDAL VELOCITIES INITIAL CONDITIONS
!
      DOUBLE PRECISION hmin_vit_ic
!
!     COEFFICIENT TO CALIBRATE THE ATMOSPHERE-WATER EXCHANGE MODEL
!     MOVED TO WAQTEL
!      DOUBLE PRECISION C_ATMOS
!
!     SETTLING VELOCITY OF SAND
!
      DOUBLE PRECISION s3d_wcs0
      TYPE(bief_obj), TARGET :: s3d_wcs
!
!     INITIAL PERCENTAGE OF NON COHESIVE SEDIMENT
!
      DOUBLE PRECISION s3d_pvsnco0
!
!     MASS-LUMPING FOR WEAK CHARACTERISTICS
!
      DOUBLE PRECISION s3d_agglow
!
!     VALUE OF THE ATMOSPHERIC PRESSURE
!     REPLACED BY CST_PATMOS IN METEO_TELEMAC MODULE)
!
!     DOUBLE PRECISION PATMOS_VALUE
!
!     THRESHOLD FOR THICKNESS OF LAYERS FOR COUPLING WITH AED2
!
      DOUBLE PRECISION, PARAMETER :: eps_aed2 = 1.d-9
!
!     MINIMUM VOLUME OF 3D ELEMENTS
!
      DOUBLE PRECISION minimum_volume
!
!     MINIMUM DISTANCE BETWEEN PLANES CLOSE TO THE BOTTOM
!
      DOUBLE PRECISION dismin_bot
!
!     MINIMUM DISTANCE BETWEEN PLANES CLOSE TO THE FREE SURFACE
!
      DOUBLE PRECISION dismin_sur
!
!     THRESHOLD HEIGHT BEFORE CRUSHED ELEMENTS
!
      DOUBLE PRECISION min_dz
!
!-----------------------------------------------------------------------
! (7) STRINGS
!-----------------------------------------------------------------------
!
!     TITLE
!
      CHARACTER(LEN=72) titcas
!
!     VARIABLES FOR GRAPHIC PRINTOUTS
!
      CHARACTER(LEN=72) sort3d,sort2d
!
!
!
      CHARACTER(LEN=72) varimp
!
!
!
      CHARACTER(LEN=72) varim3
!
!     INITIAL CONDITIONS
!
      CHARACTER(LEN=72) cdtini
!
!     CHOICE OF ELEMENT
!
      CHARACTER(LEN=72) element
!
!
!
!     CHARACTER(LEN=3)  BINGEO
!
!
!
!     CHARACTER(LEN=3)  BINRES
!
!
!
!     CHARACTER(LEN=3)  BINPRE
!
!
!
!     CHARACTER(LEN=3)  BINHYD
!
!     CHOICE OF EQUATIONS TO SOLVE
!
      CHARACTER(LEN=20), TARGET :: equa
!
!
!
      CHARACTER(LEN=32) varcla(10)
!
!     NAMES OF VARIABLES RECOGNISED FROM RESULTS AND GEOMETRY FILES
!
      CHARACTER(LEN=32) texte(maxvar)
!
!     NAMES OF VARIABLES RECOGNISED FROM PREVIOUS COMPUTATION FILE
!
      CHARACTER(LEN=32) textpr(maxvar)
!
!
!
      CHARACTER(LEN=32) text3(maxva3)
!
!
!
      CHARACTER(LEN=32) textp3(maxva3)
!
!
!
      CHARACTER(LEN=3) s3d_birsed
!
!
!
!     CHARACTER(LEN=3) BISUIS
!
!     NAMES OF TRACERS
!
      CHARACTER(LEN=32), ALLOCATABLE :: nametrac(:)
!
!     NAMES OF PRIVATE ARRAYS (GIVEN BY USER)
!
      CHARACTER(LEN=32) names_prive2d(4)
!
!     NAMES OF DIFFERENTING ARRAYS (GIVEN BY USER)
!
      CHARACTER(LEN=32) names_advar2d(maxvar)
!
!-----------------------------------------------------------------------
! (8) SLVCFG STRUCTURES
!-----------------------------------------------------------------------
!
!     SOLVER FOR DIFFUSION OF VELOCITIES
!
      TYPE(slvcfg) :: slvdvi
!
!     SOLVER FOR PROPAGATION
!
      TYPE(slvcfg) :: slvpro
!
!     SOLVER FOR DIFFUSION OF K-EPSILON
!
      TYPE(slvcfg) :: slvdke
!
!     SOLVER FOR DIFFUSION OF TRACERS
!
      TYPE(slvcfg), ALLOCATABLE :: slvdta(:)
!
!     SOLVER FOR VERTICAL VELOCITY COMPONENT
!     NO MORE SLVW SINCE V5P5
!
!     TYPE(SLVCFG) :: SLVW
!
!     SOLVER FOR DIFFUSION OF THE SEDIMENT
!
      TYPE(slvcfg) :: s3d_slvdse
!
!     SOLVER FOR PPE (NON-HYDROSTATIC)
!
      TYPE(slvcfg) :: slvpoi
!
!     SOLVER FOR PROJECTION (NON-HYDROSTATIC)
!     NO MORE SLVPRJ
!
!     TYPE(SLVCFG) :: SLVPRJ
!
!-----------------------------------------------------------------------
! (9) MESH STRUCTURE(S)
!-----------------------------------------------------------------------
! 2 SEPARATE MESHES, 2D AS USUAL AND 3D WITH SIGMA-MESH SPECIFIC
! FEATURES, SEE ALMESH.F
!
!     2D MESH WITH SIGMA-MESH SPECIFIC FEATURES
!
      TYPE(bief_mesh), TARGET :: mesh2d
!
!     3D MESH WITH SIGMA-MESH SPECIFIC FEATURES
!
      TYPE(bief_mesh), TARGET :: mesh3d
!
!-----------------------------------------------------------------------
! (10) ALIASES
!-----------------------------------------------------------------------
!
!     DECLARATION OF POINTERS FOR ALIASES
!     TARGETS ARE ALLOCATED AND POINTED TO IN POINT_TELEMAC3D
!
!     ALIASES FOR WORKING VECTORS, REAL 3D, INTEGER 3D, REAL 2D
!
!
!     BIEF_OBJ STRUCTURES FOR ARRAYS OF DIMENSION NPOIN3
!
      TYPE(bief_obj), POINTER :: t3_01,t3_02,t3_03,t3_04,t3_05
      TYPE(bief_obj), POINTER :: t3_06,t3_07,t3_08,t3_09,t3_10
      TYPE(bief_obj), POINTER :: t3_11,t3_12,t3_13,t3_14,t3_15
      TYPE(bief_obj), POINTER :: t3_16,t3_17,t3_18
      TYPE(bief_obj), POINTER :: prive1,prive2,prive3,prive4
!
!     BIEF_OBJ STRUCTURES FOR INTEGER ARRAYS
!
      TYPE(bief_obj), POINTER :: it1, it2, it3, it4
!
!     ALIASES OF BIEF_OBJ STRUCTURE FOR LOCAL WORK, DIMENSION NPOIN2
!
      TYPE(bief_obj), POINTER :: t2_01,t2_02,t2_03,t2_04,t2_05
      TYPE(bief_obj), POINTER :: t2_06,t2_07,t2_08,t2_09,t2_10
      TYPE(bief_obj), POINTER :: t2_11,t2_12,t2_13,t2_14,t2_15
      TYPE(bief_obj), POINTER :: t2_16,t2_17,t2_18,t2_19,t2_20
      TYPE(bief_obj), POINTER :: t2_21,t2_22,t2_23
      TYPE(bief_obj), POINTER :: prive2d1,prive2d2,prive2d3,prive2d4
!
!     2D NODE COORDINATES
!
      TYPE(bief_obj), POINTER :: x2,y2
!
!     3D NODE COORDINATES
!
      TYPE(bief_obj), POINTER :: x3,y3,z3
!
!     3D ELEMENT-ORIENTED WORKING FIELD ALLOCATED WITH THE MESH
!
      TYPE(bief_obj), POINTER :: w1
!
!     BASE TRIANGLE SURFACES
!
      TYPE(bief_obj), POINTER :: surfa2
!
!     TRIANGLE SURFACES, IN 3D
!
      TYPE(bief_obj), POINTER :: surfa3
!
!     LATERAL BOUNDARY NORMAL VECTORS DEFINED AT THE NODES
!
      TYPE(bief_obj), POINTER :: xnebor2
!
!     LATERAL BOUNDARY NORMAL VECTORS DEFINED AT THE NODES
!
      TYPE(bief_obj), POINTER :: ynebor2
!
!     2D NORMAL VECTORS DEFINED PER BOUNDARY SEGMENT
!
      TYPE(bief_obj), POINTER :: xsgbor2
!
!     2D NORMAL VECTORS DEFINED PER BOUNDARY SEGMENT
!
      TYPE(bief_obj), POINTER :: ysgbor2
!
!     3D NORMAL VECTORS DEFINED PER BOUNDARY ELEMENT
!
      TYPE(bief_obj), POINTER :: xsgbor3
!
!     3D NORMAL VECTORS DEFINED PER BOUNDARY ELEMENT
!
      TYPE(bief_obj), POINTER :: ysgbor3
!
!     3D NORMAL VECTORS DEFINED PER BOUNDARY ELEMENT
!
      TYPE(bief_obj), POINTER :: zsgbor3
!
!     CONNECTIVITY TABLES IN 2D  : (ELEMENT NUMBER AND LOCAL NODE NUMBER) --> GLOBAL NODE NUMBER
!
      TYPE(bief_obj), POINTER :: ikle2
!
!     CONNECTIVITY TABLES IN 3D : (ELEMENT NUMBER AND LOCAL NODE NUMBER) --> GLOBAL NODE NUMBER
!
      TYPE(bief_obj), POINTER :: ikle3
!
!     CONNECTIVITY TABLES IN 2D : (NODE BOUNDARY NUMBER) --> GLOBAL NODE NUMBER
!
      TYPE(bief_obj), POINTER :: nbor2
!
!     CONNECTIVITY TABLES IN 3D : (NODE BOUNDARY NUMBER) --> GLOBAL NODE NUMBER
!
      TYPE(bief_obj), POINTER :: nbor3
!
!     COORDINATES OF POINTS IN THE 3D MESH
!
      DOUBLE PRECISION, DIMENSION(:), POINTER :: x
!
!     COORDINATES OF POINTS IN THE 3D MESH
!
      DOUBLE PRECISION, DIMENSION(:), POINTER :: y
!
!     COORDINATES OF POINTS IN THE 3D MESH
!
      DOUBLE PRECISION, DIMENSION(:), POINTER :: z
!
!     NUMBER OF ELEMENTS IN THE 2D MESH
!
      INTEGER, POINTER:: nelem2
!
!     NUMBER OF ELEMENTS IN THE 3D MESH
!
      INTEGER, POINTER:: nelem3
!
!     MAXIMUM NUMBER OF ELEMENTS IN THE 2D MESH
!
      INTEGER, POINTER:: nelmax2
!
!     MAXIMUM NUMBER OF ELEMENTS IN THE 3D MESH
!
      INTEGER, POINTER:: nelmax3
!
!     NUMBER OF BOUNDARY POINTS IN THE 2D MESH
!
      INTEGER, POINTER:: nptfr2
!
!     NUMBER OF BOUNDARY POINTS IN THE 3D MESH
!
      INTEGER, POINTER:: nptfr3
!
!
!
      INTEGER, POINTER:: neleb, nelebx
!
!     MAXIMUM NUMBER OF BOUNDARY POINTS IN THE 2D MESH
!
      INTEGER, POINTER:: nptfrx2
!
!     MAXIMUM NUMBER OF BOUNDARY POINTS IN THE 3D MESH
!
      INTEGER, POINTER:: nptfrx3
!
!     DIMENSION OF 2D SPACE
!
      INTEGER, POINTER:: dim2
!
!     DIMENSION OF 3D SPACE
!
      INTEGER, POINTER:: dim3
!
!     TYPE OF 2D ELEMENT
!
      INTEGER, POINTER:: typelm2
!
!     TYPE OF 3D ELEMENT
!
      INTEGER, POINTER:: typelm3
!
!     NUMBER OF POINTS IN THE 2D MESH
!
      INTEGER, POINTER:: npoin2
!
!     NUMBER OF POINTS IN THE 3D MESH
!
      INTEGER, POINTER:: npoin3
!
!     MAXIMUM NUMBER OF POINTS IN THE 2D MESH
!
      INTEGER, POINTER:: npmax2
!
!     MAXIMUM NUMBER OF POINTS IN THE 3D MESH
!
      INTEGER, POINTER:: npmax3
!
!     MAXIMUM NUMBER OF POINTS NEIGHBOURS OF A POINT IN 2D
!
      INTEGER, POINTER:: mxptvs2
!
!     MAXIMUM NUMBER OF POINTS NEIGHBOURS OF A POINT IN 3D
!
      INTEGER, POINTER:: mxptvs3
!
!     MAXIMUM NUMBER OF ELEMENTS NEIGHBOURS OF A POINT IN 2D
!
      INTEGER, POINTER:: mxelvs2
!
!     MAXIMUM NUMBER OF ELEMENTS NEIGHBOURS OF A POINT IN 3D
!
      INTEGER, POINTER:: mxelvs3
!
!     VECTOR LENGTH OF THE MACHINE (LV2=LV3, FOR SIGMA MESH)
!
      INTEGER, POINTER:: lv
!
!-----------------------------------------------------------------------
! (11) CONSTANTS INITIALISED IN CSTKEP AND CREF_FREDSOE (FOR ZREF)
!-----------------------------------------------------------------------
!
!
!     FOR MONITORING K-EPSILON MODEL (SEE SUBROUTINE CSTKEP)
!
      DOUBLE PRECISION pernorm2,perprod,rimin,rimax
!
!     VON KARMAN CONSTANT
!
      DOUBLE PRECISION :: karman
!
!     K-EPSILON CONSTANT
!
      DOUBLE PRECISION :: cmu
!
!     K-EPSILON CONSTANTS
!
      DOUBLE PRECISION :: c1,c2,sigmak,sigmae,virt,sigmanu
!
!     SCHMIDT NUMBER
!
      DOUBLE PRECISION :: schmit
!
!     MINIMUM K
!
      DOUBLE PRECISION :: kmin
!
!     MAXIMUM K
!
      DOUBLE PRECISION :: kmax
!
!     MINIMUM EPSILON
!
      DOUBLE PRECISION :: emin
!
!     MAXIMUM EPSILON
!
      DOUBLE PRECISION :: emax
!
!     MINIMUM NU
!
      DOUBLE PRECISION :: numin
!
!     MAXIMUM NU
!
      DOUBLE PRECISION :: numax
!
!     PRANDTL NUMBER
!
      DOUBLE PRECISION, TARGET :: prandtl
!
! CV 7.0
!    FOND FICTIF : DZ1/FICTIF
!
      DOUBLE PRECISION :: fict
!
!     K-OMEGA MODEL CONSTANTS
!
      DOUBLE PRECISION :: alpha,beta,betas,omstar
!
!
!
      TYPE(bief_obj), TARGET :: zref
!
!-----------------------------------------------------------------------
!
!      12) TELEMAC-3D FILES
!
!-----------------------------------------------------------------------
!
!
!     MAXIMUM NUMBER OF FILES
!
      INTEGER, PARAMETER :: maxlu_t3d = 53
!
!     DATA STRUCTURE WITH DATA ON FILES
!
      TYPE(bief_file), TARGET :: t3d_files(maxlu_t3d)
!
!
!-----------------------------------------------------------------------
!
!      12) SAVED VALUES
!
!-----------------------------------------------------------------------
!
!  VARIABLES TO BE READ WHEN SUITE IS CALLED:
!  0 : DISCARD    1 : READ  (SAME NUMBERING AS IN NOMVAR)
!
!                  U V   H   ZF
      INTEGER :: alire2d(maxvar) = (/ 1,1,0,1,0,1,(0,lll=7,maxvar) /)
!                   U V   H   ZF
!      DATA ALIRE2D/1,1,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0/
!
!     IN 3D FILES
!                                                         U V
!                                             U V W       C C
!                                             C C C     D O O
!                                             O O O   D H N N
!                                         D   N N N   M H V V U V W
!                   Z U V W       K E     P   V V V   1 N C C D D D
      INTEGER :: alire3d(maxvar) = (/
     &              1,1,1,1,0,0,0,1,1,0,0,1,0,1,1,1,0,1,1,1,1,1,1,1,0,0,
     &              (0,lll=27,maxvar) /)
!      DATA ALIRE3D/1,1,1,1,0,0,0,1,1,0,0,1,0,1,1,1,0,1,1,1,1,1,1,1,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
!     &             0,0,0,0,0,0,0,0,0,0/
      LOGICAL s3d_charr,s3d_susp
      ! TODO: DO initialisation of value in point_telemac3d.f
      ! MURD3D_POS
      LOGICAL :: deja_murd3d_pos = .false.
!     INDIC_MURD3D_POS WILL BE A LIST OF SEGMENTS WITH NON ZERO FLUXES
      INTEGER, ALLOCATABLE :: indic_murd3d_pos (:)
      INTEGER sizeindic_murd3d_pos
!     TRISOU WAVE DRIVEN CURRENT
      LOGICAL :: dejalu_trisou = .false.
!
      SAVE
!
      END MODULE declarations_telemac3d