- 积分
- 6749
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-3-7
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
计算出来的风应力旋度太大了,有没有大神能帮忙看看问题出在哪里,帮忙改一下?
clear;clc;
u=ncread('D:\biyesheji\ECMWF\uv107.nc','u10');
v=ncread('D:\biyesheji\ECMWF\uv107.nc','v10');
lat=ncread('D:\biyesheji\ECMWF\uv107.nc','latitude');
lon=ncread('D:\biyesheji\ECMWF\uv107.nc','longitude');
roh=1.26; % kg/m^3, air density
% Computation of Wind Stresses
time=1;
[lt, ln]=size(u(:,:,time));
Tx=NaN(lt, ln);
Ty=NaN(lt, ln);
for ii=1:lt
for jj=1:ln
U=sqrt(u(ii, jj,time)^2+v(ii, jj,time)^2); % Wind speed
Cd=(0.73+0.69*U)*0.001;
tx(ii, jj,time)=Cd*roh*U*u(ii, jj,time); % kg/m^3*m/s*m/s= N/m^2
ty(ii, jj,time)=Cd*roh*U*v(ii, jj,time);
end% endfor : jj
end% endfor: ii
w=360/(24*60*60);
f=2*w*sin(lat*3.14/180);
e=1020*f;
Tx=tx./e';Ty=ty./e';
[curlz,cav]=curl(lon,lat,Tx',Ty');
以下是部分数据风应力旋度数值:
0.0182572267949581 0.0182269983924925 0.0181967061944306 0.0181827046908438 0.0181712741032243 0.0181481866165996
0.0154243018478155 0.0154117662459612 0.0153989419341087 0.0153685966506600 0.0153694767504931 0.0153604303486645
0.0248317848891020 0.0247041592374444 0.0246273297816515 0.0245415698736906 0.0244666561484337 0.0243965256959200
0.0350043345242739 0.0347553621977568 0.0345738036558032 0.0343706710264087 0.0341351181268692 0.0339316744357348
0.0183938816189766 0.0181415984407067 0.0178442383185029 0.0175463585183024 0.0172517960891128 0.0169928316026926
-0.00631732866168022 -0.00656576268374920 -0.00682262610644102 -0.00702888052910566 -0.00725689250975847 -0.00748418923467398
-0.0233717169612646 -0.0234852423891425 -0.0235829651355743 -0.0236471877433360 -0.0237413165159524 -0.0238710944540799
-0.0358034828677774 -0.0357474037446082 -0.0356631401227787 -0.0355864280136302 -0.0355055633699521 -0.0354109063046053
-0.0374405410839245 -0.0371644741389900 -0.0368553367443383 -0.0365596762858331 -0.0362594716716558 -0.0359505359083414
-0.0280086565762758 -0.0276412987150252 -0.0272768270224333 -0.0269213747233152 -0.0265586813911796 -0.0261903186328709
-0.0182936834171414 -0.0180409485474229 -0.0177897945977747 -0.0175436786375940 -0.0172975352033973 -0.0170420845970511
-0.0101613746955991 -0.00997817842289805 -0.00979009224101901 -0.00961867347359657 -0.00942353997379541 -0.00924936542287469
-0.00495310872793198 -0.00478754378855228 -0.00463476497679949 -0.00448378501459956 -0.00434258766472340 -0.00419385125860572
-0.00368631258606911 -0.00355175882577896 -0.00346622848883271 -0.00335739972069860 -0.00327302562072873 -0.00318362098187208
-0.00408481340855360 -0.00402128929272294 -0.00399416917935014 -0.00399588793516159 -0.00398076139390469 -0.00398522615432739
-0.00478352326899767 -0.00474439887329936 -0.00475260056555271 -0.00478023756295443 -0.00482332892715931 -0.00486710201948881
-0.00516830757260323 -0.00520224496722221 -0.00528797972947359 -0.00537004508078098 -0.00547436624765396 -0.00560979638248682
-0.00541491014882922 -0.00557201355695725 -0.00574649777263403 -0.00594194000586867 -0.00614172825589776 -0.00635611871257424
|
|