Sources Esope | Cast3M

coupem
C COUPEM    SOURCE    PV        20/03/24    21:16:16     10554          
      SUBROUTINE COUPEM(XYZ1,XYZ2,DIS1,DIS2,
     >                  NUM1,NUM2,MELEME,IELEM,JCOLOR)
C
C     Esclave de COUPE
C     Ajout d'un element compose des noeuds NUM1 et NUM2 et
C           de l'interpolation entre XYZ1 et XYZ2 avec les poids
C           DIS1 et DIS2
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
-INC SMELEME
C
      DIMENSION XYZ1(4),XYZ2(4)
C
      segact mcoord*mod
      NBPTI=nbpts
      NBPTS=NBPTI+1
      SEGADJ,MCOORD
C
      COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS2))
      NBPTI=NBPTI+1
      IREF=(NBPTI-1)*4
      DO IE1=1,4
        XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ2(IE1)-XYZ1(IE1))
      ENDDO
C
      IELEM=IELEM+1
      NUM(1,IELEM)=NUM1
      NUM(2,IELEM)=NUM2
      NUM(3,IELEM)=NBPTS
      ICOLOR(IELEM)=JCOLOR
C
      RETURN
      END