gs运行时发生闪退，求支招！

schliezer

各位大神好，小白正在算相对螺旋杜SRH，输入gs按回车后发生闪退，检查半天不知何故，请各位大神帮忙看看，附gs：
'reinit'
'open D:\gradsdata\ncepfnl_1_6h\fnl_20120719_00_00_c.ctl'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lat 37 40'
'set lon 102 107'
'define pi=3.1415'
tt=1
while(tt<=13)
'c'
'set t 'tt' '
zz=1
while(zz<=17)
'set z 'zz''
'define alfa=atan2(ugrdprs,vgrdprs)'
'define V=mag(ugrdprs,vgrdprs)'
'define beta=ave(alfa,lev=850,lev=300)'
'define vave=ave(V,lev=850,lev=300)'
'define cu=vave*sin(beta+2*pi/9)'
'define cv=vave*cos(beta+2*pi/9)'
zz=zz+1
endwhile
（到这里基本都是对的，但加上以下的循环就出这问题）
'define SRH=0'
zz=1
while(zz<=16)
'set z 'zz''
'define H=(ugrdprs(z+1)-cu)*(vgrdprs(z)-cv)-(ugrdprs(z)-cu)*(vgrdprs(z+1)-cv)'
'define SRH=SRH+H'
zz=zz+1
endwhile

'set gxout contour'
'd SRH'
'printim D:\gradsdata\srh'%tt%'.gif white'
tt=tt+1
endwhile
;
另外，附SRH的计算公式：
  

schliezer

自顶，别沉啊。

idfall0803

想知道 +1 @@