c
aniprop - program to calculate propagating modes of anisotropic layer
c writes a file of dispersion curves at evenly-space freq points
c
c reads a layered model like the following file (ignore the "c "s)
c
c K&H SoCal model, deep crust horizontal anisotropy
TITLE
c 3
# OF LAYERS OVER HSPACE
c 45 45
THETA,PHI ORIENTATION ANGLES FOR SYMMETRY AXIS
c 4000 5500 0.06 0.00 3175 0.03 2600
c 0 0
c 27400 6300 0.00 0.00 3637 0.00 2800
DEPTH (M), VP (M/S), "B", "C",VS (M/S), "E"
c 90 45
c 32400 6800 0.04 0.00 3925 0.02 2900
B,C,E ARE ANISOTROPIC PARAMETERS
c 0 0
c 60000 7800 0.00 0.00 4500 0.00 3200
NOTE: HSPACE MUST BE ISOTROPIC
c
c revised to calculate dispersion curves: 7/7/95
c revised to make rootfinder more focussed: 7/11/95
c revised to make computations dimensionless: started 8/3/95,
finished 8/14
c bugfix for tilted axis of symmetry: 8/16/95
c bugfix for group velocity: 8/22/95
c bugfix for defective matrices: 8/22/95
c bugfix for transition to evanescence in top layer: 9/10/95
c multiple bugfixes in the rootfinder and in the test for
c horizontal waves in the surface layer
2/5/00
c bugfix for underflow in routine grvel 3/8/00
c
c
c compile sequence
c eislib is the Eispack library
c f77 -o aniprop -fast -native -O5 aniprop.f /data/d4/park/Ritz/eislib
c
c compile sequence with Park-specific plotting programs
for debugging
c plotlib&jlib link with a complaint-that-I-ignore in Solaris
4.x
c xf77 -o aniprop -O aniprop.f /data/d4/park/Plotxy/plotlib /data/d4/park/Ritz/eislib
/data/d4/park/Ritz/jlib
c
c for hexagonally symmetric media
c reads fast axis orientation, constants A,B,C,D,E from file
animodel
c calculate quadratic eigenvalue problem based on the Christoffel
matrix
c see appendix of P. Shearer's thesis