kerbou.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\tomawac\kerbou.f
00002 !
00056                      DOUBLE PRECISION FUNCTION KERBOU
00057 !                    ********************************
00058 !
00059      &( XK1   , XK2   , FREQ1  , FREQ2  , DEPTH , TETA1 , TETA2 )
00060 !
00061 !***********************************************************************
00062 ! TOMAWAC   V6P3                                   20/06/2011
00063 !***********************************************************************
00064 !
00065 !
00066 !
00067 !
00068 !
00069 !
00070 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00071 !| DEPTH          |-->| WATER DEPTH
00072 !| FREQ1          |-->| FREQUENCY OF COMPONENT 1
00073 !| FREQ2          |-->| FREQUENCE OF COMPONENT 2
00074 !| TETA1          |-->| DIRECTION OF COMPONENT 1
00075 !| TETA2          |-->| DIRECTION OF COMPONENT 2
00076 !| XK1            |-->| WAVE NUMBER OF COMPONENT 1
00077 !| XK2            |-->| WAVE NUMBER OF COMPONENT 2
00078 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00079 !
00080       USE DECLARATIONS_TOMAWAC, ONLY : DEUPI,GRAVIT
00081 !
00082       IMPLICIT NONE
00083 !
00084 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00085 !
00086       DOUBLE PRECISION, INTENT(IN) :: XK1,XK2,FREQ1,FREQ2,TETA1,TETA2
00087       DOUBLE PRECISION, INTENT(IN) :: DEPTH
00088 !
00089 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00090 !
00091       DOUBLE PRECISION VAR1,VAR2,VAR3,VAR4,DANG
00092 !
00093       INTRINSIC COS
00094 !
00095 !-----------------------------------------------------------------------
00096 !
00097       VAR1  = XK1**2
00098       VAR2  = XK2**2
00099       VAR3  = XK1*XK2
00100       VAR4  = DEUPI**2*FREQ1*FREQ2
00101       DANG  = COS(TETA1-TETA2)
00102 !
00103       KERBOU = GRAVIT*0.5D0*(VAR1+VAR2+2.D0*VAR3*DANG) +
00104      &         (VAR4/VAR3)*((VAR1+VAR2)*DANG+VAR3*(1.D0+DANG**2))/DEPTH
00105 !
00106 !-----------------------------------------------------------------------
00107 !
00108       RETURN
00109       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0