ifabtom.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE ifabtom
!                   ******************
!
     &(ifabor,nelem2,netage)
!
!***********************************************************************
! TELEMAC3D   V6P3                                   21/08/2010
!***********************************************************************
!
!brief    FILLS 'IFABOR' OF MESH3D FOR TOMAWAC
!
!
!history  J-M HERVOUET (LNHE)
!+        28/09/2012
!+        V6P3
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IFABOR         |<->| CORRESPONDENCE BOUNDARY FACE - 2D ELEMENT 2D
!| NELEM2         |-->| NUMBER OF ELEMENTS IN 2D
!| NETAGE         |-->| NUMBER OF PLANES - 1
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
!
      USE declarations_special
      USE interface_tomawac, ex_ifabtom => ifabtom
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)          :: NELEM2,NETAGE
      INTEGER, INTENT(INOUT)       :: IFABOR(nelem2,5,netage)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM2,IETAGE,LOC
!
!=======================================================================
!  IFABOR ARRAYS
!=======================================================================
!
!     SEE SCHAR41 FOR THE MEANING OF LOC
!     1 IS : DO NOT RECOMPUTE VELOCITIES WHEN CROSSING A MESH PLANE
!     2 IS : RECOMPUTE VELOCITIES WHEN CROSSING A MESH PLANE
!
      loc = 1
!
!  IFABOR ON THE VERTICAL FACES OF THE ELEMENTS OF THE FIRST LAYER
!
      DO ielem2 = 1,nelem2
!       ALREADY DONE BY 2D
!       IFABOR(IELEM2,1,IETAGE) = IFABOR(IELEM2,1,1)
!       IFABOR(IELEM2,2,IETAGE) = IFABOR(IELEM2,2,1)
!       IFABOR(IELEM2,3,IETAGE) = IFABOR(IELEM2,3,1)
        ifabor(ielem2,4,1) = loc
        ifabor(ielem2,5,1) = loc
      ENDDO
      IF(netage.GE.2) THEN
        DO ietage = 2,netage
          DO ielem2 = 1,nelem2
            ifabor(ielem2,1,ietage) = ifabor(ielem2,1,1)
            ifabor(ielem2,2,ietage) = ifabor(ielem2,2,1)
            ifabor(ielem2,3,ietage) = ifabor(ielem2,3,1)
            ifabor(ielem2,4,ietage) = loc
            ifabor(ielem2,5,ietage) = loc
          ENDDO
        ENDDO
      ENDIF
!
!-----------------------------------------------------------------------
!
!     IFABOR ON THE BOTTOM AND SURFACE
!
      DO ielem2 = 1,nelem2
!       -1 : SOLID
!        0 : LIQUID
!        1 : PERIODICITY (TREATED AS INTERNAL BOUNDARY)
        ifabor(ielem2,4,1)      = loc
        ifabor(ielem2,5,netage) = loc
!
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END