fast marching method 计算内波相速度

在计算内孤立波传播轨迹的时候，可以应用(Jackson,2009)(J09)提出的经验公式

C(x,y)=Cmaxtanh[B1+H(x,y)B2]−−−−−−−−−−−−−−−−√

这里就使得相速度只与水深和经验参数Cmax,B_1,B_2,有关。
根据J09，经验参数为：
Cmax=2.971 m/s
B1=0.003
B2=1390.758 m
起始点源为：
Lat=20.529∘N
Lon=124.879∘E
得出相速度之后，可以解eikonal equation（EE）来确定波峰线的传播轨迹。

(∂T∂x)2+(∂T∂y)2=1C(x,y)2(EE)

解EE的方法为fast marching method(FMM).FMM是一种解最短路径的方法，与传统的解最短路径问题Dijkstra’s Method(DM)算法的不同在于，DM的路径一定要沿各点前进，而FMM的路径可以为不沿格点的斜线。
[Level Set Methods and Fast Marching Methods p94-95


计算结果如下