cpikle2.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE cpikle2
!                   ******************
!
     &(ikle3,klei3,ikles,nelem2,nelmax2,npoin2,nplan)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    EXTENDS THE CONNECTIVITY TABLE.
!+                CASE OF A MESH OF PRISMS, IKLE BUILT FROM THE
!+                CONNECTIVITY OF THE TRIANGLE MESH.
!
!history  J-M HERVOUET (LNH)
!+        23/08/99
!+        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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE3          |<->| 3D CONNECTIVITY TABLE
!| IKLES          |-->| 2D CONNECTIVITY TABLE WITH DIMENSION (3,NELEM2)
!| KLEI3          |<--| LIKE IKLE3 BUT WITH SWAPPED DIMENSIONS
!| NELEM2         |-->| NUMBER OF ELEMENTS IN 2D
!| NELMAX2        |-->| MAXIMUM NUMBER OF ELEMENTS IN 2D
!| NPLAN          |-->| NUMBER OF PLANES
!| NPOIN2         |-->| NUMBER OF POINTS IN 2D
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NELEM2,NELMAX2,NPOIN2,NPLAN
      INTEGER, INTENT(IN)    :: IKLES(3,nelem2)
      INTEGER, INTENT(INOUT) :: IKLE3(nelmax2,nplan-1,6)
      INTEGER, INTENT(INOUT) :: KLEI3(6,nelmax2,nplan-1)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM,I
!
!-----------------------------------------------------------------------
!
!     BOTTOM AND TOP OF ALL LAYERS
!
      IF(nplan.GE.2) THEN
        DO i = 1,nplan-1
          DO ielem = 1,nelem2
            ikle3(ielem,i,1) = ikles(1,ielem) + (i-1)*npoin2
            ikle3(ielem,i,2) = ikles(2,ielem) + (i-1)*npoin2
            ikle3(ielem,i,3) = ikles(3,ielem) + (i-1)*npoin2
            ikle3(ielem,i,4) = ikles(1,ielem) +  i   *npoin2
            ikle3(ielem,i,5) = ikles(2,ielem) +  i   *npoin2
            ikle3(ielem,i,6) = ikles(3,ielem) +  i   *npoin2
            klei3(1,ielem,i) = ikles(1,ielem) + (i-1)*npoin2
            klei3(2,ielem,i) = ikles(2,ielem) + (i-1)*npoin2
            klei3(3,ielem,i) = ikles(3,ielem) + (i-1)*npoin2
            klei3(4,ielem,i) = ikles(1,ielem) +  i   *npoin2
            klei3(5,ielem,i) = ikles(2,ielem) +  i   *npoin2
            klei3(6,ielem,i) = ikles(3,ielem) +  i   *npoin2
          ENDDO
        ENDDO
      ELSE
        WRITE(lu,*) 'CPIKLE2 : MINIMUM OF 2 PLANES NEEDED'
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END