高斯热源问题

不死局

我在加载高斯热源时候出现这样的问题应该怎样修改，着急用，做毕业课题的，大家帮帮忙，程序如下

/filname,laster
/title,laster
/prep7
et,1,solid70
mptemp,1,20,200,400,800,1200
mp,dens,1,7850
mpdata,c,1,1,461,517,613,895,652
mpdata,kxx,1,1,41.87,38.16,31.82,21.43,29.64
mp,dens,2,7825
mpdata,c,2,1,472,559,550,569,605
mpdata,kxx,2,1,44.17,41.88,36.45,22.17,25.66
block,0,0.0015,0,-0.006,0,0.03
block,0.0015,0.02,0,-0.006,0,0.03
cylind,0.0015,,0,0.03,0,90
vglue,1,2,3
lesize,11,,,0.006/0.00025
lesize,22,,,0.006/0.00025
lesize,7,,,0.0015/0.00025
lesize,37,,,0.0015/0.00025
lesize,8,,,0.006/0.0005,5
lesize,1,,,0.006/0.0005,5
lesize,18,,, 0.006/0.0005,1/5
lesize,15,,, 0.006/0.0005,1/5
lesize,6,,, 0.006/0.0005,1/5
lesize,3,,, 0.006/0.0005,1/5
lesize,38,,, 0.0185/0.0002,5
lesize,39,,, 0.0185/0.0002,5
lesize,40,,, 0.0185/0.0002,5
lesize,41,,, 0.0185/0.0002,5
type,1
mat,1
meshkey,1
mshkey,1
mshape,0,3D
vmesh,1
vmesh,5
type,1
mat,2
mshape,0,3D
mshkey,1
vmesh,4
finish
/solu
antype,trans
trnopt,full
lumpm,off
nlgeom,off
outpr,basic,all
outres,all,all
toffst,273
tunif,20
nsel,s,loc,x,0
dsym,symm,x
asel,s,area,,2
!asel,a,area,,19
asel,a,area,,23
asel,a,area,,12
asel,a,area,,25
asel,a,area,,24
asel,a,area,,3
asel,a,area,,22
asel,a,area,,1
!asel,a,area,,18
nsla,s,1
sf,all,conv,65,20
length=0.03
width=0.04
height=0.006
rcoating=0.0015
pi=3.1415926
rlaser=0.0015
vlaser=0.004
plaser=300
absor=0.3
esizeunit=0.00025
tinc=0.0625
tm=length/vlaser
maxx=rcoating/esizeunit+1
maz=length/esizeunit+1
maxtime=tm/tinc+1
!maxflux=3*plaser*absor/(pi*rlaser*rlaser)
*dim,laserflux,table,maxx,maz,maxtime,x,z,time
*do,k,1,maxtime,1   !按时间进行循环
    *do,i,1,maxx,1
    laserflux(i,0,k)=(i-1)*esizeunit !定义第零列的值
    *enddo
    *do,j,1,maz,1
    laserflux(0,j,k)=(j-1)*esizeunit
    *enddo
     *do,i,1,maxx,1
     *do,j,1,maz,1
      xcenter=0
      zcenter=vlaser*(k-1)*tinc+esizeunit/2
      distance=sqrt(((i-1)*esizeunit-xcenter)**2+((j-1)*esizeunit-zcenter)**2)
*if,distancem,le,rlaser,then
   laserflux(i,j,k)=0.3*plaser/(pi*rlaser**2)/exp(2*distance**2/(Rlaser**2))
!maxflux*exp(3*distance**2)/(rlaser**2))  !absor*plaser/(pi*rlaser**2)
    *else
   laserflux(i,j,k)=0
   *endif
   *enddo
  *enddo
laserflux(0,0,k)=(k-1)*tinc  !定义每个时间步的值
*enddo
*do,t,1e-6,tm,tinc
nsel,s,loc,z,0,length            
nsel,r,loc,x,0,rcoating
nsel,r,loc,y,0,rcoating
sf,all,hflux,%laserflux%
nsel,all
esel,all
time,t
autots,off
nsubst,2,5,1
kbc,1
tsres,erase
outres,all,all
solve
*enddo
outres,all,all
finish