wrfout数据计算实测风涡度

大气太阳茄

使用wrf模拟结果计算涡度，模拟中dx=dy=5km。

请问，Fortran程序中计算涡度：
vor(i,j,iz)=1./(dd*2)*( (v(i+1,j,iz)-v(i-1,j,iz))-(u(i,j+1,iz)-u(i,j-1,iz)) )

还是用
vor(i,j,iz)=mf(i,j)/(dd*2)*( (v(i+1,j,iz)-v(i-1,j,iz))-(u(i,j+1,iz)-u(i,j-1,iz)) )


（时间循环写在最外面，mf是地球放大系数，dd取了5000）

我试了以上的公式，结果用grads画图出现的结果特别不正常，量级直接达到10的14次方。请大家帮忙看下公式有什么问题？

lichiqin4

每个风速下面还要除以地球放大系数