/title, Plasmawelding
/filname, thermal 
/prep7 
blc4,0.01,0,0.065,0.035,0.0006 
et,1,solid70
uimp,1,dens,,,3960, 
uimp,1,kxx,,,35, 
uimp,1,c,,,850,
type,1 
mat,1 
esize,0.001 
vmesh,all 
finish

/solu 
/psf,hflux,,2
l=0.001
v=0.050
t1=l/v 
r=0.002
tumg=20 
max_time=85
max_x=85
max_y=35

*dim,flux2,table,max_x,max_y,max_time*max_y/(2*r*100),x,y,time 
*do,y,1,max_y/2,2
*do,k,1,max_time,1 
*do,i,1,max_x,1 
flux2(i,0,k)=(i-1)*l 
*enddo 
*do,j,1,max_y,1 
flux2(0,j,k)=(j-1)*l 
*enddo 
*do,i,1,max_x,1 
*do,j,1,max_y-2,1 
xcenter=v*(k-1)*t1 
ycenter=(y-1)*r
distance=sqrt(((i-1)*l-xcenter)**2+((j-1)*l-ycenter)**2) 
*if,distance,le,r,then 
flux2(i,j,(y-1)*max_time+k)=8000000
flux2(i,j+2,(y+1)*max_time-k)=8000000 
*else 
flux2(i,j,(y-1)*max_time+k)=0
flux2(i,j+2,(y+1)*max_time-k)=0 
*endif 
*enddo 
*enddo
flux2(0,0,(y-1)*max_time+k)=((y-1)*max_time+k-1)*t1 
flux2(0,0,(y+1)*max_time-k)=((y+1)*max_time-k-1)*t1
*enddo
*enddo

antype,trans 
tunif,tumg 
*do,tm,0.02,22.6,t1 
time,tm 
deltim,0.01,0.01,0.01 
tsres,erase 
asel,s,loc,z,0.0006,0.0006 
sfa,all,1,hflux,%flux2%
solve 
*enddo 
finish