get_bnd_ipobo.f

Go to the documentation of this file.

!                   ************************
                    SUBROUTINE get_bnd_ipobo
!                   ************************
!
     &(fformat,fid,npoin,nelebd,typ_bnd_elem,ipobo,ierr)
!
!***********************************************************************
! HERMES   V7P0                                               01/05/2014
!***********************************************************************
!
!brief    Returns an array containing
!+        1 if a point is a boundary point 0 otherwise
!
!history  Y AUDOUIN (LNHE)
!+        24/03/2014
!+        V7P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FFORMAT        |-->| FORMAT OF THE FILE
!| FID            |-->| FILE DESCRIPTOR
!| NPOIN          |-->| TOTAL NUMBER OF NODES
!| NELEBD         |-->| TOTAL NUMBER OF BOUNDARY ELEMENTS
!| TYP_BND_ELEM   |-->| TYPE OF THE BOUNDARY ELEMENT
!| IPOBO          |<->| AN ARRAY CONTAINING
!|                |   | 1 IF A POINT IS A BOUNDARY POINT 0 OTHERWISE
!| 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, NPOIN, TYP_BND_ELEM
      INTEGER,          INTENT(IN)  :: NELEBD
      INTEGER,          INTENT(INOUT) :: IPOBO(npoin)
      INTEGER,          INTENT(OUT) :: IERR
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      IF((typ_bnd_elem.EQ.type_null).OR.(nelebd.EQ.0)) THEN
        ipobo(:) = 0
        RETURN
      ENDIF
!
      SELECT CASE (fformat(1:7))
        CASE ('SERAFIN')
          CALL get_bnd_ipobo_srf(fid, npoin,ipobo, ierr)
        CASE ('MED    ')
          CALL get_bnd_ipobo_med(fid,typ_bnd_elem,npoin,ipobo,ierr)
        CASE ('CGNS   ')
          CALL get_bnd_ipobo_cgns(fid, typ_bnd_elem,npoin,ipobo, ierr)
        CASE DEFAULT
          ierr = hermes_unknown_file_format_err
          WRITE(error_message,*)
     &        'GET_BND_IPOBO: BAD FILE FORMAT: ',fformat
          RETURN
      END SELECT
!
!-----------------------------------------------------------------------
!
      RETURN
      END