Télécharger coupei.eso

Retour à la liste

Numérotation des lignes :

coupei
  1. C COUPEI SOURCE PV 20/03/24 21:16:10 10554
  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.  
  14. -INC PPARAM
  15. -INC CCOPTIO
  16. -INC SMCOORD
  17. -INC SMELEME
  18. C
  19. DIMENSION XYZ0(4),XYZ1(4),XYZ2(4),XYZ3(4)
  20. C
  21. segact mcoord*mod
  22. NBPTI=nbpts
  23. NBPTS=NBPTI+3
  24. SEGADJ,MCOORD
  25. C
  26. COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS1))
  27. NBPTI=NBPTI+1
  28. IREF=(NBPTI-1)*4
  29. DO IE1=1,4
  30. XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ1(IE1)-XYZ0(IE1))
  31. ENDDO
  32. C
  33. COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS2))
  34. NBPTI=NBPTI+1
  35. IREF=(NBPTI-1)*4
  36. DO IE1=1,4
  37. XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ2(IE1)-XYZ0(IE1))
  38. ENDDO
  39. C
  40. COEFF=ABS(DIS0)/(ABS(DIS0)+ABS(DIS3))
  41. NBPTI=NBPTI+1
  42. IREF=(NBPTI-1)*4
  43. DO IE1=1,4
  44. XCOOR(IREF+IE1)=XYZ0(IE1)+COEFF*(XYZ3(IE1)-XYZ0(IE1))
  45. ENDDO
  46. C
  47. NBPTI=NBPTI-3
  48. IELEM=IELEM+1
  49. DO IE1=1,3
  50. NUM(IE1,IELEM)=NBPTI+IE1
  51. ENDDO
  52. ICOLOR(IELEM)=JCOLOR
  53. C
  54. RETURN
  55. END
  56.  
  57.  
  58.  
  59.  

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