The TELEMAC-MASCARET system  trunk
taub_waqtel.f
Go to the documentation of this file.
1 ! **********************
2  SUBROUTINE taub_waqtel
3 ! **********************
4 !
5  &(cf,density,taub,npoin,un,vn)
6 !
7 !***********************************************************************
8 ! TELEMAC2D V7P1
9 !***********************************************************************
10 !
11 !brief COMPUTES BED SHEAR STRESS FOR WAQTEL.
12 !+ SEE THE USE OF TOB_SISYPHE
13 !
14 !
15 !
16 !history R. ATA (EDF LAB, LNHE)
17 !+ 17/03/2016
18 !+ V7P2
19 !+ First version.
20 !
21 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
22 !| CF |-->| FRICTION COEFFICIENT
23 !| DENSITY |-->| DENSITY
24 !| NPOIN |-->| TOTAL NUMBER OF MESH NODES
25 !| TAUB |<--| BED SHEAR STRESS
26 !| UN,VN |-->| VELOCITY COMPONENTS
27 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
28 !
29  USE bief
31  USE interface_waqtel, ex_taub_waqtel => taub_waqtel
32 !
33  IMPLICIT NONE
34 !
35 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
36 !
37  INTEGER , INTENT(IN) :: NPOIN
38  DOUBLE PRECISION , INTENT(IN) :: DENSITY
39  TYPE(bief_obj) , INTENT(IN) :: CF,UN,VN
40  TYPE(bief_obj) , INTENT(INOUT) :: TAUB
41 !
42 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
43 !
44  INTEGER I
45  DOUBLE PRECISION CC
46 !
47  INTRINSIC sqrt
48 !
49 !-----------------------------------------------------------------------
50 !
51 ! TOB=0.5*RHO*CF*U^2
52 !
53  cc=0.5d0*density
54 !
55  CALL cpstvc(cf,taub)
56 !
57  DO i=1,npoin
58  taub%R(i)=cc*cf%R(i)*(un%R(i)**2+vn%R(i)**2)
59  ENDDO
60 !
61 !-----------------------------------------------------------------------
62 !
63  RETURN
64  END
65 
subroutine taub_waqtel(CF, DENSITY, TAUB, NPOIN, UN, VN)
Definition: taub_waqtel.f:7
subroutine cpstvc(X, Y)
Definition: cpstvc.f:7
Definition: bief.f:3