C COUPEI    SOURCE    PV        20/03/24    21:16:10     10554          
      SUBROUTINE COUPEI(XYZ0,XYZ1,XYZ2,XYZ3,DIS0,DIS1,DIS2,DIS3,
     >                  MELEME,IELEM,JCOLOR)
C
C     Esclave de COUPE
C     Ajout d'un element correspondant a la coupe par la pointe
C           du tetraedre par un plan XYZ0 et DIS0 sont les ccordonnees
C           de la pointe et la distance au plan. Les autres correspondent
C           aux autres points
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)

-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
-INC SMELEME
C
      DIMENSION XYZ0(4),XYZ1(4),XYZ2(4),XYZ3(4)
C
      segact mcoord*mod
      NBPTI=nbpts
      NBPTS=NBPTI+3
      SEGADJ,MCOORD
C
      COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS1))
      NBPTI=NBPTI+1
      IREF=(NBPTI-1)*4
      DO IE1=1,4
        XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ1(IE1)-XYZ0(IE1))
      ENDDO
C
      COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS2))
      NBPTI=NBPTI+1
      IREF=(NBPTI-1)*4
      DO IE1=1,4
        XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ2(IE1)-XYZ0(IE1))
      ENDDO
C
      COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS3))
      NBPTI=NBPTI+1
      IREF=(NBPTI-1)*4
      DO IE1=1,4
        XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ3(IE1)-XYZ0(IE1))
      ENDDO
C
      NBPTI=NBPTI-3
      IELEM=IELEM+1
      DO IE1=1,3
        NUM(IE1,IELEM)=NBPTI+IE1
      ENDDO
      ICOLOR(IELEM)=JCOLOR
C
      RETURN
      END