filpol.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE filpol
!                   *****************
!
     &( f , c , xsom , ysom , nsom , mesh )
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    THE POINTS OF VECTOR F, WHICH ARE INSIDE OF THE
!+                POLYGON DEFINED BY VERTICES XSOM AND YSOM, ARE
!+                INITIALISED TO CONSTANT C.
!
!note     THE POLYGON IS READ ANTI-CLOCKWISE.
!
!history  J-M HERVOUET (LNH)
!+        08/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              |-->| CONSTANT TO BE PRESCRIBED INSIDE THE POLYGON.
!| F              |<->| THE RESULT (VECTOR)
!| MESH           |-->| MESH STRUCTURE
!| NSOM           |-->| NUMBER OF VERTICES IN THE POLYGON
!| XSOM           |-->| ABSCISSAE OF VERTICES IN THE POLYGON
!| YSOM           |-->| ORDINATES OF VERTICES IN THE POLYGON
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_filpol => filpol
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NSOM
      DOUBLE PRECISION, INTENT(IN) :: C,XSOM(*),YSOM(*)
      TYPE(bief_obj), INTENT(INOUT) :: F
      TYPE(bief_mesh), INTENT(IN) :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER NPOIN,IELM
!
!-----------------------------------------------------------------------
!
!  TYPE OF ELEMENT FOR VELOCITY
!
      ielm = f%ELM
!
!-----------------------------------------------------------------------
!
      IF(ielm.EQ.10) THEN
!
        CALL filp10( f%R , c , xsom , ysom , nsom ,
     &       mesh%X%R,mesh%Y%R,mesh%NELEM,mesh%NELMAX,mesh%IKLE%I)
!
      ELSEIF(ielm.EQ.11) THEN
!
        npoin = f%DIM1
        CALL filp11(f%R,c,xsom,ysom,nsom,mesh%X%R,mesh%Y%R,npoin)
!
      ELSEIF(ielm.EQ.12) THEN
!
        npoin = f%DIM1 - mesh%NELEM
        CALL filp12( f%R , c , xsom , ysom , nsom , mesh%X%R , mesh%Y%R,
     &               npoin , mesh%NELEM , mesh%NELMAX , mesh%IKLE%I)
!
      ELSE
!
!       ELM NOT IMPLEMENTED: ERROR
!
        WRITE(lu,101) ielm
101     FORMAT(1x,'FILPOL (BIEF): IELM = ',1i6,' ELEMENT NOT AVAILABLE')
        CALL plante(1)
        stop
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END