C KDOM3     SOURCE    CHAT      05/01/13    00:54:06     5004           
      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
C**** Variables de COOPTIO
C
      IMPLICIT INTEGER(I-N)
      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
     &     ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
     &     ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
     &     ,IECHO,   IIMPI,  IOSPI
     &     ,IDIM,    IFICLE, IPREFI     
     &     ,MCOORD
     &     ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
     &     ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
     &     ,NORINC,NORVAL,NORIND,NORVAD
     &     ,NUCROU, IPSAUV, IREFOR, ISAFOR
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 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