The TELEMAC-MASCARET system  trunk
tobcw_gaia.f
Go to the documentation of this file.
1 ! *********************
2  SUBROUTINE tobcw_gaia
3 ! *********************
4 !
5  & (tob, tobw, thetac, thetaw, tobcw_mean, tobcw_max, npoin)
6 !
7 !***********************************************************************
8 ! GAIA
9 !***********************************************************************
10 !
12 !
13 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
21 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
22 !
23  USE interface_gaia, ex_tobcw_gaia=>tobcw_gaia
25  IMPLICIT NONE
26 !
27 !
28 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
29 !
30  INTEGER, INTENT(IN) :: NPOIN
31 !
32  DOUBLE PRECISION, INTENT(IN) :: THETAC(npoin), THETAW(npoin)
33  DOUBLE PRECISION, INTENT(IN) :: TOB(npoin),TOBW(npoin)
34  DOUBLE PRECISION, INTENT(INOUT) :: TOBCW_MEAN(npoin)
35  DOUBLE PRECISION, INTENT(INOUT) :: TOBCW_MAX(npoin)
36  !
37 !!-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
38 !
39  INTEGER I
40  DOUBLE PRECISION DIFF_THETA, PI
41 !
42 !-----------------------------------------------------------------------
43 !
44  pi = 4.d0 * atan( 1.d0 )
45 !
46 !-----------------------------------------------------------------------
47 !
48  DO i=1,npoin
49 ! COMPUTE MEAN WAVE + CURRENT SHEAR STRESS SOULSBY (1997)
50  IF(tob(i)+tobw(i).GT.1e-9) THEN
51  tobcw_mean(i)= tob(i)*(1.d0 + 1.2d0 *
52  & (tobw(i)/(tob(i)+tobw(i)))**3.2d0)
53  ELSE
54  tobcw_mean(i) = 0.d0
55  ENDIF
56 ! COMPUTE MAXIMUM WAVE + CURRENT SHEAR STRESS SOULSBY (1997)
57  diff_theta = thetaw(i)-thetac(i)
58  tobcw_max(i) = dsqrt((tobcw_mean(i)+
59  & tobw(i)*dcos(diff_theta*pi/180.d0))**2 +
60  & (tobw(i)*dsin(diff_theta*pi/180.d0))**2)
61  ENDDO
62 !
63 !-----------------------------------------------------------------------
64 !
65  RETURN
66  END SUBROUTINE tobcw_gaia
subroutine tobcw_gaia(TOB, TOBW, THETAC, THETAW, TOBCW_MEAN, TOBCW_MAX, NPOIN)
Definition: tobcw_gaia.f:7