calcue.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE calcue
!                   *****************
!
!
!***********************************************************************
! ARTEMIS   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPUTES AN EFFECTIVE SPEED UE FOR THE ESTIMATION
!+                OF THE FRICTION DISSIPATION COEFFICIENT FW.
!
!history  J-M HERVOUET (LNH)
!+
!+
!+   LINKED TO BIEF 5.0
!
!history  D. AELBRECHT (LNH) ; D. PAUGAM ( PLACEMENT)
!+        02/06/1999
!+        V5P1
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE declarations_artemis
!
      USE declarations_special
      IMPLICIT NONE
!
      INTRINSIC abs,exp
!
      CALL vector(u0 , '=' , 'GRADF          X' , ielm ,
     &            1.d0 , phir , sbid , sbid , sbid , sbid , sbid ,
     &            mesh , msk  , maskel )
!
      CALL vector(v0 , '=' , 'GRADF          Y' , ielm ,
     &            1.d0 , phir , c , c , c , c , c ,
     &            mesh , msk  , maskel )
!
!
!      THE OLD VARIABLE U1 IS STORED IN T1
!
      CALL vector(t1, '=' , 'GRADF          X' , ielm ,
     &            1.d0 , phii , c , c , c , c , c ,
     &            mesh , msk  , maskel )
!
!
!      THE OLD VARIABLE V1 IS STORED IN T2
!
      CALL vector(t2, '=' , 'GRADF          Y' , ielm ,
     &            1.d0 , phii , c , c , c , c , c ,
     &            mesh , msk  , maskel )
!
      CALL vector(t4 , '=' , 'MASBAS          ' , ielm ,
     &            1.d0 , c , c , c , c , c , c ,
     &            mesh , msk  , maskel )
!
      CALL os('X=Y/Z   ', x=u0, y=u0, z=t4)
      CALL os('X=Y/Z   ', x=v0, y=v0, z=t4)
      CALL os('X=Y/Z   ', x=t1, y=t1, z=t4)
      CALL os('X=Y/Z   ', x=t2, y=t2, z=t4)
!
!--------------------------------------------------------------
!             COMPUTES UE
!--------------------------------------------------------------
!
      CALL os('X=C     ', x=t4 , c=0.d0)
      CALL os('X=YZ    ', x=t4 , y=u0 , z=u0)
      CALL os('X=X+YZ  ', x=t4 , y=v0 , z=v0)
      CALL os('X=X+YZ  ', x=t4 , y=t1 , z=t1)
      CALL os('X=X+YZ  ', x=t4 , y=t2 , z=t2)
!
      CALL os('X=CX    ', x=t4 , c=0.5d0)
      CALL os('X=SQR(Y)', x=t1 , y=t4)
      CALL os('X=CY    ', x=t4 , y=t1, c=1.2d0 )
!
      RETURN
      END