Esope source files

   1 : C COUPEL    SOURCE    CHAT      05/01/12    22:24:48     5004           
   2 :       SUBROUTINE COUPEL(XYZ1,XYZ2,XYZ3,DIS1,DIS2,DIS3,
   3 :      >                  NUM0,MELEME,IELEM,JCOLOR)
   4 : C     Esclave de COUPEK
   5 : C     Ajout d'un element compose du noeud NUM0
   6 : C           de l'interpolation entre XYZ1 et XYZ2 avec les poids
   7 : C           DIS1 et DIS2 et de de l'interpolation entre XYZ1 et XYZ3
   8 : C           avec les poids DIS1 et DIS3
   9 : C
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)
  12 : -INC CCOPTIO
  13 : -INC SMCOORD
  14 : -INC SMELEME
  15 : C
  16 :       DIMENSION XYZ1(4),XYZ2(4),XYZ3(4)
  17 : C
  18 :       NBPTI=XCOOR(/1)/4
  19 :       NBPTS=NBPTI+2
  20 :       SEGADJ,MCOORD
  21 : C
  22 :       COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS2))
  23 :       NBPTI=NBPTI+1
  24 :       IREF=(NBPTI-1)*4
  25 :       DO IE1=1,4
  26 :         XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ2(IE1)-XYZ1(IE1))
  27 :       ENDDO
  28 : C
  29 :       COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS3))
  30 :       NBPTI=NBPTI+1
  31 :       IREF=(NBPTI-1)*4
  32 :       DO IE1=1,4
  33 :         XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ3(IE1)-XYZ1(IE1))
  34 :       ENDDO
  35 : C
  36 :       NBPTI=NBPTI-2
  37 :       IELEM=IELEM+1
  38 :       NUM(1,IELEM)=NUM0
  39 :       NUM(2,IELEM)=NBPTI+1
  40 :       NUM(3,IELEM)=NBPTI+2
  41 :       ICOLOR(IELEM)=JCOLOR             
  42 : C
  43 :       RETURN
  44 :       END
  45 :  
  46 :