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

C Open Earth model file.
OPEN(UNIT=1,FILE=EMODEL,STATUS='OLD')
C Open fundamental solutions files.
OPEN(UNIT=2,FILE='fsolns.dat',STATUS='UNKNOWN')
OPEN(UNIT=3,FILE='fs1.dat',STATUS='UNKNOWN')
OPEN(UNIT=4,FILE='fs2.dat',STATUS='UNKNOWN')
OPEN(UNIT=7,FILE='fs3.dat',STATUS='UNKNOWN')
C Set angular frequency of Earth's rotation (WGS 84).
WE=7.292115D-5
C Calculate square of angular frequency of Earth's rotation.
WES=WE*WE
C Set value of pi.
PI=3.141592653589793D0
C Set value of universal constant of gravitation (CODATA 2006).
G=6.67428D-11
C Enter period, azimuthal number, inertial and Coriolis switches.
WRITE(6,12)
12 FORMAT(1X,'Enter period in hours, azimuthal number and',1X,
1 'inertial and Coriolis switches'/'(1 for in, 0 for out).')
READ(5,*)TPER,MAZ,INERT,ICOR
C Convert period to seconds.
TPER=3600.D0*TPER
C Calculate dimensionless angular frequency.
SIG=PI/(TPER*WE)
C Calculate square of dimensionless angular frequency.
SIGS=SIG*SIG
C Compute square of angular frequency.
ANGS=DFLOAT(INERT)
ANGS=2.D0*PI*ANGS/TPER
ANGS=ANGS*ANGS
C Compute Coriolis term.
COR=2.D0*DFLOAT(MAZ)*DFLOAT(ICOR)
COR=2.D0*PI*COR*WE/TPER
C Set maximum relative error tolerance for integrations.
EPS=1.D-5
C Set minimum starting radius for variable Earth properties.
RMIN=1200.D0
C Set maximum dimensions for interpolation.
M1=100
M2=198
M3=98
C Read in and write out inner core values for Earth model.
C Read in and write out headers.
C Read in Earth model name.
READ(1,13)(ENAME(I),I=1,10)
13 FORMAT(10A7)
WRITE(2,14)(ENAME(I),I=1,10)
WRITE(6,14)(ENAME(I),I=1,10)
14 FORMAT(5X,10A7)
C Read additional headers.
READ(1,15)NN
15 FORMAT(I10)
C Read in number of model points in each region NM(I)
C and number of integration steps NI(I).
READ(1,16)(NM(I),NI(I),I=1,4)
16 FORMAT(8I10)
C Write out Earth properties.
WRITE(6,17)

17 FORMAT(/,/,5X,'Radius',6X,'Rho',4X,'Lambda',7X,'Mu',6X,'Gzero')

WRITE(6,18)