Esope source files

   1 : C COUPEI    SOURCE    CHAT      05/01/12    22:24:37     5004           
   2 :       SUBROUTINE COUPEI(XYZ0,XYZ1,XYZ2,XYZ3,DIS0,DIS1,DIS2,DIS3,
   3 :      >                  MELEME,IELEM,JCOLOR)
   4 : C
   5 : C     Esclave de COUPE
   6 : C     Ajout d'un element correspondant a la coupe par la pointe
   7 : C           du tetraedre par un plan XYZ0 et DIS0 sont les ccordonnees
   8 : C           de la pointe et la distance au plan. Les autres correspondent
   9 : C           aux autres points
  10 : C
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 : -INC CCOPTIO
  14 : -INC SMCOORD
  15 : -INC SMELEME
  16 : C
  17 :       DIMENSION XYZ0(4),XYZ1(4),XYZ2(4),XYZ3(4)
  18 : C
  19 :       NBPTI=XCOOR(/1)/4
  20 :       NBPTS=NBPTI+3
  21 :       SEGADJ,MCOORD
  22 : C
  23 :       COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS1))
  24 :       NBPTI=NBPTI+1
  25 :       IREF=(NBPTI-1)*4
  26 :       DO IE1=1,4
  27 :         XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ1(IE1)-XYZ0(IE1))
  28 :       ENDDO
  29 : C
  30 :       COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS2))
  31 :       NBPTI=NBPTI+1
  32 :       IREF=(NBPTI-1)*4
  33 :       DO IE1=1,4
  34 :         XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ2(IE1)-XYZ0(IE1))
  35 :       ENDDO
  36 : C
  37 :       COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS3))
  38 :       NBPTI=NBPTI+1
  39 :       IREF=(NBPTI-1)*4
  40 :       DO IE1=1,4
  41 :         XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ3(IE1)-XYZ0(IE1))
  42 :       ENDDO
  43 : C
  44 :       NBPTI=NBPTI-3
  45 :       IELEM=IELEM+1
  46 :       DO IE1=1,3
  47 :         NUM(IE1,IELEM)=NBPTI+IE1
  48 :       ENDDO
  49 :       ICOLOR(IELEM)=JCOLOR             
  50 : C
  51 :       RETURN
  52 :       END
  53 :  
  54 :