lecsellim.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE lecsellim
!                   *****************
!
     &(nlim,lihbor,liubor,livbor,hbor,ubor,vbor,
     & chbord,nbor,npmax,nptfr,ncolor)
!
!***********************************************************************
! STBTEL
!***********************************************************************
!
!brief    READS THE BOUNDARY CONDITIONS FILE AND
!+                STORES IN ARRAYS THE DATA READ.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| NLIM           |-->| LOGICAL UNIT OF BOUNDARY CONDITIONS FILE
!| 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
!| NBORD          |<--| BOUNDARY NUMBERING
!| NPTFR          |-->| NUMBER OF BOUNDARY POINTS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      USE declarations_stbtel, ONLY: fformat, typ_bnd_elem
      USE interface_hermes
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NLIM
      INTEGER, INTENT(IN)    :: NPMAX
      INTEGER, INTENT(INOUT) :: LIUBOR(npmax),LIVBOR(npmax)
      INTEGER, INTENT(INOUT) :: LIHBOR(npmax)
      INTEGER, INTENT(INOUT) :: NBOR(npmax)
      INTEGER, INTENT(OUT) :: NPTFR
      DOUBLE PRECISION,  INTENT(INOUT) :: UBOR(npmax,2),VBOR(npmax,2)
      DOUBLE PRECISION,  INTENT(INOUT) :: HBOR(npmax),CHBORD(npmax)
      INTEGER, INTENT(INOUT) :: NCOLOR(npmax)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER :: I, IERR
!
      CALL get_bnd_npoin(fformat,nlim,typ_bnd_elem,nptfr,ierr)

      CALL get_bnd_value(fformat, nlim, typ_bnd_elem, nptfr, lihbor,
     &                   liubor, livbor, hbor, ubor(:,1), vbor(:,1),
     &                   chbord, .false., lihbor, hbor, hbor, hbor,
     &                   nptfr, ierr)
      WRITE(lu,*) 'NPTFR FROM LECSELLIM ', nptfr
      CALL check_call(ierr, 'LECSELIM:GET_BND_VALUE')

      CALL get_bnd_numbering(fformat, nlim, typ_bnd_elem, nptfr,
     &                       nbor, ierr)
      CALL check_call(ierr, 'LECSELIM:GET_BND_NUMBERING')

      DO i=1,nptfr
        ncolor(i) = i
      ENDDO
      ubor(:,2) = ubor(:, 1)
      vbor(:,2) = vbor(:, 1)

      END SUBROUTINE lecsellim