NCL和grads计算涡度和散度的函数差异很大

yanzi

uv2dv_cfd是不是需要固定格点才能计算？

yanzi

uv2dv_cfd是不是需要固定格点才能计算？

3711313330

楼主自己编程中用的散度公式不对，球坐标下的散度公式不是这样写的

sun晒

yuxiao19910 发表于 2016-4-11 17:41
NCL计算散度的公式为：   
dv = Dv/Dy + Du/Dx -(v/a)*tan(phi)
dv(j,i) = (v(j+1,i)-v(j-1,i))/dy2(j)+ ...

请问这个ncl的计算过程在哪找的？

irish

xo2001 发表于 2013-6-20 22:21
最近想用NCL直接计算就省得编程了，所以就想把几个常用的函数了解了解~谁知道刚做个散度涡度就有问题...

有心了~很棒~

QQtangxian

太感谢啦，谢谢啦

兔子爱吃萝卜

sun晒 发表于 2016-10-14 10:19
请问这个ncl的计算过程在哪找的？

https://www.ncl.ucar.edu/Documen ... -in/uv2dv_cfd.shtml，函数介绍里有。

兔子爱吃萝卜

sun晒 发表于 2016-10-14 10:19
请问这个ncl的计算过程在哪找的？

https://www.ncl.ucar.edu/Documen ... -in/uv2dv_cfd.shtml，函数介绍里有。

vicky510475733

google翻译的ncl散度计算官网，或许可以仔细看看下面提到的文章，应该会详细说明计算方法的。
使用居中有限差分计算纬度 - 经度网格的发散度。网格不必是全局的，并且允许 缺少值。

根据HB Bluestein的说法 [ Midlatitudes中的Synoptic-Dynamic Meteorology，1992，Oxford Univ。按p113-114]，让D代表偏导数，a代表地球半径，phi代表纬度，dx2 / dy2代表适当的纵向和纬度间距。然后，让j为纬度y-subscript，我是经度x-subscript：

    dv = Dv / Dy + Du / Dx  - （v / a）* tan（phi）

    dv（j，i）=（v（j + 1，i）-v（j-1，i））/ dy2（j）
              +（u（j，i + 1）-u（j，i-1））/ dx2（j）
              - （v（j，i）/ a）* tan（phi（j））
最后的术语说明了球体上经络的收敛。网格可以是全局的或代表有限的区域。如果 存在缺失值，则函数会在所有可能的点执行计算，但会在无法使用的坐标处返回 缺失值。

如果要在边界处估计发散（boundOpt = 2或3），则使用垂直于适当边界的单侧差分方案来确定边界值。

伽蓝鸟

球坐标下散度的公式不是你第三个自己编程的那样。