Earth’s figure and gravitation - Earth Dynamics: Deformations and Oscillations of the Rotating Earth - page 384

WRITE(6,21)

FORMAT(11X,'(km)',9X,'(gm/cc)',11X,'(10 − 11 21 ',4X,

1'(10 30 ',3X,'(10 21 kg)')

WRITE(2,22)
WRITE(6,22)
22 FORMAT(47X,'/m)',3X,'kg.m 2 )')
C Scale and write out values.
DO 23 I=1,NS
RAD=RI(I)/1000.D0
DEN=RHOI(I)/1000.D0
RECIPF=1.D0/FI(I)
FP=FPI(I)*1.D11
CMAS=CMA(I)*1.D-30
MASSC=MASS(I)*1.D-21
WRITE(2,24)RAD,DEN,RECIPF,FP,CMAS,MASSC
WRITE(6,24)RAD,DEN,RECIPF,FP,CMAS,MASSC
24 FORMAT(1X,F15.1,F15.3,F10.2,F8.3,2F10.2)
23 CONTINUE
C Write out gravity torque coefficient GAMMA.
WRITE(2,25)GAMMA
WRITE(6,25)GAMMA
25 FORMAT(//5X,'GAMMA=',1PD8.2,1X,'N.m')
C Write out equivalent mass of the inner core as though density is that
C at bottom of the outer core.
WRITE(2,26)MASSIP
WRITE(6,26)MASSIP
26 FORMAT(//5X,'MASSIP=',1PD8.2,1X,'kg')
C Write out density at the bottom of the outer core.
RHOBOC=RHOBOC/1000.D0
WRITE(2,27)RHOBOC
WRITE(6,27)RHOBOC

27 FORMAT(//5X,'density at bottom of outer core=',2PD10.4,1X,'gm/cc')

END

Tabl e 5 . 4 Profiles of parameters for the inner core as computed
for Earth model Cal8.

f

Radius Density C − A M i

cm 3 ) 10 − 11 m − 1 ) 10 30 kgm 2 ) 10 21 kg)

(km) (g / 1/ f ( × ( × ( ×

0.0 13.580 422.17 0.0 0.0 0.0

110.5 13.587 422.23 − 0.266 0.0 0.08

221.1 13.592 422.27 − 0.210 0.03 0.62

331.6 13.594 422.31 − 0.132 0.22 2.08

442.2 13.592 422.33 − 0.032 0.91 4.92

552.7 13.584 422.32 0.089 2.78 9.61

663.3 13.570 422.29 0.233 6.92 16.60

773.8 13.549 422.23 0.398 14.94 26.35

884.4 13.520 422.13 0.580 29.10 39.29

994.9 13.481 422.00 0.776 52.36 55.87

1105.5 13.426 421.82 1.029 88.50 76.50

1216.0 13.340 421.59 1.361 142.15 101.57