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

C Find remaining distance to the ICB.
H=R(N1)-XV
C Compute fraction of radius.
RATIO=H/R(N1)
C Continue integration only if remaining distance is more than
C 1.D-9 of radius.
IF(RATIO.GT.1.D-9) GO TO 51
C Convert solutions back to y-variables.
C Set value of alpha for first two solutions.
ANM2=AN-2.D0
C Convert three solutions.
DO 61 J=1,3
C Set value of alpha for third solution.
IF(J.EQ.3) ANM2=AN
C Find constant equal to the radius of the ICB to the power of 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)

61 CONTINUE

END

As an illustration of output from the programme ICFS, we show the free solu-
tions, regular at the geocentre, for degree n 1, at a period

of 4 hours. They are generated by the free constants A 1,0 , A 6,0 and A 4,1 , and are
shown in Figures 3.1, 3.2 and 3.3. In each case, the plots are shown for unit value
of the respective free constant.

= 1, azimuthal number m =

4.0

3.5

y 2 × 2 × 10 − 5

−

3.0

y 3 × 2

y 5 × 2 × 10 − 1

2.5

2.0

y 1

1.5

1.0

− y 6 × 2 × 10 − 6

y 4 × 5 × 10 − 5

0.5

0.0 0

174 347 521 695

Radius (km)

869 1042 1216

Figure 3.1 Free solution generated by the free constant A 1,0 = 1.