请问关于GRADS画那个湿位涡图的问题(提问更新~)

wtf2ftw

直接套用了版主的帖子内容http://bbs.06climate.com/forum.php?mod=viewthread&tid=9884&highlight=%CA%AA%CE%BB%CE%D0。用的是FNL资料  想绘制湿位涡图，可是总是出错，求指点~  图都是像下图所示，自己找了半天没找到可能的原因

（问题更新）没人啊？~泪奔ing~

                               
登录/注册后可看大图

目前的情况是：根据版主帖子内容分两步，1是导出算好的假相当位温.dat。2是据此算出MPV
其中第1步我之前算错了，先已改好，并确定第1步生成的dat已经正确，就是在第2步出问题，总也画不对,画的图都如下图

算mpv的GS我贴其中的一段（就是版主帖子里的，基本没改），出图没用GS，直接在GRADS里输入的d mpv1，请大家帮我看看
算MPV的gs：
'reinit'
'set gxout fwrite'
'set fwrite e:\data1\weiwo\mpv02.dat'
'open E:\data1\weiwo\thse02.ctl'
**************
i=1
while(i<=8)
'set lon 110 130'
'set lat 25 40'
'set t 'i''
'set lev 950'
'define vo=hcurl(uwnd,vwnd)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(1000-900)'
'define dthse=thse(lev=1000)-thse(lev=900)'
'define du=uwnd(lev=1000)-uwnd(lev=900)'
'define dv=vwnd(lev=1000)-vwnd(lev=900)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(thse,x)'
'define dty=cdiff(thse,y)'
'define mpv1=-g*(vo+f)*dthse/dp'  
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define mpv=mpv1+mpv2'
'd mpv1*1000000'
*****下面是算另外的高度，内容重复，省略****
..............

i=i+1
endwhile
'disable fwrite'
;
上述GS用到之前的thse02的dat的描述文件如下
DSET E:\data1\weiwo\thse02.dat
TITLE Upper Data
undef -9.99E33
xdef 21 linear 110 1.000000
ydef 16 linear 25 1.0
zdef 16 levels 1000 950 900 850 800 750 700 650 600 550 500 450 400 350 300 250
tdef 8 linear 18z10feb2003 6hr
vars 3
thse 16 99 temps
uwnd 16 99 u-wind
vwnd 16 99 v-wind
Endvars

（问题再更新）
又在论坛里搜了一遍，仔细看了一下之前看过却没认真看的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=10190
结果一试就OK了，囧，这个已经满足我的需要了
不过我问的这个为什么不行，老出错，我做完了事再慢慢想吧，要是有前辈指点那不胜感激，总之，谢谢气象家园~

north23

帮顶了