ielbor.f

Go to the documentation of this file.

!                   ***********************
                    INTEGER FUNCTION ielbor
!                   ***********************
!
     &( ielm , i )
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    GIVES THE BOUNDARY ELEMENT TYPE CORRESPONDING TO
!+                A GIVEN ELEMENT TYPE IN THE DOMAIN.
!+
!+            WHEN THERE ARE SEVERAL TYPES (AS IS THE CASE FOR THE
!+                PRISMS FOR EXAMPLE) USES INDEX I TO DISTINGUISH THEM.
!
!history  J-M HERVOUET
!+        06/02/08
!+        V5P9
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| I              |-->| CASE OF SEVERAL BOUNDARY ELEMENTS
!| IELM           |-->| TYPE OF ELEMENT IN THE DOMAIN
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
      INTEGER IELM,I
!
!-----------------------------------------------------------------------
!
      IF(ielm.EQ.10.OR.ielm.EQ.20) THEN
        ielbor = 0
      ELSEIF(ielm.EQ.11.OR.ielm.EQ.12.OR.ielm.EQ.21) THEN
        ielbor = 1
      ELSEIF(ielm.EQ.13) THEN
        ielbor = 2
      ELSEIF(ielm.EQ.30) THEN
        ielbor = 80
      ELSEIF(ielm.EQ.31) THEN
        ielbor = 81
      ELSEIF(ielm.EQ.51.AND.i.EQ.1) THEN
        ielbor = 11
      ELSEIF(ielm.EQ.51.AND.i.EQ.2) THEN
        ielbor = 61
      ELSEIF(ielm.EQ.50.AND.i.EQ.1) THEN
        ielbor = 10
      ELSEIF(ielm.EQ.50.AND.i.EQ.2) THEN
        ielbor = 60
      ELSEIF(ielm.EQ.40.AND.i.EQ.1) THEN
        ielbor = 10
      ELSEIF(ielm.EQ.40.AND.i.EQ.2) THEN
        ielbor = 70
      ELSEIF(ielm.EQ.41.AND.i.EQ.1) THEN
        ielbor = 11
      ELSEIF(ielm.EQ.41.AND.i.EQ.2) THEN
        ielbor = 71
      ELSE
        WRITE(lu,101) ielm
101     FORMAT(1x,'IELBOR (BIEF) : ',1i6,' ELEMENT NOT IMPLEMENTED')
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END