编写了一个matlab脚本计算地转流的流速,可供参考。其中地转流流速的公式为
SLA为海表面高度数据,g是重力加速度,f是科里奥利力频率。
画出来的结果如下图
所用的脚本代码如下
% using sea surface height data calculate geostrophic current
% author:Davidli date:2022.1.5
clear
% 1. load sla ,lat, lon
ncname='D:\sla\TOPEX AVISO\sla_2008.nc';
sla8=ncread(ncname,'sla');
sla=sla8(:,:,9);
lat=ncread(ncname,'latitude');
lon=ncread(ncname,'longitude');
[lon,lat]=meshgrid