The subseismic equation and boundary conditions - Earth Dynamics: Deformations and Oscillations of the Rotating Earth - page 442

CM01(3,I)=0.D0
CM01(4,I)=0.D0
87 CONTINUE
CM01(3,1)=GZEROM
CM01(3,2)=-1.D0/RHOBM
CM01(3,3)=-1.D0
CM01(4,1)=4.D0*PI*G*RHOBM
CM01(4,2)=0.D0
CM01(4,4)=1.D0
C Set coefficient matrices equal.
DO 88 I=1,4
DO 89 J=1,4
CM02(I,J)=CM01(I,J)
CM03(I,J)=CM02(I,J)
89 CONTINUE
88 CONTINUE
C Set values of constant vectors.
BV01(1)=0.D0
BV01(2)=0.D0
BV01(3)=1.D0
BV01(4)=0.D0
DO 90 I=1,3
BV02(I)=0.D0
90 CONTINUE
BV02(4)=-1.D0
C Store first constant vector.
DO 91 I=1,4
BV03(I)=BV01(I)
91 CONTINUE
C Solve for linear combination coefficients.
CALL LINSOL(CM01,BV01,4,CAUGS0,DET,4,9)
BV1(1)=BV01(1)
CALL LINSOL(CM02,BV02,4,CAUGS0,DET,4,9)
BV2(1)=BV02(1)
C Alter coefficient matrix to apply subseismic condition.
CM03(4,1)=GZEROM
CM03(4,2)=-1.D0/RHOBM
CM03(4,4)=0.D0
C Solve for subseismic coefficients.
CALL LINSOL(CM03,BV03,4,CAUGS0,DET,4,9)
BV3(1)=BV03(1)
BV3(4)=BV03(4)
C Calculate potential condition parameter.
BN=(AN+1.D0)/(4.D0*PI*G*RHOBM*BV3(1)+BV3(4))
C Scale Love numbers to dimensionless values.
84 CONTINUE
BV1(1)=GZEROM*BV1(1)
BV2(1)=GZEROM*BV2(1)/RCMB
BV3(1)=GZEROM*BV3(1)
C Write out internal Love numbers.
WRITE(6,69)N,BV1(1),BV2(1)
WRITE(2,70)BV3(1),BN
C End loop over degree.
71 CONTINUE
C Write out dimensionless angular frequency coefficients of degree-one
C reciprocal Love number.
WRITE(2,92)ALPHA,BETA,GAMA

92 FORMAT(3D23.15)

END