Télécharger la plaquette Cast3M


fr en

Télécharger kdom3.eso

Retour à la liste

Numérotation des lignes : 

kdom3

  1. C KDOM3     SOURCE    KK2000    14/04/10    21:15:12     8032
  2.       SUBROUTINE KDOM3(IP1,IP2,IP3,IP4,IP5,IP6,XCEN,VOL)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  KDOM3
  9. C
  10. C DESCRIPTION       :  Subroutine called by KDOM2
  11. C                      Given a pyramid with base (IP1,IP2,IP3,IP4),
  12. C                      IP5 = center of the base
  13. C                      IP6 = vertex of the pyramid
  14. C                      this subroutine compute the volume of the 4
  15. C                      tetrahedra (IP1,IP2,IP5,IP6),(IP2,IP3,IP5,IP6),
  16. C                                 (IP3,IP4,IP5,IP6),(IP4,IP1,IP5,IP6)
  17. C                      and its mass center.
  18. C                      VOL can be positive or negative depending on
  19. C                      the orientation of the vectors
  20. C                      (see KDOM4)
  21. C
  22. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  23. C
  24. C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
  25. C
  26. C************************************************************************
  27. C
  28. C
  29.  
  30.       IMPLICIT INTEGER(I-N)
  31.       IMPLICIT REAL*8(A-H,O-Z)
  32.  
  33. C
  34.       INTEGER IP1,IP2,IP3,IP4,IP5,IP6, I1
  35.       REAL*8 XCEN(3), VOL
  36.      &     ,X1(3),X2(3),X3(3),X4(3),VOL1,VOL2,VOL3,VOL4
  37.  
  38. -INC PPARAM
  39. -INC CCOPTIO
  40. C     (IP1,IP2,IP5,IP6)
  41.       CALL KDOM4(IP1,IP2,IP5,IP6,X1,VOL1)
  42. C     (IP2,IP3,IP5,IP6)
  43.       CALL KDOM4(IP2,IP3,IP5,IP6,X2,VOL2)
  44. C     (IP3,IP4,IP5,IP6)
  45.       CALL KDOM4(IP3,IP4,IP5,IP6,X3,VOL3)
  46. C     (IP4,IP1,IP5,IP6)
  47.       CALL KDOM4(IP4,IP1,IP5,IP6,X4,VOL4)
  48. C
  49.       VOL=VOL1+VOL2+VOL3+VOL4
  50.       DO I1=1,3,1
  51.          XCEN(I1)=((VOL1*X1(I1))+(VOL2*X2(I1))+(VOL3*X3(I1))
  52.      &        +(VOL4*X4(I1)))/VOL
  53.       ENDDO
  54. C
  55.       RETURN
  56. C
  57.       END
  58.  
  59.  
  60.  
  61.  

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