帮帮忙，看看这个宏命令哪里错了

ST大豆

diameter=0.002
length=0.01
velocity=0.04
time=diameter/velocity
csys,1
WPCSYS,,1
/solu
!nropt,full,,on
j=length/diameter
Q=0
wpstyle,,,,,,,,1
wpstyle,0.005,0.005,-1,1,0.003,1,2,,5
csys,4
solcontrol,on
timint,on
tunif,300
zint=0
*do,i,1,j
nsel,s,loc,x,diameter/2
nsel,r,loc,y,0,360
nsel,r,loc,z,zint,zint+0.002
sf,all,hflux,41101273
allsel,all
eplot
!*
time,Q+time
autots,1
nsubst,10,,,1
kbc,1
!*
tsres,erase
!*
lswrite,i,
nsel,s,loc,x,diameter/2
nsel,r,loc,y,0,360
nsel,r,loc,z,zint,zint+0.002
sfdele,all,hflux
allsel,all
zint=zint+0.002
Q=Q+time
*enddo
!
!solve

ST大豆

模拟激光加热工件的，加载进去整体都是一个温度