KRAKEN计算群速度（弹性海底）

KRAKEN输出的.prt文件中包含的群速度是根据液态海底的一个公式计算的，因此考虑弹性海底时群速度结果有误，此时可根据定义计算群速度。特征值可从.mod文件中读取，，简单的数值差分就可以得到正确的群速度了。

下图为一个弹性海底的仿真示例，左图结果为.prt文件中的群速度，右图为真实群速度。可见无剪切波时左右两图的群速度曲线一致，有剪切波时左右两图的群速度曲线存在显著差异。

左图结果为.prt文件中的群速度，右图为真实群速度

KRAKEN作者Michael B. Porter的解释:

Are you taking the group speed out of the print file? If so, then the answers will be wrong because the group speed is calculated using a formula from our book that is only valid for purely acoustic problems. The alternative way you can do this is to read the wavenumbers out of the mode file and calculate the group speed as a numerical derivative dk/ d omega or something like that.