Télécharger kdom3.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  38. C (IP1,IP2,IP5,IP6)
  39. CALL KDOM4(IP1,IP2,IP5,IP6,X1,VOL1)
  40. C (IP2,IP3,IP5,IP6)
  41. CALL KDOM4(IP2,IP3,IP5,IP6,X2,VOL2)
  42. C (IP3,IP4,IP5,IP6)
  43. CALL KDOM4(IP3,IP4,IP5,IP6,X3,VOL3)
  44. C (IP4,IP1,IP5,IP6)
  45. CALL KDOM4(IP4,IP1,IP5,IP6,X4,VOL4)
  46. C
  47. VOL=VOL1+VOL2+VOL3+VOL4
  48. DO I1=1,3,1
  49. XCEN(I1)=((VOL1*X1(I1))+(VOL2*X2(I1))+(VOL3*X3(I1))
  50. & +(VOL4*X4(I1)))/VOL
  51. ENDDO
  52. C
  53. RETURN
  54. C
  55. END
  56.  
  57.  
  58.  
  59.  

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