ptel41.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE ptel41
!                   *****************
!
     &(xel,x,ikle,nelmax,nelem)
!
!***********************************************************************
! BIEF   V6P3                                          03/01/2013
!***********************************************************************
!
!brief    GOES FROM A VECTOR BY POINTS TO A VECTOR BY ELEMENTS.
!+                CASE OF A LINEAR PRISM.
!
!history  J-M HERVOUET (LNH)
!+        03/01/2013
!+        V6P3
!+    First version.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE           |-->| CONNECTIVITY TABLE
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
!| X              |-->| VECTOR DEFINED PER POINT
!| XEL            |<--| SAME VECTOR DEFINED PER ELEMENT
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NELEM,NELMAX,IKLE(nelmax,6)
!
!-----------------------------------------------------------------------
!
!     VECTORS STRUCTURES
!
      DOUBLE PRECISION, INTENT(IN)    :: X(*)
      DOUBLE PRECISION, INTENT(INOUT) :: XEL(nelmax,6)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM
!
!-----------------------------------------------------------------------
!
      DO ielem = 1,nelem
!
        xel(ielem,1)=x(ikle(ielem,1))
        xel(ielem,2)=x(ikle(ielem,2))
        xel(ielem,3)=x(ikle(ielem,3))
        xel(ielem,4)=x(ikle(ielem,4))
        xel(ielem,5)=x(ikle(ielem,5))
        xel(ielem,6)=x(ikle(ielem,6))
!
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END