# Test of MEAM potential for SiC system

units		lj
boundary	p p f

atom_style	atomic

read_data	cipolla
mass		1 1.0
mass		2 1.0

pair_style	lj/cut 6.17

group		top type 2
group		bottom type 1
pair_coeff	1 1 0.25  2.67269615442102
pair_coeff	1 2 0.83 1.78179743628068
pair_coeff	2 2 0.25 2.67269615442102

neighbor	1.2 bin
neigh_modify	delay 1

fix             coord all coord/original


fix		1 top rigid single 
fix		2 bottom nve
fix 		3 bottom langevin 0.0 0.0 0.001 2134

		
thermo		100
timestep	0.001

compute         disp all displace/atom coord 
compute         disp1 all reduce sum c_disp[1]

thermo_style    custom step temp pe ke etotal vol c_disp1


dump		1 all custom 100000 dump.meam tag type x y z

run		100000
write_restart   osvaldo