pttoel.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE pttoel
!                   *****************
!
     &(xel,x,mesh)
!
!***********************************************************************
! BIEF   V6P3                                   21/08/2010
!***********************************************************************
!
!brief    GOES FROM A VECTOR BY POINTS TO A VECTOR BY ELEMENTS.
!
!history  J-M HERVOUET (LNH)
!+        10/01/95
!+        V5P6
!+
!
!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 (LNH)
!+        03/01/2013
!+        V6P3
!+    Adding 3D cases.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| MESH           |-->| MESH STRUCTURE
!| X              |-->| VECTOR DEFINED PER POINT
!| XEL            |<--| SAME VECTOR DEFINED PER ELEMENT
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_pttoel => pttoel
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!  VECTORS STRUCTURES
!
      TYPE(bief_obj), INTENT(IN)    :: X
      TYPE(bief_obj), INTENT(INOUT) :: XEL
!
!-----------------------------------------------------------------------
!
!  MESH STRUCTURE
!
      TYPE(bief_mesh), INTENT(IN) :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      IF(x%ELM.EQ.11) THEN
!
        CALL ptel11(xel%R,x%R,mesh%IKLE%I,mesh%NELMAX,mesh%NELEM)
!
      ELSEIF(x%ELM.EQ.31.OR.x%ELM.EQ.51) THEN
!
        CALL ptel31(xel%R,x%R,mesh%IKLE%I,mesh%NELMAX,mesh%NELEM)
!
      ELSEIF(x%ELM.EQ.41) THEN
!
        CALL ptel41(xel%R,x%R,mesh%IKLE%I,mesh%NELMAX,mesh%NELEM)
!
      ELSE
        WRITE(lu,101) x%ELM
101   FORMAT(1x,'PTTOEL (BIEF) : IELM = ',1i6,' ELEMENT NOT AVAILABLE')
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END