Bellhop射线追踪程序概述

G
笛卡尔坐标系中的几何帽形波束（默认）
OPTION(2:2)
g
射线中心坐标系中的几何帽形波束
C
采用笛卡尔坐标系Cerveny高斯波束
R
采用中心射线坐标系Cerveny高斯波束
B
几何高斯波束
设定波束移位效果等选项；可以空着
*
读取声源波束模式*.sbp（角度-幅度）文件
OPTION(3:3)
O
不读取（默认）
• SPARC 在时域直接运算 的谱积分程序
BELLHOP .ENV 文件
BELLHOP .ENV 文件
6.1 水平海底波导
flatwav_A.env / Arrives
行号
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
• 3、对于每个作业的主程序文件一句一句地消化， 深化理解每份作业的内涵，算是提高。耗时估计3-7 天。
• 4、本文档适合大三以上和非水声物理专业同志们 的快速参学。
声学工具箱的结构
• KRAKEN normal mode approach
• KRAKEN 简正波方法
• BELLHOP Gaussian beam tracing and ray tracing code
Bellhop射线追踪程序概述
招式说明：
• 1、在64位windows平台上，将本程序包的文件夹 添加在matlab“设置路径”的“添加并包含子文件 夹”。
• 2、读完文档，跑完本程序包所有招式（作业）， 射线追踪算法的运行就算初步通关，对于海洋传播 声场的二维射线表达方法基本上建立感性认识。单 单溜达作业（六路招式）耗时半天。
''
不包括波束移位效果（默认）；
S
包括波束移位效果
设置阵列类型；可以空着
OPTION(4:4)
R
点源（圆柱形坐标系）（默认）
X
线源（笛卡尔坐标系）
设置阵列类型；可以空着
OPTION(5:5)
R
直线网格，接收器置于rr(:)×rd(:)点（默认）
I
不规则网格，接收器在rr(:)，rd(:)。
Matlab 脚本
内容
'Munk / Flat / Arrives' 50.0 1 'SVW' 51 0.0 5000.0 0.0 1548.52 / 200.0 1530.29 / 250.0 1526.69 / 400.0 1517.78 / 600.0 1509.49 / 800.0 1504.30 / 1000.0 1501.38 / 1200.0 1500.14 / 1400.0 1500.12 / 1600.0 1501.02 / 1800.0 1502.57 / 2000.0 1504.62 / 2200.0 1507.02 / 2400.0 1509.69 / 2600.0 1512.55 / 2800.0 1515.56 / 3000.0 1518.67 / 3200.0 1521.85 / 3400.0 1525.10 / 3600.0 1528.38 / 3800.0 1531.70 / 4000.0 1535.04 / 4200.0 1538.39 / 4400.0 1541.76 / 4600.0 1545.14 / 4800.0 1548.52 / 5000.0 1551.91 / 'A' 0.0 5000.0 1600.00 0.0 1.8 .0 .0 / 1 1000.0 / 1 1000.0 / 1 101.0 / 'A' 101 -14.0 14.0 / 100.0 5500.0 102.0
• BELLHOP 高斯波束追踪 与射线追踪程序
• SCOOTER spectral integral (also know as wavenumber integration or the reflectivity method)
• SCOOTER 谱积分(也称为 波数积分或反射率法)
• SPARC A spectral integral code that operates directly in the time domain.
OPTIONS3（选项3） 语 法 5个字符的字符串，用于设定输出选项。
说明
设定将写入输出文件的信息类型。
R
生成声线文件
E
生成本征声线文件
OPTION(1:1)
A
生成幅度-延时文件（ascii码）
Leabharlann Baidu
a
生成幅度-延时文件（二进制）
C
相干TL计算
I
非相干TL计算
S
半相干TL计算（Lloyd镜声源模式）
设置声压计算采用的近似方法；可以空着
R
E
C\I\S
A
6.2 高斯海山(变化的海底)
• 应用高斯函数生成一个理想化的海山。
行号
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
说明
TITLE FREQ (Hz) NMEDIA SSP-TOP-WATER-OPT NMESH SIGMA Z(NSSP) Z() CP() CS() RHO() AP() AS() .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . .. . . . . BOTOPT SIGMA ZB CPB CSB RHOB APB ASB NSD SD(1:NSD) (m) NRD RD(1:NRD) (m) NRR RR(1:NRR ) (km) OPTION: 'R/E/C/A/I/S' NBEAMS ISINGLE ALPHA(1:NBEAMS) (°) STEP (m) ZBOX (m) RBOX (km)