我错在哪里，为什么运行不出来，急求答案

*dim,flux2,tab4,max_z,max_x,max_y,max_time,z,x,y,time
*do,k,1,max_time,1
*do,j,1,max_y,1
*do,h,1,max_z,1
flux2(h,0,0,k)=(h-1)*lsize
*enddo
*do,i,1,max_x,1
flux2(0,i,0,k)=(i-1)*lsize
*enddo
flux2(0,0,j,k)=(k-1)*tinc
*enddo
*do,j,1,max_y,1
*do,h,1,max_z,1
*do,i,1,max_x,1
zcenter=v*(k-1)*tinc
xcenter=lx/2
distance=sqrt(((h-1)*lsize-zcenter)**2+((i-1)*lsize-xcenter)**2)
*if,distance,le,g,then
flux2(h,i,j,k)=Qmax/exp(3*distance**2/(r**2))
*else
flux2(h,i,j,k)=0
*endif
*enddo
*enddo
*enddo
flux2(0,0,0,k)=(k-1)*tinc
*enddo