MATLAB处理WOA18浮标数据

1、利用woa18数据，查找浮标坐标点、并对坐标点处绘制声速剖面。
2、利用Krakenc计算传播损失差。
环境：
Matlab 2015b以上版本
ACTup水声学综合工具箱
数据来源：
World Ocean Atlas 2018（WOA18）
整理：
读取nc数据，分别读取经纬度、深度、温度、盐度。
根据浮标经纬度，获取坐标点所在的行和列，并以此读取浮标所在位置的温度、盐度。

MATLAB代码实现如下：

clear;
clc;
%% 提取nc数据
tem_ncid=('woa18_decav81B0_t15_04.nc');% 夏季温度数据文件名赋值，冬季同理
sal_ncid=('woa18_decav81B0_s15_04.nc');% 夏季盐度数据文件名赋值，冬季同理
lat=ncread(tem_ncid,'lat');%读取经度
lon=ncread(tem_ncid,'lon');%读取纬度
depth=ncread(tem_ncid,'depth');%读取深度
tem=ncread(tem_ncid,'t_an');%读取温度
sal=ncread(sal_ncid,'s_an');%读取盐度
%% 找浮标经纬度所在行和列
num_lon=find(lon==-1.4512500e+02);
num_lat=find(lat==50.875);
%% 提取Papa所在的温度、盐度数据
a=tem(num_lon,num_lat,:);