风应力旋度散度 求指导

QQT

1.div.mat文件是我所求得风应力散度的数值，量级是10^4，div.mat文件中已经乘了10^4.所选的区域是太平洋赤道附近-10S~10S.
但是imagesc显示时有错误得不到想要的结果  如图。

，求指点。
2.求散度的公式dx方向考虑地球曲率的影响；dy方向不需要考虑曲率影响。
我的计算公式：
dx=longitud(a)-longitude(b);
dx=dx.*(pi./180).*R.*cos(latitude*pi./180)；
dy=latitude(a)-latitude(b);
dy=dy.*(pi./180).*R
du=u(a)-u(b);
dv=v(a)-v(b);
div=du/dx+dv/dy;
curl=dv/dx-du/dy;
不知道是否正确求指点