Matlab计算风应力旋度-程序分享

Aiiky · 发表于 2018-10-25 09:15:16

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 Aiiky 于 2021-12-30 21:51 编辑

最近需要计算风应力旋度curlz，update下很早以前写的matlab程序，和大家分享下~
z方向的风应力旋度curlz的计算公式：

curlz = dv/dx - du/dy
d这里应该是偏导符号（不知道如何插入公式就直接这样打了~），u=taux，v=tauy分别为x和y方向的风应力。
程序默认的地球半径为6371 km。在Matlab 2018a上测试没有问题，程序见附件，说明请doc curlz_myself。

希望对大家有帮助~
2021-12-30
更新了下程序，这次输入的taux和tauy格式为lon*lat，而非lat*lon了，结果也为lon*lat格式，这样更加方便。
请查阅说明后使用，doc curlz_myself。

Aiiky · 发表于 2020-2-17 14:49:12

不清楚怎么才能修改帖子。最近更新了下curlz_myself程序，和之前的程序相比，程序没有变化，只是增加了更多的解释和说明。需要的同学请下载这个最新的程序~

minmin · 发表于 2018-10-25 10:18:37

{:eb502:}有用

Aiiky · 发表于 2018-11-1 09:34:45

minmin 发表于 2018-10-25 10:18
有用

嘿嘿

Givx · 发表于 2020-9-30 17:54:24

向大佬学习

MengXC · 发表于 2020-10-11 12:02:28

您好，想向您请教一下计算风应力旋度后，是不是可以计算ekman抽吸速度的大小啊！是否可以进一步讲一下ekman抽吸速度呀

Aiiky · 发表于 2021-12-30 21:57:40

MengXC 发表于 2020-10-11 12:02
您好，想向您请教一下计算风应力旋度后，是不是可以计算ekman抽吸速度的大小啊！是否可以进一步讲一下ekman ...

Ekman曾底（D_E）的速度为：

w(D_E) = curlz( tau / ( rho * f ) )
公式可以参考董老师《物理海洋学导论》103页。

因此，你可以在我提供的程序中，将taux/( rho * f )和tauy/( rho * f)带入输入，去计算curlz即可。
注意，rho你可以取常数，1024，f则是一个lon*lat的矩阵。注意单位。

Zahor · 发表于 2022-3-30 20:02:16

ganxieganxie

yylxbsd · 发表于 2024-2-18 22:41:07

很有用，非常感谢

13406040474 · 发表于 2024-3-26 21:09:30

谢谢博主很有用

		自动登录	找回密码
密码			立即注册

[源程序] Matlab计算风应力旋度-程序分享

登录后查看更多精彩内容~