rc = 3.6 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) = 4.290 re(1,1) = 2.4800 alpha(1,1) = 5.15716154 repuls(1,1) = 0.05 attrac(1,1) = 0.05 Cmin(1,1,1) = 0.36 Cmax(1,1,1) = 2.80