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)) 15 : CALL PROSY(3,3,H,E,T) 16 : AT(1,1)=H(1,1)-H(1,2)-H(1,3)+H(2,2)-H(2,3)+H(3,3) 17 : AT(1,2)=-H(1,1)-H(1,2)-H(1,3)+2.*H(2,3)+H(3,3) 18 : AT(1,3)=-H(1,2)+2.*H(1,3)-H(2,2)-H(2,3)+H(3,3) 19 : 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) 20 : AT(2,3)=-5.*H(1,2)-2.*H(1,3)-2.*H(2,3)-H(3,3) 21 : 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) 22 : AD=HM*HM*HM*RAI/(288.*A*A*A) 23 : DO3I=1,3 24 : DO3J=I,9 25 : 3 REQ(I,J)=0. 26 : DO4I=7,9 27 : II=I-6 28 : DO4J=I,9 29 : 4 REQ(I,J)=AT(II,J-6)*AD 30 : AE=1.5*AD 31 : AT(1,1)=H(3,3) 32 : AT(1,2)=-2.*H(2,3)-H(3,3) 33 : AT(1,3)=-2.*H(1,3)-H(3,3) 34 : AT(2,2)=4.*(H(2,2)+H(2,3))+H(3,3) 35 : AT(2,3)=4.*H(1,2)+2.*(H(1,3)+H(2,3))+H(3,3) 36 : AT(3,3)=4.*(H(1,1)+H(1,3))+H(3,3) 37 : DO 5 I=1,3 38 : II=I+3 39 : DO 5 J=I,3 40 : 5 REQ(II,J+3)=AT(I,J)*AE 41 : AT(1,1)=-H(1,3)+H(2,3) 42 : AT(2,1)=2.*(H(1,2)-H(2,2))+H(1,3)-H(2,3) 43 : AT(3,1)=2.*(H(1,1)-H(1,2))+H(1,3)-H(2,3) 44 : AT(1,2)=H(1,3)+H(3,3) 45 : AT(2,2)=-2.*(H(1,2)+H(2,3))-H(1,3)-H(3,3) 46 : AT(3,2)=-2.*H(1,1)-3.*H(1,3)-H(3,3) 47 : AT(1,3)=-H(2,3)-H(3,3) 48 : AT(2,3)=2.*H(2,2)+3.*H(2,3)+H(3,3) 49 : AT(3,3)=2.*(H(1,2)+H(1,3))+H(2,3)+H(3,3) 50 : DO 6 I=1,3 51 : II=I+3 52 : DO 6 J=1,3 53 : 6 REQ(II,J+6)=AT(I,J)*AD 54 : DO9I=1,9 55 : DO9J=1,9 56 : 9 PL(I,J)=0. 57 : PL(1,7)=1. 58 : PL(2,1)=1. 59 : PL(2,7)=-1. 60 : PL(3,4)=1. 61 : PL(3,7)=-1. 62 : PL(4,2)=0. 63 : PL(4,3)=0.5*(X13-X23) 64 : PL(4,5)=-PL(4,2) 65 : PL(4,6)=-PL(4,3) 66 : PL(5,5)=-0.5*Y13 67 : PL(5,6)=0.5*X23 68 : PL(5,8)=-PL(5,5) 69 : PL(5,9)=-PL(5,6) 70 : PL(6,2)=-0.5*Y13 71 : PL(6,3)=0.5*X13 72 : PL(6,8)=-PL(6,2) 73 : PL(6,9)=-PL(6,3) 74 : PL(7,1)=-1. 75 : PL(7,2)=PL(4,5) 76 : PL(7,3)=PL(4,6) 77 : PL(7,4)=1. 78 : PL(7,5)=PL(7,2) 79 : PL(7,6)=PL(7,3) 80 : PL(8,4)=-1. 81 : PL(8,5)=PL(5,8) 82 : PL(8,6)=PL(5,9) 83 : PL(8,7)=1. 84 : PL(8,8)=PL(8,5) 85 : PL(8,9)=PL(8,6) 86 : PL(9,1)=1. 87 : PL(9,2)=PL(6,2) 88 : PL(9,3)=PL(6,3) 89 : PL(9,7)=-1. 90 : PL(9,8)=PL(9,2) 91 : PL(9,9)=PL(9,3) 92 : DO11I=1,9 93 : DO11J=I,9 94 : 11 REQ(J,I)=REQ(I,J) 94 : CALL PROSY(9,9,RE,REQ,PL) 95 : RETURN 96 : END 97 :
© Cast3M 2003 - All rights reserved.
Disclaimer