Télécharger coupei.eso

Retour à la liste

Numérotation des lignes :

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

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