% ****************************************************** % PE % ****************************************************** % *** set up Munk SSP *** zs = 1000.0; % source depth f = 50; omega = 2 * pi * f % frequency in Hertz eps = 0.00737; c0 = 1500; d = 5000; % bottom depth nz = 250; % number of finite-difference points h = d / nz; h2 = h * h; % mesh spacing z = linspace( 0, d, nz ); % grid coordinates x = 2 * ( z - 1300 ) / 1300; c = c0 * ( 1 + eps * ( x - 1 + exp( -x ) ) ); plot( z, c ); view( 90, 90 ); xlabel( '
简单PE模型(Munk声速剖面)
最新推荐文章于 2024-10-28 13:35:00 发布