matlab 基于温度盐度的地转流。

kdx123

clc; clear; close all;
load('G:\生成数据\HYCOM\temp.mat')
load('G:\生成数据\HYCOM\sal.mat')
%%
t=temp(:,:,:,1);
s=sal(:,:,:,1);
v=s;
v(:)=nan;
u=v;
% %%
for i=1:length(lon)
    for j=1:length(lat)
depth1(i,j,:)=sw_pres(depth,lat(j));
    end
end
%%
ro=sw_dens(s,t,depth1);
rom=nanmean(ro(:));
rox=ro(2:end,:,:)-ro(1:end-1,:,:);
roy=ro(:,2:end,:)-ro(:,1:end-1,:);
for i=1:length(lon)-1
    for j=1:length(lat)-1
[s(i,j),al(i,j)]=sw_dist([lat(j),lat(j+1)],[lon(i),lon(i+1)],'km');
dx(i,j)=s(i,j)*cosd(al(i,j))*1e3;
dy(i,j)=s(i,j)*sind(al(i,j))*1e3;
    end
end
%%
f=sw_f(lat);
for i=1:length(lon)-1
for j=1:length(lat)-1
    v(i,j,:)=1/rom/f(i).*rox(i,j,:)/dx(i,j);
    u(i,j,:)=-1/rom/f(i).*roy(i,j,:)/dy(i,j);
end
end
%%
ro0=ro(:,:,1);
u0=u(:,:,1);
v0=u(:,:,1);
s0=sqrt(u0.^2+v0.^2);
contourf(lon,lat,s0','linecolor','none')
colorbar
colormap jet

%%

萌新

这是什么图