interface_tomawac.f

Go to the documentation of this file.

!
!  INTERFACES OF TOMAWAC SUBROUTINES
!
      MODULE interface_tomawac
!
!=======================================================================
!
!     DEFINITION OF STRUCTURES
!
      USE bief_def
!
!-----------------------------------------------------------------------
!
!     DEFINITION OF INTERFACES
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE anacos
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE anamar
      IMPLICIT NONE

        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE anaven
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE angles(XLAMD , DTPLUS, DTMOIN)
      IMPLICIT NONE
      DOUBLE PRECISION, INTENT(IN)    :: XLAMD
      DOUBLE PRECISION, INTENT(INOUT) :: DTPLUS, DTMOIN
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE bornes(B , N , A , XM , X0 , X1 )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             ::  N
      DOUBLE PRECISION, INTENT(IN)    :: B , A , XM
      DOUBLE PRECISION, INTENT(INOUT) :: X0 , X1
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE condiw ( PART, U_TEL, V_TEL, H_TEL )
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: PART
      TYPE(bief_obj), INTENT(IN)      :: U_TEL,V_TEL,H_TEL
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE conw4d
     &(cx, cy, ct, cf, xk, cg, npoin2, ndire, jf, nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2,JF
      DOUBLE PRECISION, INTENT(INOUT) :: CX(npoin2,ndire,jf)
      DOUBLE PRECISION, INTENT(INOUT) :: CY(npoin2,ndire,jf)
      DOUBLE PRECISION, INTENT(INOUT) :: CT(npoin2,ndire,jf)
      DOUBLE PRECISION, INTENT(INOUT) :: CF(npoin2,ndire,jf)
      DOUBLE PRECISION, INTENT(IN)    :: CG(npoin2,nf),XK(npoin2,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE conwac
     &( cx    , cy    , ct    , xk    , cg    , npoin2, ndire , jf    ,
     &  nf    )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2,JF
      DOUBLE PRECISION, INTENT(IN)    :: CG(npoin2,nf),XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: CX(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: CY(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: CT(npoin2,ndire)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE cormar( PART, U_TEL, V_TEL, H_TEL )
      USE bief_def
      IMPLICIT NONE
      INTEGER,          INTENT(IN) :: PART
      TYPE(bief_obj),   INTENT(IN) :: U_TEL,V_TEL,H_TEL
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION couple
     &( xk1   , yk1   , xk2   , yk2   , xk3   , yk3   , xk4   , yk4   )
      IMPLICIT NONE
      DOUBLE PRECISION , INTENT(IN)    :: XK1   , YK1   , XK2   , YK2
      DOUBLE PRECISION , INTENT(IN)    :: XK3   , YK3   , XK4   , YK4
      DOUBLE PRECISION COUPLE
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE couuti
     &(ncou,fmtcou)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NCOU
      CHARACTER(LEN=8), INTENT(IN)    :: FMTCOU
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_current
     &(ncou,fmtcou)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NCOU
      CHARACTER(LEN=8), INTENT(IN)    :: FMTCOU
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE cqueue
     &( jfre  , jbis  , coef1 )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: JFRE
      DOUBLE PRECISION, INTENT(INOUT) :: COEF1
      INTEGER, INTENT(INOUT)          :: JBIS
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION delfra( SS    )
      IMPLICIT NONE
      DOUBLE PRECISION DELFRA, SS
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE diffrac
     &( cx    , cy    , ct    , xk    , cg    , npoin2,
     &  ndire , iff   , nf    , f     , rx    , ry    ,
     &  rxx   , ryy   , neigb )
        USE bief_def
        USE declarations_tomawac, ONLY : maxnsp
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NF, NDIRE, NPOIN2, IFF
      INTEGER, INTENT(IN) :: NEIGB(npoin2,maxnsp)
      DOUBLE PRECISION, INTENT(INOUT) :: CX(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: CY(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: CT(npoin2,ndire)
      DOUBLE PRECISION, INTENT(IN)    :: CG(npoin2,nf),XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: RX(maxnsp,npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: RY(maxnsp,npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: RXX(maxnsp,npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: RYY(maxnsp,npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE dump2d(  XF1 , NP1 )
      IMPLICIT NONE
      INTEGER, INTENT(IN)          :: NP1
      DOUBLE PRECISION, INTENT(IN) :: XF1(np1)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE ecrete
     &( f     , depth , npoin2, ndire , nf    , promin)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPOIN2 , NDIRE, NF
      DOUBLE PRECISION, INTENT(IN)    :: DEPTH(npoin2),PROMIN
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE ecrspe
     &( f     , ndire , nf    , npoin2, lt, auxil , noleo ,
     &  nleo  , debres,  date , time  , knolg , mesh )
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPOIN2,NLEO,NF,NDIRE,LT
      INTEGER, INTENT(IN)             :: KNOLG(*),NOLEO(nleo)
      INTEGER, INTENT(IN)             :: DATE(3),TIME(3)
      DOUBLE PRECISION, INTENT(INOUT) :: AUXIL(ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      LOGICAL, INTENT(IN)             :: DEBRES
      TYPE(bief_mesh), INTENT(INOUT)  :: MESH
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE f1f1f1
     &( f1sf  , nf1   , iq_om1)
      IMPLICIT NONE
      INTEGER,          INTENT(IN)    :: IQ_OM1
      INTEGER,          INTENT(INOUT) :: NF1
      DOUBLE PRECISION, INTENT(INOUT) :: F1SF(*)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fbott3d
     &  (fbx, fby, fs, npoin2, xk, ndire, nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: NPOIN2, NDIRE, NF
      DOUBLE PRECISION, INTENT(IN)    :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT)    :: FBX(npoin2)
      DOUBLE PRECISION, INTENT(INOUT)    :: FBY(npoin2)  
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fdiss3d
     &  (fdx, fdy, npoin2, xk, ndire, fs, nf)
      IMPLICIT NONE
        INTEGER, INTENT(IN) :: NPOIN2, NDIRE, NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2, ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2, nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FDX(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: FDY(npoin2)  
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
          SUBROUTINE fdissk    
     &        (fdk, npoin2, ndire, fs, ztel, nz, hsmjt, fznorm, nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN)  :: NZ  
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE, NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2, ndire, nf)
      DOUBLE PRECISION, INTENT(IN) :: HSMJT(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: FDK(npoin2, nz)
      DOUBLE PRECISION, INTENT(IN) :: ZTEL(npoin2, nz)
      DOUBLE PRECISION, INTENT(INOUT) :: FZNORM(npoin2)
          END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE filt_sa
        IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION foncro
     &( x     , b     , n     , a     , xm    )
      IMPLICIT NONE
      INTEGER,          INTENT(IN)    :: N
      DOUBLE PRECISION, INTENT(IN)    :: X      , B     , A     , XM
      DOUBLE PRECISION FONCRO
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fpread( FREAD, F, NF, NDIRE, NPOIN2, EXPO )
      IMPLICIT NONE
      INTEGER,          INTENT(IN) :: NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN) :: EXPO  , F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FREAD(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE frem01( FM01 , F, NF, NDIRE, NPOIN2)
      IMPLICIT NONE
      INTEGER,          INTENT(IN) ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FM01(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE frem02( FM02, F, NF, NDIRE, NPOIN2)
      IMPLICIT NONE
      INTEGER,          INTENT(IN) ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FM02(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fremoy( FMOY, F, NF, NDIRE, NPOIN2)
      IMPLICIT NONE
      INTEGER,         INTENT(IN) ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION,INTENT(IN) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION,INTENT(INOUT) :: FMOY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE frepic( FPIC, F, NF, NDIRE, NPOIN2)
      IMPLICIT NONE
      INTEGER,INTENT(IN)             :: NF    , NDIRE , NPOIN2
      DOUBLE PRECISION,INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION,INTENT(INOUT) :: FPIC(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fric3d
     &  (cfwc, npoin2, dirhou, u_tel, v_tel,  uwbm)
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2
      TYPE(bief_obj),   INTENT(IN) :: U_TEL,V_TEL
      TYPE(bief_obj),   INTENT(INOUT)   :: CFWC
      DOUBLE PRECISION, INTENT(IN) :: UWBM(npoin2)
      DOUBLE PRECISION, INTENT(IN) :: DIRHOU(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P2 New subroutine
      INTERFACE
        SUBROUTINE frmset ( NEIGB , NPOIN2, NELEM2,
     &  IKLE  , RK    , RX    , RY    , RXX   , RYY )
        USE declarations_tomawac, ONLY : maxnsp
        IMPLICIT NONE
!      INTEGER,INTENT(IN)            :: NPOIN2, MAXNSP, NELEM2
      INTEGER,INTENT(IN)            :: NPOIN2, NELEM2
      INTEGER,INTENT(INOUT)         :: NEIGB(npoin2,maxnsp)
      INTEGER,INTENT(IN)            :: IKLE(nelem2,3)
      DOUBLE PRECISION,INTENT(INOUT):: RK(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(INOUT):: RX(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(INOUT):: RY(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(INOUT):: RXX(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(INOUT):: RYY(maxnsp,npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fsprd1
     &( fra   , ndire , spred1, teta1 , spred2, teta2 , xlamda)
      IMPLICIT NONE
      INTEGER,INTENT(IN)             :: NDIRE
      DOUBLE PRECISION,INTENT(IN)    :: SPRED1, TETA1 , SPRED2, TETA2
      DOUBLE PRECISION,INTENT(IN)    :: XLAMDA
      DOUBLE PRECISION,INTENT(INOUT) :: FRA(ndire)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fsprd2
     &( fra   , ndire , spred1, teta1 , spred2, teta2 , xlamda)
      IMPLICIT NONE
      INTEGER,INTENT(IN)             :: NDIRE
      DOUBLE PRECISION,INTENT(IN)    :: SPRED1, TETA1 , SPRED2, TETA2
      DOUBLE PRECISION,INTENT(IN)    :: XLAMDA
      DOUBLE PRECISION,INTENT(INOUT) :: FRA(ndire)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE fsprd3
     &( fra   , ndire , spred1, teta1 , spred2, teta2 , xlamda)
      IMPLICIT NONE
      INTEGER,INTENT(IN)             :: NDIRE
      DOUBLE PRECISION,INTENT(IN)    :: SPRED1, TETA1 , SPRED2, TETA2
      DOUBLE PRECISION,INTENT(IN)    :: XLAMDA
      DOUBLE PRECISION,INTENT(INOUT) :: FRA(ndire)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION gammln
     &( xx    , deupi )
      IMPLICIT NONE
      DOUBLE PRECISION GAMMLN
      DOUBLE PRECISION,INTENT(IN)    :: XX    , DEUPI
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE gauleg
     &( w_leg , x_leg , npoin )
      IMPLICIT NONE
      INTEGER ,INTENT(IN)             :: NPOIN
      DOUBLE PRECISION ,INTENT(INOUT) :: W_LEG(npoin) , X_LEG(npoin)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE ifabtom(IFABOR,NELEM2,NETAGE)
      IMPLICIT NONE
      INTEGER, INTENT(IN)          :: NELEM2,NETAGE
      INTEGER, INTENT(INOUT)       :: IFABOR(nelem2,5,netage)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE impr(LISPRD,LT,AT,ISITS,ICOD)
      IMPLICIT NONE
      INTEGER,INTENT(IN)           :: LT,ICOD,LISPRD,ISITS
      DOUBLE PRECISION,INTENT(IN)  :: AT
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE iniphy
     &( xk    , cg    , b     ,  npoin2, nf )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF    , NPOIN2
      DOUBLE PRECISION, INTENT(INOUT) :: B(npoin2,nf)  , XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: CG(npoin2,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE initab(IBOR1,IFABOR1,NELEM2_DIM,PART)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: PART,NELEM2_DIM
      INTEGER, INTENT(IN)    :: IFABOR1(nelem2_dim,3)
      INTEGER, INTENT(INOUT) :: IBOR1(nelem2_dim,7)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE intang
     &( lavant, lapres, idire , ndire , deltad)
      IMPLICIT NONE
      INTEGER, INTENT(IN)          :: NDIRE , IDIRE
      DOUBLE PRECISION, INTENT(IN) :: DELTAD
      INTEGER, INTENT(INOUT)       :: LAVANT, LAPRES
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE invert( RN    , N     , NP    )
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: N, NP
      DOUBLE PRECISION, INTENT(INOUT) :: RN(np,np)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        DOUBLE PRECISION FUNCTION kerbou
     &(xk1,xk2,freq1,freq2,depth,teta1,teta2)
      IMPLICIT NONE
      DOUBLE PRECISION, INTENT(IN) :: XK1,XK2,FREQ1,FREQ2,TETA1,TETA2
      DOUBLE PRECISION, INTENT(IN) :: DEPTH
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE kmoye2
     &( xkmoy , xk    , f     , nf    , ndire , npoin2, 
     &  aux1  , aux2  , aux3  )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: AUX1(npoin2),AUX2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: AUX3(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: XKMOY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE kmoyen
     &( xkmoy , xk    , f     , nf    , ndire ,
     &  npoin2, aux1  , aux2  , aux3  )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: AUX1(npoin2),AUX2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: AUX3(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: XKMOY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE lecdoi
     &(f1,name1,mode1,f2,name2,mode2, f3,name3,mode3,
     & npoin2,ndon,fmtdon,
     & at,tv1,tv2,f11,f12,f21,f22,f31,f32,indic,chdon,nvar,
     & texte,trouve,unitime,phastime)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NDON,NPOIN2,INDIC
      INTEGER, INTENT(IN)             :: MODE1,MODE2,MODE3
      INTEGER, INTENT(INOUT)          :: NVAR
      DOUBLE PRECISION, INTENT(INOUT) :: F1(npoin2),F2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F3(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F11(npoin2),F12(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F21(npoin2),F22(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F31(npoin2),F32(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: AT,UNITIME,PHASTIME
      DOUBLE PRECISION, INTENT(INOUT) :: TV1,TV2
      CHARACTER(LEN=8), INTENT(IN)    :: FMTDON
      CHARACTER(LEN=7), INTENT(IN)    :: CHDON
      CHARACTER(LEN=32),INTENT(IN)    :: NAME1,NAME2,NAME3
      CHARACTER(LEN=32),INTENT(INOUT) :: TEXTE(30)
      LOGICAL, INTENT(INOUT)          :: TROUVE(3)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE lecdon
     &(f1,name1,mode1,f2,name2,mode2, f3,name3,mode3,
     & npoin2,ndon,fformat,indic,chdon,texte,trouve)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NDON,NPOIN2,INDIC
      INTEGER, INTENT(IN)             :: MODE1,MODE2,MODE3
      DOUBLE PRECISION, INTENT(INOUT) :: F1(npoin2),F2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F3(npoin2)
      CHARACTER(LEN=8), INTENT(IN)    :: FFORMAT
      CHARACTER(LEN=7), INTENT(IN)    :: CHDON
      CHARACTER(LEN=32),INTENT(IN)    :: NAME1,NAME2,NAME3
      CHARACTER(LEN=32),INTENT(INOUT) :: TEXTE(30)
      LOGICAL, INTENT(INOUT)          :: TROUVE(3)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE lecdon_tomawac(FILE_DESC,PATH,NCAR,
     &                            CAS_FILE,DICO_FILE)
        USE declarations_special
      IMPLICIT NONE
      CHARACTER(LEN=PATH_LEN), INTENT(INOUT) :: FILE_DESC(4,maxkeyword)
      INTEGER, INTENT(IN)               :: NCAR
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: PATH
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: DICO_FILE
      CHARACTER(LEN=PATH_LEN), INTENT(IN)    :: CAS_FILE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE lecsui
     &(f,ndire,nf,npoin2,vent,couran,npre,fmtpre, maree,tra01)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPRE,NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TRA01(npoin2*ndire)
      LOGICAL, INTENT(IN)             :: COURAN,VENT,MAREE
      CHARACTER(LEN=8), INTENT(IN)    :: FMTPRE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE limite( F , FREQ  , NPOIN2, NDIRE , NF    )
        IMPLICIT NONE
        INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
        DOUBLE PRECISION, INTENT(IN)    :: FREQ(nf)
        DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE limwac
     &(f     , fbor  , nptfr , ndire , nf    , npoin2,
     & kent  , prive , npriv , imp_file)
      USE bief_def, ONLY : bief_file
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPTFR,NDIRE,NF,NPOIN2,NPRIV
      INTEGER, INTENT(IN)             :: KENT
      DOUBLE PRECISION, INTENT(IN):: PRIVE(npoin2,npriv)
      TYPE(bief_file), INTENT(IN)    :: IMP_FILE
      DOUBLE PRECISION, INTENT(INOUT):: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT):: FBOR(nptfr,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE maruti
     &(nmar,fmtmar)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NMAR
      CHARACTER(LEN=8), INTENT(IN)    :: FMTMAR
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE moudiss1
     &  (fwx, fwy, npoin2, xk, ndire, fs,nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FWX(npoin2), FWY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE moudiss2
     &  (fwx, fwy, npoin2, xk, ndire, fs,nf, taux1, f_int)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FWX(npoin2), FWY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F_INT(npoin2),TAUX1(npoin2)
      END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_tide
     &(nmar,fmtmar)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NMAR
      CHARACTER(LEN=8), INTENT(IN)    :: FMTMAR
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE nomvar_tomawac
     &(texte,mnemo,maxvar)
      IMPLICIT NONE
      INTEGER, INTENT(IN)         ::      MAXVAR
      CHARACTER(LEN=32), INTENT(INOUT) :: TEXTE(maxvar)
      CHARACTER(LEN=8) , INTENT(INOUT) :: MNEMO(maxvar)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE noudon
     &(f1,name1,mode1,f2,name2,mode2,f3,name3,mode3,
     & npoin,ndon,fformat,
     & at,tv1,tv2,f11,f12,f21,f22,f31,f32,indic,chdon,nvar,texte,
     & trouve,unitime,phastime)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NDON,NPOIN,INDIC
      INTEGER, INTENT(INOUT)          :: NVAR
      INTEGER, INTENT(IN)             :: MODE1,MODE2,MODE3
      DOUBLE PRECISION, INTENT(INOUT) :: F1(npoin),F2(npoin),F3(npoin)
      DOUBLE PRECISION, INTENT(INOUT) :: F11(npoin),F21(npoin)
      DOUBLE PRECISION, INTENT(INOUT) :: F12(npoin),F22(npoin)
      DOUBLE PRECISION, INTENT(INOUT) :: F31(npoin),F32(npoin)
      DOUBLE PRECISION, INTENT(IN)    :: AT,UNITIME,PHASTIME
      DOUBLE PRECISION, INTENT(INOUT) :: TV1,TV2
      CHARACTER(LEN=8), INTENT(IN)    :: FFORMAT
      CHARACTER(LEN=7), INTENT(IN)    :: CHDON
      CHARACTER(LEN=32),INTENT(IN)    :: NAME1,NAME2,NAME3
      CHARACTER(LEN=32),INTENT(IN)    :: TEXTE(30)
      LOGICAL, INTENT(INOUT)          :: TROUVE(3)
!
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE predif
     &(cx    , cy    , ikle2 , ifabor, elt   , eta   ,
     & xk    , cg    , itr01 , npoin3, npoin2, nelem2, ndire ,
     & nf    , couran, f     , rx    , ry    , rxx   , ryy   , neigb )
      USE bief_def
        USE declarations_tomawac, ONLY : maxnsp
      IMPLICIT NONE

      INTEGER,INTENT(IN)    :: NPOIN3,NPOIN2,NELEM2,NDIRE,NF
      INTEGER,INTENT(IN)    :: IKLE2(nelem2,3)
      INTEGER,INTENT(IN)    :: NEIGB(npoin2,maxnsp)
      INTEGER,INTENT(INOUT) :: IFABOR(nelem2,7)
      INTEGER,INTENT(INOUT) :: ELT(npoin3,nf), ETA(npoin3,nf)
      INTEGER,INTENT(INOUT) :: ITR01(npoin3,3)
      DOUBLE PRECISION,INTENT(IN) :: RX(maxnsp,npoin2),RY(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(IN) :: RXX(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(IN) :: RYY(maxnsp,npoin2)
      DOUBLE PRECISION,INTENT(IN) :: XK(npoin2,nf),CG(npoin2,nf)
      DOUBLE PRECISION,INTENT(IN) :: F(npoin2,ndire,nf)
      LOGICAL,INTENT(IN) :: COURAN
      TYPE(bief_obj), INTENT(INOUT) :: CX,CY
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE prenl1
     &( iangnl, coefnl, ndire , nf    , raisf , xlamd )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NDIRE , NF
      INTEGER, INTENT(INOUT)          :: IANGNL(ndire,8)
      DOUBLE PRECISION, INTENT(IN)    :: RAISF , XLAMD
      DOUBLE PRECISION, INTENT(INOUT) :: COEFNL(16)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P1 New subroutine
      INTERFACE
        SUBROUTINE prenl2
     &( iangnl, coefnl, ndire , nf    , raisf , xlamd , xmu   )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NDIRE , NF
      INTEGER, INTENT(INOUT)          :: IANGNL(ndire,16)
      DOUBLE PRECISION, INTENT(IN)    :: RAISF , XLAMD , XMU
      DOUBLE PRECISION, INTENT(INOUT) :: COEFNL(32)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P1 New subroutine
      INTERFACE
        SUBROUTINE prenl3
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE prepro
     &( cx    , cy    , ikle2 , ifabor, elt   , eta   , fre   ,
     &  xk    , cg    , itr01 , npoin3, npoin2, nelem2, ndire ,
     &  nf    , couran)
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: NPOIN3,NPOIN2,NELEM2,NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf),CG(npoin2,nf)
      INTEGER, INTENT(INOUT) :: ELT(npoin3,nf),ETA(npoin3,nf), FRE(*)
      INTEGER, INTENT(IN)    :: IKLE2(nelem2,3)
      INTEGER, INTENT(INOUT) :: ITR01(npoin3,3),IFABOR(nelem2,7)
      LOGICAL, INTENT(IN)    :: COURAN
      TYPE(bief_obj), INTENT(INOUT) :: CX,CY
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE preqt2
        IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE propa
     &(f, b, elt, eta, fre, npoin3, npoin2,
     & ndire, nf, couran, tra01)
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPOIN3,NPOIN2,NDIRE,NF
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: B(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TRA01(npoin3,8)
      INTEGER, INTENT(INOUT) :: ELT(npoin3,nf),ETA(npoin3,nf)
      INTEGER, INTENT(INOUT) :: FRE(*)
      LOGICAL, INTENT(IN)    :: COURAN
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION qbbj78( B     , IQBBJ )
      IMPLICIT NONE
      DOUBLE PRECISION QBBJ78
      DOUBLE PRECISION, INTENT(IN) ::  B
      INTEGER, INTENT(IN)     ::  IQBBJ
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qbrek1
     &( tstot , f     , fcar  , varian, nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)            :: NF, NDIRE, NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)   :: VARIAN(npoin2),FCAR(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qbrek2
     &( tstot , f     , fcar  , varian, nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)            :: NF, NDIRE, NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)   :: VARIAN(npoin2), FCAR(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qbrek3
     & ( tstot , f     , fcar  , varian, nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)   ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)   :: VARIAN(npoin2), FCAR(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qbrek4
     & ( tstot , f     , fcar  , varian, nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)   ::          NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)   :: VARIAN(npoin2),FCAR(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qdscur
     &( tstot , tsder , f     , cf    , xk    , usold , usnew ,
     &  nf    , ndire , npoin2, f_int , betoto, betotn)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: USNEW(npoin2),USOLD(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: F_INT(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: BETOTO(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: BETOTN(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: CF(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qfrot1
     &( tstot , tsder , f     , xk    ,  nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        FUNCTION qgauss
     &( b     , n     , a     , xm    )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  N
      DOUBLE PRECISION QGAUSS
      DOUBLE PRECISION, INTENT(IN)    :: B , A , XM
        END FUNCTION
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qmout1
     &( tstot , tsder , f     , xk    , enrj  , fmoy  , xkmoy ,
     &  nf  , ndire , npoin2, taux1   )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: XKMOY(npoin2),ENRJ(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX1(npoin2),FMOY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qmout2
     &( tstot , tsder , f     , xk    , enrj  , fmoy  , xkmoy , usold ,
     &  usnew , nf    , ndire , npoin2, taux1 ,f_int  , betoto, betotn)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: USNEW(npoin2),USOLD(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: FMOY(npoin2),XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)    :: ENRJ(npoin2),XKMOY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F_INT(npoin2),TAUX1(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: BETOTO(npoin2),BETOTN(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qnlin1
     &( tstot , tsder , iangnl, nf    , ndire , npoin2, f     ,
     &  xkmoy , taux1 , taux2 , taux3 , taux4 , taux5 , dfini )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPOIN2,NDIRE,NF
      INTEGER, INTENT(IN)             :: IANGNL(ndire,8)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XKMOY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX1(npoin2),TAUX2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX3(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX4(npoin2),TAUX5(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: DFINI(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P1 New subroutine
      INTERFACE
        SUBROUTINE qnlin2
     &( tstot , tsder , iangnl,  coefnl, nf    , ndire , npoin2, f     ,
     &  xkmoy , taux1 , dfini , xcoef )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NPOIN2, NDIRE , NF
      INTEGER, INTENT(IN)             :: IANGNL(ndire,16)
      DOUBLE PRECISION, INTENT(IN)    :: XCOEF, COEFNL(32)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XKMOY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX1(npoin2), DFINI(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P1 New subroutine
      INTERFACE
        SUBROUTINE qnlin3( F1, T1TOT,T1DER, N1POIN2,N1PLAN,N1F )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: N1POIN2,N1PLAN,N1F
      DOUBLE PRECISION, INTENT(INOUT) :: T1TOT( n1poin2,n1plan,n1f)
      DOUBLE PRECISION, INTENT(INOUT) :: T1DER( n1poin2,n1plan,n1f)
      DOUBLE PRECISION, INTENT(INOUT) :: F1(n1poin2,n1plan,n1f)
        END SUBROUTINE
      END INTERFACE
      INTERFACE
!
!-----------------------------------------------------------------------
!
        SUBROUTINE qporos
     &( tstot , tsder , f , cg, lt,  xk,
     &  nf    , ndire  , npoin2   , amorp  )
!
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2,LT
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: CG(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: AMORP(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qtria1
     &( f     , xk    , nf    , ndire , npoin2, tstot , ftot  , fmoy  )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF, NDIRE, NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)    :: FTOT(npoin2) , FMOY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qtria2
     &( f     , xk    , nf    , ndire , npoin2, tstot )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF, NDIRE, NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!=======================================================================
!VBA AJOUT SUBROUTINE QVEG1 12/09/2014
!=======================================================================
!
      INTERFACE
        SUBROUTINE qveg
     &( tstot , tsder , f , varian , fmoy , xkmoy , nf    , ndire  ,
     &  npoin2   )
!
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: XKMOY(npoin2),VARIAN(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf),FMOY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!=========================================================================
!VBA FIN AJOUT SUBROUTINE QVEG1
!=========================================================================
!
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qwind1
     &( tstot , tsder , f     , xk    , usold , usnew , twold , twnew ,
     &  z0old , z0new , nf    , ndire , npoin2, told  , tnew  ,
     &  usn   , uso   , omnew , omold , betan , betao )
      IMPLICIT NONE
      INTEGER, INTENT(IN)   ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf),XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN)   :: TWOLD(npoin2), TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)   :: USOLD(npoin2), USNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)   :: Z0OLD(npoin2), Z0NEW(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TNEW(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT):: TOLD(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT):: USO(npoin2),U SN(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: OMNEW(npoin2),OMOLD(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: BETAN(npoin2), BETAO(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT):: TSDER(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qwind2
     &( tstot , tsder , f     , xk    , usold , usnew , twold , twnew ,
     &  nf    , ndire , npoin2, usn   , uso   )
      IMPLICIT NONE
      INTEGER, INTENT(IN)            :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: TWOLD(npoin2),TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)   :: USNEW(npoin2),USOLD(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: USO(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT):: USN(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT):: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf),XK(npoin2,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qwind3
     &( tstot , tsder , f     , xk    , usold , usnew , twold , twnew ,
     &  nf    , ndire , npoin2, betan , betao , dirn  , diro  )
      IMPLICIT NONE
      INTEGER, INTENT(IN)   ::    NF  , NDIRE        , NPOIN2
      DOUBLE PRECISION, INTENT(IN)   :: TWOLD(npoin2), TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)   :: USOLD(npoin2), USNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)   :: F(npoin2,ndire,nf),XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT):: DIRO(npoin2) , DIRN(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: BETAN(npoin2), BETAO(npoin2)
      DOUBLE PRECISION, INTENT(INOUT):: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT):: TSDER(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE qwindl
     &( tstot , usold , usnew , twold , twnew , nf    , ndire ,
     &  npoin2, usn   , uso   , fpmo  , fpmn )
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(INOUT) :: FPMO(npoin2),FPMN(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: TWOLD(npoin2),TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: USNEW(npoin2),USOLD(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: USO(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: USN(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE radiat
     &        ( fx    , fy    ,xk1   , fs    , cg1   ,
     &        cgsuc1, dsxxdx, dsxydx, dsxydy, dsyydy )
      USE bief_def
      USE declarations_tomawac, ONLY : ndire,nf,npoin2
      IMPLICIT NONE
      DOUBLE PRECISION, INTENT(IN)    :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: CG1(npoin2,nf),XK1(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: CGSUC1(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FX(npoin2),FY(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: DSXXDX(npoin2),DSXYDX(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: DSXYDY(npoin2),DSYYDY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE read_spectra_coords
     &(fid,np,xp,yp)
      IMPLICIT NONE
      INTEGER, INTENT(IN)                         :: FID
      INTEGER, INTENT(IN)                         :: NP
      DOUBLE PRECISION,ALLOCATABLE, INTENT(INOUT) :: XP(:),YP(:)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE rpi_intr
     &(neigb  , nb_close, rx  , ry      , rxx     , ryy ,
     & npoin2 , i       , maxnsp  , ffd , firdiv1 , firdiv2 ,
     & secdiv1, secdiv2 , secdiv3, frstdiv , scnddiv)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, MAXNSP,I
      INTEGER, INTENT(IN) :: NEIGB(npoin2,maxnsp),NB_CLOSE(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: RX(maxnsp),RY(maxnsp)
      DOUBLE PRECISION, INTENT(IN)    :: RXX(maxnsp),RYY(maxnsp)
      DOUBLE PRECISION, INTENT(INOUT) :: SECDIV1,SECDIV2,SECDIV3
      DOUBLE PRECISION, INTENT(INOUT) :: FIRDIV1,FIRDIV2
      DOUBLE PRECISION, INTENT(IN)    :: FFD(npoin2)
      LOGICAL, INTENT(IN)             :: FRSTDIV,SCNDDIV
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
!V6P2 New subroutine
      INTERFACE
        SUBROUTINE rpi_invr
     &( x     , y     , neigb , nb_close, rk_d , rx_d , ry_d  , rxx_d ,
     &  ryy_d , npoin2, i     , quo   , ac    , maxnsp, mindist )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: NPOIN2, MAXNSP, I
      INTEGER, INTENT(IN)    :: NEIGB(npoin2,maxnsp), NB_CLOSE(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: QUO, AC
      DOUBLE PRECISION, INTENT(IN)    :: X(npoin2), Y(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: MINDIST(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: RK_D(maxnsp)
      DOUBLE PRECISION, INTENT(INOUT) :: RX_D(maxnsp), RY_D(maxnsp)
      DOUBLE PRECISION, INTENT(INOUT) :: RXX_D(maxnsp), RYY_D(maxnsp)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE semimp
     &( f     ,   cf , xk   ,  nf  , ndire, npoin2,
     &  iangnl, tstot, tsder, told,  tnew,  z0new, twnew, 
     &  taux1,  taux2, taux3, taux4, taux5, taux6, taux7,
     &  mdia, ianmdi, coemdi, fbor)
        USE declarations_special
        USE declarations_tomawac, ONLY : nptfr
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2,NDIRE,NF, IANGNL(*)
      DOUBLE PRECISION, INTENT(INOUT) :: Z0NEW(npoin2), TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX1(npoin2),
     &                                   taux2(npoin2),taux3(npoin2),
     &                                   taux4(npoin2),taux5(npoin2),
     &                                   taux6(npoin2),taux7(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSDER(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TOLD(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: TNEW(npoin2,ndire)
      DOUBLE PRECISION, INTENT(INOUT) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(IN) :: CF(*)
!....MDIA method declarations
      INTEGER, INTENT(IN) :: MDIA, IANMDI(*)
      DOUBLE PRECISION, INTENT(IN) ::  COEMDI(*)
      DOUBLE PRECISION, INTENT(IN)   :: FBOR(nptfr,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE sor3d
     &(f,ndire,nf,npoin2,vent, couran,maree,titre,tra01,mesh3d)
      USE bief_def
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TRA01(npoin2*ndire)
      LOGICAL, INTENT(IN)             :: COURAN,VENT,MAREE
      CHARACTER(LEN=80), INTENT(IN)   :: TITRE
      TYPE(bief_mesh), INTENT(IN)     :: MESH3D
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE speini
     &( f     , spec  , fra   , uv    , vv    ,
     &  fremax, fetch , sigmaa, sigmab, gamma , fpic  , hm0   , alphil,
     &  teta1 , spred1, teta2 , spred2, xlamda, npoin2, ndire , nf    ,
     &  inispe, depth , frabi )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NPOIN2, NDIRE , NF    , INISPE, FRABI
      DOUBLE PRECISION, INTENT(IN)    :: FREMAX, FETCH , SIGMAA
      DOUBLE PRECISION, INTENT(IN)    :: SIGMAB, GAMMA
      DOUBLE PRECISION, INTENT(IN)    :: FPIC  , HM0   , ALPHIL, TETA1
      DOUBLE PRECISION, INTENT(IN)    :: SPRED1, TETA2
      DOUBLE PRECISION, INTENT(IN)    :: SPRED2, XLAMDA
      DOUBLE PRECISION, INTENT(IN)    :: UV(*) , VV(*)
      DOUBLE PRECISION, INTENT(IN)    :: DEPTH(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FRA(ndire), SPEC(nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE spejon
     &( spec  , nf    , al    , fp     , gamma , sigmaa, sigmab,
     &  fpmin )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF
      DOUBLE PRECISION, INTENT(IN)    :: SIGMAA, SIGMAB, GAMMA, FPMIN
      DOUBLE PRECISION, INTENT(IN)    :: FP    , AL
      DOUBLE PRECISION, INTENT(INOUT) :: SPEC(nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE spetma
     &( spec  , nf    , al    , fp     , gamma , sigmaa, sigmab,
     &  fpmin , depth  )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF
      DOUBLE PRECISION, INTENT(IN)    :: SIGMAA, SIGMAB, GAMMA, FPMIN
      DOUBLE PRECISION, INTENT(IN)    :: FP, AL, DEPTH
      DOUBLE PRECISION, INTENT(INOUT) :: SPEC(nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE dir_spread( DIRSPR, F, NDIRE, NF, NPOIN2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: DIRSPR(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE stress
     &( tauwav, tstot , f     , usnew , twnew , z0new ,
     &  npoin2, ndire , nf    , xtauw , ytauw , tauhf )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NPOIN2, NDIRE , NF
      DOUBLE PRECISION, INTENT(IN)    :: USNEW(npoin2), TWNEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: Z0NEW(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: TSTOT(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUWAV(npoin2), TAUHF(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: XTAUW(npoin2) , YTAUW(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE tautot
     &( taut  , uvent , tauw  , seuil ,
     &  itr   , itrmin, itrmax)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  ITRMIN, ITRMAX
      INTEGER, INTENT(INOUT) ::  ITR
      DOUBLE PRECISION, INTENT(IN)    :: UVENT , TAUW
      DOUBLE PRECISION, INTENT(IN)    :: SEUIL
      DOUBLE PRECISION, INTENT(INOUT) :: TAUT
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE tetmoy
     &( tetam , f     ,ndire , nf    , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TETAM(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE tom_corfon
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE totnrj( VARIAN, F, NF, NDIRE, NPOIN2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::          NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: VARIAN(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE transf
     &( fa    , fr    , xk    , knew  , newf  , newf1 , taux1 , taux2 ,
     &  npoin2, ndire , nf    )
      IMPLICIT NONE
      INTEGER, INTENT(IN)    :: NPOIN2, NDIRE, NF
      INTEGER, INTENT(INOUT) :: KNEW(npoin2),NEWF(npoin2), NEWF1(npoin2)
      DOUBLE PRECISION, INTENT(IN)    :: FR(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: TAUX1(npoin2),TAUX2(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: FA(npoin2,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE ustar1
     &( ustar , z0    , tauwav,  npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: TAUWAV(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: USTAR(npoin2), Z0(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE ustar2
     &( ustar , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::  NPOIN2
      DOUBLE PRECISION, INTENT(INOUT) :: USTAR(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE uvstokes
     &       (ust, vst, wst, fs, npoin2, xk1, zfj, ndire, ztel,nz,nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN)  :: NZ,NF
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2, ndire, nf)
      DOUBLE PRECISION, INTENT(IN) :: ZFJ(npoin2)
      DOUBLE PRECISION, INTENT(IN) :: XK1(npoin2, nf)
      DOUBLE PRECISION, INTENT(INOUT) :: UST(npoin2, nz)
      DOUBLE PRECISION, INTENT(INOUT) :: VST(npoin2, nz)
      DOUBLE PRECISION, INTENT(INOUT) :: WST(npoin2)
      DOUBLE PRECISION, INTENT(IN) :: ZTEL(npoin2, nz)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE venuti
     &(nven,fmtven)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NVEN
      CHARACTER(LEN=8), INTENT(IN)    :: FMTVEN
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_wind
     &(nven,fmtven)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NVEN
      CHARACTER(LEN=8), INTENT(IN)    :: FMTVEN
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE vitfon(VIFOND,F,XK,NF,NPOIN2,NDIRE)
      IMPLICIT NONE
      INTEGER, INTENT(IN)             :: NF,NDIRE,NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: VIFOND(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE wac
     &(part, nit_ori)
      IMPLICIT NONE
      INTEGER,           INTENT(IN)      :: PART
      INTEGER, OPTIONAL, INTENT(IN)      :: NIT_ORI
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE windiss1
     &  (fwx, fwy, npoin2, xk, ndire, fs,nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FWX(npoin2), FWY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE windiss2
     &  (fwx, fwy, npoin2, xk, ndire, fs,nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FWX(npoin2), FWY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE windiss3
     &  (fwx, fwy, npoin2, xk, ndire, fs,nf)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: NPOIN2, NDIRE,NF
      DOUBLE PRECISION, INTENT(IN) :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN) :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: FWX(npoin2), FWY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE wac_init
     & (part,impres, debres,date,time)
      IMPLICIT NONE
      LOGICAL, INTENT(INOUT) :: IMPRES, DEBRES
      INTEGER, INTENT(IN) :: PART
      INTEGER, INTENT(INOUT) :: DATE(3),TIME(3)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE wipj
     & (  wip, fs, npoin2, xk, wipdx, wipdy,ndire,nf)
      INTEGER, INTENT(IN)    :: npoin2,ndire,nf
      DOUBLE PRECISION, INTENT(IN)    :: FS(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: XK(npoin2,nf)
      DOUBLE PRECISION, INTENT(INOUT) :: WIP(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: WIPDX(npoin2)
      DOUBLE PRECISION, INTENT(INOUT) :: WIPDY(npoin2)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE wnscou
     &( ck2   , freq  , depth )
      IMPLICIT NONE
      DOUBLE PRECISION, INTENT(IN)    :: FREQ  , DEPTH
      DOUBLE PRECISION, INTENT(INOUT) :: CK2
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE wpower
     &( f     , cg    , nf    , ndire , npoin2)
      IMPLICIT NONE
      INTEGER, INTENT(IN)    ::          NF    , NDIRE , NPOIN2
      DOUBLE PRECISION, INTENT(IN)    :: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(IN)    :: CG(npoin2,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_anacos
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_anamar
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_anaven
      IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_dump2d
     &(xf1, np1)
      IMPLICIT NONE
      INTEGER, INTENT(IN)          :: NP1
      DOUBLE PRECISION, INTENT(IN) :: XF1(np1)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_limwac
     &(f     , fbor  , nptfr , ndire , nf    , npoin2,
     & kent  , prive , npriv , imp_file)
      USE bief_def, ONLY : bief_file
      IMPLICIT NONE
      INTEGER, INTENT(IN)            :: NPTFR,NDIRE,NF,NPOIN2,NPRIV
      INTEGER, INTENT(IN)            :: KENT
      DOUBLE PRECISION, INTENT(IN)   :: PRIVE(npoin2,npriv)
      TYPE(bief_file), INTENT(IN)    :: IMP_FILE
      DOUBLE PRECISION, INTENT(INOUT):: F(npoin2,ndire,nf)
      DOUBLE PRECISION, INTENT(INOUT):: FBOR(nptfr,ndire,nf)
        END SUBROUTINE
      END INTERFACE
!
!-----------------------------------------------------------------------
!
      INTERFACE
        SUBROUTINE user_tom_corfon
          IMPLICIT NONE
        END SUBROUTINE
      END INTERFACE
!
!=======================================================================
!
      END MODULE interface_tomawac