filp10.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE filp10
!                   *****************
!
     &( f , c , xsom , ysom , nsom , x , y , nelem , nelmax , ikle )
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    INITIALISES A FUNCTION TO A CONSTANT VALUE
!+                INSIDE OF A POLYGON.
!
!history  C MOULIN (LNH)
!+        06/12/94
!+        V5P1
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| C              |-->| THE CONSTANT VALUE
!| F              |<->| THE FUNCTION
!| IKLE           |-->| CONNECTIVITY TABLE.
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
!| NSOM           |-->| NUMBER OF VERTICES IN THE POLYGON
!| X              |-->| ABSCISSAE OF POINTS IN THE MESH
!| XSOM           |-->| ABSCISSAE OF POINTS IN THE POLYGON
!| Y              |-->| ORDINATES OF POINTS IN THE MESH
!| YSOM           |-->| ORDINATES OF POINTS IN THE POLYGON
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_filp10 => filp10
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN)             :: NSOM,NELEM,NELMAX
      DOUBLE PRECISION, INTENT(INOUT) :: F(*)
      DOUBLE PRECISION, INTENT(IN)    :: X(*) , Y(*)
      DOUBLE PRECISION, INTENT(IN)    :: XSOM(nsom) , YSOM(nsom) , C
      INTEGER, INTENT(IN)             :: IKLE(nelmax,3)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM , I1 , I2 , I3
!
!-----------------------------------------------------------------------
!
      DO ielem = 1 , nelem
!
        i1 = ikle(ielem,1)
        i2 = ikle(ielem,2)
        i3 = ikle(ielem,3)
        IF( inpoly(x(i1),y(i1),xsom,ysom,nsom) .AND.
     &      inpoly(x(i2),y(i2),xsom,ysom,nsom) .AND.
     &      inpoly(x(i3),y(i3),xsom,ysom,nsom) ) f(ielem) = c
!
      ENDDO ! IELEM
!
!-----------------------------------------------------------------------
!
      RETURN
      END