/prep7  

/units,mpa   

*set,x0,0  

*set,y0,0

*set,rb,0.02  

*set,Q,500000   

*set,PI,3.1415926

*set,tim,0

*set,A,0.09

QC=2*A*Q/(PI*rb*rb)        

et,1,solid70  

mp,dens,1,7.97 

mp,c,1,465  

mp,kxx,1,54  

block,0,0.52,-0.125,0.125,-0.05,0    

esize,0.01,    

mshkey,1  

vmesh,all   

save

x_center=x0  

y_center=y0

QC*exp(-2*(({X}-x_center)^2+({Y}-y_center)^2)/(rb^2))

/sol

antype,4,new  

trnopt,full

toffst,0 

tunif,298  

*do,i,1,26

tim=tim+0.0004

time,tim 

x_center=x0+(i-1)*rb 

y_center=y0

*set,laser(4,0,1),x_center 

nsel,s,loc,x,x_center-0.02,x_center+0.02

nsel,r,loc,y,y_center-0.02,y_center+0.02

nsel,r,loc,z,0

sf,all,hflux,%laser% 

allsel,all

nsubst,5

solve

*enddo