get_mesh_connectivity.f

Go to the documentation of this file.

!                   ********************************
                    SUBROUTINE get_mesh_connectivity
!                   ********************************
!
     &(fformat,fid,typ_elem,ikle,nelem,ndp,ierr)
!
!***********************************************************************
! HERMES   V7P0                                               01/05/2014
!***********************************************************************
!
!brief    Returns the connectivity table for
!+        the element of type typ_elem in the mesh
!+        will do nothing if there are no element of typ_elem in the mesh
!
!history  Y AUDOUIN (LNHE)
!+        24/03/2014
!+        V7P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FFORMAT        |-->| FORMAT OF THE FILE
!| FID            |-->| FILE DESCRIPTOR
!| TYP_ELEM       |-->| TYPE OF THE ELEMENT
!| IKLE           |<->| THE CONNECTIVITY TABLE
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NDP            |-->| NUMBER OF POINTS PER ELEMENT
!| IERR           |<--| 0 IF NO ERROR DURING THE EXECUTION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE utils_serafin
      USE utils_med
      USE utils_cgns
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      CHARACTER(LEN=8), INTENT(IN)  :: FFORMAT
      INTEGER,          INTENT(IN)  :: FID
      INTEGER,          INTENT(IN)  :: TYP_ELEM
      INTEGER,          INTENT(IN)  :: NELEM
      INTEGER,          INTENT(IN)  :: NDP
      INTEGER,          INTENT(INOUT) :: IKLE(nelem*ndp)
      INTEGER,          INTENT(OUT) :: IERR
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      SELECT CASE (fformat(1:7))
        CASE ('SERAFIN')
          CALL get_mesh_connectivity_srf(fid, typ_elem, ikle, nelem,
     &                                   ndp, ierr)
        CASE ('MED    ')
          CALL get_mesh_connectivity_med(fid, typ_elem, ikle, nelem,
     &                                   ndp, ierr)
        CASE ('CGNS   ')
          CALL get_mesh_connectivity_cgns(fid, typ_elem, ikle, nelem,
     &                                   ndp, ierr)
        CASE DEFAULT
          ierr = hermes_unknown_file_format_err
          WRITE(error_message,*)
     &       'GET_MESH_CONNECTIVITY: BAD FILE FORMAT: ', fformat
          RETURN
      END SELECT
!
!-----------------------------------------------------------------------
!
      RETURN
      END