海洋wind引起的艾克曼抽吸Ekman Pumping Veolcity （EPV）的计算

navy_yehaijun

如题所示，海洋中风引起的艾克曼抽吸的计算需要用到风向vector、风应力wind stress、wind stress cur。
然而本站中关于这些方面的程序都不是很全，也不一定是正确的。curl函数是在用于这些计算时需要注意的事项又是什么呢？
希望大家给我一些好的意见。
谢谢

下面是我目前能做到的，肯定不正确，希望有人能够帮助我。

clear;clf;clc;
ncload 2011033100_2011040100_daily-ifremer-L3-MWF-GLO-20110402103252-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.25; % kg/m^3
cd=1.2e-3; % non-dimensional drag  1.5*10-3
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; % Meridional wind stress /units="N/m^2
% curl_tau = tauy-taux; % windstress curl = N/m^3  

taux=rho_air*cd*wmag.*ux;
tauy=rho_air*cd*wmag.*vy;
% EPV=curl(tau/rho_water*f); %rho_water为海水密度1020kg/m^3；f为coriolis parameter
rho_water=1020;f=2*7.292*0.00001*sin(deg2rad(30));

[curlz,cav]=curl(longitude,latitude,taux./(rho_water*f),tauy./(rho_water*f));
% curll=curl(longitude,latitude,taux./(rho_water*f),tauy./(rho_water*f));
m_proj('miller','lon',[108 122],'lat',[10 25]);  

m_pcolor(x,y,curlz);
shading interp;
colorbar

我想想看

你处理的时候，应该经纬度对应的单位不是m吧，这一块你有处理的方法吗？

Leona-lu

不知道为什么我用curl函数出来全是NaN Inf