gtsh31.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE gtsh31
!                   *****************
!
     &(shp,elt,ikle,eltcar,npoin,nelmax)
!
!***********************************************************************
! BIEF   V6P3                                   21/09/2012
!***********************************************************************
!
!brief    Gives the starting element and the barycentric coordinates of
!+        the head of characteristics, here for tetrahedra.
!
!history  J-M HERVOUET (LNHE)
!+        21/09/2012
!+        V6P3
!+        First version, inspired from gtsh41.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ELT            |<--| ELEMENT CHOSEN FOR EVERY POINT
!| ELTCAR         |-->| STARTING ELEMENT FOR LINEAR AND QUADRATIC POINTS
!| IKLE           |-->| CONNECTIVITY TABLE
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
!| NPOIN          |-->| NUMBER OF POINTS
!| NSEG           |-->| NUMBER OF SEGMENTS
!| SHP            |<--| BARYCENTRIC COORDINATES OF NODES IN THEIR
!|                |   | ASSOCIATED ELEMENT "ELT"
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)             :: NPOIN,NELMAX
      INTEGER, INTENT(IN)             :: IKLE(nelmax,4),ELTCAR(npoin)
      INTEGER, INTENT(INOUT)          :: ELT(npoin)
      DOUBLE PRECISION, INTENT(INOUT) :: SHP(4,npoin)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I,IELEM,IPLAN
!
!-----------------------------------------------------------------------
!
      DO i=1,npoin
        ielem=eltcar(i)
        elt(i) = ielem
        IF(ielem.NE.0) THEN
          IF(ikle(ielem,1).EQ.i) THEN
            shp(1,i)=1.d0
            shp(2,i)=0.d0
            shp(3,i)=0.d0
            shp(4,i)=0.d0
          ELSEIF(ikle(ielem,2).EQ.i) THEN
            shp(1,i)=0.d0
            shp(2,i)=1.d0
            shp(3,i)=0.d0
            shp(4,i)=0.d0
          ELSEIF(ikle(ielem,3).EQ.i) THEN
            shp(1,i)=0.d0
            shp(2,i)=0.d0
            shp(3,i)=1.d0
            shp(4,i)=0.d0
          ELSEIF(ikle(ielem,4).EQ.i) THEN
            shp(1,i)=0.d0
            shp(2,i)=0.d0
            shp(3,i)=0.d0
            shp(4,i)=1.d0
          ELSE
            WRITE(lu,*) 'PROBLEM IN GTSH31'
            WRITE(lu,*) 'POINT ',i,' ELEMENT ',ielem
            WRITE(lu,*) 'NPOIN=',npoin
            CALL plante(1)
            stop
          ENDIF
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END