关于螺旋度的分析与计算，欢迎大家一起讨论

圈圈圆圆圈圈

去年计算螺旋度的时候遇到了很多问题，也看了论坛里的很多帖子，今年重新捡起来研究了一下，发现有些帖子的计算还是有些问题，在此写下一些心得，欢迎大家一起讨论。
首先明确一下螺旋度的概念：速度矢和涡度矢点积的三重积分，其中速度矢和涡度矢的点积称为局地螺旋度，很多论文都有这个概念，在此不再枚举。

而我们常常用到的都是局地螺旋度，局地螺旋度包括水平螺旋度和垂直螺旋度两部分，其公式展开如下：


根据这个定义，可以发现这个公式和空气密度是没有关系的，之所以很多程序中加了12.64，是参考了冉令坤（2009）的算法，其定义和物理意义如下：


至此，很多程序出现了另一个可能引起错误的地方，即w和omega的转换，也就是给公式乘上了-10。严格来说w=-10omg是不对的，但可能计算结果近似，在此不多讨论。但很多问题就出现在了+w和-omg表示上升运动，-w和+ome表示下降运动，所以用w（omg）计算的话，+H表示气旋涡度的向上（向下）输送，其物理意义完全相反，当然螺旋度单位和量级也不同。

我用论坛里的程序改了一个NCL程序，发现画出来暴雨区的螺旋度是正值，跟冉令坤的说法完全相反，这才找到原因，是w和omg的转换问题。因此我完全参考了螺旋度的原始公式，得到的结果也不错，和未来降水落区确实有很好的对应关系。程序很简单，我把主体计算语句放在下方，供大家讨论和参考。


另外我浅浅看了一眼论坛里计算水平螺旋度的程序，我不是很明白为什么程序只读取了u和v的数据。根据我的理解，计算水平涡度应该会用到垂直速度，还请广大网友指教。

圈圈圆圆圈圈

调整了半天不知道怎么在文字中插入图片，也调整不了图片的顺序，各位海涵