integrating to get the total volumetric transport in Sv

2 次查看（过去 30 天）
Hannah Evans
Hannah Evans 2021-3-23
I'm a little lost on how to integrate the Sverdrup transport to get the total volumetric transport. the answer should be
The zonal wind stress contribution to the total Sverdrup transport: -22.5 Sv (negative sign means southward)
meridional wind stress contribution to the total Sverdrup transport: -3.4 Sv
Both zonal and meridional wind stress contribution: -25.9 Sv.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% sverdrup.m
% Reads and plots wind stress data in the subtropical North Atlantic Ocean.
% The wind stress data are defined at 2-degree intervals in longitude and
% latitude. Originally coded by Raymond G. Najjar 2004
% Modified by PRB in 2005, 2007
clear
% Read in monthly wind stress (dyne cm^2) fields at 25N and 27N in the Atlantic
load taux25n.txt
load tauy25n.txt
load taux27n.txt
load tauy27n.txt
% Compute annual averages and convert to N m^-2. (dyne is a cgs unit for force = g*cm s^-2)
tx25 = mean(taux25n)*0.1;
ty25 = mean(tauy25n)*0.1;
tx27 = mean(taux27n)*0.1;
ty27 = mean(tauy27n)*0.1;
% The wind stress data extend from 81W to 15W; that's from the east coast of Florida
% to the northwest coast of Africa.
% Plot the stresses as a function of distance in km from Florida.
imax = size(tx25,2);
i = 1:imax;
xkm = (i-1)*222*cos(26.*pi/180);
yo = i*0;
% Part I: The integrated Ekman and Sverdrup transports
% Computing Ekman Transport
omega = (7.29).*10.^-5;
f = 2.*omega.*sin((26./180).*3.14);
tx26 = (tx25+tx27);
ty26 = (ty25+ty27);
Mx26 = (ty26./f);
My26 = (-tx26./f);
% Computing wind stress curl and Sverdrup Transport
dx = 111.*cos((26./180).*3.14);
dy = 111.*1000;
a = 6.4.*10.^6;
B = (1./a).*2.*omega.*cos(26);
for i =1:imax-1
curl(i+1) = (1./B).*((ty26(i)-ty26(i))/(2.*dx) - (tx27(i)-tx25(i))/(2.*dy));
end
figure(1)
subplot (2,1,1)
plot(xkm,tx25,'-b',xkm,ty25,'-r',xkm,tx27,'-.b',xkm,ty27,'-.r',xkm,yo,'k')
legend('\tau_x at 25\circN','\tau_y at 25\circN','\tau_x at 27\circN','\tau_y at 27\circN')
title('Annual Mean Wind Stress at 25\circN and 27\circN in the Atlantic')
xlabel('distance east of Florida (km)')
ylabel('stress (N m^-^2)')
subplot (2,1,2)
plot(xkm,My26,'-.k',xkm,curl,'k')
title('Ekman and Sverdrup mass transports at 26°N');
xlabel('distance east of Florida (km)');
ylabel('transport (kg/ms)');

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Stress and Strain 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by