declarations_telemac2d.f

Go to the documentation of this file.

!                   *****************************
                    MODULE declarations_telemac2d
!                   *****************************
!
!
!***********************************************************************
! TELEMAC2D   V8P2
!***********************************************************************
!
!brief    DECLARATION OF PRINICIPAL TELEMAC2D VARIABLES
!
!history  J-M HERVOUET (LNH)
!+
!+
!+   LINKED TO BIEF 5.0
!
!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  C.COULET (ARTELIA)
!+        30/03/2012
!+        V6P2
!+   Modification for adding "bridge" file and separation of weirs and
!+   culvert file
!
!history  C.COULET / A.MASSON (ARTELIA) / P. CHASSE (CETMEF)
!+        19/07/2012
!+        V6P2
!+   Modification for adding breaches management during simulation
!
!history  C. COULET (ARTELIA)
!+        23/04/2013
!+        V6P3
!+     Correction of a bug in BUSE.F
!
!history  C.COULET / A.REBAI / E.DAVID (ARTELIA)
!+        07/06/2013
!+        V6P3
!+   Modification for new treatment of weirs
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        17/01/2014
!+        V7P0
!+   Adding a zones data file.
!
!history  Y.B. TADESSE (TUHH, INSTITUTE OF RIVER AND COASTAL ENGINEERING)
!+        14/02/2014
!+        V6P3R2
!+   Modification for adding later breach growth option
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        30/04/2014
!+        V7P0
!+   Adding 2 blocks FTILD2 and FNCAR2 for weak characteristics.
!
!history  C. VILLARET (HRW)
!+        22/09/2014
!+        V7P0
!+   Adding a logical keyword FRICOU for wave-current interactions
!
!history  R. ATA (EDF-LNHE)
!+        10/11/2014
!+        V7P0
!+   adding of waq variables
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        03/09/2015
!+        V7P1
!+   A number of arrays now allocatable, depending on new keywords
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        23/06/2016
!+        V7P2
!+   NAMETRAC now allocatable.
!
!history  C.COULET (ARTELIA)
!+        01/09/2016
!+        V7P2
!+   Modification for // treatment of weirs (type=2)
!+   Creation of a dedicated structure for weirs
!
!history  R. ATA (EDF LAB, LNHE)
!+        23/09/2016
!+        V7P2
!+   add new turbulence model of spalart-allmaras
!
!history  L.LEBALLEUR (ACTIMAR)
!+        30/03/2017
!+        V7P3
!+   Modification of MISC_NCMX_BTM and MISC_CONSTID_BTM
!+   Adding 7 extra harmonic constituents, for FES2014 database.
!
!history  B. GLANDER (BAW)
!+        28/11/2017
!+        V7P2
!+   add coupling with Nestor
!
!history  C.COULET (ARTELIA)
!+        09/10/2020
!+        V8P2
!+   Addition of output file to print discharge of weirs (type=2)
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief_def
      USE friction_def
!
!       NOTE: THIS MODULE IS ORGANISED IN 10 PARTS
!
!      01) VECTORS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!      02) MATRICES (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!      03) BLOCKS (WILL BE DECLARED AS BIEF_OBJ STRUCTURES)
!      04) INTEGERS
!      05) LOGICAL VALUES
!      06) REALS
!      07) STRINGS
!      08) SLVCFG STRUCTURES
!      09) MESH STRUCTURE
!      10) ALIASES
!
!
!       ALL BIEF_OBJ AND BIEF_MESH STRUCTURES ARE ALLOCATED
!       IN SUBROUTINE POINT_TELEMAC2D
!
!-----------------------------------------------------------------------
!
!       1) VECTORS
!
!-----------------------------------------------------------------------
!
!     COMPONENTS OF VELOCITY
!
      TYPE(bief_obj), TARGET :: u,v
!
!     DEPTH AT NEW TIME-STEP
!
      TYPE(bief_obj), TARGET :: h
!
!     TRACERS AT NEW TIME-STEP
!
      TYPE(bief_obj), TARGET :: t
!
!     K AT NEW TIME-STEP
!
      TYPE(bief_obj), TARGET :: ak
!
!     EPSILON AT NEW TIME-STEP
!
      TYPE(bief_obj), TARGET :: ep
!
!     VISCSA (OF SPALART-ALLMARAS) AT NEW TIME-STEP
!
      TYPE(bief_obj), TARGET :: viscsa
!
!     NUSA (OF SPALART-ALLMARAS) AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: nutild
!
!     NUSA (OF SPALART-ALLMARAS) AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: nun
!
!     U AND V AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: utild,vtild
!
!     H AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: htild
!
!     T AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: ttild
!
!     AK AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: aktild
!
!     EP AFTER ADVECTION BY CHARACTERISTICS
!
      TYPE(bief_obj), TARGET :: eptild
!
!     U AND V AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: un,vn
!
!     H AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: hn
!
!     T AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: tn
!
!     AK AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: akn
!
!     EP AT OLD TIME-STEP
!
      TYPE(bief_obj), TARGET :: epn
!
!     INCREMENT OF H IN ONE TIME STEP
!
      TYPE(bief_obj), TARGET :: dh
!
!     INCREMENTS OF U AND V IN ONE TIME STEP
!
      TYPE(bief_obj), TARGET :: du,dv
!
!     INCREMENT OF HN IN ONE TIME STEP
!
      TYPE(bief_obj), TARGET :: dhn
!
!     COMPONENTS OF ADVECTING FIELD
!
      TYPE(bief_obj), TARGET :: uconv,vconv
!
!     DEPTH IN THE PROPAGATION TERM I.E. IN H DIV(U)
!
      TYPE(bief_obj), TARGET :: hprop
!
!     PRESCRIBED VALUES ON BOUNDARIES FOR U AND V
!
      TYPE(bief_obj), TARGET :: ubor,vbor
!
!     VALUE ON BOUNDARIES FOR NU IN SA MODEL
!
      TYPE(bief_obj), TARGET :: nubor
!
!     VALUES ON BOUNDARIES: WALL FRICTION COEFFICIENTS
!
      TYPE(bief_obj), TARGET :: aubor,uetuta
!
!     PRESCRIBED VALUES ON BOUNDARIES FOR DEPTH
!
      TYPE(bief_obj), TARGET :: hbor
!
!     PRESCRIBED VALUES ON BOUNDARIES FOR TRACERS
!
      TYPE(bief_obj), TARGET :: tbor
!
!     PRESCRIBED VALUES ON BOUNDARIES: K AND EPSILON
!
      TYPE(bief_obj), TARGET :: kbor,ebor
!
!     FLUX ON THE BOUNDARIES: WATER AND TRACERS
!
      TYPE(bief_obj), TARGET :: flbor,flbortra
!
!     ON BOUNDARIES: COEFFICIENTS FOR HEAT FLUXES
!
      TYPE(bief_obj), TARGET :: atbor,btbor
!
!     ADIMENSIONAL FRICTION COEFFICIENTS (BOTTOM AND BOUNDARIES)
!
      TYPE(bief_obj), TARGET :: cf,cfbor
!
!     VOID STRUCTURE
!
      TYPE(bief_obj), TARGET :: s
!
!     BOTTOM TOPOGRAPHY: PER POINT, PER ELEMENT
!
      TYPE(bief_obj), TARGET :: zf,zfe
!
!     VELOCITY DIFFUSIVITY: VELOCITY, TRACERS, SAVED
!
      TYPE(bief_obj), TARGET :: visc,visct,visc_s
!
!     SOURCE TERMS IN THE MOMENTUM EQUATION
!
      TYPE(bief_obj), TARGET :: fu,fv
!
!     WAVE STRESSES (FROM ARTEMIS OR TOMAWAC)
!
      TYPE(bief_obj), TARGET :: fxwave,fywave
!
!     INITIAL DEPTH
!
      TYPE(bief_obj), TARGET :: h0
!
!     MASK FOR TRACERS
!
      TYPE(bief_obj), TARGET :: masktr
!
!     ATMOSPHERIC PRESSURE
!
!     TYPE(BIEF_OBJ), TARGET :: PATMOS
!
!     COMPONENTS OF WIND VELOCITY
!
!     TYPE(BIEF_OBJ), TARGET :: WINDX,WINDY
!
!     DENSITY
!
      TYPE(bief_obj), TARGET :: ro
!
!     INTEGRAL OF BASES, AND THE SAME AFTER ASSEMBLING IN PARALLEL
!
      TYPE(bief_obj), TARGET :: volu2d,v2dpar
!
!     INVERSE OF INTEGRAL OF BASES
!
      TYPE(bief_obj), TARGET :: unsv2d
!
!     RIGHT HAND SIDES OF LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: cv1,cv2,cv3,cv1s
!
!     RIGHT HAND SIDE OF CONTINUITY EQUATION
!
      TYPE(bief_obj), TARGET :: smh
!
!     P0 WORKING ARRAYS
!
      TYPE(bief_obj), TARGET :: te1,te2,te3
!
!     WORKING ARRAYS FOR OPTION 3 OF TREATMENT OF TIDAL FLATS
!
      TYPE(bief_obj), TARGET :: te4,te5
!
!     WORKING ARRAY (SIZE MULTIPLE OF NELMAX)
!
      TYPE(bief_obj), TARGET :: w1
!
!     WORKING ARRAY FOR PRESCRIBED DISCHARGE (FINITE VOLUMES AND FINITE ELEMENTS)
!
      TYPE(bief_obj), TARGET :: w1deb
!
!     FOR FLUXES ACCROSS SECTIONS
!
      TYPE(bief_obj), TARGET :: msksec
!
      TYPE(bief_obj), TARGET :: zflats
!
!     FRICTION COEFFICIENT, IN TERMS OF STRICKLER, CHEZY, ETC
!
      TYPE(bief_obj), TARGET :: chestr,chestr0
!
!     TYPES OF BOUNDARY CONDITIONS ON U, V, H AND TRACERS
!
      TYPE(bief_obj), TARGET :: liubor,livbor,lihbor,litbor
!
!     TYPES OF BOUNDARY CONDITIONS ON TRACERS
!
      TYPE(bief_obj), TARGET :: limtra
!
!     TYPES OF BOUNDARY CONDITIONS ON K AND EPSILON
!
      TYPE(bief_obj), TARGET :: limkep
!
!     TYPES OF BOUNDARY CONDITIONS ON SA MODEL
!
      TYPE(bief_obj), TARGET :: limsa
!
!     TYPES OF BOUNDARY CONDITIONS FOR PROPAGATION
!
      TYPE(bief_obj), TARGET :: limpro
!
!     BOUNDARY CONDITIONS ON H, U, V (COPIES OF LIHBOR,LIUBOR,LIVBOR)
!
      TYPE(bief_obj), TARGET :: clh,clu,clv
!
!     ORIGINAL BOUNDARY NODE NUMBER (OR COLOUR GIVEN BY USER)
!
      TYPE(bief_obj), TARGET :: boundary_colour
!
!     POSITIONS OF FLOATING BODIES
!
      TYPE(bief_obj), TARGET :: xflot,yflot
!
!     FOR LAGRANGIAN DRIFTS
!
      TYPE(bief_obj), TARGET :: shpflo,xlag,ylag,shplag
!
!     TIME STEP OF INITIAL RELEASE FOR FLOATING BODIES, IDEM FOR FINAL
!
      TYPE(bief_obj), TARGET :: debflo,finflo
!
!     ELEMENT WHERE THE FLOATING BODIES (DROGUES) ARE
!
      TYPE(bief_obj), TARGET :: eltflo
!
!     FOR TAGGING THE DROGUES
!
      TYPE(bief_obj), TARGET :: tagflo
!
!     TO CLASSIFY THE DROGUES
!
      TYPE(bief_obj), TARGET :: clsflo
!
!     TIME STEP AT THE BEGINNING (END) OF LAGRANGIAN DRIFTS
!
      TYPE(bief_obj), TARGET :: deblag,finlag
!
!     ELEMENT WHERE IS THE LAGRANGIAN DRIFT
!
      TYPE(bief_obj), TARGET :: eltlag
!
!     INTEGERS WORKING ARRAYS
!
      TYPE(bief_obj), TARGET :: it1,it2,it3,it4
!
!     COMPONENTS OF DISCHARGE AT TIME N
!
      TYPE(bief_obj), TARGET :: qu,qv
!
!     STORED WATER DEPTH, AND CORRECTED STORED WATER DEPTH
!
      TYPE(bief_obj), TARGET :: hstok,hcstok
!
!     SOURCE TERMS OF TRACERS
!
      TYPE(bief_obj), TARGET :: smtr
!
!     REFERENCE OF BOUNDARY NODES
!
      TYPE(bief_obj), TARGET :: logfr
!
!     FINITE VOLUME FLUXES
!
      TYPE(bief_obj), TARGET :: flux
!
!     TRACER FLUX
!
      TYPE(bief_obj), TARGET :: fluxt
!
!     TRACER FLUX OF LAST ITERATION (FOR SECOND ORDER)
!
      TYPE(bief_obj), TARGET :: fluxt_old
!
!     FLUX OF MASS FOR TRACER
!
      TYPE(bief_obj), TARGET :: fluhtemp
!
!     FLUX OF MASS AT BOUNDARY FOR TRACER
!
      TYPE(bief_obj), TARGET :: fluhbtemp
!
!     RECONSTRUCTED DEPTH, AND CORRECTED VERSION
!
      TYPE(bief_obj), TARGET :: ht,hc
!
!     VARIATIONS OF Z OF ORDER 2
!
      TYPE(bief_obj), TARGET :: dsz
!
!     ??????
!
      TYPE(bief_obj), TARGET :: secmou , ifamas
!
!     MASK FOR POINTS (1: normal; 0:masked)
!
      TYPE(bief_obj), TARGET :: maskpt
!
!     MASK FOR ELEMENTS (1: normal; 0:masked)
!
      TYPE(bief_obj), TARGET :: maskel
!
!     LIQUID BOUNDARY NUMBERS
!
      TYPE(bief_obj), TARGET :: numliq
!
!     MAXIMUM ELEVATIONS
!
      TYPE(bief_obj), TARGET :: maxz
!
!     CORRESPONDING TIMES FOR MAXIMUM ELEVATIONS
!
      TYPE(bief_obj), TARGET :: tmaxz
!
!     MAXIMUM VELOCITIES
!
      TYPE(bief_obj), TARGET :: maxv
!
!     CORRESPONDING TIMES FOR MAXIMUM VELOCITIES
!
      TYPE(bief_obj), TARGET :: tmaxv
!
!     FOR STORING ZONE NUMBERS
!
      TYPE(bief_obj), TARGET :: zone
!
!     FOR STORING RESULTS OF FOURIER ANALYSIS (AMPLITUDE AND PHASE)
!
      TYPE(bief_obj), TARGET :: ampl,phas
!
!     FRICTION COEFFICIENT FOR BOUNDARY CONDITIONS
!
      TYPE(bief_obj),TARGET :: chbord
!
!     COMPATIBLE COMPONENTS OF VELOCITY FOR DELWAQ
!
      TYPE(bief_obj),TARGET :: udel,vdel
!
!     VELOCITY WILL BE UCONV + DM1*GRAD(ZCONV)
!
      TYPE(bief_obj),TARGET :: dm1,zconv
!
!     COMPATIBLE FLUXES FOR DELWAQ
!
      TYPE(bief_obj),TARGET :: flodel
!
!     COEFFICIENTS OF LIMITATION OF FLUXES (SEGMENT, EBE)
!
      TYPE(bief_obj),TARGET :: flulim,flulimebe
!
!     EXPLICIT SOURCE TERMS FOR TRACERS
!
      TYPE(bief_obj),TARGET :: texp
!
!     TERMS DUE TO SOURCES FOR TRACERS
!
      TYPE(bief_obj),TARGET :: tscexp
!
!     IMPLICIT SOURCE TERMS FOR TRACERS
!
      TYPE(bief_obj),TARGET :: timp
!
!     FLUXES FOR FINITE VOLUMES
!
      TYPE(bief_obj),TARGET :: flux_old
!
!     FOR TIDAL BOUNDARY CONDITIONS
!
      TYPE(bief_obj),TARGET :: hbtide,ubtide,vbtide,numtide
!
!     FOR RAIN-EVAPORATION
!
      TYPE(bief_obj),TARGET :: pluie
!
!     FOR RAIN-EVAPORATION OPTION WITH RUNOFF-RAINFALL PREDICTION
!
      TYPE(bief_obj),TARGET :: accrof,cn,zfslop,accr,potmaxret,iabst
!
!     FOR CULVERTS/SIPHONS MANAGEMENT
!
      TYPE(bief_obj),TARGET :: entbus, sorbus, lrgbus, haubus, altbus
      TYPE(bief_obj),TARGET :: csbus , cebus , delbus, angbus, lbus
      TYPE(bief_obj),TARGET :: ubus, vbus, dbus, tbus, clpbus
      TYPE(bief_obj),TARGET :: secbus, longbus, aabus
      TYPE(bief_obj),TARGET :: cv, c56, cv5, c5, ctrash, fricbus, circ
!
!     NEIGHBORS OF SEGMENT (FOR WAF SCHEME)
!
      TYPE(bief_obj),TARGET :: neiseg
!
!     FOR BREACHES MANAGEMENT
!
      TYPE(bief_obj),TARGET :: optnbr, tdecbr, durbr, zfinbr, zdecbr
      TYPE(bief_obj),TARGET :: numpsd, nbndbr, indbr, zcrbr, polwdt
      TYPE(bief_obj),TARGET :: optero, dkaxcr, dkaycr, pondsb, nponbr
      TYPE(bief_obj),TARGET :: curbrw, finbrw, inibrw
!
!     FOR WEIR MANAGEMENT (TYPE WEIR = 1)
!
      TYPE(bief_obj),TARGET :: npsing
      TYPE(bief_obj),TARGET :: ndga1, ndgb1
      TYPE(bief_obj),TARGET :: zdig, phidig
!
!     FOR WEIR MANAGEMENT (TYPE WEIR = 2) SEE TYPE WEIR HEREUNDER
!
!     VARIABLES FOR SECONDARY CURRENTS (SECCURRENTS)
!
      TYPE(bief_obj),TARGET :: sec_tau, sec_r
!
!     WALL DISTANCES
!
      TYPE(bief_obj),TARGET :: wdist
!
!     REFERENCE LEVEL FOR NESTOR
!
      TYPE(bief_obj),TARGET :: zrl
!
!     CORRECTION FOR SECOND ORDER IN FV
!
      TYPE(bief_obj),TARGET :: corr_i,corr_j,corr_zl,corr_zr
      TYPE(bief_obj),TARGET :: corr_hl,corr_hr,corr_ul,corr_ur
      TYPE(bief_obj),TARGET :: corr_vl,corr_vr
!
!-----------------------------------------------------------------------
!
!       2) MATRICES
!
!-----------------------------------------------------------------------
!
!       MATRICES
!
!       SYSTEM SOLVED IN PROPAG WILL BE :
!
!     ( AM1  BM1  BM2 )   (DH)      ( CV1 )
!     (               )   (  )      (     )
!     ( CM1  AM2  A23 )   (U )  =   ( CV2 )
!     (               )   (  )      (     )
!     ( CM2  A32  AM3 )   (V )      ( CV3 )
!
!
!     MATRICES OF THE FINAL LINEAR SYSTEM OF THE PRIMITIVE EQUATIONS OPTION
!
      TYPE(bief_obj), TARGET :: am1,am2,am3,bm1,bm2,cm1,cm2,a23,a32
!
!     BOUNDARY MATRIX
!
      TYPE(bief_obj), TARGET :: mbor
!
!     BM1 AND BM2 SAVED
!
      TYPE(bief_obj), TARGET :: bm1s,bm2s
!
!-----------------------------------------------------------------------
!
!       3) BLOCKS
!
!-----------------------------------------------------------------------
!
!     DIFFUSION MATRIX
!
      TYPE(bief_obj), TARGET :: tm1
!
!     BLOCK FOR DIRICHLET VALUES HBOR, UBOR AND VBOR
!
      TYPE(bief_obj), TARGET :: dirbor
!
!     BLOCK OF MASKS
!
      TYPE(bief_obj), TARGET :: mask
!
!     BLOCKS OF WORKING ARRAYS
!
      TYPE(bief_obj), TARGET :: tb,tb2
!
!     BLOCK OF PRIVATE VECTORS
!
      TYPE(bief_obj), TARGET :: prive
!
!     BLOCK OF DIFFERENTIATED VECTORS
!
      TYPE(bief_obj), TARGET :: advar
!
!     BLOCKS OF VARIABLES AT TIME N+1, AT TIME N
!
      TYPE(bief_obj), TARGET :: f,fn
!
!     BLOCKS OF VARIABLES AFTER ADVECTION BY CHARACTERISTICS (STRONG AND WEAK)
!
      TYPE(bief_obj), TARGET :: ftild,ftild2
!
!     BLOCK OF VARIABLES TO BE ADVECTED BY CHARACTERISTICS (STRONG AND WEAK)
!
      TYPE(bief_obj), TARGET :: fncar,fncar2
!
!     BLOCKS OF MATRICES FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: mat
!
!     BLOCKS OF DIFFERENTIATED MATRICES FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: admat
!
!     BLOCKS OF RIGHT HAND SIDES FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: rhs
!
!     BLOCKS OF DIFFERENTIATED RIGHT HAND SIDES FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: adrhs
!
!     BLOCKS OF UNKNOWNS FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: unk
!
!     BLOCKS OF DIFFERENTIATED UNKNOWNS FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: adunk
!
!     BLOCK OF CLANDESTINE VARIABLES
!
      TYPE(bief_obj), TARGET :: varcl
!
!     BLOCK OF VARIABLES FOR OUTPUT
!
      TYPE(bief_obj), TARGET :: varsor
!                                                       WAC2
!
!     BLOCK OF VARIABLES FOR SENDING TO TOMAWAC
!
      TYPE(bief_obj), TARGET :: tel2tom
!
!     BLOCK OF VARIABLES FOR RECEIVING FROM TOMAWAC
!
      TYPE(bief_obj), TARGET :: tom2tel
!
!     NR OF VARIBLES SEND TO TOMAWAC
!
      INTEGER :: nvartel2tom
!
!     NR OF VARIBLES RECEIVED FROM TOMAWAC
!
      INTEGER :: nvartom2tel
!
!-----------------------------------------------------------------------
!
!       4) INTEGERS
!
!-----------------------------------------------------------------------
!
!     MAXIMUM NUMBER OF SOURCES
!
      INTEGER maxsce
!
!     MAXIMUM NUMBER OF BOUNDARIES
!
      INTEGER maxfro
!
!     MAXIMUM NUMBER OF OUTPUT VARIABLES
!
      INTEGER, PARAMETER :: maxvar = 100
      INTEGER :: nvar_t2d
!
!     MAXIMUM NUMBER OF TRACERS
!
      INTEGER maxtra
!
!     ADDRESS OF TRACERS IN VARSOR
!
      INTEGER :: adr_trac
!
!     MAXIMUM NUMBER OF POINTS TO DEFINE SOURCES REGIONS
!
      INTEGER maxptsce
!
!     NUMBER OF SOURCE REGIONS
!
      INTEGER nreg
!
!     NUMBER OF POINTS TO DEFINE SOURCE REGIONS
!
      INTEGER, ALLOCATABLE :: pts_reg(:)
!
!     NUMBER OF POINTS OF THE MESH CONTAINED IN SOURCE REGIONS
!
      INTEGER, ALLOCATABLE :: tnp(:)
!
!     NODES NUMBER OF POINTS CONTAINED IN THE SOURCE REGIONS
!
      INTEGER, ALLOCATABLE :: pt_in_poly(:,:)
!
!     NUMBER OF NODES FOR THE CONTROL SECTIONS (TWO TIMES THE NUMBER OF CONTROL SECTIONS)
!
      INTEGER ncp
!
!     ARRAY CONTAINING THE GLOBAL NUMBER OF POINTS IN THE CONTROL SECTIONS
!
      INTEGER, ALLOCATABLE :: ctrlsc(:)
!
!     MAXIMUM RANK OF LOGICAL UNITS AS DECLARED IN SUBMIT STRINGS IN THE DICTIONARY
!     06/10/2015: 45
!     15/03/2016: 46
!     30/08/2016: 48
!     28/11/2017: 52   coupling with Nestor
!     09/10/2020: 53   output weirs discharge
!
      INTEGER, PARAMETER :: maxlu_t2d = 53
!
!     GEOMETRY FILE NUMBER
!
      INTEGER, TARGET :: t2dgeo
!
!     BOUNDARY CONDITIONS FILE NUMBER
!
      INTEGER, TARGET :: t2dcli
!
!     PREVIOUS COMPUTATION FILE NUMBER
!
      INTEGER, TARGET :: t2dpre
!
!     RESULTS FILE NUMBER
!
      INTEGER, TARGET :: t2dres
!
!     BOTTOM TOPOGRAPHY FILE NUMBER
!
      INTEGER t2dfon
!
!     BINARY DATA FILE 1,2
!
      INTEGER t2dbi1,t2dbi2
!
!     ASCII ATMOSPHERIC FILE
!
      INTEGER t2atma
!
!     BINARY ATMOSPHERIC FILE
!
      INTEGER, TARGET :: t2atmb
!
!     FORMATTED DATA FILE 1,2
!
      INTEGER, TARGET :: t2dfo1,t2dfo2
!
!     BINARY RESULTS FILE NUMBER
!
      INTEGER t2drbi
!
!     FORMATTED RESULTS FILE NUMBER
!
      INTEGER t2drfo
!
!     REFERENCE FILE NUMBER
!
      INTEGER t2dref
!
!     LIQUID BOUNDARIES FILE NUMBER
!
      INTEGER, TARGET :: t2dimp
!
!     FRICTION DATA FILE NUMBER
!
      INTEGER t2dcof
!
!     DELWAQ FILES NUMBERS
!
      INTEGER t2ddl1,t2ddl2,t2ddl3,t2ddl4,t2ddl5,t2ddl6,t2ddl7,t2ddl8
      INTEGER t2ddl9,t2dl10,t2dl11
!
!     STAGE-DISCHARGE CURVES FILE NUMBER
!
      INTEGER t2dmab
!
!     SOURCES FILE NUMBER
!
      INTEGER t2dvef
!
!     DATA SOURCES FILE NUMBER
!
      INTEGER t2dsdn
!
!     SECTIONS INPUT FILE NUMBER
!
      INTEGER t2dsec
!
!     SECTIONS OUTPUT FILE NUMBER
!
      INTEGER t2dseo
!
!     MIGRHYCAR STEERING FILE NUMBER
!
      INTEGER t2dmig
!
!     HARMONIC CONSTANTS FILE NUMBER
!
      INTEGER t2dhar
!
!     TIDAL MODEL FILE NUMBER
!
      INTEGER t2dtid
!
!     ASCII TIDAL MODEL DATABASE FILE NUMBER
!
      INTEGER t2dbdd
!
!     BINARY TIDAL MODEL DATABASES 1 AND 2 FILE NUMBER
!
      INTEGER t2dbb1,t2dbb2
!
!     WEIRS DATA FILE NUMBER
!
      INTEGER t2dseu
!
!     WEIRS OUTPUT FILE NUMBER
!
      INTEGER t2dwop
!
!     CULVERTS/TUBES/BRIDGES DATA FILE NUMBER
!
      INTEGER t2dbus
!
!     BREACHES DATA FILE NUMBER
!
      INTEGER t2dbrc
!
!     ASCII AND BINARY DROGUES FILE
!
      INTEGER t2dflo,t2dblo
!
!     PREVIOUS DROGUES FILE
!
      INTEGER t2dplo
!
!     DATA DROGUES POSITIONING ASCII FILE NUMBER
!
      INTEGER t2dply
!
!     ZONES FILE
!
      INTEGER t2dzfi
!
!     WAQ STEERING FILE (INTERNAL WAQ NOT THOSE OF DELWAQ)
!
!      INTEGER T2DWAQ
!
!     WAQ DICO (INTERN WAQ NOT THOSE OF DELWAQ)
!
!      INTEGER T2DWQD
!
!     KHIONE STEERING FILE
!
!      INTEGER T2DICE
!
!     KHIONE DICO
!
!      INTEGER T2DICD
!
!     FLUXLINE FILE NUMBER (NEW METHOD FOR CROSS SECTIONS)
!
      INTEGER t2dflx
!
!     NESTOR ACTION FILE
!
      INTEGER t2nact
!
!     NESTOR POLYGON FILE
!
      INTEGER t2npol
!
!     NESTOR REFERENCE LEVEL FILE
!
      INTEGER t2nref
!
!     NESTOR RESTART FILE
!
      INTEGER t2nrst
!
!     LAMBERT 93 CONVERSION FILE
!
      INTEGER t2dl93
!
!     GRAPHIC PRINTOUT PERIOD
!
      INTEGER, TARGET ::  leoprd
!
!     LISTING PRINTOUT PERIOD
!
      INTEGER, TARGET :: lisprd
!
!     NUMBER OF TIME STEPS
!
      INTEGER, TARGET :: nit
!
!     TYPE OF ADVECTION (1:u and v, 2:h, 3:tracers, 4:k and epsilon)
!
      INTEGER iconvf(4)
!
!     SCHEME FOR ADVECTION OF TRACERS
!
      INTEGER, ALLOCATABLE :: iconvft(:)
!
!     TURBULENCE MODEL
!
      INTEGER, TARGET :: iturb
!
!     LAW OF BOTTOM FRICTION
!
      INTEGER kfrot
!
!     NUMBER OF SUB-ITERATIONS FOR NON-LINEARITIES
!
      INTEGER nsoui
!
!     NOT USED (PROPAGATION OPTION)
!
!     INTEGER OPTPRO
!
!     INITIAL GUESS FOR H
!
      INTEGER iordrh
!
!     INITIAL GUESS FOR U
!
      INTEGER iordru
!
!     NUMBER OF SUB-ITERATIONS FOR NON-LINEARITIES
!
      INTEGER nsousi
!
!     NUMBER OF FIRST TIME STEP FOR GRAPHIC PRINTOUTS
!
      INTEGER , TARGET :: ptinig
!
!     NUMBER OF FIRST TIME STEP FOR LISTING PRINTOUTS
!
      INTEGER ptinil
!
!     TURBULENCE MODEL FOR SOLID BOUNDARIES
!
      INTEGER lisrug
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED DISCHARGE
!
      INTEGER ndebit
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED ELEVATION
!
      INTEGER ncote
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED VELOCITY
!
      INTEGER nvites
!
!     NUMBER OF BOUNDARIES WITH PRESCRIBED TRACER
!
      INTEGER ntrace
!
!     NUMBER OF CLANDESTINE VARIABLES
!
      INTEGER nvarcl
!
!     VECTOR LENGTH
!
      INTEGER lvmac
!
!     OPTION FOR THE TREATMENT OF TIDAL FLATS
!
      INTEGER optban
!
!     NUMBER OF DROGUES, MAXIMUM NUMBER OF DROGUES (FOR DIMENSIONING ARRAYS)
!
      INTEGER nflot,nflot_max
!
!     PRINTOUT PERIOD FOR DROGUES
!
      INTEGER floprd
!
!     NUMBER OF LAGRANGIAN DRIFTS
!
      INTEGER nlag
!
!     BOTTOM SMOOTHINGS
!
      INTEGER lisfon
!
!     DISCRETIZATIONS IN SPACE
!
      INTEGER discre(5)
!
!     NUMBER OF SINK/SOURCES
!
      INTEGER nrejet
!
!     NUMBER OF SINK/SOURCES GIVEN BY THEIR NODE NUMBERS
!
      INTEGER nrejeti
!
!     NUMBER OF SINK/SOURCES WITH GIVEN VELOCITY
!
      INTEGER nrejeu
!
!     NUMBER OF SINK/SOURCES WITH GIVEN TRACER
!
      INTEGER nrejtr
!
!     SUPG OPTION
!
      INTEGER optsup(4)
!
!     ORIGINAL DATE OF TIME
!
      INTEGER, TARGET :: mardat(3)
!
!     ORIGINAL HOUR OF TIME
!
      INTEGER, TARGET :: martim(3)
!
!     ORIGINAL TIME OFFSET in SECONDS
!
!     DOUBLE PRECISION TEL_OFFSET
!
!     VARIOUS DISCRETISATION TIMES
!
      INTEGER ielm,ielm0,ielm1,ielmh,ielmu,ielmt,ielmk,
     &        ielme,ielmnu
!
!     MATRIX STORAGE
!
      INTEGER optass
!
!     MATRIX-VECTOR PRODUCT
!
      INTEGER produc
!
!     NUMBER OF WEIRS
!
      INTEGER nweirs
!
!     NUMBER OF NEIGHBOURS PROCESSORS IN CASE OF // (FOR WEIRS)
!
      INTEGER n_nghb_weirs
!
!     NUMBER OF NODES USED FOR COMPUTATION,
!
      INTEGER nweirs_nodes
!
!     NUMBER OF NEIGHBOURS PROCESSORS IN CASE OF // (FOR WEIRS NODES)
!
      INTEGER n_nghb_w_nodes
!
!     NUMBER OF NODES SEND IN CASE OF //
!
      INTEGER n_wn_send
!
!     NUMBER OF NEIGHBOURS PROCESSORS IN CASE OF // (FOR WEIRS NODES SEND)
!
      INTEGER n_nghb_wn_send
!
!     NUMBER OF CULVERTS/SIPHONS
!
      INTEGER nbuse
!
!     NUMBER OF BOUNDARIES WITH GIVEN OPTION FOR LIQUID BOUNDARIES
!
      INTEGER ntypfr
!
!     OPTION FOR LIQUID BOUNDARIES
!
      INTEGER, ALLOCATABLE :: frtype(:)
!
!     TREATMENT OF THE LINEAR SYSTEM
!
      INTEGER, TARGET ::  solsys
!
!     VELOCITY PROFILES
!
      INTEGER, ALLOCATABLE :: provel(:)
!
!     NUMBER OF LIQUID BOUNDARIES
!
      INTEGER nfrliq
!
!     NEAREST POINTS OF SOURCES
!
      INTEGER, ALLOCATABLE :: isce(:)
!
!     OPTION FOR THE DIFFUSION OF TRACERS
!
      INTEGER opdtra
!
!     OPTION FOR THE DIFFUSION OF VELOCITIES
!
      INTEGER opdvit
!
!     TYPE OF SOURCES
!
      INTEGER optsou
!
!     NUMBER OF PRIVATE ARRAYS, NUMBER OF PRIVATE ARRAYS WITH GIVEN NAME
!
      INTEGER npriv,n_names_priv
!
!     NUMBER OF DIFFERENTIATING ARRAYS, AND THOSE WITH A GIVEN NAME
!
      INTEGER nadvar,n_names_advar
!
!     NUMBER OF DIRECTIONS FOR DIFFERENTIATING in VECTOR MODES
!
      INTEGER ad_numofdir
!
!     DELWAQ PRINTOUT PERIOD
!
      INTEGER waqprd
!
!     TREATMENT OF NEGATIVE DEPTHS
!
      INTEGER opt_hneg
!
!     CURRENT ITERATION NUMBER
!
      INTEGER, TARGET :: lt
!
!     RECORD NUMBER IN THE WAVE DRIVEN CURRENTS FILE
!
      INTEGER npth
!
!     NUMBER OF USER-DEFINED ZONES (E.G. FOR SETTING FRICTION) IN SUBROUTINE DEF_ZONES
!
      INTEGER nzone
!
!     NUMBER OF PERIODS FOR FOURIER ANALYSIS
!
      INTEGER nperiaf
!
!     OPTION FOR PARAMETER ESTIMATION, IDENTIFICATION METHOD
!
      INTEGER optid
!
!     NUMBER OF REMARKABLE POINTS
!
      INTEGER npts
!
!     LIST OF REMARKABLE POINTS
!
      INTEGER list_pts(100)
!
!     CHOICE OF COST-FUNCTION
!
      INTEGER optcost
!
!     MAXIMUM NUMBER OF ITERATIONS FOR IDENTIFICATION
!
      INTEGER maxest
!
!     COUPLING PERIOD (WITH SISYPHE)
!
      INTEGER, TARGET ::  percou
!
!     COUPLING PERIOD FOR TOMAWAC
!
      INTEGER, TARGET :: percou_wac
!
!     FINITE VOLUME SCHEME
!
      INTEGER optvf
!
!     ORIGINS OF THE COORDINATE SYSTEM
!
      INTEGER i_orig,j_orig
!
!     NUMBER OF TRACERS
!
      INTEGER, TARGET :: ntrac
!
!     NUMBER TOTAL OF TRACERS
!
      INTEGER ntrtot
!
!     INDEX FOR TEMPERATURE, FOR SALINITY
!     IND_FT,IND_FC ?
      INTEGER ind_t,ind_s,ind_sec
!
!     INDEX FOR SUSPENDED SEDIMENTS
!
      INTEGER ind_sed
!
!     NUMBER OF POINTS GIVEN FOR EACH DISCHARGE-ELEVATIONS CURVES
!
      INTEGER, ALLOCATABLE :: pts_curves(:)
!
!     OPTION FOR TREATING STAGE-DISCHARGE CURVES
!
      INTEGER, ALLOCATABLE :: sta_dis_curves(:)
!
!     DEBUGGER
!
      INTEGER, TARGET :: debug
!
!     DEPTH IN FRICTION TERMS
!
      INTEGER hfrot
!
!     LAW OF FRICTION ON LATERAL BOUNDARIES
!
      INTEGER kfrotl
!
!     TREATMENT OF FLUXES AT THE BOUNDARIES
!
      INTEGER dirflu
!
!     OPTION FOR TIDAL BOUNDARY CONDITIONS
!
      INTEGER, ALLOCATABLE, TARGET :: bnd_tide(:)
      INTEGER tidaltype
!
!     OPTION FOR TSUNAMI GENERATION
!
      INTEGER opttsunami
!
!     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
!
!     RECORD NUMBER FOR RESTART
!
      INTEGER, TARGET :: start_record
!
!     NUMBER OF BREACHES
!
      INTEGER nbrech
!
!     SPATIAL PROJECTION TYPE, AND A SAVE OF ITS ORIGINAL VALUE
!
      INTEGER protyp,protyp0
!
!     OPTION FOR CHARACTERISTICS,
!     NUMBER OF GAUSS POINTS FOR WEAK CHARACTERISTICS
!
      INTEGER optcha,ngauss
!
!     MAXIMUM NUMBER OF ITERATIONS FOR ADVECTION SCHEMES
!
      INTEGER maxadv
!
!     STOCHASTIC DIFFUSION MODEL (FOR PARTICLES ALGAE, OIL SPILLS...)
!
      INTEGER stocha
!
!     OPTION FOR WEIR TREATMENT
!
      INTEGER typseuil
!
!     STORAGE OF MAXIMUM NUMBER OF ELEMENTS ON ALL WEIRS
!
      INTEGER maxnps
!
!     PSI SCHEME OPTION
!
!     INTEGER OPTPSI
!
!     ADVECTION SCHEME OPTIONS FOR VELOCITY AND K-EPSILON AND SA
!
      INTEGER optadv_vi,optadv_ke,optadv_sa
!
!     ADVECTION SCHEME OPTIONS FOR TRACERS
!
      INTEGER, ALLOCATABLE :: optadv_tr(:)
!
!     OPTION FOR WIND
!
      INTEGER optwind
!
!     NUMBER OF CORRECTIONS FOR DISTRIBUTIVE SCHEMES
!     NUMBER OF SUB-STEPS FOR DISTRIBUTIVE SCHEMES
!
      INTEGER nco_dist,nsp_dist
!
!     COUNTER FOR GRAPHICAL AND LISTING OUTPUTS
!
      INTEGER, TARGET :: compleo,complis
!
!     OPTION FOR THE TREATMENT OF CULVERTS
!
      INTEGER optbuse
!
!     OPTION FOR HYDROSTATIC RECONSTRUCTION
!
      INTEGER hropt
!
!     OPTION FOR RAINFALL-RUNOFF MODEL
!
      INTEGER runoffopt
!
!     ANTECEDENT MOISTURE CONDITION (FOR SCS_CN RAINFALL-RUNOFF MODEL)
!
      INTEGER amc
!
!     RATIO FOR INITIAL ABSTRACTION (FOR SCS_CN RAINFALL-RUNOFF MODEL)
!
      INTEGER iascnopt
!
!     NUMBER OF REFINEMENT LEVELS FOR THE CONVERGENCE STUDY
!     CURRENT LEVEL
!
      INTEGER rlevels,rlevel
!
!     CHOICE OF THE NUMERICAL SCHEME FOR FINITE VOLUME
!
      INTEGER icin
!
!     CHOICE OF SPACE AND TIME ORDER FOR FINITE VOLUME
!
      INTEGER sorder, torder
!
!     CHOICE OF FLUX LIMITOR FOR SECOND ORDER IN SPACE
!
      INTEGER ilimhz, ilimuv, ilimt
!
!     WEAK OR STRONG IMPOSITION OF KINETIC SCHEME BND CONDITION
!
      INTEGER bndcin
!
!-----------------------------------------------------------------------
!
!       5) LOGICAL VALUES
!
!-----------------------------------------------------------------------
!
!     USED IN DEBSCE
!
      LOGICAL, ALLOCATABLE :: okdebsce(:)
!
!     USED IN TRSCE
!
      LOGICAL, ALLOCATABLE :: oktrsce(:,:)
!
!     USED IN Q
!
      LOGICAL, ALLOCATABLE :: okq(:)
!
!     USED IN SL
!
      LOGICAL, ALLOCATABLE :: oksl(:)
!
!     USED IN TR
!
      LOGICAL, ALLOCATABLE :: oktr(:,:)
!
!     USED IN VIT
!
      LOGICAL, ALLOCATABLE :: okvit(:)
!
!     IF YES, WATER DISCHARGE OF SOURCES IS GIVEN
!
      LOGICAL okrejet
!
!     IF YES, COMPUTATION CONTINUED
!
      LOGICAL debu
!
!     IF YES, LISTING PRINTOUT IS REQUIRED
!
      LOGICAL listin
!
!     IF YES, DIFFUSION OF VELOCITY
!
      LOGICAL difvit
!
!     IF YES, AIR PRESSURE IS TAKEN INTO ACCOUNT
!
      LOGICAL atmos
!
!     IF YES, ADVECTION
!
      LOGICAL conv
!
!     IF YES, ADVECTION OF U AND V (1), H (2), TRACERS (3), K AND EPSILON (4)
!
      LOGICAL convv(4)
!
!     IF YES, DIFFUSION OF TRACERS
!
      LOGICAL dift
!
!     IF YES, CORIOLIS
!
      LOGICAL coriol
!
!     IF YES, INFORMATION PRINTED ON LISTING
!
      LOGICAL infogr
!
!     IF YES, MASS-BALANCE
!
      LOGICAL bilmas
!
!     IF YES, H CLIPPING
!
      LOGICAL cliph
!
!     IF YES, WIND IS TAKEN INTO ACCOUNT
!
      LOGICAL vent
!
!     IF YES, PROPAGATION
!
      LOGICAL propa
!
!     IF YES, SPHERICAL COORDINATES
!
      LOGICAL spheri
!
!     IF YES, LINEARISED PROPAGATION
!
      LOGICAL prolin
!
!     IF YES, VALIDATION
!
      LOGICAL valid
!
!     IF YES, VERIFICATION OF LIMITS
!
      LOGICAL verlim
!
!     IF YES, CONTINUITY CORRECTION
!
      LOGICAL corcon
!
!     IF YES, INITIAL TIME SET TO ZERO
!
      LOGICAL, TARGET :: raztim
!
!     IF YES, INFORMATION ABOUT K-EPSILON MODEL PRINTED ON LISTING
!
      LOGICAL infoke
!
!     IF YES, INFORMATION ABOUT SPALART-ALLMARAS MODEL PRINTED ON LISTING
!
      LOGICAL infonu
!
!     IF YES, THERE ARE MASKED ELEMENTS
!
      LOGICAL msk
!
!     IF YES, STOP WHEN A STEADY STATE IS REACHED
!
      LOGICAL stoper
!
!     IF YES, ELEMENTS MASKED BY USER
!
      LOGICAL mskuse
!
!     TO KNOW IF A VARIABLE WILL BE EXITED ON FILE, ON LISTING
!
      LOGICAL sorleo(maxvar),sorimp(maxvar)
!
!     VARIABLE DENSITY. IF YES, DENSITY EFFECTS INCLUDED
!
      LOGICAL rovar
!
!     TIDE GENERATING FORCE
!
      LOGICAL maree
!
!     IF YES, THERE IS A TREATMENT OF BOUNDARY CONDITIONS WITH
!     THOMPSON METHOD
!
      LOGICAL thomfr
!
!     IF YES, C-U PRECONDITIONING
!
      LOGICAL preccu
!
!     IF YES, VARIABLE TIME-STEP
!
      LOGICAL dtvari
!
!     IF YES, TIDAL FLATS
!
      LOGICAL bandec
!
!     IF YES, PRINT INFORMATIONS ON LISTING
!
      LOGICAL entet
!
!     IF YES, OIL SPILL MODEL
!
      LOGICAL spill_model
!
!     IF YES, WAVE DRIVEN CURRENTS ARE TAKEN INTO ACCOUNT
!
      LOGICAL courou
!
!     IF YES, WAVE FRICTION ENHANCEMENT IS TAKEN INTO ACCOUNT
!
      LOGICAL fricou
!
!     IF YES, VERTICAL STRUCTURES ARE TAKEN INTO ACCOUNT
!
      LOGICAL vertic
!
!     IF YES, THE USER WANTS TO DEFINE ZONES
!
      LOGICAL defzon
!
!     IF YES, COMPATIBLE COMPUTATION OF FLUXES (RELATED TO FLUXES THROUGH SECTIONS)
!
      LOGICAL comflu
!
!     IF YES, PRINTS CUMULATED FLOWRATES
!
      LOGICAL cumflo
!
!     IF YES, OUTPUT OF SALINITY FOR DELWAQ
!
      LOGICAL sali_del
!
!     IF YES, OUTPUT OF TEMPERATURE FOR DELWAQ
!
      LOGICAL temp_del
!
!     IF YES, OUTPUT OF VELOCITY FOR DELWAQ
!
      LOGICAL velo_del
!
!     IF YES, OUTPUT OF DIFFUSIVITY FOR DELWAQ
!
      LOGICAL diff_del
!
!     IF YES, RAIN OR EVAPORATION TAKEN INTO ACCOUNT
!
      LOGICAL rain
!
!     IF YES, INFERENCE OF MINOR CONSTITUENTS
!
      LOGICAL intmicon
!
!     IF YES, THE USER WANTS TO DEFINE BREACHES
!
      LOGICAL breche
!
!     IF YES, THE PARTICLES ARE ALGAE
!
      LOGICAL algae
!
!     SECONDARY CURRENTS
!
      LOGICAL seccurrents
!
!     IF YES, CROSS SECTIONS WITH FLUXLINE METHOD
!
      LOGICAL doflux
!
!     OPTION FOR CONVERGENCE STUDIES
!
      LOGICAL convergence
!
!     IF YES KEEP LONG LAT IN RESULT FILE
!
      LOGICAL keep_lonlat
!
!     IF YES, USE NESTOR
!
      LOGICAL nestor
!
!     IF YES, COEFFICIENT OF WIND INFLUENCE VARIES WITH WIND SPEED
!
      LOGICAL fairaccu
!
!     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 STEP
!
      DOUBLE PRECISION, TARGET :: dt, dtcas,dtn
!
!     END TIME
!
      DOUBLE PRECISION, TARGET :: tmax
!
!     INITIAL TIME STEP
!
      DOUBLE PRECISION dtini
!
!     GRAVITY ACCELERATION
!
      DOUBLE PRECISION grav
!
!     FRICTION COEFFICIENT
!
      DOUBLE PRECISION ffon
!
!     COEFFICIENT OF CORIOLIS
!
      DOUBLE PRECISION fcor
!
!     COEFFICIENT OF WIND INFLUENCE
!
      DOUBLE PRECISION, TARGET :: fair
!
!     COMPONENTS OF WIND VELOCITY,ITS NORM AND DIRECTION
!     REPLACED BY CST_WINDX,CST_WINDY,CST_WINDS,CST_WINDD
!     IN METEO_TELEMAC MODULE
!
!     DOUBLE PRECISION FUAIR,FVAIR,WIND_SPD(2)
      DOUBLE PRECISION wind_spd(2)
!
!     INITIAL VALUES OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE :: trac0(:)
!
!     COEFFICIENT FOR DIFFUSION OF TRACERS
!
      DOUBLE PRECISION, ALLOCATABLE :: difnu(:) 
!
!     IMPLICITATION FOR TRACERS
!
      DOUBLE PRECISION tetat
!
!     VELOCITY DIFFUSIVITY
!
      DOUBLE PRECISION, TARGET ::  propnu
!
!     THRESHOLD VALUES FOR STOPPING THE COMPUTATION
!
      DOUBLE PRECISION bornes(8)
!
!     IMPLICITATION FOR DEPTH
!
      DOUBLE PRECISION tetac
!
!     IMPLICITATION FOR VELOCITY
!
      DOUBLE PRECISION tetau
!
!     IMPLICITATION FOR DIFFUSION OF VELOCITY
!
      DOUBLE PRECISION tetad
!
!     MASS-LUMPING COEFFICIENT ON H
!
      DOUBLE PRECISION aggloc
!
!     MASS-LUMPING COEFFICIENT ON VELOCITY
!
      DOUBLE PRECISION agglou
!
!     MASS-LUMPING FOR WEAK CHARACTERISTICS
!
      DOUBLE PRECISION agglow
!
!     MINIMUM VALUE OF DEPTH
!
      DOUBLE PRECISION hmin
!
!     TIME STEP REDUCTION FOR K-EPSILON MODEL
!
!     DOUBLE PRECISION REDUC
!
!     MEAN DEPTH FOR LINEARIZATION
!
      DOUBLE PRECISION haulin
!
!     WATER DENSITY
!
      DOUBLE PRECISION roeau
!
!     LATITUDE OF ORIGIN POINT
!
      DOUBLE PRECISION lambd0
!
!     ROUGHNESS COEFFICIENT OF BOUNDARIES
!
      DOUBLE PRECISION sb
!
!     INITIAL ELEVATION
!
      DOUBLE PRECISION cotini
!
!     INITIAL DEPTH
!
      DOUBLE PRECISION hautin
!
!     PRESCRIBED FLOWRATES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: debit(:)
!
!     PRESCRIBED ELEVATIONS
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: cote(:)
!
!     PRESCRIBED VELOCITIES
!
      DOUBLE PRECISION, ALLOCATABLE ::  vites(:)
!
!     PRESCRIBED TRACERS VALUES
!
      DOUBLE PRECISION, ALLOCATABLE :: tracer(:)
!
!     FLUXES AT BOUNDARIES
!
      DOUBLE PRECISION, ALLOCATABLE, TARGET :: flux_boundaries(:)
!
!     ABSCISSAE, ORDINATES OF SOURCES
!
      DOUBLE PRECISION, ALLOCATABLE :: xsce(:),ysce(:)
!
!     ABSCISSAE, ORDINATES OF POINTS FOR SOURCES REGIONS
!
      DOUBLE PRECISION, ALLOCATABLE :: xcoo(:,:),ycoo(:,:)
!
!     AREA OF POINTS INSIDE SOURCES REGIONS
!
      DOUBLE PRECISION, ALLOCATABLE :: area_p(:)
!
!     SOURCES WATER DISCHARGES TAKEN FROM THE STEERING FILE
!
      DOUBLE PRECISION, ALLOCATABLE ::  dsce(:)
!
!     SOURCES DISCHARGES WITH VARIATION IN TIME
!
      DOUBLE PRECISION, ALLOCATABLE :: dsce2(:)
!
!     SOURCES TRACER DISCHARGES TAKEN FROM THE STERING FILE
!
      DOUBLE PRECISION, ALLOCATABLE :: tsce(:,:)
!
!     SOURCES OF TRACERS WITH VARIATIONS IN TIME
!
      DOUBLE PRECISION, ALLOCATABLE :: tsce2(:,:)
!
!     VELOCITIES OF THE SOURCES ALONG X,Y
!
      DOUBLE PRECISION, ALLOCATABLE :: usce(:),vsce(:)
!
!     UPWIND COEFFICIENTS FOR SUPG (1:u and v, 2:h, 3:tracers, 4:k and epsilon)
!
!     DOUBLE PRECISION COSUPG(4)
!
!     STOP CRITERIA (ORDER: H,U AND V,T)
!
      DOUBLE PRECISION criper(3)
!
!     REFERENCE TEMPERATURE FOR COMPUTING DENSITY
!
      DOUBLE PRECISION tmoy
!
!     DIRECTION OF NORTH, COUNTER-CLOCK-WISE, STARTING FROM VERTICAL AXIS
!
      DOUBLE PRECISION nord
!
!     NON-DIMENSIONAL DISPERSION COEFFICIENTS
!
      DOUBLE PRECISION elder(2)
!
!     LONGITUDE OF ORIGIN POINT
!
      DOUBLE PRECISION phi0
!
!     MASS-LUMPING ON TRACERS
!
      DOUBLE PRECISION agglot
!
!     DURATION OF COMPUTATION
!
      DOUBLE PRECISION duree
!
!     DESIRED COURANT NUMBER
!
      DOUBLE PRECISION cflwtd
!
!     MINIMUM DEPTH TO TAKE WIND INTO ACCOUNT
!
      DOUBLE PRECISION hwind
!
!     THRESHOLD FOR NEGATIVE DEPTHS
!
      DOUBLE PRECISION hneg
!
!     FREE SURFACE GRADIENT COMPATIBILITY
!
      DOUBLE PRECISION tetazcomp
!
!     TIME RANGE FOR FOURIER ANALYSIS
!
      DOUBLE PRECISION tafbgn,tafend
!
!     CURRENT TIME, INITIAL TIME
!
      DOUBLE PRECISION, TARGET ::  at,at0
!
!     FOURIER ANALYSIS PERIODS
!
      DOUBLE PRECISION periaf(50)
!
!     ARRAY OF REALS TO READ INTO SELAFIN FILES
!
      REAL, ALLOCATABLE :: w(:)
!
!     TOLERANCES FOR IDENTIFICATION
!
      DOUBLE PRECISION tolest(4)
!
!     NEWMARK TIME INTEGRATION COEFFICIENT
!
      DOUBLE PRECISION gamma
!
!     ARRAY WITH STAGE-DISCHARGE CURVES
!
      DOUBLE PRECISION, ALLOCATABLE :: qz(:,:,:)
!
!     RAIN OR EVAPORATION IN MM PER DAY
!     REPLACED BY CST_RAINFALL IN METEO_TELEMAC MODULE)
!
!     DOUBLE PRECISION RAIN_MMPD
!
!     PHYSICAL CHARACTERISTICS OF THE TSUNAMI
!
      DOUBLE PRECISION coetsunami(10)
!
!     VALUES OF TRACERS IN THE RAIN
!
      DOUBLE PRECISION, ALLOCATABLE :: train(:)
!
!     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
!
!     THRESHOLD DEPTH FOR RECEDING PROCEDURE
!
      DOUBLE PRECISION hrec
!
!     DIAMETRE OF ALGAE (PER CLASS)
!
      DOUBLE PRECISION, ALLOCATABLE :: dalgae(:)
!
!     DENSITY OF ALGAE (PER CLASS)
!
      DOUBLE PRECISION, ALLOCATABLE :: ralgae(:)
!
!     THICKNESS OF ALGAE (PER CLASS)
!
      DOUBLE PRECISION, ALLOCATABLE :: ealgae(:)
!
!     TIME OF RELEASE OF ALGAE (PER CLASS)
!
      DOUBLE PRECISION, ALLOCATABLE :: talgae(:)
!
!     ALGAE TYPES (PER CLASS)
!
      INTEGER         , ALLOCATABLE :: yalgae(:)
!
!     ALGAE RELEASE TYPES (PER CLASS)
!
      INTEGER         , ALLOCATABLE :: rel_algae(:)
!
!     WAVE SHEAR STRESS FOR ALGAE DISLODGEMENT 1
!
      DOUBLE PRECISION, ALLOCATABLE :: tw1_algae(:)
!
!     WAVE SHEAR STRESS FOR ALGAE DISLODGEMENT 2
!
      DOUBLE PRECISION, ALLOCATABLE :: tw2_algae(:)
!
!     RATE OF DEGRADATION CONSTANT FOR ALGAE
!
      DOUBLE PRECISION, ALLOCATABLE :: a_algae(:)
!
!     VALUE OF THE ATMOSPHERIC PRESSURE
!     REPLACED BY CST_PATMOS IN METEO_TELEMAC MODULE)
!
!     DOUBLE PRECISION PATMOS_VALUE
!
!     PRODUCTION COEFFICIENT FOR SECONDARY CURRENTS
!     DISSIPATION COEFFICIENT FOR SECONDARY CURRENTS
!
      DOUBLE PRECISION sec_as,sec_ds
!
!     VON KARMAN CONSTANT
!
      DOUBLE PRECISION karman
!
!     DURATION OF RAIN OR EVAPORATION IN HOURS
!
      DOUBLE PRECISION rain_hdur
!
!     TOLERANCE FOR FV SCHEME WETTING AND DRYING
!
      DOUBLE PRECISION eps_fv
!
!     ALLOCATABLE ARRAYS WITH DIMENSION MAXTRA, FOR TRACERS
!
!     MASSOU: MASS CREATED BY SOURCE TERM DURING THE TIME STEP
!     MASTR0: INITIAL MASS
!     MASTR2: CURRENT MASS
!     MASTEN: MASS ENTERED THROUGH BOUNDARIES
!     MASTOU: TOTAL MASS CREATED BY SOURCE TERM
!     MASTRAIN: TOTAL MASS BROUGHT BY RAIN OR EVAPORATION
!     FLUTENT AND FLUTSOR: FOR FINITE VOLUMES
!
      DOUBLE PRECISION, ALLOCATABLE :: massou(:)  ,mastr0(:)
      DOUBLE PRECISION, ALLOCATABLE :: mastr2(:)  ,flutsor(:)
      DOUBLE PRECISION, ALLOCATABLE :: masten(:)  ,mastou(:)
      DOUBLE PRECISION, ALLOCATABLE :: mastrain(:),flutent(:)
!
!     SOME PARAMETERS FOR SA TURBULENCE MODEL
!
      DOUBLE PRECISION :: sigmanu = 2.d0/3.d0
      DOUBLE PRECISION :: numin = 1.d-8
      DOUBLE PRECISION :: numax = 1.d10
!
!     CALIBRATION COEFFICIENTS FOR THE 2D MIXING LENGTH TURBULENCE MODEL
!
      DOUBLE PRECISION calmixlength(2)
!
!-----------------------------------------------------------------------
!
!       7) STRINGS
!
!-----------------------------------------------------------------------
!
!     TITLE OF STUDY
!
      CHARACTER(LEN=72) titcas
!
!     VARIABLES FOR GRAPHIC PRINTOUT
!
      CHARACTER(LEN=72) vardes
!
!     INITIAL CONDITIONS
!
      CHARACTER(LEN=72) cdtini
!
!     VARIABLES TO BE PRINTED
!
      CHARACTER(LEN=72) varimp
!
!     EQUATIONS
!
      CHARACTER(LEN=20), TARGET :: equa
!
!     NAMES OF CLANDESTINE VARIABLES
!
      CHARACTER(LEN=32) varcla(10)
!
!     NAMES AND UNITS OF VARIABLES
!
      CHARACTER(LEN=32) texte(maxvar),textpr(maxvar)
!
!     NAMES OF REMARKABLE POINTS
!
      CHARACTER(LEN=32) name_pts(100)
!
!     NAMES OF TRACERS
!
      CHARACTER(LEN=32), ALLOCATABLE :: nametrac(:)
!
!     NAMES OF PRIVATE ARRAYS (GIVEN BY USER)
!
      CHARACTER(LEN=32) names_prive(4)
!
!     NAMES OF DIFFERENTING ARRAYS (GIVEN BY USER)
!
      CHARACTER(LEN=32) names_advar(maxvar)
!
!-----------------------------------------------------------------------
!
!       8) SLVCFG STRUCTURES
!
!-----------------------------------------------------------------------
!
!     STRUCTURE WITH SOLVER FOR PROPAGATION
!
      TYPE(slvcfg) :: slvpro
!
!     STRUCTURE WITH SOLVER OPTIONS FOR K
!
      TYPE(slvcfg) :: slvk
!
!     STRUCTURE WITH SOLVER OPTIONS FOR E
!
      TYPE(slvcfg) :: slvep
!
!     STRUCTURE WITH SOLVER OPTIONS FOR SA
!
      TYPE(slvcfg) :: slvnu
!
!     SOLVER FOR DIFFUSION OF TRACERS
!
      TYPE(slvcfg), ALLOCATABLE :: slvtra(:)
!
!-----------------------------------------------------------------------
!
!       9) MESH STRUCTURE
!
!-----------------------------------------------------------------------
!
!     MESH STRUCTURE
!
      TYPE(bief_mesh), TARGET :: mesh
!
!     FINE MESH STRUCTURE IN CASE OF CONVERGENCE STUDY
!
      TYPE(bief_mesh), TARGET :: finemesh
!
!-----------------------------------------------------------------------
!
!      10) ALIASES
!
!-----------------------------------------------------------------------
!
!     DECLARATION OF POINTERS FOR ALIASES.
!     TARGETS ARE DEFINED IN POINT_TELEMAC2D.
!
!     ALIASES FOR WORKING VECTORS IN TB
!
      TYPE(bief_obj),POINTER :: t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11
      TYPE(bief_obj),POINTER :: t12,t13,t14,t15,t16,t17,t18,t19,t20
      TYPE(bief_obj),POINTER :: t21,t22
!
!     ALIASES FOR VECTOR IN BLOCK PRIVE
!
      DOUBLE PRECISION,POINTER :: prive1(:),prive2(:)
      DOUBLE PRECISION,POINTER :: prive3(:),prive4(:)
!
!     CONNECTIVITY TABLE
!
      TYPE(bief_obj), POINTER :: ikle
!
!     COORDINATES OF POINTS IN THE MESH
!
      DOUBLE PRECISION, DIMENSION(:), POINTER :: x,y
!
!     NUMBER OF ELEMENTS IN THE MESH
!
      INTEGER, POINTER:: nelem
!
!     MAXIMUM NUMBER OF ELEMENTS IN THE MESH
!
      INTEGER, POINTER:: nelmax
!
!     NUMBER OF BOUNDARY POINTS IN THE MESH
!
      INTEGER, POINTER:: nptfr
      INTEGER nptfr2
!
!     MAXIMUM NUMBER OF BOUNDARY POINTS IN THE MESH
!
      INTEGER, POINTER:: nptfrx
!
!     TYPE OF ELEMENT
!
      INTEGER, POINTER:: typelm
!
!     NUMBER OF POINTS IN THE MESH
!
      INTEGER, POINTER:: npoin
!
!     MAXIMUM NUMBER OF POINTS IN THE MESH
!
      INTEGER, POINTER:: npmax
!
!     MAXIMUM NUMBER OF POINTS NEIGHBOURS OF A POINT
!
      INTEGER, POINTER:: mxptvs
!
!     MAXIMUM NUMBER OF ELEMENTS NEIGHBOURS OF A POINT
!
      INTEGER, POINTER:: mxelvs
!
!     VECTOR LENGTH OF THE MACHINE
!
      INTEGER, POINTER:: lv
!
!     NUMBER OF SEGMENT IN A MESH
!
      INTEGER, POINTER:: nseg
!
!=======================================================================
!
!       ADDED DECLARATIONS FOR THE IDENTIFICATION OF PARAMETERS
!
!
!-----------------------------------------------------------------------
!
!       1) VECTORS
!
!-----------------------------------------------------------------------
!
!     ADJOINT VARIABLES
!
      TYPE(bief_obj), TARGET :: pp,qq,rr
!
!     MEASUREMENTS
!
      TYPE(bief_obj), TARGET :: ud,vd,hd
!
!     DIRECT VARIABLES AT ITERATION IT
!
      TYPE(bief_obj), TARGET :: uu,vv,hh
!
!     DIRECT VARIABLES AT ITERATION IT+1
!
      TYPE(bief_obj), TARGET :: uit1,vit1,hit1
!
!     BOUNDARY VALUES FOR ADJOINT VARIABLE (ONLY DIRICHLET)
!
      TYPE(bief_obj), TARGET :: pbor,qbor,rbor
!
!     VECTOR USED TO CHANGE THE SET OF STRICKLERS
!
      TYPE(bief_obj) :: desc
!
!     GRADIENT OF COST FUNCTION (ITERATION K)
!
      TYPE(bief_obj) :: gradj
!
!     GRADIENT OF COST FUNCTION (ITERATION K-1)
!
      TYPE(bief_obj) :: gradjn
!
!     SET OF STRICKLERS (ZONES): NEW
!
      TYPE(bief_obj) :: setstr
!
!     SET OF STRICKLERS' (ZONES): OLD
!
      type(bief_obj) :: setstr2
!
!     GRADIENTS
!
      TYPE(bief_obj) :: alpha1,alpha2,alpha3
!
!     VARIABLES TRANSMITTED FROM TOMAWAC TO SISYPHE
!
      TYPE(bief_obj), TARGET :: dirmoy,hm0,tpr5,orbvel,cosdir,sindir
!
!     VARIABLES TRANSMITTED FROM SISYPHE TO TELEMAC
!
      TYPE(api_cpl), TARGET :: sis_cpl
!
!     VARIABLES TRANSMITTED FROM GAIA TO TELEMAC
!
      TYPE(api_cpl), TARGET :: gai_cpl
!
!-----------------------------------------------------------------------
!
!       2) MATRICES
!
!-----------------------------------------------------------------------
!
!       MATRICES
!
!       SYSTEM SOLVED WILL BE :
!
!     ( TAM1  TCM1  TCM2 )   (PP )      ( CV1 )
!     (                  )   (   )      (     )
!     ( TBM1  TAM2  A23  )   (QQ )  =   ( CV2 )
!     (                  )   (   )      (     )
!     ( TBM2  A32   TAM3 )   (RR )      ( CV3 )
!
!
!     MATRICES FOR ADJOINT SYSTEM
!
      TYPE(bief_obj), TARGET :: tam1,tam2,tam3,tbm1,tbm2,tcm1,tcm2
!
!-----------------------------------------------------------------------
!
!       3) BLOCKS
!
!-----------------------------------------------------------------------
!
!     BLOCKS OF MATRICES, RIGHT HAND SIDES AND UNKNOWNS
!     FOR LINEAR SYSTEMS
!
!     BLOCKS OF MATRICES, RIGHT HAND SIDES FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: matadj
!
!     BLOCKS OF UNKNOWNS FOR LINEAR SYSTEMS
!
      TYPE(bief_obj), TARGET :: unkadj
!
!     BLOCK OF ADJOINT VARIABLES FOR OUTPUT
!
      TYPE(bief_obj), TARGET :: varsora
!
!     BLOCK OF DIRICHLET CONDITION FOR ADJOINT VARIABLES
!
      TYPE(bief_obj), TARGET :: adjdir
!
!-----------------------------------------------------------------------
!
!     4) INTEGERS
!
!     NUMBER OF VARIABLES FOR LITENR (READ BY SKIPGEO)
!
      INTEGER nvarres
!
!     NUMBER OF ITERATIONS FOR PARAMETERS CALIBRATION
!
      INTEGER nitera
!
!-----------------------------------------------------------------------
!
!     5) LOGICALS
!
!     IF YES: ACTIVATES THE ADJOINT MODE IN TELEMAC2D/PROPAG
!
      LOGICAL adjo
!
!     WHETHER THE ADJOINT SYSTEM RESULTS ARE PRINTED; NB 100=MAXVAR
!
      LOGICAL sorleoa(100)
!
!     WHETHER ...; NB 100=MAXVAR
!
      LOGICAL sorimpa(100)
!
!     FOR TRACERS IMPLICIT SOURCE TERMS
!
      LOGICAL, ALLOCATABLE :: yasmi(:)
!
!-----------------------------------------------------------------------
!
!       6) STRINGS
!
!-----------------------------------------------------------------------
!
!     NAMES OF VARIABLES IN RES FILE
!
      CHARACTER(LEN=32) texres(100)
!
!     NAMES OF VARIABLES IN RBI FILE
!
      CHARACTER(LEN=32) texrbi(100)
!
!     NAMES OF VARIABLES IN REF FILE
!
      CHARACTER(LEN=32) texref(100)
!
!     KEYWORD "PARAMETER ESTIMATION"
!
      CHARACTER(LEN=72) estime
!
!=======================================================================
!
!     COPY OF SUBMIT STRINGS IN THE DICTIONARY
!
      CHARACTER(LEN=PATH_LEN) submit(4,300)
!
!-----------------------------------------------------------------------
!
!       7) FRICTION DEFINITION BY ZONE
!
!-----------------------------------------------------------------------
!
!     ???????
!
      TYPE(friction_obj), TARGET :: frtab
!
!     ???????
!
      TYPE(bief_obj), TARGET :: kfropt, nkfrot
!
!     DEFAULT MANNING COEFFICIENT
!
      TYPE(bief_obj), TARGET :: ndefma
!
!     COEFFICIENTS FOR VEGETATION LAW AND VEGETATION LAW NUMBER
!
      TYPE(bief_obj), TARGET :: vcoeff, veglaw
!
!     ???????
!
      TYPE(bief_obj), TARGET :: ndef_b, kfro_b
!
!     NUMBER OF ZONES
!
      INTEGER:: nzones
!
!     MAXIMUM NUMBER OF FRICTION DOMAINS
!
      INTEGER:: nzonmx
!
!     IF YES, VEGETATION FRICTION
!
      LOGICAL:: vegetation
!
!     IF YES, FRICTION DATA
!
      LOGICAL:: frictb
!
!     DEFAULT MANNING VALUE (FOR COLEBROOK-WHITE LAW)
!
      DOUBLE PRECISION :: ndef
!
!-----------------------------------------------------------------------
!
!      8) TELEMAC-2D FILES
!
!-----------------------------------------------------------------------
!
!     TELEMAC-2D FILES
!
      TYPE(bief_file), TARGET :: t2d_files(maxlu_t2d)
!
!-----------------------------------------------------------------------
!
!      9) SECTIONS
!
!-----------------------------------------------------------------------
!
!     FOR CONTROL SECTIONS
!     ALLOCATED AND INITIALISED IN READ_SECTIONS_TELEMAC2D
!
      type(chain_type), ALLOCATABLE :: chain(:)
!
!
!-----------------------------------------------------------------------
!
!      10) WEIRS
!
!-----------------------------------------------------------------------
!
!     FOR WEIRS IF TYPE = 2
!     ALLOCATED AND INITIALISED IN LECSNG2
!
      type(weir_element)     , ALLOCATABLE :: weirs(:)
      type(weir_element_proc), ALLOCATABLE :: weirs_proc(:)
      type(weir_nodes)       , ALLOCATABLE :: wnodes(:)
      type(weir_nodes_proc)  , ALLOCATABLE :: wnodes_proc(:)
      type(weir_nodes)       , ALLOCATABLE :: wn_send(:)
      type(weir_nodes_proc)  , ALLOCATABLE :: wn_send_proc(:)
!
      INTEGER :: len_w_buf, len_wt_buf, len_wh_buf
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: w_buf_recv
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: w_buf_send
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: wt_buf_recv
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: wt_buf_send
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: wh_buf_recv
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: wh_buf_send
!
      INTEGER, DIMENSION(:), ALLOCATABLE :: weir_breach_start
      INTEGER, DIMENSION(:), ALLOCATABLE :: weir_breach_end
!
!     SAVED VARIABLES
!
      ! majtrac
      DOUBLE PRECISION, ALLOCATABLE :: dst_mt(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: dsp_mt(:),dsm_mt(:),corrt_mt(:)
      DOUBLE PRECISION,ALLOCATABLE  :: gradi_mt(:),gradj_mt(:)
      DOUBLE PRECISION,ALLOCATABLE  :: gradij_mt(:),gradji_mt(:)
      LOGICAL :: deja_mt = .false.
      ! read_fic_polygone and read_fic_points
      INTEGER nply
      INTEGER, PARAMETER :: max_ply = 21
      INTEGER, PARAMETER :: maxpoin_ply = 2100
      INTEGER            :: npoin_ply(max_ply)
      DOUBLE PRECISION   :: value_ply(maxpoin_ply)
      DOUBLE PRECISION   :: x_ply(maxpoin_ply),y_ply(maxpoin_ply)
      ! read_fic_source
      INTEGER nvalue_rfs,nlig_rfs,il1_rfs,il2_rfs
      INTEGER, PARAMETER :: maxvalue_rfs=2100
      CHARACTER(LEN=9) choix_rfs(maxvalue_rfs),lastwhat_rfs
      DOUBLE PRECISION tl1_rfs,tl2_rfs,lastat_rfs
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: infic_rfs
      DOUBLE PRECISION, DIMENSION(:)  , ALLOCATABLE :: time_rfs
      INTEGER src_ref_date(6)
      LOGICAL :: deja_rfs=.false.
      DOUBLE PRECISION :: src_tel_offset
      ! bilan.f
      DOUBLE PRECISION flux1_old,masse0,masse1,masse2,masent,masset
      DOUBLE PRECISION fluxn_old
      ! fluxpr_telemac2d
      DOUBLE PRECISION, ALLOCATABLE :: work_fpr(:)
      LOGICAL :: old_method_fpr=.false.
      LOGICAL :: init_fpr=.true.
      INTEGER :: nseo_fpr
      ! testeur
!     FORMER ARRAY !!!
      DOUBLE PRECISION, ALLOCATABLE :: flux_testeur(:)
      LOGICAL :: deja_testeur = .false.
      ! nodalf_schureman
      INTEGER year,nday,hour,minute,second
      ! flucin
!     AUTOMATIC EX ARRAYS!!!!!!
      DOUBLE PRECISION, ALLOCATABLE :: dsh_fc(:,:),dsu_fc(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: dsv_fc(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: dsp_fc(:),dsm_fc(:),dsz_fc(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: corr_fc(:),dtll_fc(:)
      DOUBLE PRECISION, ALLOCATABLE :: gradi_fc(:,:),gradj_fc(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: gradij_fc(:,:),gradji_fc(:,:)
      LOGICAL :: deja_fc = .false.
      ! nodalf_pugh
      INTEGER year_np,nday_np,hour_np,minute_np,second_np
      ! read_fic_frliq
      LOGICAL :: deja_rff = .false.
      DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: infic_rff
      DOUBLE PRECISION, DIMENSION(:)  , ALLOCATABLE :: time_rff
      INTEGER, PARAMETER :: maxvalue_rff = 2100
      CHARACTER(LEN=9) choix_rff(maxvalue_rff),lastwhat_rff
      INTEGER :: il1_rff,il2_rff,nvalue_rff,nlig_rff
      DOUBLE PRECISION :: tl1_rff,tl2_rff,lastat_rff
      INTEGER liq_ref_date(6)
      DOUBLE PRECISION :: liq_tel_offset
      ! spectre
      INTEGER :: afbgn,afend
      DOUBLE PRECISION, ALLOCATABLE :: am(:,:),bm(:,:),ha(:)
      LOGICAL :: deja = .false.
      ! bord_tidal_bc
      LOGICAL :: deja_tbc = .false.
      INTEGER, ALLOCATABLE :: firsttide(:),lasttide(:)
      ! flusec_telemac2d
      LOGICAL :: deja_fs = .false.
      LOGICAL :: old_method_fs=.false.
      DOUBLE PRECISION, ALLOCATABLE :: vlx_fs(:),volneg_fs(:)
      DOUBLE PRECISION, ALLOCATABLE :: volpos_fs(:)
      INTEGER, ALLOCATABLE :: nseg_fs(:),liste_fs(:,:,:)
      ! bilant1
      DOUBLE PRECISION mastr0_bt1(100),mastr1_bt1(100)
      DOUBLE PRECISION mastr2_bt1(100),masten_bt1(100)
      DOUBLE PRECISION mastou_bt1(100),dirtot_bt1(100)
      ! preres_telemac2d
      LOGICAL :: deja1_preres = .false.
      LOGICAL :: deja2_preres = .false.
      LOGICAL :: deja3_preres = .false.
      ! breach
      LOGICAL :: dejalu_breach = .false.
      ! thomps
      LOGICAL :: init_thomps = .true.
      TYPE(bief_obj) :: fncar1,ftild1
      ! bord_tide
      LOGICAL :: deja_bt = .false.
      INTEGER nptfrl,ntide
      INTEGER, ALLOCATABLE :: firsttide_bt(:),lasttide_bt(:)
      INTEGER, ALLOCATABLE :: shifttide_bt(:)
      DOUBLE PRECISION, ALLOCATABLE :: ahm2(:),phm2(:)
      DOUBLE PRECISION, ALLOCATABLE :: ahs2(:),phs2(:)
      DOUBLE PRECISION, ALLOCATABLE :: ahn2(:),phn2(:)
      DOUBLE PRECISION, ALLOCATABLE :: ahm4(:),phm4(:)
      DOUBLE PRECISION, ALLOCATABLE :: aum2(:),pum2(:)
      DOUBLE PRECISION, ALLOCATABLE :: aus2(:),pus2(:)
      DOUBLE PRECISION, ALLOCATABLE :: aun2(:),pun2(:)
      DOUBLE PRECISION, ALLOCATABLE :: aum4(:),pum4(:)
      DOUBLE PRECISION, ALLOCATABLE :: avm2(:),pvm2(:)
      DOUBLE PRECISION, ALLOCATABLE :: avs2(:),pvs2(:)
      DOUBLE PRECISION, ALLOCATABLE :: avn2(:),pvn2(:)
      DOUBLE PRECISION, ALLOCATABLE :: avm4(:),pvm4(:)
      DOUBLE PRECISION, ALLOCATABLE :: fhxm2(:),fhym2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fhxs2(:),fhys2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fhxn2(:),fhyn2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fhxm4(:),fhym4(:)
      DOUBLE PRECISION, ALLOCATABLE :: fuxm2(:),fuym2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fuxs2(:),fuys2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fuxn2(:),fuyn2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fuxm4(:),fuym4(:)
      DOUBLE PRECISION, ALLOCATABLE :: fvxm2(:),fvym2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fvxs2(:),fvys2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fvxn2(:),fvyn2(:)
      DOUBLE PRECISION, ALLOCATABLE :: fvxm4(:),fvym4(:)
      DOUBLE PRECISION upvm2,upvn2,upvs2,ffmn2,ffm4
      DOUBLE PRECISION fm2,fs2,fn2,fm4
      ! bord_tide_misc
      INTEGER, PARAMETER   :: misc_ncmx_btm = 61
      CHARACTER(LEN=4) :: misc_constid_btm(misc_ncmx_btm) =
     &                 (/'2MK6','2MN6','2MS6','2N2 ','2Q1 ',
     &                   '2SM2','2SM6','E2  ','Tta1','J1  ',
     &                   'K1  ','K2  ','KJ2 ','KQ1 ','L2  ',
     &                   'La2 ','M1  ','M2  ','M3  ','M4  ',
     &                   'M6  ','M8  ','Mf  ','MK3 ','MK4 ',
     &                   'MKS2','Mm  ','MN4 ','MO3 ','MP1 ',
     &                   'MS4 ','MSf ','MSK6','MSN2','MSN6',
     &                   'MSqm','Mtm ','Mu2 ','N2  ','N4  ',
     &                   'Nu2 ','O1  ','OO1 ','P1  ','Pi1 ',
     &                   'Q1  ','R2  ','Ro1 ','S1  ','S2  ',
     &                   'S4  ','Sa  ','Sig1','SK4 ','SN4 ',
     &                   'Ssa ','T2  ','Phi1','Ki1 ','Psi1',
     &                   'Z0  '/)
      LOGICAL :: deja_btm = .false.
      INTEGER :: ntide_btm, nptfrl_btm,nwaves_btm,nwaves2_btm
      INTEGER, ALLOCATABLE :: firsttide_btm(:),lasttide_btm(:)
      INTEGER, ALLOCATABLE :: shifttide_btm(:)
      INTEGER, ALLOCATABLE :: indw_btm(:),indw2_btm(:),indw3_btm(:)
      CHARACTER(LEN=4), ALLOCATABLE :: namewave_btm(:)
      DOUBLE PRECISION, ALLOCATABLE :: ah_btm(:,:),ph_btm(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: au_btm(:,:),pu_btm(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: av_btm(:,:),pv_btm(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: lon_btm(:,:),lat_btm(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: upv_btm(:),ff_btm(:)
      DOUBLE PRECISION, ALLOCATABLE :: omega_btm(:),phcalhw_btm(:)
      ! flusec_t2d.f
      TYPE fluxline
        INTEGER, ALLOCATABLE :: sectionids(:)
        INTEGER, ALLOCATABLE :: direction(:)
        INTEGER              :: nofsections
      END TYPE fluxline
      LOGICAL :: deja_flusect2d = .false.
      INTEGER :: numberoflines_flusect2d
      DOUBLE PRECISION :: time_flusect2d
      TYPE(fluxline), ALLOCATABLE :: fluxlinedata_flusect2d(:)
      DOUBLE PRECISION, ALLOCATABLE :: flx_flusect2d(:,:)
      DOUBLE PRECISION, ALLOCATABLE :: volflux_flusect2d(:,:)
      ! collect_values.f
      INTEGER :: msg_tag = 5000
      ! collect_values_trac.f
      INTEGER :: msg_tag_tra = 5000
      ! collect_values_wave.f
      INTEGER :: msg_tag_w = 5000
      ! prosou

      LOGICAL :: prosou_dejalu = .false.
      LOGICAL :: prosou_fv_dejalu = .false.
!
      ! telemac2d.F
      DOUBLE PRECISION kmax
      DOUBLE PRECISION emax,sigmae,sigmak,c1
      ! YAFLULIM: If yes, flulim will be applied to segment fluxes when
      ! calling cvdftr
      ! YAFLODEL If yes, flodel has been computed
      LOGICAL :: yaflulim, yaflulimebe, yaflodel
      INTEGER leoprd_charr
      INTEGER ref_t2d
      INTEGER trouve(maxvar+10)

      DOUBLE PRECISION kmin,flusor,mass_rain,fluent
      DOUBLE PRECISION masses,relaxs,relaxb
      DOUBLE PRECISION emin,estar,schmit,c2,cmu

      LOGICAL akep,sa
      LOGICAL imp
      LOGICAL,TARGET::leo,leo_trac

      CHARACTER(LEN=24), PARAMETER :: code1='TELEMAC2D               '

      INTEGER ioptan_t2d

      INTEGER jjjj
      INTEGER :: alire(maxvar) = (/
     &            1,1,0,1,1,1,0,0,0,1,1,1,(0,jjjj=13,maxvar) /)
!     ADVECTION FIELD USED FOR SISYPHE CALL
!
      TYPE(bief_obj), POINTER :: usis,vsis
!     Number of iterations asked in the case file for api
      INTEGER :: total_iter
!
      TYPE(bief_obj), TARGET :: nes_dz
!
      DOUBLE PRECISION, ALLOCATABLE :: avail(:,:,:)
!
      ! Criteria to stop the computation
      LOGICAL :: arret1, arret2, arret3
!
      SAVE
!
!     ============================================
!
      END MODULE declarations_telemac2d