point_telemac3d.f

Go to the documentation of this file.

!                   **************************
                    SUBROUTINE point_telemac3d
!                   **************************
!
!
!***********************************************************************
! TELEMAC3D   V8P2
!***********************************************************************
!
!brief    ALLOCATES TELEMAC3D STRUCTURES.
!
!history  JACEK A. JANKOWSKI PINXIT
!+        **/03/1999
!+
!+   FORTRAN 95 VERSION
!
!history  J-M HERVOUET (LNHE)
!+        19/10/2009
!+        V6P0
!+
!
!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  J-M HERVOUET (LNHE)
!+        02/04/2012
!+        V6P2
!+   New variables for graphic printouts added, for a clean restart.
!
!history  J-M HERVOUET (EDF R&D, LNHE)
!+        18/06/2013
!+        V6P3
!+   Size of IT1,2,3,4 modified in case of weak characteristics.
!
!history  C. VILLARET & T. BENSON & D. KELLY (HR-WALLINGFORD)
!+        27/02/2014
!+        V7P0
!+   New developments in sediment merged on 25/02/2014.
!
!history  J-M HERVOUET (EDF R&D, LNHE)
!+        14/05/2014
!+        V7P0
!+   ZCONV and ZFLATS now declared as type 15, discontinuous linear.
!
!history  G. ANTOINE (EDF LAB, LNHE)
!+        19/09/2014
!+        V7P0
!+   Adding variables for mixed sediment, see IF(S3D_MIXTE).
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        28/07/2015
!+        V7P1
!+   Adding aliases and 2D private arrays.
!
!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
!+   Adding ZN.
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        22/08/2016
!+        V7P2
!+   Adding TB2, a block of work arrays for the LIPS advection scheme.
!
!history  C.-T. PHAM (EDF, LNHE)
!+        01/03/2017
!+        V7P2
!+   Allowing k-epsilon model on a direction and not on the other.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE declarations_telemac3d
      USE meteo_telemac !, ONLY: WINDX,WINDY,PATMOS,TAIR,POINT_METEO
      USE declarations_waqtel, ONLY: waqprocess
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER CFG(2),CFG2D(2),CFGMURD(2),CFGBOR2D(2),CFGMURD_TF(2)
      INTEGER ITRAC,ITAB,IELM,IELV,IELH,STATUT,NTR,I,NSEG,SIZ,K
      LOGICAL YESWEAK,YESLIPS
      CHARACTER(LEN=1) TYPDIA, TYPEXT
!
!-----------------------------------------------------------------------
!
      IF(listin) THEN
        WRITE(lu,21)
      ENDIF
 21   FORMAT(1x,/,1x,'POINT_TELEMAC3D: MEMORY ALLOCATION',/)
!
!-----------------------------------------------------------------------
!
!     WILL THERE BE WEAK CHARACTERISTICS?
!
      yesweak=.false.
      IF(optcha.GT.1) THEN
        IF(schcvi.EQ.adv_car) yesweak=.true.
        IF(schcke.EQ.adv_car) yesweak=.true.
        DO itrac=1,ntrac
          IF(schcta(itrac).EQ.adv_car) yesweak=.true.
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
!     WILL THERE BE A LIPS SCHEME?
!
      yeslips=.false.
!     LOOKING AT VELOCITIES
      IF((schcvi.EQ.adv_nsc.OR.schcvi.EQ.adv_psi).AND.
     &   optadv_vi.EQ.4) THEN
        yeslips=.true.
      ENDIF
!     LOOKING AT K-EPSILON
      IF((schcke.EQ.adv_nsc.OR.schcke.EQ.adv_psi).AND.
     &   optadv_ke.EQ.4) THEN
        yeslips=.true.
      ENDIF
!     LOOKING AT TRACERS
      IF(ntrac.GT.0) THEN
        DO itrac=1,ntrac
          IF((schcta(itrac).EQ.adv_nsc.OR.schcta(itrac).EQ.adv_psi).AND.
     &        optadv_tr(itrac).EQ.4) THEN
            yeslips=.true.
          ENDIF
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
! DECLARES DISCRETISATION TYPES HERE
!
      ielm0 = 10*(ielmh/10) ! FOR TELEMAC2D
      ielm1 = ielbor(ielmh,1) ! BOUNDARY DISCRET. FOR TELEMAC2D
!
! IELMU IS 12, WHEN QUASI-BUBBLE FREE SURFACE REQUIRED (IN LECDON)
! IF NOT, IELMU=IELMH
!
      ielmx=max(ielmu,ielm2h,ielmh) ! IT WILL BE MAX. DISCR. IN 2D
!
! STORAGE TYPE AND MATRIX-VECTOR PRODUCT TYPE
!
      cfg(1) = optass
      cfg(2) = produc   ! PRODUC=1 HARD IN LECDON
      cfg2d(1) = optass2d
      cfg2d(2) = produc   ! PRODUC=1 HARD IN LECDON
!     MURD MATRIX WITH EDGE-BASED STORAGE FOR TIDAL FLATS
      cfgmurd_tf(1)=3
      cfgmurd_tf(2)=1
!     NORMAL MURD MATRIX WITH EBE STORAGE BECAUSE (SEE MURD3D.F) PSI
!     SCHEME IS DONE IN EBE STORAGE
      cfgmurd(1)=1
      cfgmurd(2)=1
!
      cfgbor2d(1)=1
      cfgbor2d(2)=1
!
!=======================================================================
!
!                     *********************
!                     *  MESH - GEOMETRY  *
!                     *********************
!
! TWO MESHES ARE ALLOCATED: (1) 2D BASE MESH, (2) 3D SIGMA-MESH
!
! ALLOCATES THE 2D MESH STRUCTURE FOR TELEMAC2D
! DISCRETISATION IELMH GIVEN IN LECDON
! IELMX = IELMU IF QUASI-BUBBLE ELEMENT REQUIRED, OTHERWISE IELMH
!
      equa = 'NO_EQUATION_IS_GIVEN'
!
      CALL almesh(mesh2d,'MESH2D',ielmx,spheri,cfg2d,
     &            t3d_files(t3dgeo)%FMT,
     &            t3d_files(t3dgeo)%LU,equa,0,nplan=1,
     &            projection=protyp,lati0=latit,longi0=longit)
      nseg=mesh2d%NSEG
!
! ALIASES FOR CERTAIN COMPONENTS OF THE 2D MESH STRUCTURE
!
      x2      => mesh2d%X
      y2      => mesh2d%Y
!
      surfa2  => mesh2d%SURFAC
      xnebor2 => mesh2d%XNEBOR
      ynebor2 => mesh2d%YNEBOR
      xsgbor2 => mesh2d%XSGBOR
      ysgbor2 => mesh2d%YSGBOR
      ikle2   => mesh2d%IKLE
      nbor2   => mesh2d%NBOR   ! PREVIOUSLY SIMPLY NBOR
!
      nelem2  => mesh2d%NELEM
      nelmax2 => mesh2d%NELMAX  ! PREVIOUSLY NELMA2 (ADAPTIVITY OUTLOOK)
      nptfr2  => mesh2d%NPTFR   ! PREVIOUSLY SIMPLY NPTFR
      nptfrx2 => mesh2d%NPTFRX
      dim2    => mesh2d%DIM1
      typelm2 => mesh2d%TYPELM
      npoin2  => mesh2d%NPOIN
      npmax2  => mesh2d%NPMAX
      mxptvs2 => mesh2d%MXPTVS
      mxelvs2 => mesh2d%MXELVS
      lv      => mesh2d%LV      ! MESH-CHECKED? 2D=3D FOR SIGMA MESH
!
!-----------------------------------------------------------------------
! ALLOCATES THE 3D MESH STRUCTURE (EQUA=EMPTY) (READ AGAIN?)
!
      equa = 'NO_EQUATION_IS_GIVEN'
!
      CALL almesh(mesh3d,'MESH3D',ielm3,spheri,cfg,
     &            t3d_files(t3dgeo)%FMT,
     &            t3d_files(t3dgeo)%LU,equa,0,nplan=nplan,
     &            projection=protyp,lati0=latit,longi0=longit)
!
! ALIAS FOR CERTAIN COMPONENTS OF THE 3D MESH STRUCTURE
! THEY ARE DEFINED IN DECLARATIONS
!
      x       => mesh3d%X%R
      y       => mesh3d%Y%R
      z       => mesh3d%Z%R
      x3      => mesh3d%X    ! POINTERS
      y3      => mesh3d%Y
      z3      => mesh3d%Z
      surfa3  => mesh3d%SURFAC
      xsgbor3 => mesh3d%XSGBOR
      ysgbor3 => mesh3d%YSGBOR
      zsgbor3 => mesh3d%ZSGBOR
      ikle3   => mesh3d%IKLE
      nbor3   => mesh3d%NBOR
!
      w1      => mesh3d%W
!
      nelem3  => mesh3d%NELEM
      nelmax3 => mesh3d%NELMAX   ! PREVIOUSLY NELMA3 (ADAPTIVITY?)
      neleb   => mesh3d%NELEB
      nelebx  => mesh3d%NELEBX
      nptfr3  => mesh3d%NPTFR
      nptfrx3 => mesh3d%NPTFRX
      dim3    => mesh3d%DIM1
      typelm3 => mesh3d%TYPELM
      npoin3  => mesh3d%NPOIN
      npmax3  => mesh3d%NPMAX
      mxptvs3 => mesh3d%MXPTVS
      mxelvs3 => mesh3d%MXELVS
!
!-----------------------------------------------------------------------
! VARIOUS MESH PARAMETER FIX
!
      netage = nplan - 1
!
!     NUMBER OF ADVECTED VARIABLES (3 FOR VELOCITY, 2 FOR K-EPSILON
!                                   THEN TRACERS)
      nvbil  = 5 + ntrac
!
      IF (listin) THEN
        WRITE(lu,32)
     &             typelm2,npoin2,nelem2,nptfr2,typelm3,npoin3,nelem3,
     &             nplan,neleb,nptfr3+2*npoin2,nptfr3,npoin2,npoin2
      ENDIF
!
 32   FORMAT(/,' 2D MESH',/,
     &         ' -------',//,
     &         ' 2D ELEMENT TYPE                : ',i8,/,
     &         ' NUMBER OF 2D NODES             : ',i8,/,
     &         ' NUMBER OF 2D ELEMENTS          : ',i8,/,
     &         ' NUMBER OF 2D BOUNDARY NODES    : ',i8,///,
     &         ' 3D MESH',/,
     &         ' -------',//,
     &         ' 3D ELEMENT TYPE                : ',i8,/,
     &         ' NUMBER OF 3D NODES             : ',i8,/,
     &         ' NUMBER OF 3D ELEMENTS          : ',i8,/,
     &         ' NUMBER OF LEVELS               : ',i8,/,
     &         ' NUMBER OF BOUNDARY ELEMENTS    : ',i8,/,
     &         ' TOTAL NUMBER OF BOUNDARY NODES : ',i8,/,
     &         ' INCLUDING   LATERAL BOUNDARIES : ',i8,/,
     &         '                        SURFACE : ',i8,/,
     &         '                         BOTTOM : ',i8,/)
!
!-----------------------------------------------------------------------
! FIELDS CONNECTED STRONGLY WITH GEOMETRY AND MESH
! REAL
!
      CALL bief_allvec(1, zprop,  'ZPROP ', ielm3 , 1, 1,mesh3d)
!     USED ONLY WITH SUPG ADVECTION
      CALL bief_allvec(1, zt   ,  'ZT    ', ielm3 , 1, 1,mesh3d)
!     USED ONLY WITH DISTRIBUTIVE SCHEMES WITH UPWIND OF DERIVATIVE
!     IN TIME
      CALL bief_allvec(1, zn   ,  'ZN    ', ielm3 , 1, 1,mesh3d)
!
      CALL allblo(gradzf, 'GRADZF')
!     23/11/2010 NOW GRADZF CONTAINS GRADIENTS OF ALL PLANES
      CALL bief_allvec_in_block(gradzf,2,1,'GRAZF ',ielm3 ,1,1,mesh3d)
!
      CALL allblo(gradzs, 'GRADZS')
      CALL bief_allvec_in_block(gradzs,2,1,'GRAZS ',ielm2h,1,1,mesh2d)
      CALL allblo(gradzn, 'GRADZN')
      CALL bief_allvec_in_block(gradzn,2,1,'GRAZN ',ielm2h,1,1,mesh2d)
!
      CALL bief_allvec(1, dssudt, 'DSSUDT', ielm2h, 1, 1,mesh2d)
!
!     DESCRIBES THE MESH ON THE VERTICAL
!
      CALL bief_allvec(1, zstar       , 'ZSTAR ' , nplan ,1,0,mesh3d)
      CALL bief_allvec(1, zplane      , 'ZPLANE' , nplan ,1,0,mesh3d)
      CALL bief_allvec(1, zchar       , 'ZCHAR ' , nplan ,1,0,mesh3d)
      CALL bief_allvec(2, transf_plane, 'TRAPLA' , nplan ,1,0,mesh3d)
!
! MESH AND GEOMETRY, INTEGERS
!
      CALL bief_allvec(2,lihbor, 'LIHBOR',         ielm1,1,1,mesh2d)
      CALL bief_allvec(2,numliq, 'NUMLIQ',        ielm2v,1,1,mesh3d)
      CALL bief_allvec(2,numliq_elm, 'NLQELM',10*(ielm2v/10),1,1,mesh3d)
      CALL bief_allvec(2,boundary_colour ,'BNDCOL',ielm1,1,1,mesh2d)
      CALL bief_allvec(2,limpro, 'LIMPRO',         ielm1,6,1,mesh2d)
!
! OPEN BED BOUNDARIES VARIABLES
!
      IF(bedbou)THEN
        CALL bief_allvec(2,nliqbed, 'NLIQBE', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1,bedflu , 'BEDFLU', ielm2h, 1, 1,mesh2d)
        DO i=1,nbedflo
          nliqbed%I(i)=0
        ENDDO
      ELSE
        CALL bief_allvec(2,nliqbed, 'NLIQBE', 0     , 1, 0,mesh2d)
        CALL bief_allvec(1,bedflu , 'BEDFLU', 0     , 1, 0,mesh2d)
      ENDIF
!
!-----------------------------------------------------------------------
!
!  BLOCK OF BOUNDARY CONDITIONS MASKS FOR PROPAGATION
!
!  BLOCK OF 9 VECTORS FOR MASKING (2D AND 3D)
!
      CALL allblo(mask  , 'MASK  ')
      CALL bief_allvec_in_block(mask,9,1,'MASK  ',ielm1,1,1,mesh2d)
      CALL allblo(mask_3d, 'MASK3D')
      CALL bief_allvec_in_block(mask_3d,9,1,'MSK3D ',
     &                          10*(ielm2v/10),1,1,mesh3d)
!
! MASKING
!
      IF (msk.OR.(optban.EQ.1.AND.nonhyd)) THEN
        CALL bief_allvec(1, maskel, 'MASKEL',10*(ielm3/10),1,1,mesh3d)
        CALL bief_allvec(1, maskpt, 'MASKPT', ielm3,       1,1,mesh3d)
      ELSE
        CALL bief_allvec(1, maskel, 'MASKEL', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, maskpt, 'MASKPT', 0, 1, 0,mesh3d)
      ENDIF
!
!    MASKS FOR LATERAL BOUNDARY ELEMENTS...
!
!     USED WHEN ELEMENTS ARE MASKED
      CALL bief_allvec(1, maskbr, 'MASKBR',10*(ielm2v/10), 1, 1,mesh3d)
!     USED AS WORK ARRAY IN DEBIMP_3D
      CALL bief_allvec(1, masktr, 'MASKTR',10*(ielm2v/10), 1, 1,mesh3d)
!
!=======================================================================
!                     ********************
!                     *    VARIABLES     *
!                     ********************
!
!-----------------------------------------------------------------------
! HORIZONTAL VELOCITY (U,V)
!
! U AND V VELOCITY COMPONENTS
!
      CALL bief_allvec(1, un,       'UN    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, vn,       'VN    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, uc,       'UC    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, vc,       'VC    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, ud,       'UD    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, vd,       'VD    ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, u,        'U     ',ielm3,  1,2,mesh3d)
      CALL bief_allvec(1, v,        'V     ',ielm3,  1,2,mesh3d)
      CALL bief_allvec(1, s0u,      'S0U   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, s0v,      'S0V   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, s1u,      'S1U   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, s1v,      'S1V   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, smu,      'SMU   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, smv,      'SMV   ',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, uborf,    'UBORF ',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, vborf,    'VBORF ',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, uborl,    'UBORL ',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, vborl,    'VBORL ',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, ubors,    'UBORS ',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, vbors,    'VBORS ',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, auborf,   'AUBORF',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, avborf,   'AVBORF',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, auborl,   'AUBORL',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, uetcal,   'UETCAL',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, avborl,   'AVBORL',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, aubors,   'AUBORS',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, avbors,   'AVBORS',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, buborf,   'BUBORF',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, buborl,   'BUBORL',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, bubors,   'BUBORS',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, bvborf,   'BVBORF',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, bvborl,   'BVBORL',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, bvbors,   'BVBORS',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, uborsave, 'UBSAVE',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, vborsave, 'VBSAVE',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, wborsave, 'WBSAVE',ielm2v, 1,1,mesh3d)
      CALL bief_allvec(1, wdist  ,  'WDIST ',ielm2h, 1,1,mesh2d)
      CALL bief_allvec(1, strain,   'STRAIN',ielm3,  1,1,mesh3d)
      CALL bief_allvec(1, rotan ,   'ROTAN ',ielm3,  1,1,mesh3d)
!
      IF(tidaltype.EQ.0) THEN
        CALL bief_allvec(1,hbtide ,'HBTIDE',0,1,0,mesh2d)
        CALL bief_allvec(1,ubtide ,'UBTIDE',0,1,0,mesh2d)
        CALL bief_allvec(1,vbtide ,'VBTIDE',0,1,0,mesh2d)
        CALL bief_allvec(2,numtide,'NUMTID',0,1,0,mesh2d)
      ELSE
        CALL bief_allvec(1,hbtide ,'HBTIDE',ielm1 ,1,1,mesh2d)
        CALL bief_allvec(1,ubtide ,'UBTIDE',ielm1 ,1,1,mesh2d)
        CALL bief_allvec(1,vbtide ,'VBTIDE',ielm1 ,1,1,mesh2d)
        CALL bief_allvec(2,numtide,'NUMTID',ielm1 ,1,1,mesh2d)
      ENDIF
!
! FRICTION VELOCITY **2 ON THE BOTTOM
!
      CALL bief_allvec(1, uetcar, 'UETCAR', ielm2h, 1,1,mesh2d)
!
! PLANE ON THE BOTTOM (THE FIRST ONE WITH A REAL ELEMENT HEIGHT ABOVE)
!
      CALL bief_allvec(2, ipbot,  'IPBOT ', ielm2h, 1,1,mesh2d)
!
! BOUNDARY CONDITION TYPES/
! ATTRIBUTES FOR THE HORIZONTAL VELOCITY COMPONENTS
!
      CALL bief_allvec(2, liubof,     'LIUBOF', ielm2h, 1,1,mesh2d)
! LIUBOL, LIVBOL, LIWBOL, LIEBOL, LIKBOL, LITABL:
! VALUES SAVED IN SECOND DIMENSION BECAUSE ADVECTION SCHEMES (CVDF3D)
! MAY CHANGE THE VALUES
      CALL bief_allvec(2, liubol,     'LIUBOL', ielm2v, 2,1,mesh3d)
      CALL bief_allvec(2, liubos,     'LIUBOS', ielm2h, 1,1,mesh2d)
!
      CALL bief_allvec(2, livbof,     'LIVBOF', ielm2h, 1,1,mesh2d)
      CALL bief_allvec(2, livbol,     'LIVBOL', ielm2v, 2,1,mesh3d)
      CALL bief_allvec(2, livbos,     'LIVBOS', ielm2h, 1,1,mesh2d)
!
! BOUNDARY CONDITION SAVED WHEN THOMPSON CONDITIONS USED
!
      IF(thomfr) THEN
        CALL bief_allvec(2,lihbor_user, 'HB_USE',ielm1 ,1,1,mesh2d)
        CALL bief_allvec(2,liubol_user, 'UB_USE',ielm2v,1,1,mesh3d)
        CALL bief_allvec(2,livbol_user, 'UV_USE',ielm2v,1,1,mesh3d)
      ENDIF
!
!-----------------------------------------------------------------------
! W VERTICAL VELOCITY COMPONENT (REAL VALUES)
! TREATED DIFFERENTLY DEPENDING ON THE (NON)HYDROSTATIC OPTION
!
      IF (nonhyd) THEN
        CALL bief_allvec(1, wn,     'WN    ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, wc,     'WC    ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, wd,     'WD    ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, wconv , 'WCONV ', ielm3, 1, 1,mesh3d)
      ELSE ! DUMMIES
        CALL bief_allvec(1, wn,     'WN    ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, wc,     'WC    ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, wd,     'WD    ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, wconv , 'WCONV ', 0, 1, 0,mesh3d)
      ENDIF
      CALL bief_allvec(1, w,      'W     ', ielm3,  1, 1,mesh3d)
      IF (nonhyd) THEN
        CALL bief_allvec(1, s0w,    'S0W   ', ielm3,  1, 1,mesh3d)
        CALL bief_allvec(1, s1w,    'S1W   ', ielm3,  1, 1,mesh3d)
      ELSE ! DUMMIES
        CALL bief_allvec(1, s0w,    'S0W   ', 0,  1, 0,mesh3d)
        CALL bief_allvec(1, s1w,    'S1W   ', 0,  1, 0,mesh3d)
      ENDIF
      CALL bief_allvec(1, wborf,  'WBORF ', ielm2h, 1, 1,mesh2d)
      CALL bief_allvec(1, wborl,  'WBORL ', ielm2v, 1, 1,mesh3d)
      CALL bief_allvec(1, wbors,  'WBORS ', ielm2h, 1, 1,mesh2d)
      IF (nonhyd) THEN
        CALL bief_allvec(1, awborf, 'AWBORF', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, awborl, 'AWBORL', ielm2v, 1, 1,mesh3d)
        CALL bief_allvec(1, awbors, 'AWBORS', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, bwborf, 'BWBORF', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, bwborl, 'BWBORL', ielm2v, 1, 1,mesh3d)
        CALL bief_allvec(1, bwbors, 'BWBORS', ielm2h, 1, 1,mesh2d)
      ELSE ! DUMMIES
        CALL bief_allvec(1, awborf, 'AWBORF', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, awborl, 'AWBORL', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, awbors, 'AWBORS', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, bwborf, 'BWBORF', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, bwborl, 'BWBORL', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, bwbors, 'BWBORS', 0, 1, 0,mesh2d)
      ENDIF
!
      CALL bief_allvec(2, liwbof,    'LIWBOF', ielm2h, 1, 1,mesh2d)
      CALL bief_allvec(2, liwbol,    'LIWBOL', ielm2v, 2, 1,mesh3d)
      CALL bief_allvec(2, liwbos,    'LIWBOS', ielm2h, 1, 1,mesh2d)
!
! SIGMA-TRANSFORMED VALUES / NO BLOCK
!
      CALL bief_allvec(1,ws,'WS    ',ielm3,1,1,mesh3d)
!
! VECTOR OF SIZE NUMBER OF POINTS IN THE ELEMENT * NELEM
!
      CALL bief_allvec(1,wel,'WEL   ',10*(ielm3/10),
     &                 bief_nbpel(ielm3,mesh3d),1,mesh3d)
!
!=======================================================================
! HYDRODYNAMIC PRESSURE SPECIFIC FOR THE NON-HYDROSTATIC OPTION
!
      IF (nonhyd) THEN
        CALL bief_allvec(1, dp,     'DP    ', ielm3,  1, 2,mesh3d)
        CALL bief_allvec(1, dpn,    'DPN   ', ielm3,  1, 1,mesh3d)
        CALL bief_allvec(1, ph,     'PH    ', ielm3,  1, 1,mesh3d)
        CALL bief_allvec(1, pborf,  'PBORF ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, pborl,  'PBORL ', ielm2v, 1, 1,mesh3d)
        CALL bief_allvec(1, pbors,  'PBORS ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(2, lipbof, 'LIPBOF', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(2, lipbol, 'LIPBOL', ielm2v, 1, 1,mesh3d)
        CALL bief_allvec(2, lipbos, 'LIPBOS', ielm2h, 1, 1,mesh2d)
      ELSE
        CALL bief_allvec(1, dp,     'DP    ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, dpn,    'DPN   ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, ph,     'PH    ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, pborf,  'PBORF ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, pborl,  'PBORL ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, pbors,  'PBORS ', 0, 1, 0,mesh2d)
        CALL bief_allvec(2, lipbof, 'LIPBOF', 0, 1, 0,mesh2d)
        CALL bief_allvec(2, lipbol, 'LIPBOL', 0, 1, 0,mesh3d)
        CALL bief_allvec(2, lipbos, 'LIPBOS', 0, 1, 0,mesh2d)
      ENDIF
!
!=======================================================================
!
      CALL allblo(trn      ,'TRN   ')
      CALL allblo(tac      ,'TAC   ')
      CALL allblo(ta       ,'TA    ')
      CALL allblo(s0ta     ,'S0TA  ')
      CALL allblo(s1ta     ,'S1TA  ')
      CALL allblo(taborf   ,'TABORF')
      CALL allblo(taborl   ,'TABORL')
      CALL allblo(tabors   ,'TABORS')
      CALL allblo(atabof   ,'ATABOF')
      CALL allblo(atabol   ,'ATABOL')
      CALL allblo(atabos   ,'ATABOS')
      CALL allblo(btabof   ,'BTABOF')
      CALL allblo(btabol   ,'BTABOL')
      CALL allblo(btabos   ,'BTABOS')
      CALL allblo(litabf   ,'LITABF')
      CALL allblo(litabl   ,'LITABL')
      CALL allblo(litabs   ,'LITABS')
      CALL allblo(trborsave,'TBSAVE')
      CALL allblo(ta_sce   ,'TA_SCE')
      CALL allblo(wchu     ,'WCHU  ')
!
      IF(ntrac.NE.0) THEN
        ielm   = ielm3
        ielh   = ielm2h
        ielv   = ielm2v
        statut = 1
      ELSE
        ielm   = 0
        ielh   = 0
        ielv   = 0
        statut = 0
      ENDIF
!
! GENERIC NAMES SHORTENED IN ORDER TO PROVIDE THE AUTOMATIC
! NUMBERING MECHANISM!
!
! MAX(NTRAC,1) BELOW : TO HAVE AT LEAST ONE ARRAY, EVEN EMPTY
!                      TO PUT IN THE CALL TO CONLIM
!                      AND CALL TO TEL4DEL
!
      CALL bief_allvec_in_block(trn,    ntrac,
     &                          1, 'TRN   ', ielm, 1, statut,mesh3d)
      CALL bief_allvec_in_block(tac,    ntrac,
     &                          1, 'TAC   ', ielm, 1, statut,mesh3d)
      CALL bief_allvec_in_block(ta,     max(ntrac,1),
     &                          1, 'TA    ', ielm, 1, statut,mesh3d)
      CALL bief_allvec_in_block(s0ta,   ntrac,
     &                          1, 'S0TA  ', ielm, 1, statut,mesh3d)
      CALL bief_allvec_in_block(s1ta,   ntrac,
     &                          1, 'S1TA  ', ielm, 1, statut,mesh3d)
      CALL bief_allvec_in_block(taborf, ntrac,
     &                          1, 'TABF  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(taborl, max(ntrac,1) ,
     &                          1, 'TABL  ', ielv, 1, statut,mesh3d)
      CALL bief_allvec_in_block(tabors, ntrac,
     &                          1, 'TABS  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(atabof, ntrac,
     &                          1, 'ATAF  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(atabol, max(ntrac,1),
     &                          1, 'ATAL  ', ielv, 1, statut,mesh3d)
      CALL bief_allvec_in_block(atabos, ntrac,
     &                          1, 'ATAS  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(btabof, ntrac,
     &                          1, 'BTAF  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(btabol, max(ntrac,1),
     &                          1, 'BTAL  ', ielv, 1, statut,mesh3d)
      CALL bief_allvec_in_block(btabos, ntrac,
     &                          1, 'BTAS  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(litabf, ntrac,
     &                          2, 'LTAF  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(litabl, max(ntrac,1) ,
     &                          2, 'LTAL  ', ielv, 2, statut,mesh3d)
      CALL bief_allvec_in_block(litabs, ntrac,
     &                          2, 'LTAS  ', ielh, 1, statut,mesh2d)
      CALL bief_allvec_in_block(trborsave,ntrac,
     &                          1, 'TBSA  ', ielv, 1, statut,mesh3d)
      CALL bief_allvec_in_block(ta_sce,ntrac,
     &                          1, 'TSCE  ', nsce, 1, 0     ,mesh3d)
!                                                     SIZE NSCE ALWAYS
      CALL bief_allvec_in_block(wchu,     max(ntrac,1),
     &                          1, 'WCHU  ', ielm, 1, statut,mesh3d)
!
!=======================================================================
! K-EPSILON MODEL
!
      IF(iturbv.EQ.3.OR.iturbh.EQ.3.OR.iturbv.EQ.7.OR.iturbh.EQ.7) THEN
        ielm   = ielm3
        ielh   = ielm2h
        ielv   = ielm2v
        statut = 1
      ELSE
        ielm   = 0
        ielh   = 0
        ielv   = 0
        statut = 0
      ENDIF
!
      CALL bief_allvec(1, akn,      'AKN   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, akc,      'AKC   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s0ak,     'S0AK  ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s1ak,     'S1AK  ',ielm,1, statut,mesh3d )
!
      CALL bief_allvec(1, epn,      'EPN   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, epc,      'EPC   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s0ep,     'S0EP  ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s1ep,     'S1EP  ',ielm,1, statut,mesh3d )
!
      CALL bief_allvec(1, kborf,    'KBORF ',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, kborl,    'KBORL ',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, kborsave, 'KBSAVE',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, kbors,    'KBORS ',ielh,1, statut,mesh2d )
!
      CALL bief_allvec(1, eborf,    'EBORF ',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, eborl,    'EBORL ',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, eborsave, 'EBSAVE',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, ebors,    'EBORS ',ielh,1, statut,mesh2d )
!
      CALL bief_allvec(1, akborf, 'AKBORF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(1, akborl, 'AKBORL', ielv, 1, statut,mesh3d )
      CALL bief_allvec(1, akbors, 'AKBORS', ielh, 1, statut,mesh2d )
!
      CALL bief_allvec(1, aeborf, 'AEBORF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(1, aeborl, 'AEBORL', ielv, 1, statut,mesh3d )
      CALL bief_allvec(1, aebors, 'AEBORS', ielh, 1, statut,mesh2d )
!
      CALL bief_allvec(1, bkborf, 'BKBORF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(1, bkborl, 'BKBORL', ielv, 1, statut,mesh3d )
      CALL bief_allvec(1, bkbors, 'BKBORS', ielh, 1, statut,mesh2d )
!
      CALL bief_allvec(1, beborf, 'BEBORF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(1, beborl, 'BEBORL', ielv, 1, statut,mesh3d )
      CALL bief_allvec(1, bebors, 'BEBORS', ielh, 1, statut,mesh2d )
!
      CALL bief_allvec(2, likbof, 'LIKBOF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(2, likbol, 'LIKBOL', ielv, 2, statut,mesh3d )
      CALL bief_allvec(2, likbos, 'LIKBOS', ielh, 1, statut,mesh2d )
!
      CALL bief_allvec(2, liebof, 'LIEBOF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(2, liebol, 'LIEBOL', ielv, 2, statut,mesh3d )
      CALL bief_allvec(2, liebos, 'LIEBOS', ielh, 1, statut,mesh2d )
!
      IF(    iturbv.EQ.3.OR.iturbh.EQ.3.OR.iturbv.EQ.7.OR.iturbh.EQ.7
     &   .OR.iturbv.EQ.5.OR.iturbh.EQ.5.OR.iturbv.EQ.9.OR.iturbh.EQ.9
     &   .OR.iturbv.EQ.6) THEN
        ielm   = ielm3
        ielh   = ielm2h
        ielv   = ielm2v
        statut = 1
      ELSE
        ielm   = 0
        ielh   = 0
        ielv   = 0
        statut = 0
      ENDIF
!
!     AK AND EP ARE USED WHEN WE RECONSTRUCT KINETIC ENERGY AND
!     DISSIPATION WITH THE SPALART TURBULENT VISCOSITY OR GOTM
!
      CALL bief_allvec(1, ak,       'AK    ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, ep,       'EP    ',ielm,1, statut,mesh3d )
!
      IF(iturbv.EQ.5.OR.iturbh.EQ.5.OR.iturbv.EQ.9.OR.iturbh.EQ.9
     &   .OR.iturbv.EQ.6) THEN
        ielm   = ielm3
        ielh   = ielm2h
        ielv   = ielm2v
        statut = 1
      ELSE
        ielm   = 0
        ielh   = 0
        ielv   = 0
        statut = 0
      ENDIF
!
      CALL bief_allvec(1, nun,      'NUN   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, nuc,      'NUC   ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, nu,       'NU    ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s0nu,     'S0NU  ',ielm,1, statut,mesh3d )
      CALL bief_allvec(1, s1nu,     'S1NU  ',ielm,1, statut,mesh3d )
!
      CALL bief_allvec(1, nuborf,   'NUBORF',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, nuborl,   'NUBORL',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, nubors,   'NUBORS',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, nuborsave,'NBSAVE',ielv,1, statut,mesh3d )

      CALL bief_allvec(1, anuborf,  'ANUBORF',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, anuborl,  'ANUBORL',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, anubors,  'ANUBORS',ielh,1, statut,mesh2d )

      CALL bief_allvec(1, bnuborf,  'BNUBORF',ielh,1, statut,mesh2d )
      CALL bief_allvec(1, bnuborl,  'BNUBORL',ielv,1, statut,mesh3d )
      CALL bief_allvec(1, bnubors,  'BNUBORS',ielh,1, statut,mesh2d )
!
      CALL bief_allvec(2, linubof, 'LINUBOF', ielh, 1, statut,mesh2d )
      CALL bief_allvec(2, linubol, 'LINUBOL', ielv, 2, statut,mesh3d )
      CALL bief_allvec(2, linubos, 'LINUBOS', ielh, 1, statut,mesh2d )
!
!=======================================================================
!
!                     *******************
!                     * OTHER VARIABLES *
!                     *******************
!
!  VARIOUS VELOCITY COMPONENTS 3D, 2D, BOUNDARY CONDITION VALUES
!
      CALL bief_allvec(1, u2d,        'U2D   ', ielmu,  1, 1,mesh2d)
      CALL bief_allvec(1, un2d,       'UN2D  ', ielmu,  1, 1,mesh2d)
      CALL bief_allvec(1, fu,         'FU    ', ielmu,  1, 1,mesh2d)
      CALL bief_allvec(1, v2d,        'V2D   ', ielmu,  1, 1,mesh2d)
      CALL bief_allvec(1, vn2d,       'VN2D  ', ielmu,  1, 1,mesh2d)
      CALL bief_allvec(1, fv,         'FV    ', ielmu,  1, 1,mesh2d)
!     DIMENSION 2 : SEE IN LECLIM, USED IN BORD3D
      CALL bief_allvec(1, ubor2d,     'UBOR2D', ielm1, 2, 1,mesh2d)
      CALL bief_allvec(1, vbor2d,     'VBOR2D', ielm1, 2, 1,mesh2d)
!
      CALL bief_allvec(1, flbor ,     'FLBOR ', ielm1, 1, 1,mesh2d)
      CALL bief_allvec(1, flblim,     'FLBLIM', ielm1, 1, 1,mesh2d)
!
      CALL bief_allvec(1, uconv ,      'UCONV ', ielm3, 1, 1,mesh3d)
      CALL bief_allvec(1, vconv ,      'VCONV ', ielm3, 1, 1,mesh3d)
      CALL bief_allvec(1, vconvc,      'VCONVC', ielm3, 1, 1,mesh3d)
      CALL bief_allvec(1, uconvc,      'UCONVC', ielm3, 1, 1,mesh3d)
      CALL bief_allvec(1, wsconv,      'WSCONV', ielm3, 1, 1,mesh3d)
      CALL bief_allvec(1, dm1   ,      'DM1   ', ielm3, 1, 1,mesh3d)
!
!-----------------------------------------------------------------------
! WATER DEPTH AND VARIABLES DERIVED FROM IT,
! PRINCIPALLY T2D DISCRETISATION OF DEPTH IELMH = IELM2H
!
      CALL bief_allvec(1, h,      'H     ', ielmh, 1, 1,mesh2d)
      CALL bief_allvec(1, hn,     'HN    ', ielmh, 1, 1,mesh2d)
      CALL bief_allvec(1, hprop,  'HPROP ', ielmh, 1, 1,mesh2d)
      CALL bief_allvec(1, nuwave, 'NUWAVE', 10   , 1, 1,mesh2d)
      CALL bief_allvec(1, hbor,   'HBOR  ', ielm1, 1, 1,mesh2d)
      CALL bief_allvec(1, dh,     'DH    ', ielmh, 1, 2,mesh2d)
      CALL bief_allvec(1, smh,    'SMH   ', ielmh, 1, 1,mesh2d)
!
      IF(rain) THEN
        CALL bief_allvec(1,pluie,'PLUIE ',ielmh,1,2,mesh2d)
        IF(ncsize.GT.1) THEN
          CALL bief_allvec(1,parapluie,'PARAPL',ielmh,1,2,mesh2d)
        ELSE
          CALL bief_allvec(1,parapluie,'PARAPL',0    ,1,0,mesh2d)
          DEALLOCATE(parapluie%R)
          parapluie%R=>pluie%R
        ENDIF
      ELSE
        CALL bief_allvec(1,pluie    ,'PLUIE ',0    ,1,0,mesh2d)
        CALL bief_allvec(1,parapluie,'PARAPL',0    ,1,0,mesh2d)
      ENDIF
!
      IF(sorg2d(35)) THEN
        CALL bief_allvec(1,maxz,'MAXZ  ',ielmh,1,1 ,mesh2d)
      ELSE
        CALL bief_allvec(1,maxz,'MAXZ  ',0    ,1,0 ,mesh2d)
      ENDIF
      IF(sorg2d(36)) THEN
        CALL bief_allvec(1,tmaxz,'TMAXZ ',ielmh,1,1 ,mesh2d)
      ELSE
        CALL bief_allvec(1,tmaxz,'TMAXZ ',0    ,1,0 ,mesh2d)
      ENDIF
!
!-----------------------------------------------------------------------
!
      IF(courou.OR.inclus(coupling,'TOMAWAC')) THEN
        CALL bief_allvec(1,fxh,'FXH   ',ielmu,1,1,mesh2d)
        CALL bief_allvec(1,fyh,'FYH   ',ielmu,1,1,mesh2d)
      ELSE
        CALL bief_allvec(1,fxh,'FXH   ',0    ,1,0,mesh2d)
        CALL bief_allvec(1,fyh,'FYH   ',0    ,1,0,mesh2d)
      ENDIF
!
!     VARIABLES TRANSMITTED FROM TOMAWAC TO SISYPHE OR GAIA
!
      IF(inclus(coupling,'TOMAWAC')) THEN
        CALL bief_allvec(1,dirmoy,'DIRMOY',ielmh,1,1,mesh2d)
        CALL bief_allvec(1,hm0   ,'HM0   ',ielmh,1,1,mesh2d)
        CALL bief_allvec(1,tpr5  ,'TPR5  ',ielmh,1,1,mesh2d)
        CALL bief_allvec(1,orbvel,'ORBVEL',ielmh,1,1,mesh2d)
      ELSE
        CALL bief_allvec(1,dirmoy,'DIRMOY',0    ,1,0,mesh2d)
        CALL bief_allvec(1,hm0   ,'HM0   ',0    ,1,0,mesh2d)
        CALL bief_allvec(1,tpr5  ,'TPR5  ',0    ,1,0,mesh2d)
        CALL bief_allvec(1,orbvel,'ORBVEL',0    ,1,0,mesh2d)
      ENDIF
!
!-----------------------------------------------------------------------
! NUMERICAL VARIABLES
!
      CALL bief_allvec(1, volu  , 'VOLU  ', ielm3 , 1, 1,mesh3d)
      CALL bief_allvec(1, volun , 'VOLUN ', ielm3 , 1, 1,mesh3d)
      CALL bief_allvec(1, volu3d, 'VOLU3D', ielm3 , 1, 1,mesh3d)
      IF(ncsize.GT.1) THEN
        CALL bief_allvec(1,volupar,'VOLUPA',ielm3,1,1,mesh3d)
        CALL bief_allvec(1,volunpar,'VLNPAR',ielm3,1,1,mesh3d)
        CALL bief_allvec(1,volu3dpar,'VLTDPA',ielm3,1,1,mesh3d)
      ELSE
        CALL bief_allvec(1,volupar,'VOLUPA',0    ,1,0,mesh3d)
        DEALLOCATE(volupar%R)
        volupar%R=>volu%R
        CALL bief_allvec(1,volunpar,'VLNPAR',0    ,1,0,mesh3d)
        DEALLOCATE(volunpar%R)
        volunpar%R=>volun%R
        CALL bief_allvec(1,volu3dpar,'VLTDPA',0    ,1,0,mesh3d)
        DEALLOCATE(volu3dpar%R)
        volu3dpar%R=>volu3d%R
      ENDIF
      CALL bief_allvec(1, volut , 'VOLUT ', ielm3 , 1, 1,mesh3d)
      CALL bief_allvec(1, volu2d, 'VOLU2D', ielm2h, 1, 1,mesh2d)
      CALL bief_allvec(1, v2dpar, 'V2DPAR', ielm2h, 1, 1,mesh2d)
      CALL bief_allvec(1, unsv2d, 'UNSV2D', ielm2h, 1, 1,mesh2d)
      CALL bief_allvec(1, unsv3d, 'UNSV3D', ielm3 , 1, 1,mesh3d)
!
      CALL bief_allvec(1, fluint,    'FLUINT', ielm3 , 1, 1,mesh3d)
      CALL bief_allvec(1, fluext,    'FLUEXT', ielm3 , 1, 1,mesh3d)
      IF(ncsize.GT.1) THEN
        CALL bief_allvec(1,fluextpar,'FLXTPA',ielm3,1,1,mesh3d)
      ELSE
        CALL bief_allvec(1,fluextpar,'FLXTPA',1    ,1,0,mesh3d)
        DEALLOCATE(fluextpar%R)
        fluextpar%R=>fluext%R
      ENDIF
      CALL bief_allvec(1, flint2,    'FLINT2', ielm2h, 1, 1,mesh2d)
!
!-----------------------------------------------------------------------
! PHYSICAL VARIABLES
!
! (COMPONENTS OF VISCOSITY AS SEPARATE VECTORS ARE REQUIRED!)
!
      CALL allblo(viscvi, 'VISCVI')
      CALL allblo(viscnu, 'VISCNU')
      CALL bief_allvec_in_block(viscvi,3,1,'VISC  ',ielm3,1,2,mesh3d)
      CALL bief_allvec_in_block(viscnu,3,1,'VISCN ',ielm3,1,2,mesh3d)
!
! ADDRESSING IS AWKWARD, E.G. REAL PART OF THE DIFFUSIVITY STRUCTURE
! OF THE THIRD TRACER IN THE Y-DIRECTION IS VISCTA%ADR(3)%P%ADR(2)%P%R
! DEAR ME!
!
      CALL allblo(viscta, 'VISCTA')
      IF(ntrac.GT.0) THEN
        CALL allblo_in_block(viscta, ntrac, 'VBTA  ')
        DO itrac = 1,ntrac
          CALL bief_allvec_in_block ( viscta%ADR(itrac)%P, 3 , 1,
     &                              'VITA  ',ielm3,1,2,mesh3d)
        ENDDO
      ENDIF
!
      CALL bief_allvec(1, rotat,  'ROTAT ', ielm3,  1, 1,mesh3d)
!
      CALL bief_allvec(1, deltar, 'DELTAR', ielm3,  1, 1,mesh3d)
      CALL bief_allvec(1, rho   , 'RHO   ', ielm3,  1, 1,mesh3d)
      IF(iturbv.EQ.6) THEN
        CALL bief_allvec(1, rhopot, 'RHOPOT', ielm3,1, 1,mesh3d)
      ELSE
        CALL bief_allvec(1, rhopot, 'RHOPOT', 0,1,0,mesh3d)
      ENDIF
      CALL bief_allvec(1, ri,     'RI    ', ielm3,  1, 1,mesh3d)
!
      CALL bief_allvec(1, rugof,  'RUGOF ', ielmu , 1, 1,mesh2d)
      CALL bief_allvec(1, rugol,  'RUGOL ', ielm2v, 1, 1,mesh3d)
      CALL bief_allvec(1, cf   ,  'CF    ', ielmu , 1, 1,mesh2d)
!
!     WINDX, WINDY, PATMOS NOW ALLOCATED IN METEO_TELEMAC
!
!     CALL ALLBLO(WIND, 'WIND  ')
!     IF (VENT.OR.INCLUS(COUPLING,'WAQTEL')) THEN
!       CALL BIEF_ALLVEC_IN_BLOCK(WIND, 2, 1,'WIND  ',IELM2H,1,1,MESH2D)
!     ELSE
!       CALL BIEF_ALLVEC_IN_BLOCK(WIND, 2, 1,'WIND  ',0,1,0,MESH2D)
!     ENDIF
!
      IF( (inclus(coupling,'WAQTEL') ).AND.
     &    (13*int(waqprocess/13).EQ.waqprocess) ) THEN
        CALL bief_allvec(1, windspd, 'WNDSPD', ielm2h, 1,1,mesh2d)
        CALL bief_allvec(1, rho4aed2, 'RHOAED', ielm3,  1,1,mesh3d)
        CALL bief_allvec(1, temp4aed2, 'TEMPAE', ielm3,  1,1,mesh3d)
        CALL bief_allvec(1, sali4aed2, 'SALIAE', ielm3,  1,1,mesh3d)
        CALL bief_allvec(1, thick4aed2, 'THICKA', ielm3,  1,1,mesh3d)
        CALL bief_allvec(2, benth4aed2, 'BENTHA', ielm2h,  1,1,mesh2d)
      ELSE
        CALL bief_allvec(1, windspd, 'WNDSPD', 0,1,0,mesh2d)
        CALL bief_allvec(1, rho4aed2, 'RHOAED', 0,1,0,mesh3d)
        CALL bief_allvec(1, temp4aed2, 'TEMPAE', 0,1,0,mesh3d)
        CALL bief_allvec(1, sali4aed2, 'SALIAE', 0,1,0,mesh3d)
        CALL bief_allvec(1, thick4aed2, 'THICKA', 0,1,0,mesh3d)
        CALL bief_allvec(2, benth4aed2, 'BENTHA', 0,1,0,mesh2d)
      ENDIF
!
!     IF (ATMOS) THEN
!       CALL BIEF_ALLVEC(1, PATMOS, 'PATMOS', IELM2H, 1,1,MESH2D)
!     ELSE
!       CALL BIEF_ALLVEC(1, PATMOS, 'PATMOS', 0,      1,0,MESH2D)
!     ENDIF
!
!-----------------------------------------------------------------------
!     VARIABLES USED IN MASS BALANCE
!     THEY WILL CORRESPOND TO THE NUMBERING OF ADVECTED VARIABLES
!
      CALL bief_allvec(1, masini, 'MASINI', nvbil,1, 0,mesh3d)
      CALL bief_allvec(1, masse,  'MASSE ', nvbil,1, 0,mesh3d)
      CALL bief_allvec(1, massen, 'MASSEN', nvbil,1, 0,mesh3d)
      CALL bief_allvec(1, flucum, 'FLUCUM', nvbil,1, 0,mesh3d)
      CALL bief_allvec(1, flux,   'FLUX  ', nvbil,1, 0,mesh3d)
      CALL bief_allvec(1, fluxb,  'FLUXB ', nvbil,maxfro+nsce+1,
     &                 0,mesh3d)
!
! BOTTOM AS GIVEN FROM THE GEOMETRY FILE
!
      CALL bief_allvec(1, zf, 'ZF    ', ielm2h, 1, 1,mesh2d)
!
! BOTTOM GEOMETRY PER 2D-ELEMENT FOR TIDAL FLATS TREATMENT
!
      IF(msk) THEN
        CALL bief_allvec(1, zfe, 'ZFE   ', 10, 1, 1,mesh2d)
      ELSE
        CALL bief_allvec(1, zfe, 'ZFE   ',  0, 1, 0,mesh2d)
      ENDIF
!
! 3D ARRAY TO STORE DH AND HN
!
      CALL bief_allvec(1, dhhn,  'DHHN  ', ielm3,  1, 1,mesh3d)
!
!-----------------------------------------------------------------------
! DROGUES (FLOATS ...EHM, TRACERS...)
!
      CALL bief_allvec(1, xflot , 'XFLOT ',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(1, yflot , 'YFLOT ',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(1, zflot , 'ZFLOT ',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(1, shpflo, 'SHPFLO', 3*nflot_max,1, 0,mesh3d)
      CALL bief_allvec(1, shzflo, 'SHZFLO',   nflot_max,1, 0,mesh3d)
!
      CALL bief_allvec(2, tagflo, 'TAGFLO',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(2, clsflo, 'CLSFLO',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(2, eltflo, 'ELTFLO',   nflot_max,1, 0,mesh3d)
      CALL bief_allvec(2, etaflo, 'ETAFLO',   nflot_max,1, 0,mesh3d)
!
!-----------------------------------------------------------------------
!
!     VALUES AT SOURCES OF ADVECTED VARIABLES U,V,W,AK AND EP
!     ALL ARE ALLOCATED, EVEN IF NOT USED
!
      CALL bief_allvec(1,u_sce ,'U_SCE ',nsce,1,0,mesh3d)
      CALL bief_allvec(1,v_sce ,'V_SCE ',nsce,1,0,mesh3d)
      CALL bief_allvec(1,w_sce ,'W_SCE ',nsce,1,0,mesh3d)
      CALL bief_allvec(1,ak_sce,'AK_SCE',nsce,1,0,mesh3d)
      CALL bief_allvec(1,ep_sce,'EP_SCE',nsce,1,0,mesh3d)
      CALL bief_allvec(1,nu_sce,'NU_SCE',nsce,1,0,mesh3d)
!
!     COPYING USCE, VSCE AND WSCE INTO U_SCE, V_SCE AND W_SCE
!     PROVISIONNALY CANCELLING AK_SCE AND EP_SCE
!
      IF(nsce.GT.0) THEN
        DO i=1,nsce
          u_sce%R(i) =usce(i)
          v_sce%R(i) =vsce(i)
          w_sce%R(i) =wsce(i)
          ak_sce%R(i)=0.d0
          ep_sce%R(i)=0.d0
          nu_sce%R(i)=0.d0
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
! VARIOUS TABLES, AS W1, ITRAV3 , PRIVE
! W1 ALLOCATED BY ALMESH
!
      CALL allblo (prive,'PRIVE ')
      CALL bief_allvec_in_block(prive,max(4,npriv),
     &                          1,'PRIV  ',ielm3,1,2,mesh3d)
!     ALIASES
      prive1=>prive%ADR(1)%P
      prive2=>prive%ADR(2)%P
      prive3=>prive%ADR(3)%P
      prive4=>prive%ADR(4)%P
!
      CALL allblo (prive2d,'PRIVE2')
      CALL bief_allvec_in_block(prive2d,max(4,npriv2d),
     &                          1,'PRI2  ',ielm2h,1,2,mesh2d)
!     ALIASES
      prive2d1=>prive2d%ADR(1)%P
      prive2d2=>prive2d%ADR(2)%P
      prive2d3=>prive2d%ADR(3)%P
      prive2d4=>prive2d%ADR(4)%P
!
! INTEGER WORK FIELDS
!
      CALL allblo(itrav3,'ITRAV3')
      siz=bief_nbmpts(ielm3,mesh3d)
!                                         40 OR 50
      IF(yesweak) siz=max(siz,bief_nbmpts(10*(ielm3/10),mesh3d)*ngauss)
      CALL bief_allvec_in_block(itrav3,4,2,'ITR3V ',siz,1,0,mesh3d)
!
! POINTERS TO 4 INTEGER WORK FIELDS
!
      it1 => itrav3%ADR(1)%P
      it2 => itrav3%ADR(2)%P
      it3 => itrav3%ADR(3)%P
      it4 => itrav3%ADR(4)%P
!
!     BLOCKS OF ADVECTED VARIABLES AND VARIOUS BIEF_OBJ STRUCTURES
!     BUILT FOR A COLLECTIVE ADVECTION
!
      CALL allblo(bl_fc ,'BL_FC ')
      CALL addblo(bl_fc,uc)
      CALL addblo(bl_fc,vc)
      CALL addblo(bl_fc,wc)
      CALL addblo(bl_fc,akc)
      CALL addblo(bl_fc,epc)
      CALL addblo(bl_fc,nuc)
!
      CALL allblo(bl_fn ,'BL_FN ')
!     ADVECTION OF VELOCITIES IS DONE FOR THE NEXT TIME-STEP
      CALL addblo(bl_fn,u)
      CALL addblo(bl_fn,v)
      CALL addblo(bl_fn,w)
      CALL addblo(bl_fn,akn)
      CALL addblo(bl_fn,epn)
      CALL addblo(bl_fn,nun)
!
      CALL allblo(bl_s0f,'BL_S0F')
      CALL addblo(bl_s0f,s0u)
      CALL addblo(bl_s0f,s0v)
      CALL addblo(bl_s0f,s0w)
      CALL addblo(bl_s0f,s0ak)
      CALL addblo(bl_s0f,s0ep)
      CALL addblo(bl_s0f,s0nu)
!
      CALL allblo(bl_fsc,'BL_FSC')
      CALL addblo(bl_fsc,u_sce)
      CALL addblo(bl_fsc,v_sce)
      CALL addblo(bl_fsc,w_sce)
      CALL addblo(bl_fsc,ak_sce)
      CALL addblo(bl_fsc,ep_sce)
      CALL addblo(bl_fsc,nu_sce)
!
      CALL allblo(bl_bol,'BL_BOL')
      CALL addblo(bl_bol,liubol)
      CALL addblo(bl_bol,livbol)
      CALL addblo(bl_bol,liwbol)
      CALL addblo(bl_bol,likbol)
      CALL addblo(bl_bol,liebol)
      CALL addblo(bl_bol,linubol)
!
      CALL allblo(bl_borl,'BL_BOR')
      CALL addblo(bl_borl,uborl)
      CALL addblo(bl_borl,vborl)
      CALL addblo(bl_borl,wborl)
      CALL addblo(bl_borl,akborl)
      CALL addblo(bl_borl,eborl)
      CALL addblo(bl_borl,nuborl)
!
      IF(ntrac.GT.0) THEN
        DO itrac=1,ntrac
          CALL addblo(bl_fc ,    tac%ADR(itrac)%P)
          CALL addblo(bl_fn ,    trn%ADR(itrac)%P)
          CALL addblo(bl_s0f,   s0ta%ADR(itrac)%P)
          CALL addblo(bl_fsc, ta_sce%ADR(itrac)%P)
          CALL addblo(bl_bol, litabl%ADR(itrac)%P)
          CALL addblo(bl_borl,taborl%ADR(itrac)%P)
        ENDDO
      ENDIF
!
!     BLOCKS OF ADVECTED VARIABLES WITH CHARACTERISTICS
!     MODULE STREAMLINE MUST BE CHANGED TO DEAL WITH A
!     LIST IN A BLOCK AS DONE ABOVE
!
      CALL allblo(fc3d, 'FC3D  ')
      CALL allblo(fn3d, 'FN3D  ')
!     ADVECTION OF VELOCITIES IS DONE FOR THE NEXT STEP
!     SO HERE U AND NOT UN
      IF(schcvi.EQ.adv_car) THEN
        CALL addblo(fn3d,u )
        CALL addblo(fc3d,uc)
        CALL addblo(fn3d,v )
        CALL addblo(fc3d,vc)
        IF(nonhyd) THEN
          CALL addblo(fn3d,w )
          CALL addblo(fc3d,wc)
        ENDIF
      ENDIF
      IF(schcke.EQ.adv_car.AND.(iturbh.EQ.3.OR.
     &                          iturbh.EQ.7.OR.
     &                          iturbv.EQ.3.OR.
     &                          iturbv.EQ.7)) THEN
        CALL addblo(fn3d,akn)
        CALL addblo(fc3d,akc)
        CALL addblo(fn3d,epn)
        CALL addblo(fc3d,epc)
!     FOR SPALART
      ELSEIF(schcke.EQ.adv_car.AND.(iturbv.EQ.5.OR.
     &                              iturbh.EQ.5.OR.
     &                              iturbv.EQ.9.OR.
     &                              iturbh.EQ.9)) THEN
        CALL addblo(fn3d,nun)
        CALL addblo(fc3d,nuc)
      ENDIF
      IF(ntrac.GT.0) THEN
        DO itrac=1,ntrac
          IF(schcta(itrac).EQ.adv_car) THEN
            CALL addblo(fn3d,trn%ADR(itrac)%P)
            CALL addblo(fc3d,tac%ADR(itrac)%P)
          ENDIF
        ENDDO
      ENDIF
!
!=======================================================================
!
!                   *********************************
!                   * STRUCTURES FOR THE RESOLUTION *
!                   *  OF LINEAR EQUATION SYSTEMS   *
!                   *          F I R S T            *
!                   *        M A T R I C E S        *
!                   *********************************
!
! BEWARE : 2D AND ESPECIALLY 3D MATRICES OCCUPY
! ======   A LARGE CHUNK OF MEMORY
!
!-----------------------------------------------------------------------
! 3D MATRICES
!
! S.U.P.G.
!
      typdia = '0'
      typext = '0'
      IF(n_adv(adv_sup).GT.0) THEN
        typdia = 'Q'
        typext = 'Q'
      ENDIF
      CALL bief_allmat(msupg,'SUPG  ',
     &                 ielm3,ielm3,cfg,typdia,typext,mesh3d)
!
! M.U.R.D.
!
      typdia = '0'
      typext = '0'
      IF(n_adv(adv_nsc).GT.0.OR.n_adv(adv_psi).GT.0) THEN
        typdia = 'Q'
        typext = 'Q'
      ENDIF
      CALL bief_allmat(mmurd,'MURD  ',
     &                 ielm3,ielm3,cfgmurd,typdia,typext,mesh3d)
!
! M.U.R.D. (EDGE-BASED FOR TIDAL FLATS)
!
      typdia = '0'
      typext = '0'
      IF(n_adv(adv_nsc_tf).GT.0) THEN
        typdia = 'Q'
        typext = 'Q'
      ENDIF
      CALL bief_allmat(murd_tf,'MURDTF',ielm3,ielm3,cfgmurd_tf,
     &                 typdia,typext,mesh3d)
!
! DIFFUSION
!
      typdia = 'Q'
      typext = 'S'
      CALL bief_allmat(mdiff,'DIFF  ',
     &                 ielm3,ielm3,cfg,typdia,typext,mesh3d)
!
! THE 3D WORK MATRICES (ALWAYS ALLOCATED AS NON SYMMETRICAL)
!                       SEE USE OF MTRA2%X IN WAVE_EQUATION
!
      typdia = 'Q'
      typext = 'Q'
      CALL bief_allmat (mtra2, 'MTRA2 ',
     &                  ielm3, ielm3, cfg   , typdia, typext,mesh3d)
      CALL bief_allmat (mtra1, 'MTRA1 ',
     &                  ielm3, ielm3, cfg   , typdia, typext,mesh3d)
!
!-----------------------------------------------------------------------
! 3 2D MATRICES (IELM2H) - (EACH OF THEM ALLOCATED IN A
! SEPARATE BLOCK), ALL IN BLOCK MAT2D
! E.G. WORK MATRIX 3 IS: MAT2D%ADR(3)%P
! THEY ALL GET NAMED MAT2D...............
!
      CALL allblo (mat2d, 'MAT2D ')
      CALL allblo_in_block(mat2d,3,'MAT2D ')
      DO itab = 1,3
        CALL bief_deallobj(mat2d%ADR(itab)%P)
        DEALLOCATE(mat2d%ADR(itab)%P)
        CALL first_all_biefobj(mat2d%ADR(itab)%P)
        CALL bief_allmat(mat2d%ADR(itab)%P, 'MAT2D ',
     &                   ielmu, ielmu, cfg2d, 'Q', 'Q',mesh2d)
        mat2d%ADR(itab)%P%FATHER='MAT2D '
      END DO
!
! ANOTHER WORK MATRIX IS BUILT USING ALMESH FOR THE 2D MESH: MTRA2
!
! BOUNDARY MATRIX FOR 2D
!
      CALL bief_allmat(mbor2d, 'MBOR2D',
     &                 ielm1, ielm1, cfgbor2d, 'Q','Q',mesh2d)
!
! 2D WORK MATRIX FOR IELMU.
!
      CALL bief_allmat(matr2h,'MATR2H',ielmu ,ielmu ,cfg,'Q','Q',mesh2d)
!
!=======================================================================
!
!                   *********************************
!                   * STRUCTURES FOR THE RESOLUTION *
!                   *  OF LINEAR EQUATION SYSTEMS   *
!                   *  S E C O N D   M E M B E R S  *
!                   *         V E C T O R S         *
!                   *********************************
!
!-----------------------------------------------------------------------
! FOR 3D PART
! COMPUTES THE NUMBER OF 3D WORK FIELDS AS A FUNCTION OF
! SOLVER AND PRECONDITIONING TYPE AND NUMBER OF VARIABLES
! CVDF3D NEEDS AT LEAST 10.
      ntr = 10
!
      IF(slvdvi%SLV.EQ.7) ntr = max(ntr,2+2*slvdvi%KRYLOV)
      IF(ntrac.GT.0) THEN
        DO itrac=1,ntrac
          IF(slvdta(itrac)%SLV.EQ.7) THEN
            ntr = max(ntr,2+2*slvdta(itrac)%KRYLOV)
          ENDIF
        ENDDO
      ENDIF
      IF (slvdke%SLV.EQ.7) ntr = max(ntr,2+2*slvdke%KRYLOV)
      IF (s3d_slvdse%SLV.EQ.7) ntr = max(ntr,2+2*s3d_slvdse%KRYLOV)
      IF (slvpoi%SLV.EQ.7) ntr = max(ntr,2+2*slvpoi%KRYLOV)
! NO MORE SLVPRJ
!     IF (SLVPRJ%SLV.EQ.7) NTR = MAX(NTR,2+2*SLVPRJ%KRYLOV)
!
      i = ntr
!
      IF (3*(slvdvi%PRECON/3).EQ.slvdvi%PRECON .OR.
     &    3*(slvdke%PRECON/3).EQ.slvdke%PRECON .OR.
     &    3*(s3d_slvdse%PRECON/3).EQ.s3d_slvdse%PRECON .OR.
     &    3*(slvpoi%PRECON/3).EQ.slvpoi%PRECON) ntr = i+2
!
      IF(ntrac.GT.0) THEN
        DO itrac=1,ntrac
          IF(3*(slvdta(itrac)%PRECON/3).EQ.slvdta(1)%PRECON) ntr = i+2
        ENDDO
      ENDIF
!
      IF (iturbv.EQ.3.OR.iturbh.EQ.3) ntr = max(ntr,12)
      IF (iturbv.EQ.6) ntr = max(ntr,12)
      IF (iturbv.EQ.7.OR.iturbh.EQ.7) ntr = max(ntr,18)
!     TRAV3 WILL BE USED IN BIEF_VALIDA WITH THIS SIZE
!     SEE ALIRE3D IN TELEMAC3D.F AND VARSO3 IN POINT_TELEMAC3D.F
!     ON 04/07/2012 ADR_TRAC=29
      ntr = max(ntr,adr_trac+ntrac-1)
!
      CALL allblo(trav3, 'TRAV3 ')
      CALL bief_allvec_in_block(trav3,ntr,1,'TRAV  ',ielm3,1,2,mesh3d)
!
!     ANOTHER SUCH BLOCK FOR THE LIPS SCHEME
!
      CALL allblo(tb2, 'TB2   ')
      IF(yeslips) THEN
!       9 IS THE VALUE CORRESPONDING TO SOLVER 7=GMRES
        CALL bief_allvec_in_block(tb2,9,1,'TB2   ',ielm3,1,2,mesh3d)
      ELSE
        CALL bief_allvec_in_block(tb2,9,1,'TB2   ',    0,1,0,mesh3d)
      ENDIF
!
! POINTERS TO 3D WORK VECTORS; FIRST 10 EXIST FOR SURE
!
      t3_01 => trav3%ADR(01)%P
      t3_02 => trav3%ADR(02)%P
      t3_03 => trav3%ADR(03)%P
      t3_04 => trav3%ADR(04)%P
      t3_05 => trav3%ADR(05)%P
      t3_06 => trav3%ADR(06)%P
      t3_07 => trav3%ADR(07)%P
      t3_08 => trav3%ADR(08)%P
      t3_09 => trav3%ADR(09)%P
      t3_10 => trav3%ADR(10)%P
      t3_11 => trav3%ADR(11)%P
      t3_12 => trav3%ADR(12)%P
      t3_13 => trav3%ADR(13)%P
      t3_14 => trav3%ADR(14)%P
      t3_15 => trav3%ADR(15)%P
      t3_16 => trav3%ADR(16)%P
      t3_17 => trav3%ADR(17)%P
      t3_18 => trav3%ADR(18)%P
!
! SECOND MEMBER 3D
!
      CALL bief_allvec(1, sem3d, 'SEM3D ', ielm3, 1, 2,mesh3d)
!
! RIGHT HAND SIDE OF CONTINUITY EQUATIONS IF SOURCES
!
      CALL allblo(sources, 'SOURCE')
      i=nsce
      IF(ncsize.GT.1) i=2*nsce
      IF(nsce.GT.0) THEN
        IF(optsou.EQ.1) THEN
          CALL bief_allvec_in_block(sources,i,1,'SCE   ',ielm3,
     &                              1,1,mesh3d)
        ELSE IF(optsou.EQ.2) THEN
          i=1
          IF(ncsize.GT.1) i=2
          CALL bief_allvec_in_block(sources,i,1,'SCE   ',ielm3,
     &                              1,1,mesh3d)
        ENDIF
      ELSE
        CALL bief_allvec_in_block(sources,i,1,'SCE   ',0    ,1,0,mesh3d)
      ENDIF
!
!-----------------------------------------------------------------------
! FOR 2D PART
!
! COMPUTES THE NUMBER OF 2D WORK FIELDS AS A FUNCTION OF
! SOLVER AND PRECONDITIONING TYPE AND NUMBER OF VARIABLES
! FIRST 23 ALLOCATED FOR SURE...
!
      ntr = 23
      IF (slvpro%SLV.EQ.7) ntr = max(ntr,6+6*slvpro%KRYLOV)
! NO MORE SLVW
!     IF (SLVW%SLV.EQ.7)   NTR = MAX(NTR,2+2*SLVW%KRYLOV)
!     IF (3*(SLVPRO%PRECON/3) .EQ. SLVPRO%PRECON .OR.
!    &    3*(SLVW%PRECON/3)   .EQ. SLVW%PRECON ) NTR = NTR + 6
      IF (3*(slvpro%PRECON/3) .EQ. slvpro%PRECON) ntr = ntr + 6
!
!     SEE VARSOR BELOW
      ntr=max(ntr,adr_trav2+2*ntrac)
!
      CALL allblo(trav2, 'TRAV2 ')
      CALL bief_allvec_in_block(trav2,ntr,1,'TR2D  ',ielmu,1,2,mesh2d)
!
! POINTERS TO FIRST 10 2D WORK VECTORS (21 EXIST FOR SURE)
!
      t2_01 => trav2%ADR( 1)%P
      t2_02 => trav2%ADR( 2)%P
      t2_03 => trav2%ADR( 3)%P
      t2_04 => trav2%ADR( 4)%P
      t2_05 => trav2%ADR( 5)%P
      t2_06 => trav2%ADR( 6)%P
      t2_07 => trav2%ADR( 7)%P
      t2_08 => trav2%ADR( 8)%P
      t2_09 => trav2%ADR( 9)%P
      t2_10 => trav2%ADR(10)%P
      t2_11 => trav2%ADR(11)%P
      t2_12 => trav2%ADR(12)%P
      t2_13 => trav2%ADR(13)%P
      t2_14 => trav2%ADR(14)%P
      t2_15 => trav2%ADR(15)%P
      t2_16 => trav2%ADR(16)%P
      t2_17 => trav2%ADR(17)%P
      t2_18 => trav2%ADR(18)%P
      t2_19 => trav2%ADR(19)%P
      t2_20 => trav2%ADR(20)%P
      t2_21 => trav2%ADR(21)%P
      t2_22 => trav2%ADR(22)%P
      t2_23 => trav2%ADR(23)%P
!
!     SECOND MEMBERS BLOCK - 3 NEEDED
!
      CALL allblo(sem2d,'SEM2D ')
      CALL bief_allvec_in_block(sem2d,3,1,'SEM2D ',ielmu,1,2,mesh2d)
!
!     WORK FIELDS WITH A DIMENSION OF THE MAX. 2D ELEMENT NUMBER
!     IELM = 10
!
      CALL bief_allvec(1,   te1, 'TE1   ', 10, 1, 1,mesh2d)
      CALL bief_allvec(1,   te2, 'TE2   ', 10, 1, 1,mesh2d)
      CALL bief_allvec(1,   te3, 'TE3   ', 10, 1, 1,mesh2d)
!
!     PIECE-WISE LINEAR FREE SURFACE
!
      CALL bief_allvec(1,zflats, 'ZFLATS', 15, 1, 1,mesh2d)
      CALL bief_allvec(1,zconv , 'ZCONV ', 15, 1, 1,mesh2d)
!
!     TWO COMPONENTS OF ZCONV HORIZONTAL GRADIENT
!
      CALL bief_allvec(1,grazco, 'GRAZCO', 10, 2, 1,mesh2d)
!
! DELWAQ FLOWS BETWEEN POINTS
!
      CALL bief_allvec(1,flodel,'FLODEL',
     &                 nseg*nplan+netage*npoin2, 1, 0,mesh3d)
!     FULL SIZE OF FLOPAR PROBABLY NOT USEFUL IF NOT PARALLEL MODE
      CALL bief_allvec(1,flopar,'FLOPAR',
     &                 nseg*nplan+netage*npoin2, 1, 0,mesh3d)
! LIMITATION OF 2D SEGMENT FLUXES
      CALL bief_allvec(1,flulim,'FLULIM',nseg,1,0,mesh3d)
!
!-----------------------------------------------------------------------
! ALLOCATES VOID STRUCTURE
!
      CALL bief_allvec(1, svide, 'SVIDE ', 0, 1, 1,mesh3d)
!
!=======================================================================
! SEDIMENT
!  IF CONSOLIDATION IS MODELLED USING THE MULTI-LAYER MODEL,
!  THE MAXIMUM NUMBER OF PLANES DISCRETISING THE BED: S3D_NPFMAX
!  EQUALS (S3D_NCOUCH+1) :
!
! ALSO ALLOCATE 1DV BLOCK (NPLANx11) FOR SET_DIF WORK ARRAYS
!     allocate trav1 for all cases (it might be used by gaia)
      CALL allblo (trav1,'TRAV1 ')
      CALL bief_allvec_in_block(trav1,11,1,'TRAV1 ',nplan,1,0,mesh3d)
      IF(s3d_sedi) THEN
        CALL bief_allvec(1, s3d_wchu,      'S3D_WCHU', ielm3, 1,1,
     &                  mesh3d)
      ELSE
        CALL bief_allvec(1, s3d_wchu,      'S3D_WCHU', 0    , 1,0,
     &                   mesh3d)
      ENDIF
!
! use ncouch + 1 directly
      IF(s3d_tasse.AND. s3d_sedi) s3d_npfmax= s3d_ncouch+ 1
!
!   CV   IF((S3D_GIBSON.OR.S3D_TASSE).AND.S3D_SEDI) THEN
!      IF(S3D_NCOUCH.GT.1) THEN
!        CALL BIEF_ALLVEC(1,S3D_EPAI,'EPAI  ',S3D_NCOUCH*NPOIN2,1,0,MESH3D)
!      ELSE
!        CALL BIEF_ALLVEC(1,S3D_EPAI,'EPAI  ',0,                1,0,MESH3D)
!      ENDIF
!
      IF(s3d_tasse.AND.s3d_sedi) THEN
        CALL bief_allvec(1, s3d_ivide, 'IVIDE ', s3d_npfmax*npoin2,
     &                   1,0,mesh3d)
      ELSE
        CALL bief_allvec(1, s3d_ivide, 'IVIDE ', 0,1,0,mesh3d)
      ENDIF
!
      IF(s3d_tasse.AND.s3d_sedi) THEN
        CALL bief_allvec(1, s3d_temp, 'TEMP  ', s3d_ncouch*npoin2,
     &                   1,0,mesh3d)
      ELSE
        CALL bief_allvec(1, s3d_temp, 'TEMP  ', 0,1,0,mesh3d)
      ENDIF
!
      IF(s3d_sedi) THEN
        CALL bief_allvec(1, s3d_dmoy,  'DMOY  ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_cref,  'CREF  ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, zref,  'ZREF  ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_hdep,  'HDEP  ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_fluer, 'FLUER ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_pdepo, 'PDEPO ', ielm2h, 1, 1,mesh2d)
! CV ...
        CALL bief_allvec(1, s3d_fludp, 'FLUDP ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_fludpt,'FLUDPT', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_zf_s,  'ZF_S  ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(1, s3d_esomt, 'ESOMT ', ielm2h, 1, 1,mesh2d)
!
!
! Pour modele multi multicouche
! ?
!      IF(S3D_NCOUCH. GT.1) THEN
!         CALL BIEF_ALLVEC(1, S3D_CONC, 'CONC  ', S3D_NCOUCH, 1, 0,MESH2D)
!         CALL BIEF_ALLVEC(1, S3D_TOCE, 'TOCE  ', S3D_NCOUCH, 1, 0,MESH2D)
!      ELSE
!        CALL BIEF_ALLVEC(1, S3D_CONC, 'CONC  ',      0, 1, 0,MESH2D)
!        CALL BIEF_ALLVEC(1, S3D_TOCE, 'TOCE  ',      0, 1, 0,MESH2D)
!      ENDIF
! 7.0 Big change
!
        ALLOCATE(s3d_conc(npoin2,s3d_ncouch))
        ALLOCATE(s3d_epai(npoin2,s3d_ncouch))  ! MUD BED LAYER THICKNESS
        ALLOCATE(s3d_toce(npoin2,s3d_ncouch))  ! CONCENTRATION OF MUD BED LAYER
!
        CALL allblo(s3d_laytoce, 'LAYTOCE') ! LAYER BED SHEAR STRESS
        CALL bief_allvec_in_block(s3d_laytoce,s3d_ncouch,1,'LAYTOCE',
     &                            0,1,0,mesh2d)
        DO k=1,s3d_ncouch
          s3d_laytoce%ADR(k)%P%R=>s3d_toce(1:npoin2,k)
          s3d_laytoce%ADR(k)%P%MAXDIM1=npoin2
          s3d_laytoce%ADR(k)%P%DIM1=npoin2
        ENDDO
!
        CALL allblo(s3d_laythi, 'LAYTHI') ! LAYER BED THICKNESS
        CALL bief_allvec_in_block(s3d_laythi,s3d_ncouch,1,'LAYTHI',
     &                            0,1,0,mesh2d)
        DO k=1,s3d_ncouch
          s3d_laythi%ADR(k)%P%R=>s3d_epai(1:npoin2,k)
          s3d_laythi%ADR(k)%P%MAXDIM1=npoin2
          s3d_laythi%ADR(k)%P%DIM1=npoin2
        ENDDO
!...CV
        CALL bief_allvec(1, zr,    'ZR    ', ielm2h, 1, 1,mesh2d)
        CALL bief_allvec(2, npf,   'NPF   ', ielm2h, 1, 1,mesh2d)
      ELSE
        CALL bief_allvec(1, s3d_dmoy,  'DMOY  ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, s3d_cref,  'CREF  ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, zref,  'ZREF  ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, s3d_hdep,  'HDEP  ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, s3d_fluer, 'FLUER ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, s3d_pdepo, 'PDEPO ',      0, 1, 0,mesh2d)
        CALL bief_allvec(1, zr,    'ZR    ',      0, 1, 0,mesh2d)
        CALL bief_allvec(2, npf,   'NPF   ',      0, 1, 0,mesh2d)
      ENDIF
!
!     NEW VARIABLES FOR MIXED SEDIMENT
!
      IF(s3d_mixte) THEN
        CALL bief_allvec(1,s3d_epaico,'EPAICO',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_epainco,'EPANCO',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_pvsco,'PVSCO ',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_pvsnco,'PVSNCO',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fluerc,'FLUERC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fluernc,'FLUENC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fludptc,'FLUDTC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fludptnc,'FLUTNC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fludpc,'FLUDPC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_fludpnc,'FLUPNC',ielm2h,1,1,mesh2d)
        CALL bief_allvec(1,s3d_wcs ,'WCS   ',ielm3 ,1,1,mesh3d)
      ELSE
        CALL bief_allvec(1,s3d_epaico,'EPAICO',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_epainco,'EPANCO',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_pvsco,'PVSCO ',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_pvsnco,'PVSNCO',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fluerc,'FLUERC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fluernc,'FLUENC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fludptc,'FLUDTC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fludptnc,'FLUTNC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fludpc,'FLUDPC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_fludpnc,'FLUPNC',0     ,1,0,mesh2d)
        CALL bief_allvec(1,s3d_wcs ,'WCS   ',0     ,1,0,mesh3d)
      ENDIF
!
!-----------------------------------------------------------------------
!
!  ARRAYS USED FOR TUBES/BRIDGES
!
      IF(nbuse.GT.0) THEN
        CALL bief_allvec(1,altbus ,'ALTBUS ',nbuse,2,0,mesh3d)
        CALL bief_allvec(1,csbus  ,'CSBUS  ',nbuse,2,0,mesh3d)
        CALL bief_allvec(1,cebus  ,'CEBUS  ',nbuse,2,0,mesh3d)
        CALL bief_allvec(1,angbus ,'ANGBUS ',nbuse,2,0,mesh3d)
        CALL bief_allvec(1,lbus   ,'LBUS   ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,dbus   ,'DBUS   ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,lrgbus ,'LRGBUS ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,haubus ,'HAUBUS ',nbuse,2,0,mesh3d)
        CALL bief_allvec(1,secbus ,'SECBUS ',nbuse,1,0,mesh3d)
        CALL bief_allvec(2,clpbus ,'CLPBUS ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,cv     ,'CV     ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,c56    ,'C56    ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,cv5    ,'CV5    ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,c5     ,'C5     ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,ctrash ,'CTRASH ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,fricbus,'FRICBUS',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,longbus,'LONGBUS',nbuse,1,0,mesh3d)
        CALL bief_allvec(2,circ   ,'CIRC   ',nbuse,1,0,mesh3d)
        CALL bief_allvec(1,delbus ,'DELBUS ',nbuse,2,0,mesh3d)
        CALL bief_allvec(2,aabus  ,'AABUS  ',nbuse,1,0,mesh3d)
        DO i=1, nbuse
          dbus%R(i)=0.d0
          secbus%R(i)=1.d0
        ENDDO
        CALL allblo(tbus ,'TBUS  ')
        IF(ntrac.GT.0) THEN
          CALL bief_allvec_in_block(tbus,ntrac,1,'TBUS  ',
     &                              nbuse,2,0,mesh3d)
        ELSE
          CALL bief_allvec_in_block(tbus,1    ,1,'TBUS  ',
     &                              nbuse,2,0,mesh3d)
        ENDIF
      ELSE
        CALL bief_allvec(1,altbus ,'ALTBUS ',0,2,0,mesh3d)
        CALL bief_allvec(1,csbus  ,'CSBUS  ',0,2,0,mesh3d)
        CALL bief_allvec(1,cebus  ,'CEBUS  ',0,2,0,mesh3d)
        CALL bief_allvec(1,angbus ,'ANGBUS ',0,2,0,mesh3d)
        CALL bief_allvec(1,lbus   ,'LBUS   ',0,1,0,mesh3d)
        CALL bief_allvec(1,dbus   ,'DBUS   ',0,1,0,mesh3d)
        CALL bief_allvec(1,lrgbus ,'LRGBUS ',0,1,0,mesh3d)
        CALL bief_allvec(1,haubus ,'HAUBUS ',0,2,0,mesh3d)
        CALL bief_allvec(1,secbus ,'SECBUS ',0,1,0,mesh3d)
        CALL bief_allvec(2,clpbus ,'CLPBUS ',0,1,0,mesh3d)
        CALL bief_allvec(1,cv     ,'CV     ',0,1,0,mesh3d)
        CALL bief_allvec(1,c56    ,'C56    ',0,1,0,mesh3d)
        CALL bief_allvec(1,cv5    ,'CV5    ',0,1,0,mesh3d)
        CALL bief_allvec(1,c5     ,'C5     ',0,1,0,mesh3d)
        CALL bief_allvec(1,ctrash ,'CTRASH ',0,1,0,mesh3d)
        CALL bief_allvec(1,fricbus,'FRICBUS',0,1,0,mesh3d)
        CALL bief_allvec(1,longbus,'LONGBUS',0,1,0,mesh3d)
        CALL bief_allvec(2,circ   ,'CIRC   ',0,1,0,mesh3d)
        CALL bief_allvec(1,delbus ,'DELBUS ',0,2,0,mesh3d)
        CALL bief_allvec(2,aabus  ,'AABUS  ',0,1,0,mesh3d)
        CALL allblo(tbus ,'TBUS  ')
        IF(ntrac.GT.0) THEN
          CALL bief_allvec_in_block(tbus,ntrac,1,'TBUS  ',
     &                              0,2,0,mesh3d)
        ELSE
          CALL bief_allvec_in_block(tbus,1    ,1,'TBUS  ',
     &                              0,2,0,mesh3d)
        ENDIF
      ENDIF
!
!-----------------------------------------------------------------------
! VARSOR BLOCK FOR 2D OUTPUT COMPATIBILITY
!
      CALL allblo(varsor ,'VARSOR')
!
      CALL addblo(varsor,u2d)             ! U  01 VELOCITY ALONG X
      CALL addblo(varsor,v2d)             ! V  02 VELOCITY ALONG Y
      CALL addblo(varsor,t2_10)           ! C  03 WAVE CELERITY
      CALL addblo(varsor,h)               ! H  04 DEPTH
      CALL addblo(varsor,t2_01)           ! S  05 FREE SURFACE
      CALL addblo(varsor,zf)              ! B  06 BOTTOM
      CALL addblo(varsor,t2_02)           ! F  07 FROUDE NUMBER
      CALL addblo(varsor,t2_03)           ! Q  08 Q SCALAR FLOW RATE
      CALL addblo(varsor,svide)           ! T  09 NOT IMPLEMENTED
      CALL addblo(varsor,svide)           ! D  10 NOT IMPLEMENTED
      CALL addblo(varsor,svide)           ! E  11 NOT IMPLEMENTED
      CALL addblo(varsor,svide)           ! K  12 NOT IMPLEMENTED
      CALL addblo(varsor,t2_04)           ! I  13 HU
      CALL addblo(varsor,t2_05)           ! J  14 HV
      CALL addblo(varsor,t2_06)           ! M  15 SQRT(U**2+V**2)
      CALL addblo(varsor,windx)           ! X  16 WIND COMPONENT X DIRECT
      CALL addblo(varsor,windy)           ! Y  17 WIND COMPONENT Y DIRECT
      CALL addblo(varsor,patmos)          ! P  18 ATMOSPHERIC PRESSURE
      CALL addblo(varsor,rugof)           ! W  19 FRICTION COEFFICIENT
      CALL addblo(varsor,svide)         !!! A  20
      CALL addblo(varsor,svide)         !!! G  21
      CALL addblo(varsor,svide)         !!! L  22
      CALL addblo(varsor,zr)              ! RB 23 RIGID BED
      CALL addblo(varsor,s3d_hdep)        ! HD 24 BED THICKNESS
      CALL addblo(varsor,s3d_fluer)       ! EF 25 EROSION FLUX
      CALL addblo(varsor,s3d_fludp)       ! DF 26 DEPOSITION FLUX
      CALL addblo(varsor,prive2d%ADR(1)%P)!    27 MNEMO PRIVE1
      CALL addblo(varsor,prive2d%ADR(2)%P)!    28 MNEMO PRIVE2
      CALL addblo(varsor,prive2d%ADR(3)%P)!    29 MNEMO PRIVE3
      CALL addblo(varsor,prive2d%ADR(4)%P)!    30 MNEMO PRIVE4
      CALL addblo(varsor,t2_07)           ! US 31 FRICTION VELOCITY
      CALL addblo(varsor,t2_11)           ! QS 32 SOLID DISCHARGE
      CALL addblo(varsor,t2_12)           ! QS 33 SOLID DISCHARGE ALONG X
      CALL addblo(varsor,t2_13)           ! QS 34 SOLID DISCHARGE ALONG Y
      CALL addblo(varsor,maxz)            ! MAXZ 35 HIGH WATER MARK
      CALL addblo(varsor,tmaxz)           ! TMAXZ 36 HIGH WATER TIME
      CALL addblo(varsor,s3d_esomt)       ! S3D_ESOMT37 BED EVOLUTION DUE TO DUSPENSION
      CALL addblo(varsor,tair)            ! TAIR 38 AIR TEMPERATURE
      CALL addblo(varsor,t2_14)           ! USURF 39 VELOCITY ALONG X AT THE FREE SURFACE
      CALL addblo(varsor,t2_15)           ! VSURF 40 VELOCITY ALONG Y AT THE FREE SURFACE
      CALL addblo(varsor,t2_16)           ! WSURF 41 VELOCITY ALONG Z AT THE FREE SURFACE
      CALL addblo(varsor,t2_17)           ! MSURF 42 VELOCITY MAGNITUDE AT THE FREE SURFACE

      IF(nvar_hyd.NE.varsor%N) THEN
        WRITE(lu,*) 'MESSAGE TO DEVELOPPERS:'
        WRITE(lu,*) 'NVAR_HYD DIFFERENT THAN VARSOR SIZE ',varsor%N
        WRITE(lu,*) 'IN POINT_TELEMAC3D'
        CALL plante(1)
        stop
      ENDIF
!
!     VARIABLES ADR_TRAC_2D TO ADR_TRAC_2D+NTRAC-1
!
      IF(ntrac.GT.0) THEN
        DO i=1,ntrac
!         SIZE OF TRAV2 MUST BE GREATER THAN ADR_TRAV2+2*NTRAC
          CALL addblo(varsor,trav2%ADR(adr_trav2+i)%P)
        ENDDO
        DO i=1,ntrac
          CALL addblo(varsor,trav2%ADR(adr_trav2+ntrac+i)%P)
        ENDDO
      ENDIF
!
!     S3D_LAYTHI: FROM ADR_TRAC_2D+NTRAC TO ADR_TRAC_2D+2*NTRAC+S3D_NCOUCH
!
      IF(s3d_ncouch.GT.0.AND.s3d_sedi) THEN
        DO i=1,s3d_ncouch
          CALL addblo(varsor,s3d_laythi%ADR(i)%P)
        ENDDO
      ENDIF
!
      CALL addblo(varsor,s3d_pvsco)       ! ADR_TRAC_2D+2*NTRAC+S3D_NCOUCH
      CALL addblo(varsor,s3d_pvsnco)      ! ADR_TRAC_2D+1+2*NTRAC+S3D_NCOUCH
!
!
! QUASI - OTHER VARIABLES, AN EMPTY BLOCK
!
      CALL allblo(varcl,'VARCL ')
!
!-----------------------------------------------------------------------
! VARSOR BLOCK FOR 3D OUTPUT
!
      CALL allblo(varso3 ,'VARSO3')
! 1
      CALL addblo(varso3,z3)
! 2
      CALL addblo(varso3,u)
! 3
      CALL addblo(varso3,v)
! 4
      CALL addblo(varso3,w)
! 5
      CALL addblo(varso3,viscvi%ADR(1)%P)
! 6
      CALL addblo(varso3,viscvi%ADR(2)%P)
! 7
      CALL addblo(varso3,viscvi%ADR(3)%P)
! 8
      CALL addblo(varso3,ak)
! 9
      CALL addblo(varso3,ep)
! 10
      CALL addblo(varso3,ri)
! 11
      CALL addblo(varso3,deltar)
! 12 : DYNAMIC PRESSURE
      CALL addblo(varso3,dp)
! 13 : HYDROSTATIC PRESSURE
      CALL addblo(varso3,ph)
! 14 : ADVECTION FIELD U
      CALL addblo(varso3,uconv)
! 15 : ADVECTION FIELD V
      CALL addblo(varso3,vconv)
! 16 : ADVECTION FIELD V
      CALL addblo(varso3,wconv)
! 17 : NOT USED !!!!!!!!!!!!!
      CALL addblo(varso3,volun)
! 18 : DM1
      CALL addblo(varso3,dm1)
! 19 : DH AND HN
      CALL addblo(varso3,dhhn)
! 20 : UCONVC
      CALL addblo(varso3,uconvc)
! 21 : VCONVC
      CALL addblo(varso3,vconvc)
! 22 : UD
      CALL addblo(varso3,ud)
! 23 : VD
      CALL addblo(varso3,vd)
! 24 : WD
      CALL addblo(varso3,wd)
! 25 : PRIVE 1
      CALL addblo(varso3,prive%ADR(1)%P)
! 26 : PRIVE 2
      CALL addblo(varso3,prive%ADR(2)%P)
! 27 : PRIVE 3
      CALL addblo(varso3,prive%ADR(3)%P)
! 28 : PRIVE 4
      CALL addblo(varso3,prive%ADR(4)%P)
! 29 : US
      CALL addblo(varso3,ustokes)
! 30 : VS
      CALL addblo(varso3,vstokes)
! 31 : WS
      CALL addblo(varso3,wstokes)
!

!     THE FOLLOWING ADDRESS FOR TRACERS IS ADR_TRAC = 32
!     THIS NUMBER IS HARDCODED IN DECLARATIONS_TELEMAC3D
!
      IF(nvar_t3d.NE.varso3%N) THEN
        WRITE(lu,*) 'MESSAGE TO DEVELOPPERS:'
        WRITE(lu,*) 'NVAR_T3D DIFFERENT THAN VARSO3 SIZE ',varso3%N
        WRITE(lu,*) 'IN POINT_TELEMAC3D'
        CALL plante(1)
        stop
      ENDIF
!
! VARIABLES NUMBER "ADR_TRAC" TO "ADR_TRAC" + NTRAC - 1
!
      IF(ntrac.GT.0) THEN
        DO i=1,ntrac
          CALL addblo(varso3,ta%ADR(i)%P)
        ENDDO
      ENDIF
!
! VARIABLES NUMBER "ADR_TRAC" + NTRAC TO "ADR_TRAC" + 4*NTRAC - 1
!
      IF(ntrac.GT.0) THEN
        DO i=1,ntrac
          CALL addblo(varso3,viscta%ADR(i)%P%ADR(1)%P)
          CALL addblo(varso3,viscta%ADR(i)%P%ADR(2)%P)
          CALL addblo(varso3,viscta%ADR(i)%P%ADR(3)%P)
        ENDDO
      ENDIF
!
      IF(inclus(coupling,'TOMAWACT3D')) THEN
        CALL bief_allvec(1, wip,     'WIP   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, wipdx,   'WIPDX ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, wipdy,   'WIPDY ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, ustokes, 'USTOKE', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, vstokes, 'VSTOKE', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, wstokes, 'WSTOKE', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, wst1,    'WST1  ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, wst2,    'WST2  ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, us2d,    'US2D  ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, vs2d,    'VS2D  ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, dus2d,   'DUS2D ', ielmh, 1, 2,mesh2d)
        CALL bief_allvec(1, dvs2d,   'DVS2D ', ielmh, 1, 2,mesh2d)
        CALL bief_allvec(1, dus2dx,  'DUS2DX', ielmh, 1, 2,mesh2d)
        CALL bief_allvec(1, dvs2dy,  'DVS2DY', ielmh, 1, 2,mesh2d)
        CALL bief_allvec(1, ustx,    'USTX  ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, vsty,    'VSTY  ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, fdx,     'FDX   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, fdy,     'FDY   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, fbx,     'FBX   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, fby,     'FBY   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, cfwc,    'CFWC  ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, fdk,     'FDK   ', ielm3, 1, 1,mesh3d)
        CALL bief_allvec(1, fwx,     'FWX   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, fwy,     'FWY   ', ielmu, 1, 1,mesh2d)
        CALL bief_allvec(1, hborcli, 'HBORCL', ielm1, 1, 1,mesh2d)
        CALL bief_allvec(1, uborcli, 'UBORCL', ielm1, 1, 1,mesh2d)
        CALL bief_allvec(1, vborcli, 'VBORCL', ielm1, 1, 1,mesh2d)
      ELSE
        CALL bief_allvec(1, wip,     'WIP   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, wipdx,   'WIPDX ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, wipdy,   'WIPDY ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, ustokes, 'USTOKE', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, vstokes, 'VSTOKE', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, wstokes, 'WSTOKE', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, wst1,    'WST1  ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, wst2,    'WST2  ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, us2d,    'US2D  ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, vs2d,    'VS2D  ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, dus2d,   'DUS2D ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, dvs2d,   'DVS2D ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, dus2dx,  'DUS2DX', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, dvs2dy,  'DVS2DY', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, ustx,    'USTX  ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, vsty,    'VSTY  ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, fdx,     'FDX   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, fdy,     'FDY   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, fbx,     'FBX   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, fby,     'FBY   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, cfwc,    'CFWC  ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, fdk,     'FDK   ', 0, 1, 0,mesh3d)
        CALL bief_allvec(1, fwx,     'FWX   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, fwy,     'FWY   ', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, hborcli, 'HBORCL', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, uborcli, 'UBORCL', 0, 1, 0,mesh2d)
        CALL bief_allvec(1, vborcli, 'VBORCL', 0, 1, 0,mesh2d)
      ENDIF
!
!     METEO MODULE
!
      CALL point_meteo(t3d_files,t3atma,t3atmb,mesh2d,ielm2h,vent,atmos,
     &                 inclus(coupling,'WAQTEL'))
!
!-----------------------------------------------------------------------
!
      IF(listin) THEN
        WRITE(lu,42)
      ENDIF
 42   FORMAT(1x,/,' END OF MEMORY ALLOCATION  ',/)
!
!-----------------------------------------------------------------------
!
      RETURN
      END