comp_ikle.f

Go to the documentation of this file.

!                   ********************
                    SUBROUTINE comp_ikle
!                   ********************
!
     &(ikle,iklbor,eltseg,nbor,nelbor,nulone,
     & ielm,nelem,nelmax,npoin,nptfr,neleb,nelebx)
!
!***********************************************************************
! BIEF   V7P0                                    28/03/2014
!***********************************************************************
!
!brief    EXTENDS THE CONNECTIVITY TABLES AND ARRAY NBOR.
!
!history  J-M HERVOUET (LNHE)
!+        20/03/2008
!+        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
!
!history  J-M HERVOUET (LNHE)
!+        16/05/2012
!+        V6P2
!+   New arguments passed to CPIK13.
!
!history  J-M HERVOUET (LNHE)
!+        28/03/2014
!+        V7P0
!+   Adaoted to new numbering of boundary segments.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ELTSEG         |-->| SEGMENT NUMBERS OF AN ELEMENT
!| IELM           |-->| TYPE OF ELEMENT
!| IKLBOR         |<->| CONNECTIVITY TABLE FOR BOUNDARY POINTS
!| IKLE           |<->| CONNECTIVITY TABLE FOR ALL POINTS
!| NBOR           |<->| GLOBAL NUMBERS OF BOUNDARY POINTS
!| NELBOR         |-->| BOUNDARY ELEMENT THAT CONTAINS SEGMENT K
!| NELEB          |-->| NUMBER OF BOUNDARY SEGMENTS
!| NELEBX         |-->| MAXIMUM NUMBER OF BOUNDARY SEGMENTS
!| NELEM          |-->| NOMBRE D'ELEMENTS
!| NELMAX         |-->| NOMBRE MAXIMUM D'ELEMENTS
!| NPOIN          |-->| NOMBRE DE SOMMETS DU MAILLAGE
!| NPTFR          |-->| NUMBER OF (LINEAR) BOUNDARY POINTS
!| NULONE         |-->| LOCAL NUMBER OF K IN ELEMENT NELBOR
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_comp_ikle => comp_ikle
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NELEM,NELMAX,IELM,NPOIN,NPTFR
      INTEGER, INTENT(IN)    :: NELEB,NELEBX
      INTEGER, INTENT(IN)    :: ELTSEG(nelmax,3)
      INTEGER, INTENT(IN)    :: NELBOR(nelebx),NULONE(nelebx)
      INTEGER, INTENT(INOUT) :: IKLE(nelmax,*),IKLBOR(nelebx,*),NBOR(*)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      IF(ielm.EQ.12) THEN
!
        CALL cpik12(ikle,nelem,nelmax,npoin)
!
      ELSEIF(ielm.EQ.13.OR.ielm.EQ.14) THEN
!
        CALL cpik13(ikle,iklbor,eltseg,nbor,nelbor,nulone,
     &              nelem,nelmax,npoin,nptfr,neleb,nelebx)
!
      ELSE
!
        WRITE(lu,11) ielm
11      FORMAT(1x,'COMP_IKLE: DISCRETIZATION NOT IMPLEMENTED:',i6)
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END