Geology Reference
In-Depth Information
Y(1,1)=1.D0
Y(2,1)=2.D0*(AN-1.D0)*MU(1)
Y(3,1)=1.D0/AN
Y(4,1)=Y(2,1)/AN
Y(5,1)=4.D0*PI*G*RHO(1)/AN
Y(6,1)=0.D0
C Set initial values of second fundamental solution.
Y(1,2)=0.D0
Y(2,2)=0.D0
Y(3,2)=0.D0
Y(4,2)=0.D0
Y(5,2)=1.D0/AN
Y(6,2)=1.D0
C Set initial values of third fundamental solution.
Y(3,3)=P2(AN,LAMBDA(1),MU(1))/P1(AN,LAMBDA(1),MU(1))
Y(2,3)=-Q1(AN,LAMBDA(1),MU(1))*Y(3,3)+Q2(AN,LAMBDA(1),MU(1))
Y(1,3)=-AN*Y(3,3)+1.D0/MU(1)
Y(4,3)=1.D0
Y(5,3)=2.D0*PI*G*RHO(1)*((AN+3.D0)*Y(1,3)-AN*(AN+1.D0)*Y(3,3))/
1 (2.D0*AN+3.D0)
Y(6,3)=(AN+2.D0)*Y(5,3)-4.D0*PI*G*RHO(1)*Y(1,3)
C Calculate coefficients of second terms in series expansions.
C Calculate coefficients of second terms
C for first fundamental solution.
YSCAL(3,1)=RHO(1)*((3.D0-AN)*GAMMA+ANGS+2.D0*WES-COR/AN)/
1 P1(AN,LAMBDA(1),MU(1))
YSCAL(2,1)=-Q1(AN,LAMBDA(1),MU(1))*YSCAL(3,1)
YSCAL(1,1)=-AN*YSCAL(3,1)
YSCAL(4,1)=0.D0
YSCAL(5,1)=2.D0*PI*G*RHO(1)*((AN+3.D0)*YSCAL(1,1)
1 -AN*(AN+1.D0)*YSCAL(3,1))/(2.D0*AN+3.D0)
YSCAL(6,1)=(AN+2.D0)*YSCAL(5,1)-4.D0*PI*G*RHO(1)*YSCAL(1,1)
C Calculate coefficients of second terms
C for second fundamental solution.
YSCAL(3,2)=RHO(1)/P1(AN,LAMBDA(1),MU(1))
YSCAL(2,2)=-Q1(AN,LAMBDA(1),MU(1))*YSCAL(3,2)
YSCAL(1,2)=-AN*YSCAL(3,2)
YSCAL(4,2)=0.D0
YSCAL(5,2)=2.D0*PI*G*RHO(1)*((AN+3.D0)*YSCAL(1,2)
1 -AN*(AN+1.D0)*YSCAL(3,2))/(2.D0*AN+3.D0)
YSCAL(6,2)=(AN+2.D0)*YSCAL(5,2)-4.D0*PI*G*RHO(1)*YSCAL(1,2)
C Generate matrix for evaluation of second terms in power series
C expansion of third fundamental solution.
CALL MATRIX(A,RHO,LAMBDA,MU,PI,G,AN,WE)
C Set right-hand side for calculation of second terms
C in third solution.
YA(1)=0.D0
YA(2)=RHO(1)*(-(4.D0*GAMMA+ANGS+2.D0*WES)*Y(1,3)
1 +(AN*(AN+1.D0)*GAMMA+COR)*Y(3,3)-Y(6,3))
YA(3)=0.D0
YA(4)=RHO(1)*((GAMMA+COR/(AN*(AN+1.D0)))*Y(1,3)
1 -(ANGS-COR/(AN*(AN+1.D0)))*Y(3,3)-Y(5,3))
YA(5)=0.D0
YA(6)=0.D0
C Store matrix A in AS.
DO 27 I=1,6
DO 28 J=1,6
AS(I,J)=A(I,J)
28
CONTINUE
27
CONTINUE
Search WWH ::
Custom Search