任意两点间的垂直剖面和时间循环

Maziy

菜鸟一个，下面这个gs是自己写的，但是运行的时候老是提示 out of buffer space ，尝试过如果将垂直剖面和时间循环分开来分别运行是可行的。实在是找不到原因，请各位大神指点
********************************************************8
'reinit'
'open **.ctl'
'set mpdset hires'
'set display color white'
'set grads off'
frame = 1
endtime=25
runscript = 1
dis_t = 1
r=1
lat1=0
lat2=0.004
lon1=0
lon2=0.004
'set x 1'
'set y 1'
'set z 1 19'


while(runscript)
'set t ' dis_t
lon=lon1
'collect 'r' free'


while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect  'r' gr2stn(dbz,'lon','lat')'
lon=lon+0.0001
endwhile
'set lon 0 0.004'
'set z 1 19'


'set gxout shaded'
'd coll2gr('r',19)'


if(ans)
'printim 'frame'.gif gif '

frame=frame+1
endif


pull dummy


if ( dis_t=endtime )
runscript=0
endif
dis_t = dis_t + 1
r=r+1


endwhile