slopes.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE slopes
!                   *****************
!
     &(coef,z,mesh)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPUTES THE COEFFICIENT 1 / COS(ALFA)
!+                WHERE ALFA IS THE SLOPE OF A TRIANGULAR ELEMENT.
!+
!+            THIS COEFFICIENT IS USED IN THE BOTTOM FRICTION
!+                TERM.
!
!warning  THE JACOBIAN MUST BE POSITIVE
!
!history  J-M HERVOUET (LNH)
!+        17/08/94
!+        V5P5
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| COEF           |<--| RESULT
!| MESH           |-->| MESH STRUCTURE
!| Z              |-->| BOTTOM ELEVATIONS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_slopes => slopes
!
      USE declarations_special
      IMPLICIT NONE
!
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  STRUCTURES OF VECTOR
!
      TYPE(bief_obj), INTENT(INOUT) :: COEF
      TYPE(bief_obj), INTENT(IN)    :: Z
!
!-----------------------------------------------------------------------
!
!  MESH STRUCTURE
!
      TYPE(bief_mesh), INTENT(IN)   :: MESH
!
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
      INTEGER NELEM,NELMAX,IELM
!
!-----------------------------------------------------------------------
!
      ielm   = mesh%TYPELM
      nelem  = mesh%NELEM
      nelmax = mesh%NELMAX
!
!-----------------------------------------------------------------------
!
      coef%ELM = ielm
      coef%DIM1= bief_nbpts(ielm,mesh)
!
!-----------------------------------------------------------------------
!
      IF(ielm.EQ.10) THEN
        CALL slop10(coef%R,mesh%XEL%R,mesh%YEL%R,
     &              z%R,mesh%IKLE%I,nelem,nelmax)
      ELSE
        WRITE(lu,301) mesh%TYPELM
301     FORMAT(1x,'SLOPES (BIEF) : UNKNOWN ELEMENT:',1i6)
        CALL plante(0)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END