Earth deformations - Earth Dynamics: Deformations and Oscillations of the Rotating Earth - page 268

1 +(AN*(AN+1.D0)*GAMMA+COR)*YSCAL(3,3)-YSCAL(6,3))
YA(3)=0.D0
YA(4)=RHO(1)*((GAMMA+COR/(AN*(AN+1.D0)))*YSCAL(1,3)
1 -(ANGS-COR/(AN*(AN+1.D0)))*YSCAL(3,3)-YSCAL(5,3))
YA(5)=0.D0
YA(6)=0.D0
C Copy stored matrix AS into matrix A and add 2.0 to diagonal of A.
DO 33 I=1,6
DO 34 J=1,6
A(I,J)=AS(I,J)
IF(I.EQ.J)A(I,J)=A(I,J)+2.D0
34 CONTINUE
33 CONTINUE
C Solve for third terms in third solution.
CALL LINSOL(A,YA,6,CAUG,DET,6,13)
YSCAL2(1,3)=YA(1)
YSCAL2(2,3)=YA(2)
YSCAL2(3,3)=YA(3)
YSCAL2(4,3)=YA(4)
YSCAL2(5,3)=YA(5)
YSCAL2(6,3)=YA(6)
GO TO 35
C Begin power series expansion of fundamental solution
C for degree N=0.
26 CONTINUE
C Set first two solutions to zero.
DO 36 I=1,6
DO 37 J=1,2
Y(I,J)=0.D0
YSCAL(I,J)=0.D0
YSCAL2(I,J)=0.D0
37 CONTINUE
36 CONTINUE
C Set values for first three terms of fundamental solution.
Y(1,3)=1.D0
Y(2,3)=3.D0*LAMBDA(1)+2.D0*MU(1)
Y(3,3)=0.D0
Y(4,3)=0.D0
Y(5,3)=2.D0*PI*G*RHO(1)
Y(6,3)=0.D0
YSCAL(1,3)=-RHO(1)*(4.D0*GAMMA+ANGS+2.D0*WES)/
1 (10.D0*(LAMBDA(1)+2.D0*MU(1)))
YSCAL(2,3)=(5.D0*LAMBDA(1)+6.D0*MU(1))*YSCAL(1,3)
YSCAL(3,3)=0.D0
YSCAL(4,3)=0.D0
YSCAL(5,3)=PI*G*RHO(1)*YSCAL(1,3)
YSCAL(6,3)=0.D0
YSCAL2(1,3)=-RHO(1)*(4.D0*GAMMA+ANGS+2.D0*WES)*YSCAL(1,3)/
1 (4.D0*(7.D0*LAMBDA(1)+14.D0*MU(1)))
YSCAL2(2,3)=(7.D0*LAMBDA(1)+10.D0*MU(1))*YSCAL2(1,3)
YSCAL2(3,3)=0.D0
YSCAL2(4,3)=0.D0
YSCAL2(5,3)=2.D0*PI*G*RHO(1)*YSCAL2(1,3)/3.D0
YSCAL2(6,3)=0.D0
35 CONTINUE
C Calculate terms in power series expansions at minimum radius RMIN.
XV=RMIN
XVS=XV*XV
XV3=VX*XVS
XV4=XVS*XVS