writesellim.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE writesellim
!                   *****************
!
     &(nlim,lihbor,liubor,livbor,hbor,ubor,vbor,
     & chbord,nbor,npmax,nptfr)
!
!***********************************************************************
! 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 declarations_special
      USE declarations_stbtel, ONLY: out_format, typ_bnd_elem
!
      IMPLICIT NONE
!
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NLIM
      INTEGER, INTENT(IN)    :: NPTFR
      INTEGER, INTENT(IN)    :: NPMAX
      INTEGER, INTENT(INOUT) :: LIUBOR(npmax),LIVBOR(npmax)
      INTEGER, INTENT(INOUT) :: LIHBOR(npmax)
      INTEGER, INTENT(INOUT) :: NBOR(npmax)
      DOUBLE PRECISION,  INTENT(INOUT) :: UBOR(npmax),VBOR(npmax)
      DOUBLE PRECISION,  INTENT(INOUT) :: HBOR(npmax),CHBORD(npmax)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I, IERR
      INTEGER, ALLOCATABLE :: NCOLOR(:)

      ALLOCATE(ncolor(nptfr))
      DO i=1,nptfr
        ncolor(i) = i
      ENDDO
!
      CALL set_bnd(out_format, nlim, typ_bnd_elem, nptfr, 1, nbor,
     &             nptfr, lihbor, liubor, livbor, hbor, ubor, vbor,
     &             chbord, lihbor, hbor, hbor, hbor, ncolor, ierr)
      CALL check_call(ierr, 'WRITESELLIM:SET_BND')

      DEALLOCATE(ncolor)
!
      END SUBROUTINE