Rotating fluids and the outer core - Earth Dynamics: Deformations and Oscillations of the Rotating Earth - page 417

10 FORMAT(1X,'Type in Earth model file name.')
READ(5,11)EMODEL
11 FORMAT(A20)
C Open Earth model file.
OPEN(UNIT=1,FILE=EMODEL,STATUS='OLD')
OPEN(UNIT=2,FILE='decomp.dat',STATUS='UNKNOWN')
C Set maximum dimensions for interpolation.
M1=100
M2=198
M3=98
C Read in and write out Earth model.
C Read in and write out headers.
READ(1,12)(ENAME(I),I=1,10)
12 FORMAT(10A8)
WRITE(6,13)(ENAME(I),I=1,10)
13 FORMAT(11X,10A8)
WRITE(2,12)(ENAME(I),I=1,10)
READ(1,14)NN
14 FORMAT(I10)
READ(1,15)(NM(I),NI(I),I=1,4)
15 FORMAT(8I10)
WRITE(6,16)
16 FORMAT(/16X,'Number of model points',1X,
1 'and number of integration steps.'/)
WRITE(6,15)(NM(I),NI(I),I=1,4)
WRITE(6,17)
17 FORMAT(//5X,'Radius',6X,'Rho',4X,'Lambda',7X,'Mu',6X,'Gzero')
WRITE(6,18)
18 FORMAT(6X,'(km)',5X,'(gm/cc)',2X,'(kbars)',4X,'(kbars)',2X,
1'(cm/sec)'/)
C Read in Earth model, inner core, outer core, mantle and crust.
K=0
DO 19 M=1,4
N1=NM(M)
READ(1,20)(R(I),RHO(I),LAMBDA(I),MU(I),GZERO(I),I=1,N1)
20 FORMAT(1X,F10.1,F10.2,F10.1,F10.1,F10.1)
DO 21 I=1,N1
C Scale Earth model to SI values and store.
J=K+I
RI(J)=R(I)*1.D3
RHOI(J)=RHO(I)*1.D3
LAMBDAI(J)=LAMBDA(I)*1.D8
MUI(J)=MU(I)*1.D8
GZEROI(J)=GZERO(I)*1.D-2
C Write Earth model out to screen.
WRITE(6,20)R(I),RHO(I),LAMBDA(I),MU(I),GZERO(I)
21 CONTINUE
NK(M)=K
K=K+N1
19 CONTINUE
C Begin calculation of decompression factor f.
C Set up interpolation for the outer core.
N1=NM(2)
N2=2*N1-2
N3=N1-2
C Put outer core values in active locations.
DO 22 I=1,N1
J=NK(2)+I
R(I)=RI(J)
RHO(I)=RHOI(J)
MU(I)=MUI(J)