ekman抽吸速度计算

荒火啊啊啊

ekman抽吸速度计算公式

荒火啊啊啊

rho_water=1025;
f=zeros(size(lat,1),1);%f只与sinfai有关，所以经度的间隔也与纬度有关
Llat=zeros(size(lat,1)-1,size(long,1));%不同纬度之间的间隔
Llon=zeros(size(lat,1),size(long,1)-1);%不同经度之间的间隔
for i=1:size(f,1)
    f(i)=sw_f(lat(i));
end
for i=1:size(Llat,1)
    for j=1:size(Llat,2)
    Llat(i,j)=sw_dist([lat(i),lat(i+1)],[long(j),long(j)],'km');
    end
end
for i=1:size(Llon,1)
    for j=1:size(Llon,2)
    Llon(i,j)=sw_dist([lat(i),lat(i)],[long(j),long(j+1)],'km');
    end
end
%%计算ekman抽吸的速度（利用三点法计算）%%
ekw=zeros(size(lat,1)-2,size(long,1)-2);
for i=1:size(lat)-2
    for j=1:size(long)-2
    ekw(i,j)=(((tauy(i+2,j+1)-tauy(i,j+1))/(f(i+1)*rho_water))/((Llat(i,j)+Llat(i+1,j))*1000)) ...
        -(((taux(i+1,j+2)-taux(i+1,j))/f(i+1)*rho_water)/((Llon(i,j)+Llon(i,j+1))*1000));   
    end
end
a=nan(size(ekw,1),1);
ekw=[a ekw a];%%重新定义数组


能力有限 上面的也不知道对不对

陶蜜云

请问楼主，Ekman抽吸速度和Ekman流速度是不一样的吧

Leona-lu

楼楼 你怎么没用到curl函数啊