1 : C COUPEM SOURCE CHAT 05/01/12 22:24:51 5004 2 : SUBROUTINE COUPEM(XYZ1,XYZ2,DIS1,DIS2, 3 : > NUM1,NUM2,MELEME,IELEM,JCOLOR) 4 : C 5 : C Esclave de COUPE 6 : C Ajout d'un element compose des noeuds NUM1 et NUM2 et 7 : C de l'interpolation entre XYZ1 et XYZ2 avec les poids 8 : C DIS1 et DIS2 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) 17 : C 18 : NBPTI=XCOOR(/1)/4 19 : NBPTS=NBPTI+1 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 : IELEM=IELEM+1 30 : NUM(1,IELEM)=NUM1 31 : NUM(2,IELEM)=NUM2 32 : NUM(3,IELEM)=NBPTS 33 : ICOLOR(IELEM)=JCOLOR 34 : C 35 : RETURN 36 : END 37 : 38 :
© Cast3M 2003 - All rights reserved.
Disclaimer