Télécharger ripp9.eso

Retour à la liste

Numérotation des lignes :

ripp9
  1. C RIPP9 SOURCE CHAT 05/01/13 03:00:48 5004
  2. SUBROUTINE RIPP9(RE,E,RAI,X13,X23,Y13,HM,A)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION E(3,1)
  6. DIMENSION RE(9,1),AT(3,3),T(3,3),H(3,3),PL(9,9),REQ(9,9)
  7. T(1,1)=Y13*Y13
  8. T(1,2)=T(1,1)
  9. T(1,3)=-2.*T(1,1)
  10. T(2,1)=X23**2
  11. T(2,2)=X13**2
  12. T(2,3)=-2.*X13*X23
  13. T(3,1)=-2.*X23*Y13
  14. T(3,2)=-2.*X13*Y13
  15. T(3,3)=-(T(3,1)+T(3,2))
  16. CALL PROSY(3,3,H,E,T)
  17. AT(1,1)=H(1,1)-H(1,2)-H(1,3)+H(2,2)-H(2,3)+H(3,3)
  18. AT(1,2)=-H(1,1)-H(1,2)-H(1,3)+2.*H(2,3)+H(3,3)
  19. AT(1,3)=-H(1,2)+2.*H(1,3)-H(2,2)-H(2,3)+H(3,3)
  20. AT(2,2)=H(1,1)+3.*H(1,2)+3.*H(1,3)+9.*H(2,2)+9.*H(2,3)+3.*H(3,3)
  21. AT(2,3)=-5.*H(1,2)-2.*H(1,3)-2.*H(2,3)-H(3,3)
  22. AT(3,3)=9.*H(1,1)+3.*H(1,2)+9.*H(1,3)+H(2,2)+3.*H(2,3)+3.*H(3,3)
  23. AD=HM*HM*HM*RAI/(288.*A*A*A)
  24. DO3I=1,3
  25. DO3J=I,9
  26. 3 REQ(I,J)=0.
  27. DO4I=7,9
  28. II=I-6
  29. DO4J=I,9
  30. 4 REQ(I,J)=AT(II,J-6)*AD
  31. AE=1.5*AD
  32. AT(1,1)=H(3,3)
  33. AT(1,2)=-2.*H(2,3)-H(3,3)
  34. AT(1,3)=-2.*H(1,3)-H(3,3)
  35. AT(2,2)=4.*(H(2,2)+H(2,3))+H(3,3)
  36. AT(2,3)=4.*H(1,2)+2.*(H(1,3)+H(2,3))+H(3,3)
  37. AT(3,3)=4.*(H(1,1)+H(1,3))+H(3,3)
  38. DO 5 I=1,3
  39. II=I+3
  40. DO 5 J=I,3
  41. 5 REQ(II,J+3)=AT(I,J)*AE
  42. AT(1,1)=-H(1,3)+H(2,3)
  43. AT(2,1)=2.*(H(1,2)-H(2,2))+H(1,3)-H(2,3)
  44. AT(3,1)=2.*(H(1,1)-H(1,2))+H(1,3)-H(2,3)
  45. AT(1,2)=H(1,3)+H(3,3)
  46. AT(2,2)=-2.*(H(1,2)+H(2,3))-H(1,3)-H(3,3)
  47. AT(3,2)=-2.*H(1,1)-3.*H(1,3)-H(3,3)
  48. AT(1,3)=-H(2,3)-H(3,3)
  49. AT(2,3)=2.*H(2,2)+3.*H(2,3)+H(3,3)
  50. AT(3,3)=2.*(H(1,2)+H(1,3))+H(2,3)+H(3,3)
  51. DO 6 I=1,3
  52. II=I+3
  53. DO 6 J=1,3
  54. 6 REQ(II,J+6)=AT(I,J)*AD
  55. DO9I=1,9
  56. DO9J=1,9
  57. 9 PL(I,J)=0.
  58. PL(1,7)=1.
  59. PL(2,1)=1.
  60. PL(2,7)=-1.
  61. PL(3,4)=1.
  62. PL(3,7)=-1.
  63. PL(4,2)=0.
  64. PL(4,3)=0.5*(X13-X23)
  65. PL(4,5)=-PL(4,2)
  66. PL(4,6)=-PL(4,3)
  67. PL(5,5)=-0.5*Y13
  68. PL(5,6)=0.5*X23
  69. PL(5,8)=-PL(5,5)
  70. PL(5,9)=-PL(5,6)
  71. PL(6,2)=-0.5*Y13
  72. PL(6,3)=0.5*X13
  73. PL(6,8)=-PL(6,2)
  74. PL(6,9)=-PL(6,3)
  75. PL(7,1)=-1.
  76. PL(7,2)=PL(4,5)
  77. PL(7,3)=PL(4,6)
  78. PL(7,4)=1.
  79. PL(7,5)=PL(7,2)
  80. PL(7,6)=PL(7,3)
  81. PL(8,4)=-1.
  82. PL(8,5)=PL(5,8)
  83. PL(8,6)=PL(5,9)
  84. PL(8,7)=1.
  85. PL(8,8)=PL(8,5)
  86. PL(8,9)=PL(8,6)
  87. PL(9,1)=1.
  88. PL(9,2)=PL(6,2)
  89. PL(9,3)=PL(6,3)
  90. PL(9,7)=-1.
  91. PL(9,8)=PL(9,2)
  92. PL(9,9)=PL(9,3)
  93. DO11I=1,9
  94. DO11J=I,9
  95. 11 REQ(J,I)=REQ(I,J)
  96. CALL PROSY(9,9,RE,REQ,PL)
  97. RETURN
  98. END
  99.  
  100.  

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