风应力wind stress 和风应力旋度 wind stress curl

navy_yehaijun

最近在做风应力wind stress 和风应力旋度 wind stress curl的计算，查看了wlzhongouc的帖子--NCEP风场求散度和涡度。并下载了其中的程序。因为数据来源不同（wlzhongouc的帖子没有给出数据），我做出来的风应力wind stress 和风应力旋度 wind stress curl总感觉不对，请各位指导一下。

数据为ASCAT。
ASCAT ocean surface winds
http://oceanwatch.pifsc.noaa.gov/las/servlets/constrain?var=76


clear;clf;clc;
ncload 2013090100_2013100100_monthly-ifremer-L3-MWF-GLO-20131028095330-01.0.nc
size=size(eastward_wind);
ux=reshape(eastward_wind(1,:,:),size(3),size(4));
ux=double(ux);
vy=reshape(northward_wind(1,:,:),size(3),size(4));
vy=double(vy);
speed=reshape(wind_speed(1,:,:),size(3),size(4));
speed=double(speed);
longitude=double(longitude);latitude=double(latitude);
LATLIMS=[min(latitude) max(latitude)];  LONLIMS=[min(longitude) max(longitude)];
[x,y]=meshgrid(longitude,latitude);


rho_air=1.23; % kg/m^3
cd=1.2e-3; % non-dimensional drag
CD = (0.73 + 0.069*ux).*10^-3;


wmag = (ux.^2+vy.^2).^.5;% wind speed
% taux=rho_air*cd*wmag*u; %Zonal Wind Stress N/m^2
% tauy=rho_air*cd*wmag*v; %wind stress curl"/units="N/m^2


taux=rho_air*cd*wmag.*ux;
tauy=rho_air*cd*wmag.*vy;
[curlz,cav]=curl(longitude,latitude,taux,tauy);
pcolor(x,y,curlz);


风应力旋度的大小一般为10-4（5），而我计算出来的约为10-3，多了一个量级。
多谢指导，欢迎交流！！

Aiiky

我以前写过一个计算 风应力旋度 的matlab程序和例子，见附件。希望对大家有帮助~
2017-5-28，更新下对程序的输入的补充：输入的风应力taux和tauy的格式是lat*lon，x轴对应lon，y轴对应lat，而非lon*lat。否则会报错~

tsl08

计算风应力旋度的话，不是应该将经纬度转化为实际距离的么？

QQT

我也在求风场旋度散度可否交流一下？

tsl08

楼主搞明白了没有啊~

沙颖凯

http://bbs.06climate.com/forum.php?mod=viewthread&tid=20665

daweilee1989

牛人  顶一个 好东西

麦田_smile

请问如何求风应力旋度呢？是否有公式？

麦田_smile

QQT 发表于 2014-6-26 19:33
我也在求风场旋度散度可否交流一下？

请问如何求风应力旋度？之前没有计算过

QQT

麦田_smile 发表于 2014-10-15 21:02
请问如何求风应力旋度？之前没有计算过

二维旋度公式

Gavin

wind stress curf 有没有数据可以下啊？