C KDOM3     SOURCE    KK2000    14/04/10    21:15:12     8032
      SUBROUTINE KDOM3(IP1,IP2,IP3,IP4,IP5,IP6,XCEN,VOL)
C
C************************************************************************
C
C PROJET            :  CASTEM 2000
C
C NOM               :  KDOM3
C
C DESCRIPTION       :  Subroutine called by KDOM2
C                      Given a pyramid with base (IP1,IP2,IP3,IP4),
C                      IP5 = center of the base
C                      IP6 = vertex of the pyramid
C                      this subroutine compute the volume of the 4
C                      tetrahedra (IP1,IP2,IP5,IP6),(IP2,IP3,IP5,IP6),
C                                 (IP3,IP4,IP5,IP6),(IP4,IP1,IP5,IP6)
C                      and its mass center.
C                      VOL can be positive or negative depending on
C                      the orientation of the vectors
C                      (see KDOM4)
C
C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
C
C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
C
C************************************************************************
C
C

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)

C
      INTEGER IP1,IP2,IP3,IP4,IP5,IP6, I1
      REAL*8 XCEN(3), VOL
     &     ,X1(3),X2(3),X3(3),X4(3),VOL1,VOL2,VOL3,VOL4

-INC PPARAM
-INC CCOPTIO
C     (IP1,IP2,IP5,IP6)
      CALL KDOM4(IP1,IP2,IP5,IP6,X1,VOL1)
C     (IP2,IP3,IP5,IP6)
      CALL KDOM4(IP2,IP3,IP5,IP6,X2,VOL2)
C     (IP3,IP4,IP5,IP6)
      CALL KDOM4(IP3,IP4,IP5,IP6,X3,VOL3)
C     (IP4,IP1,IP5,IP6)
      CALL KDOM4(IP4,IP1,IP5,IP6,X4,VOL4)
C
      VOL=VOL1+VOL2+VOL3+VOL4
      DO I1=1,3,1
         XCEN(I1)=((VOL1*X1(I1))+(VOL2*X2(I1))+(VOL3*X3(I1))
     &        +(VOL4*X4(I1)))/VOL
      ENDDO
C
      RETURN
C
      END