风应力旋度量级不对

zzw9701

计算出来的风应力旋度太大了，有没有大神能帮忙看看问题出在哪里，帮忙改一下？
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

76481786

支持一下啊啊啊啊

小涡牛

经纬度应该转换为实际的地面距离（单位：m）