Geology Reference
In-Depth Information
LAMBDA(I)=LAMBDAI(J)
GZERO(I)=GZEROI(J)
24 CONTINUE
C Calculate gravity gradient at geocentre.
GAMMA=(4.D0*PI*G*RHO(1)-2.D0*WES)/3.D0
C Set minimum radius for variable Earth properties.
R(1)=RMIN
C Set gravity at minimum radius.
GZERO(1)=GAMMA*RMIN
C Construct interpolation matrix.
CALL SPMAT(N1,N2,N3,C,R,B,M1,M2,M3)
C Begin loop over degree for inner core.
WRITE(6,25)
25 FORMAT(//11X,'Degree',14X,'Love number h1',14X,'Love number h2')
DO 26 NP1=1,NN
N=NP1-1
AN=DFLOAT(N)
C Test if N=0.
IF(N.EQ.0) GO TO 27
C Test if N=1 and set counter.
IF(N.NE.1) GO TO 28
ICOUNT=1
C Store ANGS and COR
ANGSS=ANGS
CORS=COR
C Set ANGS and COR to zero for calculation of DC value.
ANGS=0.D0
COR=0.D0
C Begin power series expansions of fundamental solutions for N>0.
C Set initial values of fundamental solutions.
C Set initial values of first fundamental solution.
28 CONTINUE
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 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)
Search WWH ::
Custom Search