taub_waqtel.f

Go to the documentation of this file.

!                   **********************
                    SUBROUTINE taub_waqtel
!                   **********************
!
     &(cf,density,taub,npoin,un,vn)
!
!***********************************************************************
! TELEMAC2D   V7P1
!***********************************************************************
!
!brief    COMPUTES BED SHEAR STRESS FOR WAQTEL.
!+        SEE THE USE OF TOB_SISYPHE
!
!
!
!history  R. ATA (EDF LAB, LNHE)
!+        17/03/2016
!+        V7P2
!+    First version.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| CF             |-->| FRICTION COEFFICIENT
!| DENSITY        |-->| DENSITY
!| NPOIN          |-->| TOTAL NUMBER OF MESH NODES
!| TAUB           |<--| BED SHEAR STRESS
!| UN,VN          |-->| VELOCITY COMPONENTS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      USE interface_waqtel, ex_taub_waqtel => taub_waqtel
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER          , INTENT(IN)    :: NPOIN
      DOUBLE PRECISION , INTENT(IN)    :: DENSITY
      TYPE(bief_obj)   , INTENT(IN)    :: CF,UN,VN
      TYPE(bief_obj)   , INTENT(INOUT) :: TAUB
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
      DOUBLE PRECISION CC
!
      INTRINSIC sqrt
!
!-----------------------------------------------------------------------
!
!     TOB=0.5*RHO*CF*U^2
!
      cc=0.5d0*density
!
      CALL cpstvc(cf,taub)
!
      DO i=1,npoin
        taub%R(i)=cc*cf%R(i)*(un%R(i)**2+vn%R(i)**2)
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END