Télécharger coupel.eso

Retour à la liste

Numérotation des lignes :

  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.  
  47.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales