- 积分
- 465
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-12
- 最后登录
- 1970-1-1
|
发表于 2015-9-30 11:15:13
|
显示全部楼层
本帖最后由 Aiiky 于 2015-9-30 11:16 编辑
目前都是参数化计算海面风应力,你可以参考下<大洋风应力和风场的计算与分析> 李克让等.然后选择以下公式.我师姐计算的matlab程序如下,
其中uv为风速大小,uwinds即东西方向速度,vwinds为南北方向速度.
for i = 1:xx
disp(i)
for j = 1:yy
for t = 1:tt
if uv(i,j,t) < 11
Cd = 1.2*10^-3;
ustress(i,j,t) = 1.2*Cd*uv(i,j,t)*uwinds(i,j,t);
vstress(i,j,t) = 1.2*Cd*uv(i,j,t)*vwinds(i,j,t);
elseif uv(i,j,t) >= 11 && uv(i,j,t) < 25
Cd = (0.49+0.065*uv(i,j,t))*10^-3;
ustress(i,j,t) = 1.2*Cd*uv(i,j,t)*uwinds(i,j,t);
vstress(i,j,t) = 1.2*Cd*uv(i,j,t)*vwinds(i,j,t);
elseif uv(i,j,t) >= 25
Cd = 2.115*10^-3;
ustress(i,j,t) = 1.2*Cd*uv(i,j,t)*uwinds(i,j,t);
vstress(i,j,t) = 1.2*Cd*uv(i,j,t)*vwinds(i,j,t);
end
clear Cd
end
end
end
你可以重点关注下参数
|
|