Motion-Stress from disper.d

The file, disper.d can be edited and disper run in a different way, computing the motion-stress vector for a given frequency and mode. For example, from the file disper.tmp, note the phase velocity for a particular mode. Pick a frequency of interest, say 32.2265605 Hz. We scroll down disper.tmp:

       31.2499981 |    259.6905012   291.4923132
       31.7382793 |    258.6609733   286.9434926
       32.2265605 |    257.4416940   282.8358906
       32.7148417 |    255.9845241   279.2012189
       33.2031230 |    254.2462253   276.0619743

We edit the pvlcty= line of disper.d to use the phase velocity. For the fundamental mode, the replacement is:

 <    pvlcty=0.0, pfreq=0.0, zend=100.0,
---
>    pvlcty=257.4416940, pfreq=32.2265605, zend=10.0,

We also change zend, choosing a more relevant depth of interest which depends on the frequency. Low frequencies extend deeper than high frequencies. We also change the computational depth interval, deltz. The replacement is:

 <   deltz=    1.0000,
---
>   deltz=    .0010,
We change the computational increment from 1 meter to .001 meters. Review the relevant lines of the disper.d file shown above in section 9.2.5 for example. Running disper with the edited disper.d file will replace disper.tmp with a new version listing the motion-stress vector computed at the new deltz interval. Figure 51 is the result.

Figure 51: DISPER: Re-running disper to compute the motion stress vectors. See section 9.2.6.1 for how to edit disper.d. The file, mat2.m created this plot. Blue is horizontal, Red is vertical.
\includegraphics[scale=.7]{FigureMS.pdf}