The subseismic equation and boundary conditions - Earth Dynamics: Deformations and Oscillations of the Rotating Earth - page 431

2 Y(6,6),YSCAL(6,3),YSCAL2(6,3),YA(6),A(6,6),AS(6,6),CAUG(6,13),
3 YIC(6,3),Y1(6,6),Y2(6,6),DELTA(6,3),CM1(5,5),CM2(5,5),
4 CM3(5,5),BV1(5),BV2(5),BV3(5),CAUGC(3,7),CMC1(3,3),CMC2(3,3),
5 CMC3(3,3),BVC1(3),BVC2(3),BVC3(3),CAUGS(5,11),CAUGS0(4,9),
6 CM01(4,4),CM02(4,4),CM03(4,4),BV01(4),BV02(4),BV03(4)
DOUBLE PRECISION MU(100),LAMBDA(100),MUI(300),LAMBDAI(300),K1(6,6)
CHARACTER*7 ENAME
CHARACTER*20 EMODEL
C Enter Earth model file name.
WRITE(6,10)
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')
C Open Love number file.
OPEN(UNIT=2,FILE='love.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 in hours, azimuthal number, inertial and Coriolis switches.
WRITE(6,12)
12 FORMAT(1X,'Enter period in hours, azimuthal number and inertial',
1 ' 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.D1
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.
C Earth model name.
READ(1,13)(ENAME(I),I=1,10)
13 FORMAT(10A7)
WRITE(2,13)(ENAME(I),I=1,10)
WRITE(6,14)(ENAME(I),I=1,10)
14 FORMAT(11X,10A7)
C Number of Love numbers to be computed.
READ(1,15)NN