comp_fac.f

Go to the documentation of this file.

!                   *******************
                    SUBROUTINE comp_fac
!                   *******************
!
     &(eltseg,oriseg,ifabor,nelem,npoin,ifac)
!
!***********************************************************************
! BIEF   V7P1
!***********************************************************************
!
!brief    COMPLETES THE ARRAY IFAC FOR QUADRATIC POINTS
!+                AT THE INTERFACE BETWEEN 2 SUBDOMAINS.
!+
!
!history  J-M HERVOUET (LNHE)
!+        24/10/08
!+        V5P9
!+   First version
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!history  J-M HERVOUET (EDF LAB, LNHE)
!+        30/07/2015
!+        V7P1
!+   FAC replaced by IFAC
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ELTSEG         |-->| GIVES THE SEGMENT NUMBER OF EDGES OF ELEMENTS
!| ORISEG         |-->| ORIENTATION OF SEGMENTS (1 OR 2).
!| FAC            |<->| COEFFICIENT FOR COMPUTING DOT PRODUCTS IN //
!| IFABOR         |-->| -2 MEANS INTERFACE WITH ANOTHER SUB-DOMAIN
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NPOIN          |-->| NUMBER OF POINTS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_comp_fac => comp_fac
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)    :: NELEM,NPOIN
      INTEGER, INTENT(IN)    :: IFABOR(nelem,3),ELTSEG(nelem,3)
      INTEGER, INTENT(IN)    :: ORISEG(nelem,3)
      TYPE(bief_obj), INTENT(INOUT) :: IFAC
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM
!
!-----------------------------------------------------------------------
!
      DO ielem=1,nelem
!
!       ALL INITIALISED WITH 1
!
        ifac%I(npoin+eltseg(ielem,1))=1
        ifac%I(npoin+eltseg(ielem,2))=1
        ifac%I(npoin+eltseg(ielem,3))=1
!
!       ON PARALLEL INTERFACES, SEGMENTS WITH ORISEG=2 ARE SET TO 0
!       THE SEGMENT ON THE OTHER SIDE WILL HAVE 1
!
        IF(ifabor(ielem,1).EQ.-2.AND.oriseg(ielem,1).EQ.2) THEN
          ifac%I(npoin+eltseg(ielem,1))=0
        ENDIF
        IF(ifabor(ielem,2).EQ.-2.AND.oriseg(ielem,2).EQ.2) THEN
          ifac%I(npoin+eltseg(ielem,2))=0
        ENDIF
        IF(ifabor(ielem,3).EQ.-2.AND.oriseg(ielem,3).EQ.2) THEN
          ifac%I(npoin+eltseg(ielem,3))=0
        ENDIF
!
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END