Télécharger la plaquette Cast3M


fr en

Télécharger coupej.eso

Retour à la liste

Numérotation des lignes : 

coupej

  1. C COUPEJ    SOURCE    PV        20/03/24    21:16:12     10554          
  2.       SUBROUTINE COUPEJ(XYZ1,XYZ2,XYZ3,XYZ4,DIS1,DIS2,DIS3,DIS4,
  3.      >                  MELEME,IELEM,JCOLOR)
  4. C
  5. C     Esclave de COUPE
  6. C     Ajout de 2 elements correspondant a la coupe par un cote
  7. C           du tetraedre par un plan. XYZ1, DIS1 et 2 sont les ccordonnees
  8. C           de la pointe et la distance au plan pour les points du premier
  9. C           segment. 3 et 4  correspondent 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 XYZ1(4),XYZ2(4),XYZ3(4),XYZ4(4)
  20. C
  21.       segact mcoord*mod
  22.       NBPTI=nbpts
  23.       NBPTS=NBPTI+4
  24.       SEGADJ,MCOORD
  25. C
  26.       COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS3))
  27.       NBPTI=NBPTI+1
  28.       IREF=(NBPTI-1)*4
  29.       DO IE1=1,4
  30.         XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ3(IE1)-XYZ1(IE1))
  31.       ENDDO
  32. C
  33.       COEFF=ABS(DIS1)/(ABS(DIS1)+ABS(DIS4))
  34.       NBPTI=NBPTI+1
  35.       IREF=(NBPTI-1)*4
  36.       DO IE1=1,4
  37.         XCOOR(IREF+IE1)=XYZ1(IE1)+COEFF*(XYZ4(IE1)-XYZ1(IE1))
  38.       ENDDO
  39. C
  40.       COEFF=ABS(DIS2)/(ABS(DIS2)+ABS(DIS3))
  41.       NBPTI=NBPTI+1
  42.       IREF=(NBPTI-1)*4
  43.       DO IE1=1,4
  44.         XCOOR(IREF+IE1)=XYZ2(IE1)+COEFF*(XYZ3(IE1)-XYZ2(IE1))
  45.       ENDDO
  46. C
  47.       COEFF=ABS(DIS2)/(ABS(DIS2)+ABS(DIS4))
  48.       NBPTI=NBPTI+1
  49.       IREF=(NBPTI-1)*4
  50.       DO IE1=1,4
  51.         XCOOR(IREF+IE1)=XYZ2(IE1)+COEFF*(XYZ4(IE1)-XYZ2(IE1))
  52.       ENDDO
  53. C
  54.       NBPTI=NBPTI-4
  55.       IELEM=IELEM+1
  56.       NUM(1,IELEM)=NBPTI+1
  57.       NUM(2,IELEM)=NBPTI+2
  58.       NUM(3,IELEM)=NBPTI+4
  59.       ICOLOR(IELEM)=JCOLOR
  60.       IELEM=IELEM+1
  61.       NUM(1,IELEM)=NBPTI+1
  62.       NUM(2,IELEM)=NBPTI+4
  63.       NUM(3,IELEM)=NBPTI+3
  64.       ICOLOR(IELEM)=JCOLOR
  65. C
  66.       RETURN
  67.       END
  68.  
  69.  
  70.  
  71.  

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