The result of running the test script DOLAND is essentially this:
Freq (hz) : 9.99999978E-03
c (m/s) : 1.00000000
seafloor : F
Model :HOMO3.5.mod
Model name: HOMO3.5.mod
LAYER H P-VEL S-VEL DENSITY
1 0.1000E-02 6.000 3.500 2.500 These are from the given velocity model
2 0.1000E-02 6.000 3.500 2.500
3 0.1000E-02 6.000 3.500 2.500
4 0.1000E-02 6.000 3.500 2.500
5 0.1000E-02 6.000 3.500 2.500
6 0.1000E-02 6.000 3.500 2.500
7 0.1000E-02 6.000 3.500 2.500
8 0.1000E-02 6.000 3.500 2.500
9 0.1000E-02 6.000 3.500 2.500
...............................................
34 0.1000E-01 6.000 3.500 2.500
35 0.1000E-01 6.000 3.500 2.500
36 0.1000E-01 6.000 3.500 2.500
COMPZ : 3.9386949257647075E-010 This is the ratio of Uz (positive down) to Pressure at surface
TILT : 6.1432600164909040E-008
ZZ I Z(mid)(m) Ka Kb Kr Kkappa Kmu Krho
XX I Z(mid)(m) Ka Kb Kr Kkappa Kmu Krho
ZZ 1 0.50000 -0.24366 0.08215 -0.08074 -0.06656 -0.01430 0.00000 Recall that this is experimental code. The ZZ
XX 1 0.50000 -0.24366 0.08215 -0.08074 -0.06655 -0.01420 0.00001 rows compute all kernels by numerical differentiation.
ZZ 2 1.50000 -0.21486 0.06819 -0.07333 -0.05868 -0.01474 0.00000
XX 2 1.50000 -0.21486 0.06819 -0.07333 -0.05869 -0.01465 0.00001 The rows marked XX compute the Ka, Kb and Kr by
ZZ 3 2.50000 -0.18945 0.05250 -0.06835 -0.05180 -0.01679 0.00000 numerical differentiation but the Kkappa, Kmu and
XX 3 2.50000 -0.18945 0.05250 -0.06835 -0.05175 -0.01673 0.00012 Krho kernels are derived from the first two using
........................................................................................... the relation from Tanimoto and Wang (2019).
ZZ 35 115.00000 0.00000 -0.00003 -0.00001 -0.00000 -0.00001 0.00000
XX 35 115.00000 0.00000 -0.00003 -0.00001 0.00000 -0.00001 -0.00000 To access only the XX results, do a
ZZ 36 125.00000 0.00000 -0.00001 -0.00001 -0.00000 -0.00001 0.00000 grep XX Land.txt
XX 36 125.00000 0.00000 -0.00001 -0.00001 0.00000 -0.00001 -0.00000
The output of the ocean-bottom computation from the script DOSEA is similar. The only difference is in the header information. The phase velocity is computed as a funciton of the water depth abd frequency. Since the phase velocities are much greater than those use for the solid-atmosphere problem, the kernels are significant to greater depths. In addition, the model applies to the material in the seafloor and does not have an initial fluid layer. The
Freq (hz) : 9.99999978E-03 c (m/s) : 92.3199005 seafloor : T
g (m/s^2) : 9.80000019
depth (m) : 1000.00000 Model name: HOMOSEA3.5.mod LAYER H P-VEL S-VEL DENSITY 1 0.1000 6.000 3.500 2.500 2 0.1000 6.000 3.500 2.500 3 0.1000 6.000 3.500 2.500 4 0.1000 6.000 3.500 2.500 5 0.1000 6.000 3.500 2.500 6 0.1000 6.000 3.500 2.500 ............................................... 69 0.1000 6.000 3.500 2.500 70 0.1000 6.000 3.500 2.500 71 0.1000 6.000 3.500 2.500 72 0.1000 6.000 3.500 2.500 COMPZ : 3.6381283763375905E-008 COMPR : 1.2383965852791334E-008 ZZ I Z(mid)(m) Ka Kb Kr Kkappa Kmu Krho XX I Z(mid)(m) Ka Kb Kr Kkappa Kmu Krho ZZ 1 50.00000 -0.00263 0.00088 -0.00087 -0.00072 -0.00015 0.00000 XX 1 50.00000 -0.00263 0.00088 -0.00087 -0.00072 -0.00015 0.00000 ZZ 2 150.00000 -0.00229 0.00072 -0.00079 -0.00063 -0.00016 0.00000 XX 2 150.00000 -0.00229 0.00072 -0.00079 -0.00063 -0.00016 0.00000 ZZ 3 250.00000 -0.00200 0.00053 -0.00073 -0.00055 -0.00019 0.00000 ........................................................................................... ZZ 71 7050.00000 -0.00000 -0.00002 -0.00001 -0.00000 -0.00001 0.00000 XX 71 7050.00000 -0.00000 -0.00002 -0.00001 -0.00000 -0.00001 -0.00000 ZZ 72 7150.00000 -0.00001 -0.00010 -0.00005 -0.00001 -0.00005 -0.00000 XX 72 7150.00000 -0.00001 -0.00010 -0.00005 -0.00000 -0.00005 0.00000