!                    *****************
                     SUBROUTINE ANAVEN
!                    *****************
!
     &(UV,VV,X,Y,NPOIN2,AT,DDC,VX_CTE,VY_CTE,WIND)
!
!***********************************************************************
! TOMAWAC   V6P3                                   08/06/2011
!***********************************************************************
!
!brief    SPECIFIES AN ANALYTICAL WIND
!+               (CAN BE VARIABLE IN TIME).
!
!warning  USER SUBROUTINE; MUST BE CODED BY THE USER; THIS IS MERELY AN EXAMPLE
!code
!+      DO IP=1,NPOIN2
!+        UV(IP)=VX_CTE
!+        VV(IP)=VY_CTE
!+      ENDDO
!
!history  M. BENOIT (LNH)
!+        07/06/95
!+        V1P0
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  G.MATTAROLO (EDF - LNHE)
!+        08/06/2011
!+        V6P1
!+   Translation of French names of the variables in argument
!
!history  J-M HERVOUET (EDF R&D, LNHE)
!+        26/11/2013
!+        V6P3
!+   Example added to copy the wind in Telemac-2D memory.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| AT             |-->| COMPUTATION TIME
!| DDC            |-->| DATE OF COMPUTATION BEGINNING
!| NPOIN2         |-->| NUMBER OF POINTS IN 2D
!| UV             |<--| WIND VELOCITY ALONG X AT THE MESH POINTS
!| VV             |<--| WIND VELOCITY ALONG Y AT THE MESH POINTS
!| VX_CTE         |-->| WIND ALONG X (CONSTANT VALUE IN STEERING FILE)
!| VY_CTE         |-->| WIND ALONG Y (CONSTANT VALUE IN STEERING FILE)
!| X              |-->| ABSCISSAE OF POINTS IN THE MESH
!| Y              |-->| ORDINATES OF POINTS IN THE MESH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      
      USE DECLARATIONS_SPECIAL
      USE INTERFACE_TOMAWAC, EX_ANAVEN => ANAVEN
      USE DECLARATIONS_TELEMAC3D, ONLY : FUAIR,FVAIR
!
!     WHEN COUPLING WITH TELEMAC-2D, THIS WILL ALLOW A COPY OF WIND
!     IN TELEMAC-2D
!
      !USE DECLARATIONS_TELEMAC2D, ONLY : WINDX,WINDY
      USE BIEF
      IMPLICIT NONE
!
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)             ::  NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: AT    , DDC   , VX_CTE, VY_CTE
      DOUBLE PRECISION, INTENT(IN)    :: X (NPOIN2)    , Y (NPOIN2)
      DOUBLE PRECISION, INTENT(INOUT)    :: UV(NPOIN2)    , VV(NPOIN2)
       !&, WIND
      TYPE (BIEF_OBJ), INTENT(INOUT)    :: wind
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER  IP
      !DOUBLE PRECISION UCONST, VCONST
!
!      UCONST= 20.D0
!      VCONST= 0.D0
!!
      !DO IP=1,NPOIN2
      !  UV(IP)=UCONST
      !  VV(IP)=VCONST
      !ENDDO
!
!     EXAMPLE IF COUPLING WITH TELEMAC-2D
!     READING TOMAWAC WIND IN TELEMAC-2D SHOULD BE POSSIBLE ALSO
!
      !DO IP=1,NPOIN2
       !UV(IP)=WINDX%R(IP)
       !VV(IP)=WINDY%R(IP)
      !ENDDO
      DO IP=1,NPOIN2
       UV(IP)=WIND%ADR(1)%P%R(IP)
       VV(IP)=WIND%ADR(2)%P%R(IP)
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END