rc = 4.0
delr = 0.1
augt1 = 0
erose_form = 2
ialloy = 2

zbl(1,1) = 0
nn2(1,1) = 1
rho0(1) = 1.000
Ec(1,1) = 8.090
re(1,1) = 2.8600
alpha(1,1) = 4.92761091
repuls(1,1) = 0.00
attrac(1,1) = 0.00
Cmin(1,1,1) = 0.25
Cmax(1,1,1) = 2.80