get_mesh_nelem.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE get_mesh_nelem
!                   *************************
!
     &(fformat,fid,typ_elem,nelem,ierr)
!
!***********************************************************************
! HERMES   V7P0                                               01/05/2014
!***********************************************************************
!
!brief    Returns the number of elements of type typ_elem in the mesh file
!
!history  Y AUDOUIN (LNHE)
!+        24/03/2014
!+        V7P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FFORMAT        |-->| FORMAT OF THE FILE
!| FID            |-->| FILE DESCRIPTOR
!| TYP_ELEM       |-->| TYPE OF THE ELEMENT
!| NELEM          |<->| THE NUMBER OF 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
      INTEGER,          INTENT(IN)  :: TYP_ELEM
      INTEGER,          INTENT(OUT) :: NELEM
      INTEGER,          INTENT(OUT) :: IERR
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      nelem = 0
      IF(typ_elem.EQ.type_null) THEN
        RETURN
      ENDIF
!
      SELECT CASE (fformat(1:7))
        CASE ('SERAFIN')
          CALL get_mesh_nelem_srf(fid, typ_elem, nelem, ierr)
        CASE ('MED    ')
          CALL get_mesh_nelem_med(fid, typ_elem, nelem, ierr)
        CASE ('CGNS   ')
          CALL get_mesh_nelem_cgns(fid, typ_elem, nelem, ierr)
        CASE DEFAULT
          ierr = hermes_unknown_file_format_err
          WRITE(error_message,*)
     &        'GET_MESH_NELEM: BAD FILE FORMAT: ',fformat
          RETURN
      END SELECT
!
!-----------------------------------------------------------------------
!
      RETURN
      END