get_bnd_connectivity.f

Go to the documentation of this file.

!                   **********************
                    SUBROUTINE get_bnd_connectivity
!                   **********************
!
     &(fformat,fid,typ_bnd_elem,nelebd,ndp,ikle_bnd,ierr)
!
!***********************************************************************
! HERMES   V7P0                                               01/05/2014
!***********************************************************************
!
!brief    Reads the connectivity of the boundary elements
!
!history  Y AUDOUIN (LNHE)
!+        24/03/2014
!+        V7P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FFORMAT        |-->| FORMAT OF THE FILE
!| FID            |-->| FILE DESCRIPTOR
!| TYP_BND_ELEM   |-->| TYPE OF THE BOUNDARY ELEMENTS
!| NELEBD         |-->| NUMBER OF BOUNDARY ELEMENTS
!| NDP            |-->| NUMBER OF POINTS PER ELEMENT
!| IKLE_BND       |<->| THE CONNECTIVITY OF THE BOUNDARY ELEMENTS
!| 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,TYP_BND_ELEM,NELEBD,NDP
      INTEGER, INTENT(INOUT) :: IKLE_BND(ndp*nelebd)
      INTEGER, INTENT(OUT) :: IERR
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      IF((typ_bnd_elem.EQ.type_null).OR.(nelebd.EQ.0)) RETURN
!
      SELECT CASE (fformat(1:7))
        CASE ('SERAFIN')
          CALL get_bnd_connectivity_srf(fid, typ_bnd_elem,nelebd,ndp,
     &                          ikle_bnd,ierr)
        CASE ('MED    ')
          CALL get_bnd_connectivity_med(fid, typ_bnd_elem,nelebd,ndp,
     &                          ikle_bnd,ierr)
        CASE ('CGNS   ')
          CALL get_bnd_connectivity_cgns(fid, typ_bnd_elem,nelebd,ndp,
     &                          ikle_bnd,ierr)
        CASE DEFAULT
          ierr = hermes_unknown_file_format_err
          WRITE(error_message,*)
     & 'GET_BND_CONNECTIVITY: BAD FILE FORMAT: ', fformat
          RETURN
      END SELECT
!
!-----------------------------------------------------------------------
!
      RETURN
      END