请各位高手帮忙改一下生成混凝土搅拌车拌筒和螺旋线的GRIP程序的错误

DUHUIDUHUI

编译和链接都没有错误，只是在载入ug中时报错说数组超出范围，错误在第48行，已用红色标记
) `  ]. t& e. P  e3 t+ cNUMBER/r(3),h(3),bt(3),i,st,stm(3),al,a3,p1,p3
8 K1 w, ~. R' ]                                           $$1-定义变量
2 _& p# ^, Q/ Z3 gENTITY/ax,ps1(1001),ps2(1001),m(3),g(3),sl1,$
      sl2,pt(4),pt0(2)
pds1:                                      $$2-生成输入参数的对话框
+ }( T; ~4 ~7 _+ W; R- T# }/ _( P    PARAM/'请输入罐体和螺旋线参数:',$
      '前锥小端半径:',r(1),$
5 N% n# ?# q9 i3 y& Q. m! u      '前锥高度:',h(1),$
      '中柱半径:',r(2),$
+ g$ o! {1 K! K6 a3 {2 `      '中柱高度:',h(2),$
      '后锥小端半径:',r(3),$
      '后锥高度:',h(3),$
      '前锥螺旋角:',bt(1),$
      '中锥螺旋角:',bt(2),$
      '后锥螺旋角:',bt(3),rsp
! `8 [8 F6 v1 C% ~" }: z7 B. b   JUMP/pds1:,stop1:,rsp
   ax=LINE/(pt0(1)=POINT/0,0,0),(pt0(2)=POINT/0,0,600)
% Q/ @$ |/ W& r8 y" q. u7 z8 L- L5 `                                           $$3-生成罐体
0 z" Y/ j. O4 x8 V9 A0 ^   g(1)=LINE/(pt(1)=POINT/r(1),0,0),(pt(2)=POINT/r(2),0,h(1))
   g(2)=LINE/pt(2),(pt(3)=POINT/r(2),0,(h(1)+h(2)))
   g(3)=LINE/PT(3),(pt(4)=POINT/r(3),0,(h(1)+h(2)+h(3)))
   m(1)=REVSRF/g(1),AXIS,ax,0,360
   m(2)=REVSRF/g(2),AXIS,ax,0,360
" ]# O# ]9 Z5 }( D6 E; ~( D   m(3)=REVSRF/g(3),AXIS,ax,0,360
   a1=ATANF((r(2)-r(1)/h(1)))          $$4-计算基本参数
; C1 c/ G" z. K: n1 G1 i2 J   a3=ATANF((r(3)-r(2)/h(3)))
   p1=(r(2)-r(1))/h(1)*COSF(a1)*COSF(bt(1))/SINF(bt(1))
   p3=(r(3)-r(2))/h(3)*COSF(a3)*COSF(bt(3))/SINF(bt(3))
   stm(1)=LOGF((r(2)-r(1))/r(1)+1)/p1
   stm(2)=stm(1)+h(2)/r(2)*SINF(bt(2))/cosf(bt(2))
& ~( R8 ~7 X; i1 f" w- v& S   stm(3)=stm(2)+LOGF((r(3)-r(2))/r(2)+1)/p3
' j1 f% f7 O0 w1 `; ]  N) G8 S/ _5 J   i=1
. \5 ~7 e0 f, G+ s$ R7 i2 W. |0 b   STR1:                              $$5-计算曲线上点的坐标
   st=(i-1)*stm(3)/1000
" y! r1 ]' X4 t! M- ?( @. {  IFTHEN/st<=stm(1)
   xc=r(1)*EXPF(p1*st)*COSF(st*180/3.14159)
   yc=r(1)*EXPF(p1*st)*SINF(st*180/3.14159)
   zc=r(1)*h(1)/(r(2)-r(1))*(EXPF(pl*st)-1)
4 ~4 y. J% v* Q/ l7 N  ELSEIF/stm(1)<st AND st<=stm(2)
1 \8 d5 W; H, k/ k4 D9 y   xc=r(2)*COSF(st*180/3.14159)
   yc=r(2)*SINF(st*180/3.14159)
+ D5 U8 C( s) G   zc=r(2)*COSF(bt(2))/SINF(bt(2))*(st-stm(1))+h(1)
  ELSEIF/stm(2)<st AND st<=stm(3)
   xc=r(2)*EXPF(p3*(st-stm(2)))*COSF(st*180/3.14159)
. \/ g8 s; C& i, [. y6 y   yc=r(2)*EXPF(p3*(st-stm(2)))*SINF(st*180/3.14159)
   zc=r(2)*h(3)/(r(3)-r(2))*(EXPF(p3*(st-stm(2)))-1)+h(1)+h(2)
$ P6 P- K) P# Y0 V  ENDIF
4 D6 G, f* }  A   ps1(i)=POINT/xc,yc,zc            $$6-生成螺旋线
   ps2(i)=POINT/-xc,-yc,zc
3 k  h5 C5 U4 j   i=i+1
: @$ N+ `) R1 t2 q, O. F1 X  IF/st<=stm(3),JUMP/STR1:
+ Q' A( e8 q  v) b  sl1=SPLINE/ps1(1..(i-1))
  sl2=SPLINE/ps2(1..(i-1))
! {( J3 h; C: d9 L BLANK/ax,g,pt,pt0,ps1(1..(i-1)),ps1(1..(i-1))
( w! H  d* m% r                                       $$结束处理
  Y+ F, _9 f' q  H) ^  stop1:
   HALT

27699885

亲自测试过这个程序，对数螺旋线grip吧，也曾出现过这问题qq137303928