filp11.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE filp11
!                   *****************
!
     &( f , c , xsom , ysom , nsom , x , y , npoin )
!
!***********************************************************************
! 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
!| NPOIN          |-->| NUMBER OF POINTS
!| 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_filp11 => filp11
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NSOM , NPOIN
      DOUBLE PRECISION, INTENT(INOUT) :: F(*)
      DOUBLE PRECISION, INTENT(IN) :: X(*) , Y(*)
      DOUBLE PRECISION, INTENT(IN) :: XSOM(nsom) , YSOM(nsom) , C
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
!
!-----------------------------------------------------------------------
!
      DO i = 1 , npoin
!
        IF(inpoly(x(i),y(i),xsom,ysom,nsom)) f(i) = c
!
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END