toe1: Theory of Everything, Rev 1 computes radial integrals of mode-mode interactions for path-integral coupling approximation. Differs from toe0 in calculating radial integrals of coupling kernels against Legendre polynomials of order 0 thru 4. This allows the anisotropic perturbation to have radial dependence. Input is a REM and its eigenfunctions, plus the orientation angles (psi,xi) of an axis of symmetry and parameters A,B,C,D,E along the path. Output of this routine is a big file of interaction scalars for modes that satisfy |l_1-l_2| = s.
toe1 reads an input file named "ph_toe1." An example:
premi0BE45_75_300
0.25 16 0.9
0
75 300 5 20
90 45
0.0 0.02828 0.0 0.0 0.01414
/data/d8/park/Modes/premiss
/data/d8/park/Modes/premissnl
/data/d8/park/Modes/premisseig
/data/d8/park/Modes/premistnl
/data/d8/park/Modes/premisteig
****EXPLANATION**********
premi0BE45_75_300 ! output file for integration scalars
0.25 16 0.9 ! f1 f2 fco --> freq range to compute, coupling half-bandwidth. This choice of parameters will compute interation scalars between a target mode and all modes that lie within 0.9 mHz of its degenerate multiplet frequency.
0 ! the maximum overtone number to include in the calculation. "0" includes fundamentals only.
75 300 5 20 ! depth range for perturbation, in km (75-300 km), number of Legendre polynomials to use (5, up to fourth-order polynomial), and maximum allowed |l1-l2| between coupling pairs. This last parameter enables the user to skip mode pairs whose angular wavenumbers differ too much to couple significantly. In practice, if s=20 is the maximum wavenumber the user intends to impose on a source-receiver great circle, there is no reason to compute interactions for |l1-l2|>20.
90 45 ! psi and xi, the polar angles of w-hat, the axis of symmetry. psi is measured from local vertical, and xi CCW from the radial component along the great-circle path.
0.0 0.02828 0.0 0.0 0.01414 ! A,B,C,D,E --> IMPORTANT: first radial function is the first Legendre polynomial, which is a constant 1/sqrt(2). If the user wishes to prescribe 1% pk2pk shear anisotropy, E=0.01414 here. For the case shown, we have prescibed 2% pk2pk anisotropic Vp & 1% pk2pk Vs, both with cos 2\xi azimuthal variation.
/data/d8/park/Modes/premiss ! reference earth model (isotropic PREM), card-deck format
/data/d8/park/Modes/premissnl ! eigenfrequency file for spheroidal modes for above model
/data/d8/park/Modes/premisseig ! eigenfunction file for spheroidal modes for above model
/data/d8/park/Modes/premistnl ! eigenfrequency file for toroidal modes for above model
/data/d8/park/Modes/premisteig ! eigenfunction file for toroidal modes for above model