Télécharger la plaquette Cast3M


fr en

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