BELLHOP-3d模型垂直剖面仿真

初学BELLHOP-3d模型，然后想看看3d模型是否能够涵盖2d垂直剖仿真的功能，以前篇文章下载声学工具箱中的munk文件为例，对传播损失进行仿真。

可以使用bellhop3d进行2d和3d仿真，运行plotshd均可得到二维垂直剖面的传播损失图（plotshdpol是画三维水平剖面的）

2d仿真：

参考slice2dHatCart

2d模式设置：

在bellhop3d运行的env文件里只需要在计算类型后边的数字标为2即设置为bellhop的2d模型

2d角度设置：

 在此基础上,我们探究不同于2d模型多出的角度设置参数问题：

Ntheta设置：1及1往上的数字均可，但是最终plotshd这个函数在提取声压的时候只会在提取第一个角度下的传播损失值；

Nbeta设置：无关，设置数值无用；

env文件角度设置：

plotshd声压提取方式：

 得到的传播损失如下图所示：

 3d仿真：

参考slice3dHatCart

3d模式设置：

和2d一样，只需要将数字写为3即可 

 

 3d角度设置：

Ntheta设置：与上述2d一致，设置1以及1往上。

Nbeta设置：需要注意，需要设置2以及2往上的数字才可以，否则会报错。

 通过这两个传播损失图的对比，可以发现结果是一样的。

加入海底对比分析

0°和180°垂直剖面如下图所示，加入海底后可以更加清晰的进行对比，呈现2d和3d结果的一致性，仿真时均可拿来使用。

2d结果：

                           （a）0°                                                              （b）180°

 3d结果：

                        （a）0°                                                              （b）180°

 海深5000m时传播损失的水平剖面如下图所示，可以看出0°时60km传播损失有所减小，之前始终处于深蓝色；180°时则是40km开始。与上面垂直剖面结果一致，说明垂直剖面的正确性。

再换个深度1500m，可以看出180°时20km左右有明显影区，与垂直剖面结果一致，再次验证其正确性。