################################
#Template
#Must contain the variables
#  initconf 
#  outtraj
#  outfile
#  interaction
#  run_length
#  samplingrate
#  particle



# VARIABLES

variable iname index test/out.xyz    
variable fname index test/out2.dcd
variable interactions index test/interactions

# Initialization
#correspond to x=y=z=1
lattice fcc 4
units		lj
#boundary	f f f
boundary	s s s
atom_style	molecular
log 		log.txt
read_data	${iname}
neighbor 2.0 multi

#########################
#Interactions

include test/interactions
group particle type 1 2

dump 10 all dcd 10 ${fname}



#####################################################
# Equilibration (Langevin dynamics at 5000 K)

velocity 	particle create 1.0 1231
fix		1 particle nve/limit 0.005
fix		addscale particle langevin 1.0 1.0  1 904297  


thermo_style	custom step temp 
thermo          10000
timestep	0.005 
run 1000

#####################################################

write_data test/out2.xyz


