label loopa clear units real atom_style charge region universe block -10 10 -10 10 -10 10 units box lattice sc 3 create_box 1 universe mass 1 12.09 create_atoms 1 box pair_style reax/c NULL checkqeq no pair_coeff * * ffield.reax.cho C run 0 run 0 jump SELF loopa