Geology Reference
In-Depth Information
IF(J.EQ.3) ANM2=AN
C Find constant equal to the radius of the ICB to the power alpha.
CONST=RIOB**ANM2
Y(2,J)=CONST*Y(2,J)
Y(4,J)=CONST*Y(4,J)
CONST=CONST*RIOB
Y(1,J)=CONST*Y(1,J)
Y(3,J)=CONST*Y(3,J)
Y(6,J)=CONST*Y(6,J)
CONST=CONST*RIOB
Y(5,J)=CONST*Y(5,J)
56 CONTINUE
C Test if N=0.
IF(N.EQ.0) GO TO 57
C Begin construction of coefficient matrix.
DO 58 I=1,3
CMC1(1,I)=Y(4,I)
CMC1(2,I)=GZEROC*Y(1,I)-Y(2,I)/RHOAC-Y(5,I)
CMC1(3,I)=4.D0*PI*G*RHOAC*Y(1,I)+Y(6,I)
58 CONTINUE
C Set coefficient matrices equal.
DO 59 I=1,3
DO 60 J=1,3
CMC2(I,J)=CMC1(I,J)
CMC3(I,J)=CMC2(I,J)
60 CONTINUE
59 CONTINUE
C Set constant vectors.
BVC1(1)=0.D0
BVC1(2)=1.D0
BVC1(3)=0.D0
BVC2(1)=0.D0
BVC2(2)=0.D0
BVC2(3)=-1.D0
C Store first constant vector.
DO 61 I=1,3
BVC3(I)=BVC1(I)
61 CONTINUE
C Solve for linear combination coefficients.
CALL LINSOL(CMC1,BVC1,3,CAUGC,DET,3,7)
CALL LINSOL(CMC2,BVC2,3,CAUGC,DET,3,7)
C Alter coefficient matrix to apply subseismic condition.
DO 62 I=1,3
CMC3(3,I)=GZEROC*Y(1,I)-Y(2,I)/RHOAC
62 CONTINUE
C Solve for subseismic coefficients.
CALL LINSOL(CMC3,BVC3,3,CAUGC,DET,3,7)
C Calculate unscaled Love numbers.
SUM1=0.D0
SUM2=0.D0
SUM3=0.D0
DO 63 I=1,3
SUM1=SUM1+BVC1(I)*Y(1,I)
SUM2=SUM2+BVC2(I)*Y(1,I)
SUM3=SUM3+BVC3(I)*Y(1,I)
63 CONTINUE
C Calculate radial derivative of V1n at inner core boundary.
SUM=0.D0
DO 64 I=1,3
SUM=SUM+(4.D0*PI*G*RHOAC*Y(1,I)+Y(6,I))*BVC3(I)
64
CONTINUE
Search WWH ::




Custom Search