编写了一个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(lon,lat);
lat=lat';
lon=lon';
% 2. calculate f
f=coriolisf(lat);
% 3. cal geo_u_v
g=9.8;
R=6400000; % Earth radius
U=zeros(52,48);
V=zeros(52,48);
for i=2:51
for j=2:47
dx=(l