Télécharger geocst.eso

Retour à la liste

Numérotation des lignes :

  1. C GEOCST SOURCE CHAT 05/01/13 00:17:47 5004
  2. SUBROUTINE GEOCST(XEL,GEOM)
  3. C=======================================================================
  4. C
  5. C CALCULE DES CONSTANTES D UN TRIANGLE POUR LE DKT
  6. C
  7. C ENTREES
  8. C XEL(3,3)= COORDONNEES DE L ELEMENT
  9. C SORTIES
  10. C GEOM(17)=DES CONSTANTES POUR LE DKT
  11. C
  12. C=======================================================================
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. DIMENSION XEL(3,*),GEOM(*)
  16. DATA TROIS/3.D0/,SIX/6.D0/
  17. X12=XEL(1,1)-XEL(1,2)
  18. X23=XEL(1,2)-XEL(1,3)
  19. X31=XEL(1,3)-XEL(1,1)
  20. Y12=XEL(2,1)-XEL(2,2)
  21. Y23=XEL(2,2)-XEL(2,3)
  22. Y31=XEL(2,3)-XEL(2,1)
  23. AIREDO=-Y31*X12+Y12*X31
  24. VJ11=Y31/AIREDO
  25. VJ12=Y12/AIREDO
  26. VJ21=-X31/AIREDO
  27. VJ22=-X12/AIREDO
  28. XSQ12=X12*X12
  29. XSQ23=X23*X23
  30. XSQ31=X31*X31
  31. YSQ12=Y12*Y12
  32. YSQ23=Y23*Y23
  33. YSQ31=Y31*Y31
  34. DIST12=XSQ12+YSQ12
  35. DIST23=XSQ23+YSQ23
  36. DIST31=XSQ31+YSQ31
  37. CL4=-SIX*X23/DIST23
  38. CL5=-SIX*X31/DIST31
  39. CL6=-SIX*X12/DIST12
  40. SL4=-SIX*Y23/DIST23
  41. SL5=-SIX*Y31/DIST31
  42. SL6=-SIX*Y12/DIST12
  43. D4=TROIS*X23*Y23/DIST23
  44. D5=TROIS*X31*Y31/DIST31
  45. D6=TROIS*X12*Y12/DIST12
  46. SU4=TROIS*YSQ23/DIST23
  47. SU5=TROIS*YSQ31/DIST31
  48. SU6=TROIS*YSQ12/DIST12
  49. GEOM(1)=CL4
  50. GEOM(2)=CL5
  51. GEOM(3)=CL6
  52. GEOM(4)=SL4
  53. GEOM(5)=SL5
  54. GEOM(6)=SL6
  55. GEOM(7)=D4
  56. GEOM(8)=D5
  57. GEOM(9)=D6
  58. GEOM(10)=SU4
  59. GEOM(11)=SU5
  60. GEOM(12)=SU6
  61. GEOM(13)=VJ11
  62. GEOM(14)=VJ12
  63. GEOM(15)=VJ21
  64. GEOM(16)=VJ22
  65. GEOM(17)=AIREDO
  66. RETURN
  67. END
  68.  
  69.  

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