ncl计算散度问题

20141170086

如图为计算200hPa高度的散度对应的ncl脚本，散度计算用了两种方法，一种是用中央差分求导进行计算(D1)，一种是直接用ncl的函数uv2dvF计算(div1)，图中分别为D1、div1和D1-div1。可以看出在中低纬度区域，两种方法计算差别不大，在中高纬度地区差异较大。特别是在极地地区，由于用中央差分对x方向求导时，dX在高纬度距离较短(见脚本第31行)，导致在极地区域求导后的值异常大，这样的方法在极地地区应该是有问题的。直接用ncl函数uv2dvF求散度得到的结果(div1)也有点奇怪，在中高纬度区域有奇怪的带状分布。在ncl安装路径下也没有找到uv2dvF函数对应的脚本，不知道它是怎么算的，所以想请问一下大家对lat和lon方向求导时是怎么求的，我脚本里的求导方法是不是有问题？谢谢大家！

墨家大宝

不会NCL。会不会因为地球其实是椭球导致你当正球算高纬处dx算得有偏差？而内置函数考虑了椭球？

伽蓝鸟

球坐标的散度公式跟直角坐标不一样。

Masterpiece

uv2dvF计算散度项的两个微分(梯度)时，用的是球面上的球谐函数(spherical harmonics)，相比中央差分的方案有更快的速度和计算精度。