请大家都来帮帮忙！

zhungki

各位大侠，小弟最近在用高斯移动热源的时候，多加了一个时间条件结果反复调试也没能成功，运行的过程总是出现"Invalid branch label(and) on *if commond The *if command is ignored ”小弟把程序发到上面请各位大侠指点帮忙！
/prep7
/triad,off
et,1,solid70
keyopt,1,2,1        !用单元表面温度来评价对流系数
et,2,surf152
r,2,5.67e-8        !设定斯蒂芬－波尔兹曼常数
keyopt,2,4,1        !无中间节点
keyopt,2,5,0        !无外部节点
keyopt,2,8,3        !用单元表面温度来评价对流系数


pi=3.1415926
ef1=0.25
ef2=0.55
I=97 !焊接电流
U=19.2 !焊接电压
r1=0.00117
r2=0.005
Qmax1=3*ef1*U*I/(pi*r1**2)
Qmax2=3*ef2*U*I/(pi*r2**2)
y_dist=0.04
x_dist=0.08
z_dist=0.005
lsize=0.001 !单元网格长度
v=0.002 !速度
tm=x_dist/v !焊接总时间
tinc=lsize/v !单步时间长度

max_time=1+tm/tinc
max_x=1+x_dist/lsize
max_y=1+y_dist/lsize

/view,1,1,1,1
block,0,0.08,0,0.005,0,0.005
block,0,0.08,0.005,0.015,0,0.005
block,0,0.08,0.015,0.025,0,0.005
block,0,0.08,0.025,0.04,0,0.005
allsel
vglue,all
numcmp,all
save

mshkey,1
vatt,1,1,1
type,1
esize,lsize
vmesh,1
esize,2*lsize
vmesh,2
esize,4*lsize
vmesh,3
esize,8*lsize
vmesh,4
eplot
/image,save,mesh,bmp

type,2
asel,s,,,1,21
asel,u,,,4
asel,u,,,7,8
amesh,all
/pnum,line,off
/pnum,defa
eplot


*dim,flux2,table,max_x,max_y,max_time,x,y,time,t
*do,k,1,max_time,1
  *do,i,1,max_x,1
   flux2(i,0,k)=(i-1)*lsize
  *enddo
   
  *do,j,1,max_y,1
   flux2(0,j,k)=(j-1)*lsize
  *enddo
   
  *do,i,1,max_x,1
    *do,j,1,max_y,1
      *do,t,1e-6,tm,1
     xcenter=v*(k-1)*tinc
     ycenter=0
     distance=sqrt((abs((i-1)*lsize-xcenter))**2+((j-1)*lsize-ycenter)**2)
    *if,distance,le,r1,and,mod(t,3),NE,0,then
        flux2(i,j,k)=Qmax1*(1-distance/r1)+Qmax2*exp(-3*distance**2/(r2**2))
    *else
    *if,distance,GE,r1,and,distance,le,r2,and,mod(t,3),EQ,0,then
        flux2(i,j,k)=Qmax2*exp(-3*(distance**2)/(r2**2))
    *endif
      *endif
    *enddo
  *enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
finish