Geology Reference
In-Depth Information
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 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 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 29 I=1,6
DO 30 J=1,6
AS(I,J)=A(I,J)
30 CONTINUE
29 CONTINUE
C Solve for second terms in third solution.
CALL LINSOL(A,YA,6,CAUG,DET,6,13)
YSCAL(1,3)=YA(1)
YSCAL(2,3)=YA(2)
YSCAL(3,3)=YA(3)
YSCAL(4,3)=YA(4)
YSCAL(5,3)=YA(5)
YSCAL(6,3)=YA(6)
C Set right-hand side for calculation of third terms in first solution.
YA(1)=0.D0
YA(2)=RHO(1)*(-(4.D0*GAMMA+ANGS+2.D0*WES)*YSCAL(1,1)
1 +(AN*(AN+1.D0)*GAMMA+COR)*YSCAL(3,1)-YSCAL(6,1))
YA(3)=0.D0
YA(4)=RHO(1)*((GAMMA+COR/(AN*(AN+1.D0)))*YSCAL(1,1)
1 -(ANGS-COR/(AN*(AN+1.D0)))*YSCAL(3,1)-YSCAL(5,1))
YA(5)=0.D0
YA(6)=0.D0
C Copy stored matrix AS into matrix A.
DO 31 I=1,6
DO 32 J=1,6
A(I,J)=AS(I,J)
32 CONTINUE
31 CONTINUE
C Solve for third terms in first solution.
CALL LINSOL(A,YA,6,CAUG,DET,6,13)
YSCAL2(1,1)=YA(1)
YSCAL2(2,1)=YA(2)
YSCAL2(3,1)=YA(3)
YSCAL2(4,1)=YA(4)
YSCAL2(5,1)=YA(5)
YSCAL2(6,1)=YA(6)
C Set right-hand side for calculation of third terms in second solution.
YA(1)=0.D0
Search WWH ::
Custom Search