求助关于hdivg函数

asu

各位大侠：
存在问题如下
1.偶然发现fortran编程得到的散度与grads中hdivg函数得到的散度不一样。
2.发现fortran编程得到的散度与利grads中利用cdiff函数得到的散度一致。
因此求问为什么二者不一样？

grads计算如下：
'reinit'
'open ../data/u.ec.ctl'
'open ../data/v.ec.ctl'
'enable print ../gmf/test.gmf'
'set x 1 94'
'set y 1 41'
'set lev 500'
'set t 20'
'define dux=cdiff(t,x)'
'define dvy=cdiff(t.2,y)'
'define dx=6370000*cos(lat*3.14159/180)*cdiff(lon,x)*3.14159/180'
'define dy=6370000*cdiff(lat,y)*3.14159/180'
'define tx=dux/dy'
'define ty=dvy/dy'
'set x 2 93'
'set y 2 40'
'set gxout shaded'
'd (tx+ty)*100000'
'set x 1 94'
'set y 1 41'
'define a=hdivg(t,t.2)'
'set gxout contour'
'set x 2 93'
'set y 2 40'
'd a*100000'
'cbarn'
'print'
'disable print'
;

图如下：

river

去看下官方文档吧http://grads.iges.org/grads/gadoc/gadocindex.html

sjysjxyf66

这个图做的好漂亮，要学习一下