calueb2.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\artemis\calueb2.f
00002 !
00050                      SUBROUTINE CALUEB2
00051 !                    ******************
00052 !
00053 !
00054 !***********************************************************************
00055 ! ARTEMIS   V7P0                                   06/2014
00056 !***********************************************************************
00057 !
00058 !
00059 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00060 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00061 !
00062       USE BIEF
00063       USE DECLARATIONS_TELEMAC
00064       USE DECLARATIONS_ARTEMIS
00065 !
00066       IMPLICIT NONE
00067       INTEGER LNG,LU
00068       COMMON/INFO/LNG,LU
00069 !
00070       INTEGER I
00071 !
00072       DOUBLE PRECISION PI,RADDEG,DEUKD,KD,UI
00073       DOUBLE PRECISION ZERO, BID
00074 !
00075       INTRINSIC SQRT, ATAN2, DMOD, ABS, COS, SIN
00076 !
00077 !-----------------------------------------------------------------------
00078 !
00079 !
00080       PARAMETER (ZERO = 1.D-10)
00081       PARAMETER (PI = 3.1415926535897932384626433D0)
00082       PARAMETER (RADDEG = 57.29577951D0)
00083 
00084 !
00085 !
00086 !=======================================================================
00087 ! CASE 1 : USE OF BOTTOM FRICTION COMPUTED FROM THE POTENTIAL
00088 !=======================================================================
00089 ! VELOCITY IN Z=-D FOR RANDOM WAVES (see TOMAWAC vitfon.f)
00090 !    UE^2 = UE^2 + 2 * Sp(f,teta) df dteta * K g / SINH(DEUKD) = UE^2 + Ai²*K g / SINH(DEUKD)
00091 !    Ue :Ai = Hi/2 = (HHO/2)
00092 !    But here we look to RMS value : Ue => Ue/sqrt(2) and we write Ai²=HHO/8
00093 !=======================================================================
00094 !
00095 !=======================================================================
00096 ! CASE 2 : USE OF BOTTOM FRICTION COMPUTED FROM THE STOKES LINEAR THEORY
00097 !          + linear hypothesis abs(U) * U = 8/3PI U
00098 !=======================================================================
00099 ! VELOCITY IN Z=-D FOR RANDOM WAVES
00100 !   Ui = H/2 * (gk/w)*1/ch(kD)
00101 !=======================================================================
00102 !
00103       IF (FORMFR.EQ.1) THEN
00104         DO I=1,NPOIN
00105           DEUKD=2D0*K%R(I)*H%R(I)
00106 !         UPDATE OF UEB WITH ACTUAL WAVE COMPONENT
00107           UEB%R(I)=UEB%R(I)+HHO%R(I)**2/8D0*K%R(I)*GRAV/SINH(DEUKD)
00108         ENDDO
00109       ELSEIF (FORMFR.EQ.2) THEN
00110         DO I=1,NPOIN
00111           KD=K%R(I)*H%R(I)
00112 !         UPDATE OF UEB WITH ACTUAL WAVE COMPONENT
00113           UI      =(8D0/(3D0*PI))*(HHO%R(I)/2D0)*GRAV*K%R(I)
00114           UI      =UI/(OMEGA*COSH(KD))
00115           UEB%R(I)=UEB%R(I)+UI**2
00116         ENDDO
00117       ELSE
00118         WRITE(6,*) 'YOUR OPTION FOR BOTTOM FRICTION IS NOT CORRECT'
00119         WRITE(6,*) 'CALUEB2 : PLEASE USE FORMFR = 1 OR FORMFR = 2 '
00120         WRITE(6,*) 'THE CODE IS GOING TO STOP.....................'
00121         STOP
00122       ENDIF
00123 !
00124       RETURN
00125       END SUBROUTINE

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