set_bnd.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE set_bnd
!                   ******************
!
     &(fformat,fid,type_bnd_elt,nelebd,ndp,ikle,nptfr,lihbor,liubor,
     & livbor,hbor,ubor,vbor,chbord,
     & litbor,tbor,atbor,btbor,color,ierr)
!
!***********************************************************************
! HERMES   V7P0                                               01/05/2014
!***********************************************************************
!
!brief    Writes the boundary information into the mesh file
!
!history  Y AUDOUIN (LNHE)
!+        24/03/2014
!+        V7P0
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| FFORMAT        |-->| FORMAT OF THE FILE
!| FID            |-->| FILE DESCRIPTOR
!| TYPE_BND_ELT   |-->| TYPE OF THE BOUNDARY ELEMENTS
!| NELEBD         |-->| NUMBER OF BOUNDARY ELEMENTS
!| NDP            |-->| NUMBER OF POINTS PER BOUNDARY ELEMENT
!| IKLE           |-->| CONNECTIVITY ARRAY FOR THE BOUNDARY ELEMENTS
!| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
!| LIHBOR         |-->| TYPE OF BOUNDARY CONDITIONS ON DEPTH
!| LIUBOR         |-->| TYPE OF BOUNDARY CONDITIONS ON U
!| LIVBOR         |-->| TYPE OF BOUNDARY CONDITIONS ON V
!| HBOR           |<--| PRESCRIBED BOUNDARY CONDITION ON DEPTH
!| UBOR           |<--| PRESCRIBED BOUNDARY CONDITION ON VELOCITY U
!| VBOR           |<--| PRESCRIBED BOUNDARY CONDITION ON VELOCITY V
!| CHBORD         |<--| FRICTION COEFFICIENT AT BOUNDARY
!| LITBOR         |-->| PHYSICAL BOUNDARY CONDITIONS FOR TRACERS
!| TBOR           |<--| PRESCRIBED BOUNDARY CONDITION ON TRACER
!| ATBOR,BTBOR    |<--| THERMAL EXCHANGE COEFFICIENTS.
!| COLOR          |<--| Boundary color of the boundary 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)  :: TYPE_BND_ELT
      INTEGER,          INTENT(IN)  :: NELEBD
      INTEGER,          INTENT(IN)  :: NDP
      INTEGER,          INTENT(IN)  :: IKLE(nelebd*ndp)
      INTEGER,          INTENT(IN)  :: NPTFR
      INTEGER,          INTENT(IN)  :: LIUBOR(nptfr),LIVBOR(nptfr)
      INTEGER,          INTENT(IN)  :: LIHBOR(nptfr),LITBOR(nptfr)
      DOUBLE PRECISION, INTENT(IN)  :: UBOR(nptfr),VBOR(nptfr)
      DOUBLE PRECISION, INTENT(IN)  :: HBOR(nptfr),CHBORD(nptfr)
      DOUBLE PRECISION, INTENT(IN)  :: TBOR(nptfr),ATBOR(nptfr)
      DOUBLE PRECISION, INTENT(IN)  :: BTBOR(nptfr)
      INTEGER,          INTENT(IN)  :: COLOR(nptfr)
      INTEGER,          INTENT(OUT) :: IERR
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      SELECT CASE (fformat(1:7))
        CASE ('SERAFIN')
          CALL set_bnd_srf(fid,type_bnd_elt,nelebd,ndp,ikle,
     &                     lihbor,liubor,
     &                     livbor,hbor,ubor,vbor,chbord,
     &                     litbor,tbor,atbor,btbor,color,ierr)
        CASE ('MED    ')
          CALL set_bnd_med(fid,type_bnd_elt,nelebd,ndp,ikle,nptfr,
     &                     lihbor,liubor,livbor,litbor,ierr)
        CASE ('CGNS')
          CALL set_bnd_cgns(fid,type_bnd_elt,nelebd,ndp,ikle,
     &                     lihbor,liubor,
     &                     livbor,hbor,ubor,vbor,chbord,
     &                     litbor,tbor,atbor,btbor,ierr)
        CASE DEFAULT
          ierr = hermes_unknown_file_format_err
          WRITE(error_message,*)
     &         'GET_SET_BND: BAD FILE FORMAT: ',fformat
          RETURN
      END SELECT
!
!-----------------------------------------------------------------------
!
      RETURN
      END